Population pharmacokinetics of mycophenolic acid and its 2 glucuronidated metabolites in kidney transplant recipients

Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation

Enterohepatic circulation (EHC) is a complex process where drugs undergo secretion and reabsorption from the intestinal lumen multiple times, resulting in pharmacokinetic profiles with multiple peaks. The impact of EHC on area under the curve (AUC) has been a topic of extensive debate, questioning the suitability of conventional AUC estimation methods. Moreover, a universal model for accurately estimating AUC in EHC scenarios is lacking. To address this gap, we conducted a simulation study evaluating five empirical models under various sampling strategies to assess their performance in AUC estimation. Our results identify the most suitable model for EHC scenarios and underscore the critical role of meal-based sampling strategies in accurate AUC estimation. Additionally, we demonstrate that while the trapezoidal method performs comparably to other models with a large number of samples, alternative models are essential when sample numbers are limited. These findings not only illuminate how EHC influences AUC but also pave the way for the application of empirical models in real-world drug studies.

A Population Pharmacokinetic-Pharmacodynamic Model of Navtemadlin, its Major Active Metabolite (M1) and Serum Macrophage Inhibitory Cykokine-1 (MIC-1)

Navtemadlin is a potent, selective, orally available inhibitor of murine double minute 2 that restores p53 activity to induce apoptosis in TP53 wild-type malignancies. Using richly sampled pharmacokinetic (PK) and pharmacodynamic (PD) data from healthy volunteers, a population PK/PD model was developed. A population PK (PPK) model described the PK characteristics of navtemadlin and its major metabolite acyl glucuronide (M1) and quantified enterohepatic recirculation (EHR). Post hoc individual PK parameters from this model were coupled with PD data for serum macrophage inhibitory cytokine-1 (MIC-1, GDF15), a cytokine biomarker of p53 activation, to construct a population PK/PD model that described plasma concentration-driven MIC-1 excursions and enabled simulation of the extent and duration of navtemadlin PD effects. The median apparent clearance (CL/F) and apparent central volume (V2/F) of navtemadlin were 36.4 L/hr and 159 L. The typical maximum stimulatory effect (S_max) was close to the median maximum MIC-1 ratio to baseline of 7.29 in observed data. Simulation revealed a dose-dependent increase of MIC-1 with steady state attained in approximately 7 days, in a 7-day-on/21-day-off dose regimen. Elevated MIC-1 concentrations persist through 17-19 days, leaving about 9-11 PD-free days in a 28-day cycle.

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.

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.

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.

Development of a joint population pharmacokinetic model of ezetimibe and its conjugated metabolite

Ezetimibe (EZE) is an extensively used antihyperlipidemic drug with an important cholesterol lowering activity. It undergoes extensive first-pass metabolism to form its active glucuronide metabolite (EZEG). Both drugs exhibit complex pharmacokinetic profiles attributed mainly to repetitive enterohepatic kinetics. The aim of the present study was the investigation of EZE and EZEG pharmacokinetics (PK), through the development of a joint population pharmacokinetic model able to characterize their kinetic processes and enterohepatic recirculation simultaneously. Concentration-time data derived from a bioequivalence study in 28 healthy subjects were used for the analysis. Population PK modeling was performed on the obtained data using nonlinear mixed effect modeling approach, where different methodologies were applied for the description of the complex metabolism and recirculation processes of the two compounds. EZE and EZEG concentrations were best described by a population PK model incorporating first-pass metabolism and an enterohepatic recirculation loop, accounting for the recycling process of the two moieties. This is the first joint population pharmacokinetic model describing the kinetics of both EZE and EZEG.

On the population pharmacokinetics and the enterohepatic recirculation of total ezetimibe

Ezetimibe is a potent cholesterol absorption inhibitor, with an erratic pharmacokinetic (PK) profile, attributed to an extensive enterohepatic recirculation (EHC). The aim of this study was to develop a population PK model able to adequately characterize the complex distribution processes of total ezetimibe. The analysis was performed on the individual concentration-time data obtained from 28 healthy subjects who participated in a bioequivalence study comparing two oral ezetimibe formulations. The population PK analysis was performed using nonlinear mixed effect modeling, where different EHC models were developed and evaluated for their performance. Total ezetimibe pharmacokinetics was best described by a four-compartment model featuring EHC through the inclusion of an additional gallbladder compartment, which was assumed to release drug at specific time-intervals consistent with food intake. The final PK model was able to adequately estimate the population pharmacokinetic parameters and to allow for a formal characterization of the pharmacokinetic profile and the secondary peaks due to enterohepatic recirculation.

Nonlinear relationship between enteric-coated mycophenolate sodium dose and mycophenolic acid exposure in Han kidney transplantation recipients

The aim of the research was to investigate the pharmacokinetics (PK) of enteric-coated mycophenolate sodium (EC-MPS) by quantification of the active metabolite of mycophenolic acid (MPA) after multiple escalating oral doses in Han kidney transplant recipients. A total of 28 Han postoperative kidney transplant recipients were given a multiple-dose of 540, 720 or 900 mg of EC-MPS two times a day in combination with tacrolimus for 6 days. Blood specimens were collected at each time point from 0 to 12 h after EC-MPS administration. MPA plasma concentrations were measured by UPLC—UV. The relationship between the EC-MPS dose and its PK parameters was assessed. In the range from 540 to 900 mg, C_max and AUC_0—12h did not increase with dose escalation. The AUC_0—12h, C_max, C₀ and T_max for the 540 720 and 900 mg doses were not significantly different, respectively (P >0.05). AUC_0—12 h and C_max were increased less than proportionally with increasing EC-MPS dose levels. Inter-individual variability in AUC_0—12h, C_max and C₀ were considerable. Nonlinear PK relationships were found from the doses of 540–900 mg of EC-MPS.

How accurate and precise are limited sampling strategies in estimating exposure to mycophenolic acid in people with autoimmune disease?

Mycophenolic acid (MPA) is a potent immunosuppressant agent, which is increasingly being used in the treatment of patients with various autoimmune diseases. Dosing to achieve a specific target MPA area under the concentration-time curve from 0 to 12 h post-dose (AUC12) is likely to lead to better treatment outcomes in patients with autoimmune disease than a standard fixed-dose strategy. This review summarizes the available published data around concentration monitoring strategies for MPA in patients with autoimmune disease and examines the accuracy and precision of methods reported to date using limited concentration-time points to estimate MPA AUC12. A total of 13 studies were identified that assessed the correlation between single time points and MPA AUC12 and/or examined the predictive performance of limited sampling strategies in estimating MPA AUC12. The majority of studies investigated mycophenolate mofetil (MMF) rather than the enteric-coated mycophenolate sodium (EC-MPS) formulation of MPA. Correlations between MPA trough concentrations and MPA AUC12 estimated by full concentration-time profiling ranged from 0.13 to 0.94 across ten studies, with the highest associations (r (2) = 0.90-0.94) observed in lupus nephritis patients. Correlations were generally higher in autoimmune disease patients compared with renal allograft recipients and higher after MMF compared with EC-MPS intake. Four studies investigated use of a limited sampling strategy to predict MPA AUC12 determined by full concentration-time profiling. Three studies used a limited sampling strategy consisting of a maximum combination of three sampling time points with the latest sample drawn 3-6 h after MMF intake, whereas the remaining study tested all combinations of sampling times. MPA AUC12 was best predicted when three samples were taken at pre-dose and at 1 and 3 h post-dose with a mean bias and imprecision of 0.8 and 22.6 % for multiple linear regression analysis and of -5.5 and 23.0 % for maximum a posteriori (MAP) Bayesian analysis. Although mean bias was less when data were analysed using multiple linear regression, MAP Bayesian analysis is preferable because of its flexibility with respect to sample timing. Estimation of MPA AUC12 following EC-MPS administration using a limited sampling strategy with samples drawn within 3 h post-dose resulted in biased and imprecise results, likely due to a longer time to reach a peak MPA concentration (t max) with this formulation and more variable pharmacokinetic profiles. Inclusion of later sampling time points that capture enterohepatic recirculation and t max improved the predictive performance of strategies to predict EC-MPS exposure. Given the considerable pharmacokinetic variability associated with mycophenolate therapy, limited sampling strategies may potentially help in individualizing patient dosing. However, a compromise needs to be made between the predictive performance of the strategy and its clinical feasibility. An opportunity exists to combine research efforts globally to create an open-source database for MPA (AUC, concentrations and outcomes) that can be used and prospectively evaluated for AUC target-controlled dosing of MPA in autoimmune diseases.

Development of population PK model with enterohepatic circulation for mycophenolic acid in patients with childhood-onset systemic lupus erythematosus

AIM

This study aimed to develop a population pharmacokinetic (PK) enterohepatic recycling model for MPA in patients with childhood-onset systemic lupus erythematosus (cSLE).

METHODS

MPA concentration-time data were from outpatients on stable oral mycophenolate mofetil (MMF) and collected under fasting conditions, with standardized meals (1 and 4 h post-dose). Sampling times were pre-dose, 20, 40 min, 1, 1.5, 2, 3, 4, 6 and 9 h, post dose. The population PK analysis simultaneously modelled MPA and 7-O-MPA-β-glucuronide (MPAG) concentrations using nonlinear mixed effect modelling.

RESULTS

PK analysis included 186 MPA and MPAG concentrations (mg l(-1)) from 19 patients. cSLE patients, age range 10-28 years, median 16.5 years were included. Mean ± SD disease duration was 3.8 ± 3.7 years. The final PK model included a gallbladder compartment for enterohepatic recycling and bile release time related to meal times, with first order absorption and single series of transit compartments. The PK estimates for MPA were CL(1) /F 25.3 l h(-1), V(3) /F 20.9 l, V(4) /F 234 l and CL(2) /F 19.8 l h(-1).

CONCLUSION

The final model fitted the complex processes of absorption and enterohepatic circulation (EHC) in those treated with MMF for cSLE and could be applied in Bayesian dose optimization algorithms.

Interaction between gentamicin and mycophenolate mofetil in experimentally induced pyelonephritis

Acute pyelonephritis (APN) is an inflammatory disease that leads to kidney malfunction. The objective of this investigation was to evaluate the impact of gentamicin (GEN) and ceftriaxone (CEF) alone and in combination with mycophenolate mofetil (MMF) on experimentally induced APN. Forty two Wistar male rats were assigned into seven groups +APN, APN +GEN, APN+CEF, APN+MMF, APN+GEN+MMF and APN+CEF+MMF. APN was induced by injecting E. coli in the left kidney. The control and +APN groups were treated with normal saline while the other APN groups received GEN, CEF, or MMF alone and/or in combination for 2 weeks. The elevated total white blood cells count and increased level of creatinine and blood urea nitrogen (BUN) in +APN groups returned to normal levels following 14 days treatment with GEN and CEF. Co-administration of GEN with MMF could not recover the APN-induced changes and resulted in a significant (P < 0.05) elevation of creatinine and BUN levels. Histopathological studies supported the biochemical findings as GEN and CEF alone could partly restore the APN-induced degeneration and leukocytic infiltration; however, the combination therapy of GEN plus MMF failed to reduce the APN-induced damages. The antibacterial susceptibility test demonstrated that the strain of E. coli used in this study was susceptible to GEN and CEF and the combination therapy did not change the antibacterial potency. These findings suggest that co-administration of GEN with MMF in APN may enhance kidney damage and the adverse effects of combination therapeutic regimen could be related partly to incompatibility of these compounds.

The evolution of population pharmacokinetic models to describe the enterohepatic recycling of mycophenolic acid in solid organ transplantation and autoimmune disease

With the increasing use of mycophenolic acid (MPA) as an immunosuppressant in solid organ transplantation and in treating autoimmune diseases such as systemic lupus erythematosus, the need for strategies to optimize therapy with this agent has become increasingly apparent. This need is largely based on MPA's significant between-subject and between-occasion (within-subject) pharmacokinetic variability. While there is a strong relationship between MPA exposure and effect, the relationship between drug dose, plasma concentration and exposure (area under the concentration-time curve [AUC]) is very complex and remains to be completely defined. Population pharmacokinetic models using various approaches have been proposed over the past 10 years to further evaluate the pharmacokinetic and pharmacodynamic behaviour of MPA. These models have evolved from simple one-compartment linear iterations to complex multi-compartment versions that try to include various factors, which may influence MPA's pharmacokinetic variability, such as enterohepatic recycling and pharmacogenetic polymorphisms. There have been major advances in the understanding of the roles transport mechanisms, metabolizing and other enzymes, drug-drug interactions and pharmacogenetic polymorphisms play in MPA's pharmacokinetic variability. Given these advances, the usefulness of empirical-based models and the limitations of nonlinear mixed-effects modelling in developing mechanism-based models need to be considered and discussed. If the goal is to individualize MPA dosing, it needs to be determined whether factors which may contribute significantly to variability can be utilized in the population pharmacokinetic models. Some pharmacokinetic models developed to date show promise in being able to describe the impact of physiological processes such as enterohepatic recycling. Most studies have historically been based on retrospective data or poorly designed studies which do not take these factors into consideration. Modelling typically has been undertaken using non-controlled therapeutic drug monitoring data, which do not have the information content to support the development of complex mechanistic models. Only a few recent modelling approaches have moved away from empiricism and have included mechanisms considered important, such as enterohepatic recycling. It is recognized that well thought-out sampling schedules allow for better evaluation of the pharmacokinetic data. It is not possible to undertake complex absorption modelling with very few samples being obtained during the absorption phase (which has often been the case). It is important to utilize robust AUC monitoring which is now being propagated in the latest consensus guideline on MPA therapeutic drug monitoring. This review aims to explore the biological factors that contribute to the clinical pharmacokinetics of MPA and how these have been introduced in the development of population pharmacokinetic models. An overview of the processes involved in the enterohepatic recycling of MPA will be provided. This will summarize the components that complicate absorption and recycling to influence MPA exposure such as biotransformation, transport, bile physiology and gut flora. Already published population pharmacokinetic models will be examined, and the evolution of these models away from empirical approaches to more mechanism-based models will be discussed.

Population pharmacokinetics of mycophenolic acid and metabolites in patients with glomerulonephritis

Mycophenolic acid (MPA) is an inosine monophosphate dehydrogenase inhibitor used for glomerulonephritis treatment. The objective of the current study was to develop a population pharmacokinetic model for MPA and metabolites in glomerulonephritis to enable appropriate design of MPA regimens in these patients with alterations in kidney structure and function. Thirty-nine patients with glomerulonephritis and receiving mycophenolate mofetil were recruited to participate in a 24-hour pharmacokinetic study. Blood was collected at times 0, 0.5, 1.0, 1.5, 2, 3, 4, 6, 8, 12, and 24 hours and urine was collected over the intervals of 0 to 6, 6 to 12, and 12 to 24 hours. Plasma and urine samples were assayed for MPA and MPA glucuronide (MPAG) by high-performance liquid chromatography and for acyl-MPA glucuronide (AcMPAG) by liquid chromatography/mass spectrometry. Population pharmacokinetic analysis and covariate model building were evaluated using Non-linear Mixed Effect Modeling software (NONMEM, Version 6.2.0; ICON Development Solutions, Ellicott City, MD). The final model for MPA and its metabolites consisted of nine discrete compartments; 1) depot gastrointestinal; 2) central MPA; 3) peripheral MPA; 4) gallbladder; 5) MPA urine; 6) MPAG central; 7) MPAG urine; 8) AcMPAG central; and 9) AcMPAG urine compartment. The MPA population mean estimates for apparent nonrenal clearance (ClNR/F) and apparent central volume of distribution were 14.3 L/hr and 21.1 L, respectively. The mean population estimate for apparent renal clearance (ClR/F) was dependent on estimated creatinine clearances (eClcr); 0.0975 L/hr for eClcr 80 mL/min or less and 0.157 L/hr for eClcr greater than 80 mL/min. Covariate analyses identified: eClcr on CLNR,MPA/F (P < 0.001), eClcr (with a cutoff value at 80 mL/min) on CLR,MPA/F (P < 0.025), serum albumin on CLNR,MPA/F (P < 0.01), eClcr on CLR,MPAG/F (P < 0.001), and eClcr on CLR,AcMPAG/F (P < 0.001). Evaluation of the final model by visual predictive check showed that most of the observed values were within the 95th percent prediction interval generated from 100 simulations of the final model. The current population pharmacokinetic model demonstrated eClcr and serum albumin influenced the renal and nonrenal components of Cl/F, suggesting patients with glomerulonephritis would have highly altered MPA exposures.