- Research article
- Open Access
Usefulness of limited sampling strategy for mycophenolic acid area under the curve considering postoperative days in living-donor renal transplant recipients with concomitant prolonged-release tacrolimus
Journal of Pharmaceutical Health Care and Sciencesvolume 3, Article number: 17 (2017)
The optimal dose of mycophenolate mofetil (MMF) in renal transplant patients has been recommended to be decided on the basis of area under the concentration-time curve (AUC0-12) of mycophenolic acid (MPA). Although meta-analysis has revealed that postoperative day (POD) is an influencing factor in MPA pharmacokinetics, there are no reports regarding a limited sampling strategy (LSS) for MPA AUC in consideration of POD. The aim of this study was to construct of an LSS considering POD that appropriately expresses the MPA AUC following renal transplantation and evaluation of the usefulness.
Serum concentration–time profiles (measured AUC0-12) comprising nine sampling points over 12 h were analyzed in 36 living-donor renal transplant recipients after MMF administration with concomitant once-daily prolonged-release tacrolimus. Two LSSs were developed by stepwise multiple regression analysis (Method A: not classified by PODs; Method B: classified by PODs into POD < 31 and POD ≥ 31). Each LSS comprised four blood-sampling points within 6 h after MMF administration. Precision and reliability were verified by using root-mean-square error (RMSE), correlation coefficient (R2), and coefficient of determination (q2) by using leave-one-out cross-validation. The absolute values of the difference between measured and estimated AUCs (delta AUC) were compared for both estimating equations.
One-hundred samples obtained from 36 recipients for AUC0-12 comprised POD < 31 (n = 39) and POD ≥ 31 (n = 61). Estimation of AUC0-12 by Method B resulted in better accuracy and reliability (Method A: RMSE = 5.5, R2 = 0.85, q2 = 0.83; Method B: POD < 31: RMSE = 5.5, R2 = 0.86, q2 = 0.83; POD ≥ 31: RMSE = 3.9, R2 = 0.92, q2 = 0.89) and significantly lower median delta AUC compared with that by Method A (delta AUC: 2.6 (0.0–11.6) v.s. 3.9 (0.1–18.1), p = 0.032).
These results suggest that LSS, classified as POD < 31 or POD > 31, would provide more accurate and reliable estimation of MPA AUC0-12 in Japanese living-donor renal transplant patients.
Mycophenolate mofetil (MMF), a prodrug of the immunosuppressant mycophenolic acid (MPA), has been widely used for the prevention of rejection in solid organ transplant patients [1,2,3]. MMF is administered to patients who have undergone renal transplantation at a dosage of 0.5–1.5 g twice daily. After oral administration, MMF is rapidly absorbed and hydrolyzed to MPA , and is then inactivated to MPA glucuronide by UDP-glucuronosyltransferase .
Numerous studies have demonstrated the relationship between area under the concentration-time curve (AUC0-12) of MPA and both risk of rejection [4, 6,7,8,9,10,11,12] and hematologic side effects [11, 13]. A target range of 30–60 mg h/L for the MPA AUC0-12 has been proposed as a guide to MMF dosage in renal transplant patients in these studies. There were large inter- and intra-individual variations in MPA AUC0-12 . However, routine measurement of full MPA AUC0-12 for 12-h dose intervals is cumbersome and cost-prohibitive. Limited sampling strategies (LSSs) have been developed in several countries for estimating MPA AUC0-12 to overcome these difficulties .
van Hest et al.  reported that MPA pharmacokinetics are affected by the patient’s renal function, serum albumin concentration, and dosage of immunosuppressants, including calcineurin inhibitors that are dependent on the post-renal transplant period. Moreover, several studies reported that oral MPA clearance is inversely proportional to postoperative days (POD), achieving gradual stability [7, 11]. Therefore, sampling points to estimate MPA AUC0-12 might vary according to POD. However, there is no report demonstrating LSS design with consideration of POD. The aim of this study was to develop LSS with consideration of POD, and to evaluation of the usefulness of these LSSs in Japanese renal transplant patients.
This study was performed on all 36 patients who underwent living-donor renal transplantation at Mie University Hospital between November 2005 and August 2015. One-hundred serum MPA concentration–time profiles were prospectively obtained between November 2012 and September 2015.
Demographic data including concomitant drug use were obtained by reviewing electronic medical records of the patients. Concomitant drugs that are documented in Lexicomp, integrated with UpToDate (version 2014; Wolters Kluwer Health, Philadelphia, PA, USA) were considered.
Assay of serum MPA concentration
Serum was separated by centrifugation at 1700 × g for 10 min by using serum separation tubes. Serum MPA concentration was determined by using a homogeneous particle-enhanced turbidimetric inhibition immunoassay (PETINIA) technique on a DIMENSION® Xpand Plus Integrated Chemistry System (Siemens Healthcare Diagnostics K.K., Tokyo, Japan).
All patients received a basic immunosuppression regimen of MMF (CellCept; Chugai Pharmaceutical Co., Ltd., Tokyo, Japan), once-daily prolonged-release tacrolimus (Graceptor; Astellas, Tokyo, Japan,), methylprednisolone, and basiliximab (Simulect i.v. injection; Novartis Pharmaceuticals, Tokyo, Japan). In addition, patients with blood type incompatibility received rituximab at a dose of 200 mg on preoperative day 4. The MMF was taken on preoperative day 4 at a fixed starting dose of 0.5 g twice daily, 1 g twice daily from POD 0, 0.75 g twice daily from POD 15, and 0.5 g twice daily from POD 60. The starting dose of tacrolimus was 0.1 mg/kg/day, adjusted based on whole-blood concentration (target concentration: 6.0–8.0 ng/mL). Methylprednisolone was started at a dose of 20 mg/day and administered at a dose of 250 mg during surgery, reduced gradually to the maintenance dose of 4 mg/day from POD 30. Basiliximab was injected intravenously at a dose of 20 mg within 2 h before the operation and on POD 4.
Determination of MPA AUC0-12
Serum MPA concentration was determined just before dosage and at 0.5, 1, 2, 3, 4, 6, 8, and 12 h after administration. The MPA AUC0-12 was calculated using a linear trapezoidal rule. The AUC0-12 of MPA was determined on POD 7, POD 14, POD 21–28, and POD ≥ 31 according to recommendation on the application of therapeutic drug monitoring to MMF therapy in transplantation .
Development of a POD-based LSS
One-hundred serum MPA concentration–time profiles were classified into four groups on the basis of POD (POD 7, POD 14, POD 21–28 and POD ≥ 31). MPA clearance (oral MMF dose per AUC0-12) was compared between these four groups by using multiple comparison test. POD-based LSSs were developed using multiple comparison analysis. Each estimating equation was developed by using stepwise multiple regression analysis, and comprised four blood-sampling points over 6 h after drug administration.
Evaluation of estimating equations
Precision was evaluated by using Spearman’s rank correlation test, root-mean-square error (RMSE), correlation coefficient (R2) by least squares method, and correlation between measured AUC and estimated AUC. Reliability was evaluated by using coefficient of determination (q2), calculated using leave-one-out cross-validation (LOOCV).
The absolute value of the difference (delta AUC) between measured and estimated AUCs was compared between each LSS by using the Wilcoxon matched pair test.
Spearman’s rank correlation test, RMSE, and least squares method were performed using JMP® Ver. 7.0 (SAS Institute, Cary, NC, USA). Wilcoxon matched pair and Kruskal-Wallis tests were performed using GraphPad Prism Ver. 5.01 (GraphPad Software, Inc., San Diego, CA, USA). LOOCV and multiple comparisons by Wilcoxon signed-rank test (Holm’s method) were performed using GNU R Ver. 3.1.0 for windows. A P value < 0.05 was considered significant.
Development of POD-based LSS
The AUCs0-12 of MPA on POD 7 (n = 13), POD 14 (n = 16), POD 21–28 (n = 10), and POD ≥ 31 (n = 61) were determined as described in the materials and methods section. Significant differences in the oral clearance of MPA were found between POD ≥ 31 and POD 7 or POD 14 by multiple comparison analysis (Fig. 1). Moreover, MPA oral clearance tended to be lower in POD ≥ 31 compared with that in POD 21–28. There were no significant differences in MPA oral clearance between POD 7, POD 14, and POD 21–28. On the basis of these results, the LSSs for PODs < 31 or ≥ 31 were defined by using Method B.
Table 1 displays the data on patient characteristics including gender, primary disease, age, body weight, creatinine clearance estimated by Cockcroft and Gault formula, alanine aminotransferase, aspartate aminotransferase, total bilirubin, concomitant drugs, and measured MPA AUC0-12. Figure 2 shows the median concentration–time profile. Median (range) C0, tmax, and Cmax were 2.9 (0.3–7.7) μg/mL, 2.3 (0.5–6.0) h, and 11.3 (3.3–27.3) μg/mL for POD < 31, and 1.9 (0.2–5.7) μg/mL, 2.2 (0.5–6.0) h, and 11.9 (3.4–43.2) μg/mL for POD ≥ 31, respectively. Median concentrations of MPA for POD ≥ 31 during 1 and 4 h after administration were mostly constant similarly as those for POD < 31 (Fig. 2). In contrast, the AUC0-12 of MPA per dose (mg) for POD < 31 was significantly lower than that for POD ≥ 31 (median (range): 0.056 (0.031–0.098) v.s. 0.078 (0.029–0.155), P < 0.0001).
Evaluation of estimating equations
The P value, RMSE, R2, and q2 of each estimating equation are shown in Table 2. The P values of all estimating equations were less than 0.001. The worst precision (RMSE and R2) and reliability (q2) were observed in MPA AUC0-12 estimation of POD < 31 by Method A. Estimation by Method B resulted in better precision (RMSE and R2) and reliability (q2) than that by Method A. The correlation of Method B with estimated and measured AUC was better than that of Method A, as demonstrated by the results where corresponding slope, intercept, and R2 values of Method B were much closer to 1, 0, and 1, respectively (Fig. 3). Wilcoxon matched pair test indicated that delta AUC estimated by Method B was significantly lower than that by Method A (Fig. 4).
In this study, two different approaches to developing LSSs for the estimation of AUC0-12 were evaluated by statistical analyses (Method A: not classified by POD, Method B: classified by POD into < 31 or ≥ 31 based on difference in MPA systemic clearance). Precision (RMSE and R2) and reliability (q2) were compared between Method A and Method B, and Method B provided better estimation of AUC0-12 compared with Method A. Moreover, delta AUC of Method B was lower than that of Method A. Therefore, these results suggested that LSSs considering POD would provide more precise and reliable estimation of MPA AUC0-12.
It has been reported that patients within 1 month post-transplant have lower MPA AUC0-12 than patients between 3 and 6 months post-transplant [7, 11]. Moreover, van Hest et al.  reported that POD was a significant factor affecting the pharmacokinetics of MPA. In our present study, the results of multiple comparisons of MPA oral clearance between four groups classified by POD proved that MPA oral clearance on POD < 31 was higher than that of POD ≥ 31 (Fig. 1). Furthermore, estimated AUC0-12 on POD < 31 did not exhibit better precision and reliability than that on POD ≥ 31 when Method A was applied (RMSE, R2, and q2 for POD < 31 or POD ≥ 31 were 6.7, 0.78, and 0.73 or 4.2, 0.90, and 0.89, respectively (data not shown)). We also compared MPA clearance between POD < 91 and POD ≥ 91. MPA clearance was lower in POD < 91 compared with that in POD ≥ 91 (17.2 [8.6–32.7] v.s. 12.5 [6.4–34.1], P = 0.0009). However, there was no significant difference in MPA clearance between POD31-90 and POD ≥ 91 (POD31-90: 14.3 [8.6–19.8] v.s. POD ≥ 91: 12.5 [6.4–34.1], p = 0.16). Moreover, LSSs classified by PODs into POD < 91 and POD ≥ 91 was not better than Method B (POD < 91: R2 = 0.86, RMSE = 5.86, q2 = 0.82; POD ≥ 91: R2 = 0.91, RMSE = 3.54, q2 = 0.88). Therefore, data from both these previous studies and our study support our opinion that estimation equations for MPA AUC0-12 with built-in consideration of POD should have better precision and reliability. This is because there is a difference in MPA pharmacokinetics between POD < 31 and POD ≥ 31.
This study has some limitations. First, measured serum MPA concentration includes serum acyl-glucuronide metabolite (AcMPAG) concentration because serum MPA concentration was measured by PETINIA, and the antibody used in PETINIA cross-reacts with AcMPAG . Therefore, the universality of LLS developed in the present study may be limited. Second, the contribution of concomitant drugs is not completely overseen by the present study. However, in LSS by Method B, we performed multiple comparison testing of delta AUC between groups classified by concomitant drug usage, proving there is no difference in delta AUC between groups (Additional files 1, 2, and 3). Therefore, LSSs devised in this study seem to be slightly influenced by the concomitant use of drugs such as quinolone, proton pump inhibitors, and valganciclovir. Third, we cannot explain the mechanism by which mycophenolate clearance changed according to POD. Because of small number of patients involved in the present study, some MPA concentration profiles has been obtained from same patients.
Pawinski et al.  reported that AUC0-12 estimation comprising three blood-sampling points (0, 0.5, and 2.0 h) may provide good prediction of MPA AUC0-12 in renal transplant patients receiving concomitant tacrolimus. However, we could not find literature on LSS that satisfied the following conditions: i) concomitant use of tacrolimus as a calcineurin inhibitor, ii) containing data within 1 month after renal transplantation, and iii) determination of MPA using the PETINIA method. We evaluated the LSS developed by Pawinski et al., which had been well analyzed and meets two conditions (concomitant with tacrolimus and containing data within 1 month after renal transplantation). However, using this equation, we could not obtain a good correlation between the estimated and measured MPA AUCs in our study population (y = 0.74x + 8.32, R2 = 0.47, Additional file 4). The reason for this might be explained as follows: i) patients in our study were concomitantly administered once daily prolonged release tacrolimus; ii) our study population included many early post-transplant patients; and iii) the frequency of UGT1A9 (a metabolic enzyme of MPA) variants (UGT1A9*1, UGT1A9*1c, and UGT1A9*3) varies in the Caucasian, African, and Asian populations .
It has been reported that the MPA AUC0-12, measured by PETINIA method, are overestimated in comparison with MPA concentrations measured by high performance liquid chromatograph (HPLC) method . Miura et al. compared LSS on the POD28 and 1 year after transplantation in Japanese kidney transplant patients. In that study, the values of MPA AUC0-12 (mean [standard deviation (SD)]) measured by HPLC method, were 63.9 [28.9] on the POD 28 and 58.1 [24.3] on 1 year after transplantation, respectively [20, 21]. On the other hand, in the present study, the values of MPA AUC0-12 (mean [SD]) were 44.6 [14.4] on the POD 21-28 (n = 10) and 45.2 [10.2] on about 1 year after transplantation (POD 345-401), respectively (data not shown). Although dose of MMF in our study was similar to that in their study, the mean of MPA AUC0-12 in our study population were lower than that by Miura et al.. In our study, patients was administered prolonged-release tacrolimus concomitantly, whereas in the study of Miura et al., tacrolimus administered to patients was not prolonged-release formulation. Although controversy remains about the interaction between MMF and tacrolimus, prolonged-release tacrolimus formulation might less effect on pharmacokinetics of MPA, since the formulation decreases Cmax of tacrolimus. Therefore, the difference in formulation of tacrolimus might be the reason for the difference of the MPA AUC0-12 between study of Miura et al. and our present study.
Yamaguchi et al.  reported the estimation equations for MPA AUC0-12 in the Japanese population. However, the usefulness of the estimation equation developed by Yamaguchi et al. was limited because the correlation between the measured and estimated AUC0-12 at 1 and 3 months after renal transplantation was not good enough and cross validation of the estimation equation was not performed. In this study, it was demonstrated that our estimation equation was better correlation than that of Yamaguchi et al. and we validated its reliability by cross validation.
In the present study, LSS consisting of four timed samples within 6 h after administration provided accurate and reliable estimation of MPA AUC0-12 and was best among the verified estimation equations (Additional file 5). It is known that MPA plasma concentration profile shows a secondary peak at around 6 h after administration because of enterohepatic circulation  and its mean contribution to the overall MPA AUC0-12 is 37% (10–61%) . Therefore, MPA blood concentration at around 6 h after administration should be important for estimating MPA AUC0-12. Measurement of MPA AUC0-12 by using 10–12 timed blood samples is a burden on both the patient and medical staff, requiring laboratory resources, considerable quantities of patients’ blood, and a minimum 12 h stay in hospital. Therefore, LSSs that developed by this study, not only Method B but also Method A, will reduce the burden on both patients and medical staff and measurement costs.
This study suggested that LSSs with consideration of POD provide more accurate and reliable estimations of MPA AUC0-12 in Japanese renal transplant patients receiving concomitant tacrolimus therapy.
Area under the concentration-time curve
Limited sampling strategies
Particle-enhanced turbidimetric inhibition immunoassay
The European Mycophenolate Mofetil Cooperative Study Group. Placebo-controlled study of mycophenolate mofetil combined with cyclosporin and corticosteroids for prevention of acute rejection. Lancet. 1995;345:1321–5.
Kobashigawa J, Miller L, Renlund D, Mentzer R, Alderman E, Bourge R, et al. A randomized active-controlled trial of mycophenolate mofetil in heart transplant recipients. Transplantation. 1998;66:507–15.
Wiesner R, Rabkin J, Klintmalm G, McDiarmid S, Langnas A, Punch J, et al. A randomized double-blind comparative study of mycophenolate mofetil and azathioprine in combination with cyclosporine and corticosteroids in primary liver transplant recipients. Liver Transpl. 2001;7:442–50.
Shaw LM, Nicholls A, Hale M, Armstrong VW, Oellerich M, Yatscoff R, et al. Therapeutic monitoring of mycophenolic acid. A consensus panel report. Clin Biochem. 1998;31:317–22.
Bullingham RE, Nicholls AJ, Kamm BR. Clinical pharmacokinetics of mycophenolate mofetil. Clin Pharmacokinet. 1998;34:429–55.
Takahashi K, Ochiai T, Uchida K, Yasumura T, Ishibashi M, Suzuki S, et al. Pilot study of mycophenolate mofetil (RS-61443) in the prevention of acute rejection following renal transplantation in Japanese patients. Transplant Proc. 1995;27:1421–4.
Hale MD, Nicholls AJ, Bullingham RE, Hené R, Hoitsma A, Squifflet JP, et al. The pharmacokinetic-pharmacodynamic relationship for mycophenolate mofetil in renal transplantation. Clin Pharmacol Ther. 1998;64:672–83.
Nichols AJ. Opportunities for therapeutic monitoring of mycophenolate mofetil dose in renal transplantation suggested by the pharmacokinetic/pharmacodynamic relationship for mycophenolic acid and suppression of rejection. Clin Biochem. 1998;31:329–33.
van Gelder T, Hilbrands LB, Vanrenterghem Y, Weimar W, de Fijter JW, Squifflet JP, et al. A randomized double-blind, multicenter plasma concentration controlled study of the safety and efficacy of oral mycophenolate mofetil for the prevention of acute rejection after kidney transplantation. Transplantation. 1999;68:261–6.
DeNofrio D, Loh E, Kao A, Korecka M, Pickering FW, Craig KA, et al. Mycophenolic acid concentrations are associated with cardiac allograft rejection. J Heart Lung Transplant. 2000;19:1071–6.
Shaw LM, Korecka M, Aradhye S, Grossman R, Bayer L, Innes C, et al. Mycophenolic acid area under the curve values in African American and Caucasian renal transplant patients are comparable. J Clin Pharmacol. 2000;40:624–33.
Shaw LM, Korecka M, DeNofrio D, Brayman KL. Pharmacokinetic, pharmacodynamic, and outcome investigations as the basis for mycophenolic acid therapeutic drug monitoring in renal transplant patients. Clin Biochem. 2001;34:17–22.
Mourad M, Malaise J, Eddour DC, Meyer MD, König J, Schepers R, et al. Correlation of mycophenolic acid pharmacokinetic parameters with side effects in kidney transplant patients treated with mycophenolate mofetil. Clin Chem. 2001;47:88–94.
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 (Orlando). 2011;25:47–57.
van Hest R, Mathot R, Vulto A, Weimar W, van Gelder T. Predicting the usefulness of therapeutic drug monitoring of mycophenolic acid: a computer simulation. Ther Drug Monit. 2005;27:163–7.
Jeong H, Kaplan B. Therapeutic monitoring of mycophenolate mofetil. Clin J Am Soc Nephrol. 2007;2:184–91.
Pawinski T, Hale M, Korecka M, Fitzsimmons WE, Shaw LM. Limited sampling strategy for the estimation of mycophenolic acid area under the curve in adult renal transplant patients treated with concomitant tacrolimus. Clin Chem. 2002;48:1497–504.
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.
Kuypers DR, Le Meur Y, Cantarovich M, Tredger MJ, Tett SE, Cattaneo D, et al. Consensus report on therapeutic drug monitoring of mycophenolic acid in solid organ transplantation. Clin J Am Soc Nepohrol. 2010;5:341–58.
Miura M, Satoh S, Niioka T, Kagaya H, Saito M, Hayakari M, et al. Limited sampling strategy for simultaneous estimation of the area under the concentration-time curve of tacrolimus and mycophenolic acid in adult renal transplant recipients. Ther Drug Monit. 2008;30:52–9.
Miura M, Satoh S, Niioka T, Kagaya H, Saito M, Hayakari M, et al. Early phase limited sampling strategy characterizing tacrolimus and mycophenolic acid pharmacokinetics adapted to the maintenance phase of renal transplant patients. Ther Drug Monit. 2009;31:467–74.
Yamaguchi K, Fukuoka N, Kimura S, Watanabe M, Tani K, Tanaka H, et al. Limited sampling strategy for the estimation of mycophenolic acid area under the concentration-time curve treated in Japanese living-related renal transplant recipients with concomitant extended-release tacrolimus. Biol Pharm Bull. 2013;36:1036–9.
Bullingham RE, Nicholls A, Hale M. Pharmacokinetics of mycophenolate mofetil (RS61443): a short review. Transplant Proc. 1996;28:925–9.
There are no funding sources for this report.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files.
TE, KN, TI and MO conceived and designed this study. YM, HK, YS and MO supervised the conduct of this study and data collection. TE provided statistical advice on study design and analyzed the data. TE drafted the manuscript, and all authors contributed substantially to its revision. MO takes responsibility for this paper as a whole. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki and its amendments, and was approved by the ethics committee of Mie University Graduate School of Medicine and Faculty of Medicine (Approval No.: 2505).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Multiple comparison of delta AUC0-12 between five groups classified according to concomitant drug usage (all patients). (PPTX 61 kb)
Multiple comparison of delta AUC0-12 between three groups classified according to concomitant drug usage (POD<31). (PPTX 53 kb)
Multiple comparison of delta AUC0-12 between five groups classified according to concomitant drug usage (POD≥31). (PPTX 56 kb)
Correlation between measured and estimated AUC0-12 estimated by using Pawinski’s estimation formula. (PPTX 56 kb)
Equations obtained using limited sampling time points for the prediction of MPA AUC0-12, Correlations with measured MPA AUC0-12 and RMSE of estimated AUC0-12. (XLSX 11 kb)
About this article
- Mycophenolate mofetil
- Limited sampling strategy
- Living-donor renal transplantation
- Mycophenolic acid
- Therapeutic drug monitoring