Population pharmacokinetics and limited sampling strategy for therapeutic drug monitoring of mycophenolate mofetil in Japanese patients with lupus nephritis
Journal of Pharmaceutical Health Care and Sciences volume 9, Article number: 1 (2023)
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.
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.
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.
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.
Mycophenolate mofetil (MMF) is a prodrug of the immunosuppressive agent mycophenolic acid (MPA). MMF is widely used in kidney and other organ transplantations and is also recommended as induction and maintenance therapy for lupus nephritis (LN), a renal manifestation of systemic lupus erythematosus . Orally administered MMF is almost completely absorbed by the gastrointestinal tract and rapidly hydrolyzed to its active form MPA, which inhibits inosine 5'-monophosphate dehydrogenase (IMPDH), an enzyme involved in the de novo synthesis of guanosine in lymphocytes . MPA is then primarily metabolized into the pharmacologically inactive mycophenolic acid glucuronide (MPAG) in the liver  and thereafter hydrolyzed back to MPA during enterohepatic recirculation (EHC), resulting in a second peak in MPA concentration .
In administering MMF, the pharmacokinetics of MPA should be taken into consideration. MPA has a narrow therapeutic window and exhibits wide inter- and intraindividual pharmacokinetic variability. The pharmacokinetics are affected by numerous factors, including renal function, liver function, serum albumin levels, concomitant drugs, and ethnic background [4,5,6,7,8,9,10,11]. Dosage reduction or drug withdrawal is sometimes necessary when adverse effects arise, such as infections or gastrointestinal symptoms [12,13,14,15]. Although a few studies have reported correlations between pharmacokinetic parameters and adverse effects [16,17,18], much remains unknown. Because of the administrative difficulties, adjusted dosage is desired rather than fixed dosage.
The 12-h area under the concentration–time curve (AUC0-12) for MPA reportedly correlates well with clinical outcomes [19, 20] and is considered to represent a valuable tool for adjusting MMF dosages. Although trough concentration could provide an alternative to AUC in general, such single-point measurements would not be appropriate in the case of MMF due to the second peak in MPA concentration [4, 21]. In patients receiving renal transplantation, a target range of 30–60 μg∙h/mL in AUC0-12 has been proposed in general, and therapeutic drug monitoring (TDM) using AUC0-12 has been suggested as a strategy for personalized treatment . One clinical challenge with AUC0-12 is that calculation requires multiple blood samples. In MMF, 8–10 blood samples within 12 h would be suggested for ideal calculation, but in reality, blood collection at such a high frequency is not only costly and time-consuming, but also places a huge burden on patients. As a solution, a limited sampling strategy (LSS), which estimates AUC0-12 from a small number of samples, may be developed. For transplant patients receiving MMF, LSSs, varying from patient group to patient group, are commonly applied in TDM using AUC0-12 .
For Asian LN patients treated with MMF, however, the pharmacokinetic properties of MPA remain largely unexplored and few, if any, population pharmacokinetic (PPK) models, which may help in the development of LSS, appear to have been built. PPK models have been reported for transplant patients with various backgrounds [24,25,26], and a few have been described for autoimmune patients of different ethnicities [10, 27]; but these may not be applicable to Asian LN patients because of potential differences in MPA pharmacokinetics. While previous studies have reported that a range of 30–45 μg∙h/mL in AUC0-12 is associated with good clinical outcomes , an LSS appears less likely to be available for Asian LN patients.
This study aimed to develop an LSS for Japanese LN patients receiving MMF, in the hope that the LSS would help facilitate safer, more effective, and individualized therapy for such patients. For this purpose, MPA pharmacokinetics after oral administration of MMF were investigated first, and then a PPK model was constructed to describe the pharmacokinetics. From the resulting model, an LSS was developed to estimate AUC0-12 from single or double blood sampling. For comparison, two distinct approaches were employed in developing an LSS: maximum a posteriori Bayesian estimation (MAP-BE) and multiple linear regression (MLR).
Patients and data collection
Participants comprised all Japanese patients ≥ 18 years old who received MMF (CellCept®; Chugai Pharmaceutical Group Co., Tokyo, Japan) for treatment of LN at Aichi Medical University Hospital (Nagakute, Japan) between March 2015 and June 2022 and for whom data were available on multiple blood samples necessary to make a pharmacokinetic profile.
From the electronic medical records, the following demographic and clinical data of patients were retrospectively collected: sex; age; body weight; laboratory data, including serum albumin, serum creatinine, urine protein/creatinine ratio (UPC), alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bilirubin (TBIL), and C-reactive protein (CRP); dose of MMF; dose of concomitant prednisolone; other concomitant medications, including tacrolimus, proton pump inhibitors (PPIs), iron/magnesium oxide, and non-steroidal anti-inflammatory drugs (NSAIDs); and adverse effects. Estimated glomerular filtration rate (eGFR) was calculated from serum creatinine levels using a new 3-variable Japanese equation , and creatinine clearance (CLcr) was calculated using the Cockcroft-Gault equation . As adverse effects of MMF, gastrointestinal involvement, infection, leukopenia, and alopecia were reported. Of these, gastrointestinal symptoms and infections were examined as major adverse effects [12,13,14,15].
MMF was administered twice a day, every 12 h. In a steady state after multiple administration, blood samples were taken at 9 time points during the dosing interval: right before administration (C0) and 0.5, 1, 2, 3, 4, 6, 8, and 12 h after administration (C0.5, C1, C2, C3, C4, C6, C8, and C12, respectively). MPA concentrations were measured using enzyme immunoassay on cobas® 6000 c501 (Roche Diagnostics K.K.; Tokyo, Japan) with a commercially available cobas® MPA Kit (Roche Diagnostics K.K.; Tokyo, Japan). The lower limit of quantification was 0.40 μg/mL.
The parameters used for pharmacokinetic analysis were maximum MPA concentration within a dosing interval of 12 h (Cmax), time to reach maximum MPA concentration (Tmax), and AUC0-12. Cmax and Tmax were derived from the MPA concentration-versus-time profiles of patients. AUC0-12 was calculated from MPA concentrations at the 9 time points of C0 to C12, using the linear trapezoidal method. Dose-normalized Cmax and dose-normalized AUC0-12 were also calculated. Correlations were investigated between AUC0-12 and MMF dose, as well as between AUC0-12 and C0 MPA concentration, using the Spearman rank correlation test, where AUC0-12 was regarded as a dependent variable and the dose or concentration as an independent variable.
From a clinical point of view, we performed group comparisons to investigate the effect of biochemical factors on the pharmacokinetic parameters of MPA, as the previous study  did. The patients were divided into two groups by eGFR (< 81 mL/min versus ≥ 81 mL/min), and also by serum albumin level (< 3.5 g/dL versus ≥ 3.5 g/dL). Similarly, for the effect of concomitant drugs, patients were split into two groups by the presence or absence of each of the following drug administrations: tacrolimus, PPI, and iron/magnesium oxide. In all group comparisons, differences in Tmax, dose-normalized Cmax, and dose-normalized AUC0-12 were examined using the Mann–Whitney U test.
All statistical analyses of pharmacokinetics were performed using easy R (EZR) version 1.54 .
PPK model development and validation
All analyses for PPK model development and validation, as well as subsequent LSS development and validation, were performed using Phoenix® NLME™ software (version 8.1; Pharsight, Mountain View, CA, USA) and EZR version 1.54 . For the estimation algorithm, the first-order conditional estimation-extended least-squares estimation method was employed. For the overall procedure, Bonate , Owen and Fiedler-Kelly , Gabrielsson and Weiner , and Guidance for Industry: Population Pharmacokinetics  were referred to.
The PPK model was developed in a step-by-step manner. First, in the search of a base model, one-, two- and three-compartment models with first-order elimination were tested. These models were parameterized in terms of absorption rate constant (Ka), central volume of distribution (V1), central compartment clearance (CL), peripheral volume of distribution (V2), and inter-compartmental clearance (Q). Because absolute bioavailability (F) was not to be assessed in cases of oral administration, V1/F, CL/F, V2/F, and Q/F were considered to correspond to V1, CL, V2, Q, respectively. Using the resulting base model, the intra-individual variability of MPA concentration was examined by comparing additive, proportional, and mixed (additive plus proportional) residual error models. Inter-individual variability was examined using the exponential error model.
As the next step, covariate models were developed and the following variables were tested as covariates: age, sex, body weight, eGFR, CLcr, serum albumin, UPC, ALT, AST, TBIL, CRP, dose of concomitant prednisolone, and concomitant drugs of PPI and iron/magnesium oxide. In addition, an EHC model was also examined in reference to a previous study . Finally, whether the MPA absorption process was well described was examined by comparing models with and without lag time. For the same purpose, transit models with different numbers of compartments (1 to 7) were also compared.
In these steps of model development, model selection decisions were made based primarily on the likelihood ratio test. In the test, objective function values (OFVs) in nested models were compared using the χ2 test. If the difference between models resulted in OFV greater than the critical value (i.e., 6.63 units when models differ by 1 degree of freedom), the more complex model was considered to be significantly better (p < 0.01). Subsequently, relative standard errors (RSE) were calculated for the evaluation of adequacy, and the condition number (i.e. the ratio of the absolute highest and lowest eigenvalues) for overparameterization. In both calculations, a smaller value was considered to be better. As another way to see adequacy, goodness-of-fit plots were performed. Goodness-of-fit plots were examined for observed concentrations versus individual predicted concentrations, observed concentrations versus population-predicted concentrations, conditional weighted residuals versus time since the dose, and conditional weighted residuals versus population-predicted concentrations [37, 38].
For validation, models in the development process were examined using the prediction-corrected visual predictive check (pcVPC, 1000 simulations) . The precision of parameter estimates was assessed using the bootstrap method (1000 bootstrap samples).
LSS development and validation
From a clinically practical point of view, we developed single and double blood sampling strategies. In the development, we employed two distinct methods: an MAP-BE approach using the PPK model and an MLR approach . In the MAP-BE approach, from one or two observed MPA concentrations, the remaining values from C0 to C12 were estimated based on the newly developed PPK model. From the observed and estimated MPA concentrations, the predicted value of AUC0-12 (AUCpredicted) was calculated. In the MLR approach, multiple linear regression analysis was performed between one or two observed values of MPA concentration (independent variables) and AUCobserved (dependent variable), which was calculated in the same way as AUC0-12. Using the resulting MLR equation, AUCpredicted was calculated from the one or two observed values.
For validation, adjusted coefficients of determination (r2) were used to examine the regression level of AUCpredicted and AUCobserved. Accuracy and precision of prediction were evaluated at each time point from C0 to C12, using mean prediction error (MPE), 95% confidence interval (CI) of MPE, and root mean squared percentage error (RMSPE). Acceptance ranges were set as ≤ ± 10% in MPE and ≤ 25% in RMSPE ; CI of MPE was considered to be acceptable when it crossed the value of 0. In the equations below, N represents the number of pairs of AUCpredicted and AUCobserved.
As major adverse effects, gastrointestinal symptoms and infections were examined. Using the Mann–Whitney U test, patients with and without gastrointestinal involvement were compared in terms of dose, dose/body weight, C0, Cmax, AUC0-12, AUC0-12/dose, and AUC0-0.5; and so were patients with and without infection. If the correlation between a parameter and either of the major adverse occurrences was suggested, receiver operating characteristic (ROC) curve analysis was performed for the parameter, to determine the optimal cutoff. Statistical analyses were performed using EZR version 1.54 .
Thirty-four Japanese ≥ 18 years old who received MMF for treating LN were enrolled. They received MMF at doses of 250–1,000 mg every 12 h twice a day, as prescribed in the package insert. Demographic and clinical characteristics of the patients are shown in Table 1.
In 33 patients, data for blood samples taken at all 9 time points (C0 to C12) were available, making full pharmacokinetic profiles. In 1 patient, data for only 5 time points (C0, C1, C2, C3, and C4) were obtained. The data set contained 302 MPA plasma concentrations obtained from 34 patients, which were used for PPK model construction. For pharmacokinetic analysis and LSS development, full profiles of the 33 patients were used.
For the 34 patients, median (first quartile—third quartile) C0 MPA concentration was 1.8 (1.1–2.9) μg/mL; median Cmax was 16.7 (10.9–26.9) μg/mL; median Tmax was 1.0 (1.0–1.0) h; and median AUC0-12 was 51.1 (41.1–72.7) μg∙h/mL. In 32 of the 33 patients with a full profile, the second peak was observed: the median second peak of MPA concentration was 3.3 (2.3–5.6) μg/mL, and the median period from administration of MMF to the second peak was 8 (6–9) h (Fig. 1A). AUC0-12 did not correlate with MMF dose (r2 = 0.14) (Additional file 1A), but did correlate with C0 MPA concentration (r2 = 0.79) (Additional file 1B).
The effects of eGFR and serum albumin on MPA pharmacokinetics were observed in dose-normalized AUC0-12: patients with eGFR < 81 mL/min (n = 18) showed significantly higher dose-normalized AUC0-12 than those with eGFR ≥ 81 (mL/min) (n = 15) (Fig. 2A), and patients with serum albumin < 3.5 g/dL (n = 12) showed significantly lower dose-normalized AUC0-12 than those with serum albumin ≥ 3.5 g/dL (n = 21) (Fig. 2B). On dose-normalized Cmax, serum albumin had significant effects: patients with serum albumin < 3.5 g/dL displayed lower dose-normalized Cmax than those with serum albumin ≥ 3.5 g/dL (Fig. 2E F). Neither eGFR nor serum albumin had significant effects on Tmax.
The effects of concomitant administration of PPI nor iron/magnesium oxide were not observed in dose-normalized AUC0-12 (Fig. 2C D), but in dose-normalized Cmax. It was significantly lower in patients with concomitant PPI (n = 21) than in those without it (n = 12) (Fig. 2G); and likewise in patients with iron/magnesium oxide (n = 5), compared to those without it (n = 28) (Fig. 2H). On Tmax, concomitant iron/magnesium oxide showed an effect: patients receiving iron/magnesium oxide (n = 5) displayed a significantly increased Tmax compared to those without such treatment (n = 28).
PPK model development and validation
The structure of the final PPK model that adequately described MPA pharmacokinetics of the 34 patients is given in Fig. 3. As a base model, a two-compartment linear model with first-order elimination was employed. Intra-individual variabilities were best described by a proportional residual error model. Then models with different covariant patterns were developed. Of these models, the one including all the variables significant in the pharmacokinetic analysis, i.e. CLcr and serum albumin for CL and concomitant PPI and iron/magnesium oxide for V1, was best in OFV. The model with the four variables, however, was suspected for overparameterization, and therefore was compared to the other candidate model, which included CLcr and serum albumin alone as covariates. The comparison found that, although the model with the four covariates had a significantly better OFV, the model with the two covariates had a substantially decreased value in other areas. RSE (%) for V1 decreased from 32.71 in the model with the four variates to 17.71 in the model with the two covariates, RSE (%) for IIV V1 from 52.92 to 28.18, and the condition number from 5128.3 to 1848.7, indicating overall improvement. In addition, the confidence interval of 1000 Bootstrap for V1 was further examined. While it was 32.41–168.41 in the model with the four covariates, it was 13.43–32.75, decreased by about 85%, in the model with the two covariates, showing that the reliability of the estimates was much improved by excluding PPI and iron/magnesium oxide from covariates (Table 2 and Additional file 2). The model with CLcr and serum albumin alone as covariates was thus adopted for the final model. In the next step, transit compartments were added one at a time, up to seven, and each model was examined. The model with lag time and the EHC model were also examined. The model with 6 compartments significantly improved the model fit and better described the absorption process. The model with lag time or the EHC model did not improve the model fit.
All goodness-of-fit plots showed the high predictive performance of the final model, with no systematic deviations observed (Fig. 1B-E). The relative errors were comparable across all concentration ranges (Fig. 1E, Additional file 3). The prediction-corrected visual predictive check is shown in Fig. 1F. The observed median and 5th and 95th percentiles were located inside the 90% CI of the simulated data. Parameter estimates by the bootstrap method closely matched the means of corresponding parameter estimates from the final model, confirming the stability of the final model (Table 2).
LSS development and validation
Table 3 shows the best four predictive performances for each of single and double blood sampling strategies under the MAP-BE approach (using Bayesian estimation), and the MLR approach (using multiple linear regression analysis). For both approaches, double sampling strategies were more accurate than single ones. In the MAP-BE approach, the MPE and RMSPE values were all within the clinically acceptable ranges of ≤ 10% and ≤ 25%, respectively; and the CIs of MPE crossed the value of 0, except for single sampling at C0, single sampling at C4, double sampling at C1 and C3, and double sampling at C1 and C4. In the MLR approach, the RMSPE values of the single sampling strategies were all larger than the acceptable range, while those of the double sampling strategies were all within the acceptable range. The MPE values were all within the acceptable range and the CIs of MPE all crossed the value of 0.
The best agreement with AUCobserved was achieved by the double blood sampling strategy of C1 and C4 by the MLR approach, where r2 was 0.885, MPE was 2.50% (95%CI, -2.74–7.73), and RMSPE was 14.75%. Of the single sampling strategies, C6 sampling strategy by the MAP-BE approach displayed the best agreement with the AUC0-12, where r2 was 0.769, MPE was -0.72% (95%CI, -7.51–6.08), and RMSPE was 18.89%. The estimation accuracy was found to be good at near the target value of AUC0-12 (30–45 μg∙h/mL) (Fig. 4).
The results for major adverse occurrences (gastrointestinal involvement and infection) are shown in Table 4. Gastrointestinal involvement occurred in 10 of the 33 patients, and diarrhea, vomiting, and epigastric pain were observed. Infections occurred in 14 patients, and upper respiratory infection, fever, pharyngitis, bronchitis, oral herpes zoster, Malassezia, and Aspergillus were observed. No case of cytomegalovirus infection was observed. For gastrointestinal involvement, correlations were suggested for three parameters: AUC0-0.5, Cmax and AUC0-12 were significantly higher in patients with gastrointestinal involvement (n = 10) than those without gastrointestinal involvement (n = 23) (p < 0.01, p < 0.05 and p < 0.05, respectively). No correlation was suggested between infection and any of the examined parameters.
ROC curve analysis was performed for each of the three parameters for which correlation with gastrointestinal involvement was suggested. The analysis indicated the optimal cutoff for each parameter: 3.880 μg∙h/mL for AUC0-0.5, 24.360 μg/mL for Cmax, and 47.597 μg∙h/mL for AUC0-12 (Additional file 4A-C).
This study retrospectively investigated the electric medical records of 34 Japanese LN patients receiving MMF. MPA pharmacokinetic analysis revealed that renal function, serum albumin, concomitant PPI, and iron/magnesium oxide had significant effects (Fig. 2). As a PPK model, a two-compartment model with first-order elimination was developed in which the absorption process was best described by fixed 6 transit compartments (Fig. 3), and intra-individual variabilities were best described by a proportional residual error model. In the model, CLcr and serum albumin were significant covariates for CL (Table 2). For an LSS to be developed, single and double blood sampling strategies were examined by MAP-BE and MLR approaches. The results showed that the double sampling strategy of C1 and C4 by MLR best agreed with AUCobserved (r2 = 0.885). Of the single sampling strategies, the one at C6 by MAP-BE agreed best (r2 = 0.769) (Table 3).
MPA pharmacokinetics are reportedly affected by renal function. Reduced renal function can affect the MPA pharmacokinetics, leading to reduced elimination of and thereby higher exposure to MPA . Because around 97% of MPA binds to albumin, hypoalbuminemia increases the free fraction of MPA, resulting in faster clearance of and thereby reduced exposure to MPA. In fact, low serum albumin concentrations have been reported to increase CL [4, 6]. In this study too, covariate analysis identified the effects of CLcr and serum albumin on CL. Although MPA pharmacokinetics could differ among populations , the estimated CL in our study was 13.15 L/h, comparable with those in the previous studies of Thai (14.5 L/h)  and Mexican (15.4 L/h)  patients with LN. Taken together, the MPA pharmacokinetics of the population in this study were similar in many ways to those in the previous studies: dose-normalized AUC0-12 increased significantly when renal function was damaged, but decreased significantly when serum albumin concentrations were low [40, 42]; and dose-normalized Cmax decreased significantly with concomitant PPI and iron/magnesium oxide (Fig. 2) [7, 43, 44]. These effects were taken into consideration in the first PPK model, in which creatinine clearance and serum albumin were covariates for MPA clearance and concomitant PPI and iron/magnesium oxide were covariates for V1 (Additional files 2, 5).
Absorption of MPA increases at low intragastric pH as the dissolubility of MMF increases [7, 45]. Although PPI was found to be a possible covariate of V1, histamine H2-receptor antagonist, which restrains gastric acid secretion, was not. A plausible speculation is that concurrent use of histamine H2-receptor antagonists does not affect pharmacokinetics overly much because these agents wear off as time passes, whereas PPIs impede irreversibly and provide lasting effects . Besides PPI, iron/magnesium oxide decrease absorption of MPA as it forms a chelate with MMF [5, 47]. These facts may explain why concomitant PPI and iron/magnesium oxide were identified as possible covariates on V1.
PPI and iron/magnesium oxide, which are not always used for treating LN, were coadministered in 22 (64%) and 5 (14.7%) out of 34 patients, respectively, in this study. They were found to have significant effects in the pharmacokinetic analysis, and the PPK model including them, in addition to CLcr and serum albumin, as covariates significantly decreased OFV. Yet, suspecting overparameterization and considering evaluation adequacy and precision, we decided to exclude them from covariates for the final model. For confirmation of the decision, we examined LSSs from Bayesian estimation using the PPK model adding covariates of PPI and iron/magnesium oxide to the final model. As with the final model, best were the LSSs of C6 in single sampling and of C6 and C8 in double sampling (Additional file 6 and Table 3). Taking into consideration the same results in Bayesian estimation and the fact that they are not always necessary for treating LN, we found it appropriate the decision to exclude concomitant PPI and iron /magnesium oxide from covariates.
The PPK model with 6 transit compartments in this study offered a good description of the absorption process. The model did not, however, improve with the addition of EHC, which causes the second peak in MPA concentration. In a previous study, EHC was reported to account for as much as 40% of MPA AUC0-12 . In this study, the second peak was seen in all but one patient, although peak values varied (Fig. 1A). Nevertheless, the EHC model did not significantly improve model fit. The failure of the EHC model may be because of the small number of patients, and also because of the lack of data on MPAG. MPA metabolizes into MPAG in the liver, and MPAG is thus involved in EHC . In this study, however, MPAG concentrations were unknown and thus not taken into consideration. If the EHC model offered improved model fit, the prediction accuracy by MPA-BE would have been better.
As we mentioned earlier, AUC0-12 of MPA is reported to be well associated with clinical outcomes [19, 20] and is considered valuable for adjusting MMF dosage. As AUC0-12 requires frequent blood sampling, possible surrogate markers such as MMF dose and trough concentration have been examined . In this study, MMF dose did not correlate with MPA AUC0-12 (r2 = 0.14). Whereas MPA concentration at C0 correlated with MPA AUC0-12 (r2 = 0.79) (Additional file 1), RMSPE was beyond the clinically acceptable range of 25% in evaluating the prediction accuracy of C0 single blood sampling strategy by MLR. These results suggest that either MMF dose or C0 MPA concentration may not work as surrogate markers for MPA AUC0-12. Further, no single sampling strategy by MLR showed RMSPE < 25%, suggesting that MPA concentration at any single time point is less likely usable as a surrogate marker. On the other hand, Bayesian estimation using the PPK model may allow a single blood sampling to predict AUC0-12 with sufficient accuracy. In this study, the best 2 strategies of C6 and C8 by MAP-BE were good enough for clinical use (Table 3).
When double blood sampling is available, the prediction accuracy of LSS naturally increases. The best 4 strategies by both MAP-BE and MLR, except for double sampling at C1 and C3, and double sampling at C1 and C4 by MAP-BE, were good enough for clinical use in this study. Although the best was the double sampling strategy of C1 and C4 by MLR, an LSS from MAP-BE may be preferable because of a few benefits. The MAP-BE approach might better estimate AUC0-0.5 and Cmax, which this study demonstrated were associated with gastrointestinal involvement (Table 4), especially when only one or two blood samples are available. An LSS by MAP-BE would thus have a better potential of facilitating dose adjustment in consideration of both effectiveness and tolerability. In addition, estimation by MAP-BE can be flexible in terms of blood sampling timing, making LSSs by MAP-BE more practical from the clinical perspective . The accuracy of MAP-BE estimation would increase if the study size is larger and the PPK model improves.
A plausible application of this study to clinical practice would be to determine the initial dosing using the PPK model and then adjust the dosage based on the Bayesian estimation using the measured values so that AUC0-12 would be within the target range of 30–45 μg∙h/mL. AUC0-0.5 < 6.880 μg∙h/mL and Cmax < 24.360 μg/mL would be desirable for minimizing gastrointestinal symptoms.
This study showed some limitations. First, this study did not fully examine all factors that may affect MPA pharmacokinetics. Prednisolone was administered in all patients and NSAIDs in only three patients: for these medications, comparisons were therefore impossible. Data on diet, which can reportedly affect MPA pharmacokinetics , were unavailable, as were data on pharmacogenetics. Secondly, this study did not take free MPA concentration into consideration, because the data were not available. In that free drug is therapeutically active, it would be desirable to evaluate free concentration in MPA. Measuring free drug concentration, though it may be still difficult in clinical practice, should be done in a future study. As another limitation, external validation of the PPK model was not conducted due to the small number of patients. Validation from a larger study is needed in the future. Especially, further examination is warranted in a study with a larger number of patients receiving PPI and iron/magnesium oxide. Finally, the study did not investigate the effectiveness of MMF, which should also be a subject for future study.
This study developed LSSs for Japanese LN patients receiving MMF so that the optimal MMF dosage could be individually determined based on estimated AUC0-12. When blood samples are taken twice, at 1 and 4 h after MMF administration, a strategy using the MLR approach would work best. When only a single blood sampling is allowed, however, a strategy using the MAP-BE approach with the PPK model developed in this study would be best with blood sampling at 6 h after administration. Although further studies are warranted, the LSS, whether single or double sampling, would facilitate safer, more effective, and individualized therapies for Japanese LN patients.
Availability of data and materials
Area under the concentration–time curve
- AUC0-12 :
12-Hour area under the concentration–time curve
Central compartment clearance
Estimated glomerular filtration rate
- FOCE ELS:
First order conditional estimation-extended least squares
Inosine 5’-monophosphate dehydrogenase
Absorption rate constant
Rate constant of transit compartment
Lower limit of the 95% confidence interval
Limited sampling strategies
Maximum a posteriori Bayesian estimation
Multiple linear regression
Mycophenolic acid glucuronide
Mean prediction error
Mean transit time
Non-steroidal anti-inflammatory drug
Objective function value
Prediction-corrected visual predictive check
Proton pump inhibitor
- PPK model:
Population pharmacokinetic model
Root mean squared percentage error
Receiver operating characteristic
Relative standard error
Systemic lupus erythematosus
Therapeutic drug monitoring
Upper limit of the 95% confidence interval
Urine protein/creatinine ratio
Central volume of distribution
Peripheral volume of distribution
Fanouriakis A, Kostopoulou M, Cheema K, Anders HJ, Aringer M, Bajema I, et al. 2019 Update of the Joint European League Against Rheumatism and European Renal Association-European Dialysis and Transplant Association (EULAR/ERA-EDTA) recommendations for the management of lupus nephritis. Ann Rheum Dis. 2020;79:713–23.
Allison AC, Eugui EM. Mycophenolate mofetil and its mechanisms of action. Immunopharmacology. 2020;47:85–118.
Shipkova M, Armstrong VW, Wieland E, Niedmann PD, Schutz E, Brenner-Weiss G, et al. Identification of glucoside and carboxyl-linked glucuronide conjugates of mycophenolic acid in plasma of transplant recipients treated with mycophenolate mofetil. Br J Pharmacol. 1999;12:1075–82.
Bullingham RE, Nicholls AJ, Kamm BR. Clinical pharmacokinetics of mycophenolate mofetil. Clin Pharmacokinet. 1998;34:429–55.
Mino Y, Naito T, Shimoyama K, Ogawa N, Kawakami J. Pharmacokinetic variability of mycophenolic acid and its glucuronide in systemic lupus erythematosus patients in remission maintenance phase. Biol Pharm Bull. 2011;34:755–9.
de Winter BC, van Gelder T, Sombogaard F, Shaw LM, van Hest RM, Mathot RA. Limited sampling strategies for therapeutic drug monitoring of mycophenolate mofetil therapy in patients with autoimmune disease. J Pharmacokinet Pharmacodyn. 2009;36:541–64.
Schaier M, Scholl C, Scharpf D, Hug F, Bonisch-Schmidt S, Dikow R, et al. Proton pump inhibitors interfere with the immunosuppressive potency of mycophenolate mofetil. Rheumatology. 2010;49:2061–7.
Vietri M, Pietrabissa A, Mosca F, Pacifici GM. Mycophenolic acid glucuronidation and its inhibition by non-steroidal anti-inflammatory drugs in human liver and kidney. Eur J Clin Pharmacol. 2000;56:659–64.
Cattaneo D, Perico N, Gaspari F, Gotti E, Remuzzi G. Glucocorticoids interfere with mycophenolate mofetil bioavailability in kidney transplantation. Kidney Int. 2002;62:1060–7.
Romano-Aguilar M, Resendiz-Galvan JE, Medellin-Garibay SE, Milan-Segovia RDC, Martinez-Martinez MU, Abud-Mendoza C, et al. Population pharmacokinetics of mycophenolic acid in Mexican patients with lupus nephritis. Lupus. 2020;29:1067–77.
Li P, Shuker N, Hesselink DA, van Schaik RH, Zhang X, van Gelder T. Do Asian renal transplant patients need another mycophenolate mofetil dose compared with Caucasian or African American patients? Transpl Int. 2014;27:994–1004.
Ginzler EM, Dooley MA, Aranow C, Kim MY, Buyon J, Merrill JT, et al. Mycophenolate mofetil or intravenous cyclophosphamide for lupus nephritis. N Engl J Med. 2005;353:2219–28.
Moore RA, Derry S. Systematic review and meta-analysis of randomised trials and cohort studies of mycophenolate mofetil in lupus nephritis. Arthritis Res Ther. 2006;8(R182):1–10.
Abe Y, Tada K, Yamaji K, Tamura N. Real-World Experience of Safety of Mycophenolate Mofetil in 119 Japanese Patients with Systemic Lupus Erythematosus: A Retrospective Single-Center Study. Biomed Res Int. 2021;8630596:1–7.
Sarangi SC, Reeta KH, Agarwal SK, Kaleekal T, Guleria S, Gupta YK. A pilot study on area under curve of mycophenolic acid as a guide for its optimal use in renal transplant recipients. Indian J Med Res. 2012;135:84–91.
Lertdumrongluk P, Somparn P, Kittanamongkolchai W, Traitanon O, Vadcharavivad S, Avihingsanon Y. Pharmacokinetics of mycophenolic acid in severe lupus nephritis. Kidney Int. 2010;78:389–95.
Zegarska J, Hryniewiecka E, Zochowska D, Tszyrsznic W, Jazwiec R, Borowiec A, et al. Mycophenolic Acid Metabolites Acyl-Glucuronide and Glucoside Affect the Occurrence of Infectious Complications and Bone Marrow Dysfunction in Liver Transplant Recipients. Ann Transplant. 2015;20:483–92.
Sugioka N, Sasaki T, Kokuhu T, Ito Y, Shibata N, Okamoto M, et al. Clinical pharmacokinetics of mycophenolate mofetil in Japanese renal transplant recipients: A retrospective cohort study in a single center. Biol Pharm Bull. 2006;29:2099–105.
Tett SE, Saint-Marcoux F, Staatz CE, Brunet M, Vinks AA, Miura M, et al. Mycophenolate, clinical pharmacokinetics, formulations, and methods for assessing drug exposure. Transplant Rev. 2011;25:47–57.
Kuypers DR, Le Meur Y, Cantarovich M, Tredger MJ, Tett SE, Cattaneo D, et al. Transplantation Society (TTS) Consensus Group on TDM of MPA. Consensus report on therapeutic drug monitoring of mycophenolic acid in solid organ transplantation. Clin J Am Soc Nephrol. 2010;5:341–58.
Barraclough KA, Isbel NM, Staatz CE. Evaluation of the mycophenolic acid exposure estimation methods used in the APOMYGERE, FDCC, and Opticept trials. Transplantation. 2010;90:44–51.
van Gelder T, Silva HT, de Fijter JW, Budde K, Kuypers DR, Tyden G, et al. Comparing mycophenolate mofetil regimens for de novo renal transplant recipients: the fixed-dose concentration-controlled trial. Transplantation. 2008;86:1043–51.
Sobiak J, Resztak M. A Systematic Review of Multiple Linear Regression-Based Limited Sampling Strategies for Mycophenolic Acid Area Under the Concentration-Time Curve Estimation. Eur J Drug Metab Pharmacokinet. 2021;46:721–42.
Shum B, Duffull SB, Taylor PJ, Tett SE. Population pharmacokinetic analysis of mycophenolic acid in renal transplant recipients following oral administration of mycophenolate mofetil. Br J Clin Pharmacol. 2003;56:188–97.
Colom H, Andreu F, van Gelder T, Hesselink DA, de Winter BCM, Bestard O, et al. Prediction of Free from Total Mycophenolic Acid Concentrations in Stable Renal Transplant Patients: A Population-Based Approach. Clin Pharmacokinet. 2018;57:877–93.
Sheng C, Zhao Q, Niu W, Qiu X, Zhang M, Jiao Z. Effect of Protein Binding on Exposure of Unbound and Total Mycophenolic Acid: A Population Pharmacokinetic Analysis in Chinese Adult Kidney Transplant Recipients. Front Pharmacol. 2020;11:340.
de Winter BC, Neumann I, van Hest RM, van Gelder T, Mathot RA. Limited sampling strategies for therapeutic drug monitoring of mycophenolate mofetil therapy in patients with autoimmune disease. Ther Drug Monit. 2009;31:382–90.
Luszczynska P, Pawinski T. Therapeutic Drug Monitoring of Mycophenolic Acid in Lupus Nephritis: A Review of Current Literature. Ther Drug Monit. 2015;37:711–7.
Matsuo S, Imai E, Horio M, Yasuda Y, Tomita K, Nitta K, et al. Collaborators developing the Japanese equation for estimated GFR. Revised equations for estimated GFR from serum creatinine in Japan. Am J Kidney Dis. 2009;53:982–92.
Cockcroft DW, Gault MH. Prediction of creatinine clearance from serum creatinine. Nephron. 1976;16:31–41.
Luszczynska P, Pawinski T, Kunicki PK, Durlik M, Augustyniak-Bartosik H, Hurkacz M. Pharmacokinetics of free and total mycophenolic acid in adult lupus nephritis patients-implications for therapeutic drug monitoring. Eur J Clin Pharmacol. 2019;75:371–9.
Kanda Y. Investigation of the freely available easy-to-use software “EZR” for medical statistics. Bone Marrow Transplant. 2013;48:452–8.
Bonate PL. Pharmacokinetic-Pharmacodynamic Modeling and Simulation. 2nd ed. New York: Springer; 2011.
Owen JS, Fiedler-Kelly J. Introduction to Population Pharmacokinetic/Pharmacodynamic Analysis with Nonlinear Mixed Effects Models. New York: John Wiley & Sons, Inc; 2014.
Gabrielsson J, Weiner D. Pharmacokinetic and Pharmacodynamic Data Analysis: Concepts and Applications. 5th ed. Stockholm: Apotekarsocieteten; 2017.
Center for Drug Evaluation and Research, Center for Biologics Evaluation and Research. Guidance for Industry: Population Pharmacokinetics. Rockville, MD: U.S. Department of Health and Human Services, Food and Drug Administration, Center for Drug Evaluation and Research, Center for Biologics Evaluation and Research. 2022. https://www.fda.gov/media/128793/download. Accessed 1 Oct 2022.
Hooker AC, Staatz CE, Karlsson MO. Conditional weighted residuals (CWRES): a model diagnostic for the FOCE method. Pharm Res. 2007;24:2187–97.
Karlsson MO, Savic RM. Diagnosing model diagnostics. Clin Pharmacol Ther. 2007;82:17–20.
Bergstrand M, Hooker AC, Wallin JE, Karlsson MO. Prediction-corrected visual predictive checks for diagnosing nonlinear mixed-effects models. AAPS J. 2011;13:143–51.
Bergan S, Brunet M, Hesselink DA, Johnson-Davis KL, Kunicki PK, Lemaitre F, et al. Personalized Therapy for Mycophenolate: Consensus Report by the International Association of Therapeutic Drug Monitoring and Clinical Toxicology. Ther Drug Monit. 2021;43:150–200.
Punyawudho B, Lertdumrongluk P, Somparn P, Kittanamongkolchai W, Traitanon O, Avihingsanon Y, et al. Population pharmacokinetics of mycophenolate mofetil in Thai lupus nephritis patients. Int J Clin Pharmacol Ther. 2012;50:272–80.
Naesens M, de Loor H, Vanrenterghem Y, Kuypers DR. The impact of renal allograft function on exposure and elimination of mycophenolic acid (MPA) and its metabolite MPA 7-O-glucuronide. Transplantation. 2007;84:362–73.
Miura M, Satoh S, Inoue K, Kagaya H, Saito M, Suzuki T, et al. Influence of lansoprazole and rabeprazole on mycophenolic acid pharmacokinetics one year after renal transplantation. Ther Drug Monit. 2008;30:46–51.
Morii M, Ueno K, Ogawa A, Kato R, Yoshimura H, Wada K, et al. Impairment of mycophenolate mofetil absorption by iron ion. Clin Pharmacol Ther. 2000;68:613–61.
Kees MG, Steinke T, Moritz S, Rupprecht K, Paulus EM, Kees F, et al. Omeprazole impairs the absorption of mycophenolate mofetil but not of enteric-coated mycophenolate sodium in healthy volunteers. J Clin Pharmacol. 2012;52:1265–2127.
Bruno G, Zaccari P, Rocco G, Scalese G, Panetta C, Porowska B, et al. Proton pump inhibitors and dysbiosis: Current knowledge and aspects to be clarified. World J Gastroenterol. 2019;25:2706–19.
Bullingham R, Shah J, Goldblum R, Schiff M. Effects of food and antacid on the pharmacokinetics of single doses of mycophenolate mofetil in rheumatoid arthritis patients. Br J Clin Pharmacol. 1996;41:513–6.
We would like express our gratitude to Drs. Hirokazu Sugiyama, Mayumi Ito, Akimasa Asai, Fumiya Kitamura and Keisuke Kamiya for patients’ care.
Part of this work was supported by a research grant from the Aichi Kidney Foundation (grant number 2020–7 to T.M. and Y.I.).
Ethical approval and consent to participate
This study protocol was approved by the hospital ethics committee at Aichi Medical University (approval no. 2018-H258). Written informed consent was obtained from all patients.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlations of MPA AUC0-12 with Dosage (A) and C0 (B).
Parameter estimates and bootstrap results of the PPK model that included PPI and iron/magnesium oxide in covariates.
Scatter plots on logarithmic scale of observed versus predicted MPA concentrations.
Receiver operating characteristic curve (ROC) analysis of associations between gastrointestinal involvement and pharmacokinetic parameters.
Concentration versus time, goodness-of-fit plots of the PPK final model, and a visual predictive check for the final model.
Best performing single and double sampling strategies for estimation of AUC0-12 of MPA by Bayesian estimation using the PPK model that included PPI and iron/magnesium oxide in covariates.
About this article
Cite this article
Mizaki, T., Nobata, H., Banno, S. et al. Population pharmacokinetics and limited sampling strategy for therapeutic drug monitoring of mycophenolate mofetil in Japanese patients with lupus nephritis. J Pharm Health Care Sci 9, 1 (2023). https://doi.org/10.1186/s40780-022-00271-w
- Mycophenolate mofetil
- Mycophenolic acid
- Lupus nephritis
- Limited sampling strategy
- Population pharmacokinetics