Pharmacokinetics of Mycophenolate Mofetil and Development of Limited Sampling Strategy in Early Kidney Transplant Recipients

The mycophenolate mofetil (MMF) dose management for optimization of post-transplant treatment especially the early postoperative phase has been well recognized. MMF is a pro-drug of mycophenolic acid (MPA) and is widely used in Chinese renal transplant patients. Until now, the pharmacokinetic (PK) characteristics and model for the area under the concentration–time curve for the 12-h (h) of exposure (AUC_{0-12 h}) of MPA (MPA-AUC_{0-12 h}) estimation were lacking for the new formulation of MMF dispersible tablet in renal transplant patients. The aims of the study were to investigate the PK characteristics of MMF dispersible tablet by detecting the active metabolite of MPA and to establish an accuracy and precision equation for calculating MPA-AUC_{0-12 h} by limited sampling strategy (LSS) in Chinese kidney transplant patients. A total of 60 postoperative kidney transplant recipients were given a multiple-dose of MMF dispersible tablet twice daily combination with tacrolimus (Tac) and steroids. On the 5th day post-transplantation, blood specimens were collected before drug administration and up to 12 h after MMF dispersible tablet administration. Non-compartmental PK analysis was used to determine the data obtained from individual patients. Multivariate stepwise regression analysis was used to develop models for predicting MPA-AUC_{0-12 h}. The 3- and 4-point sampling models using 2 h, 4 h, 8 h and 1 h, 2 h, 4 h and 8 h, respectively, allowed accurate estimation of MPA-AUC_{0-12 h}. PK parameters of MMF dispersible tablet were obtained and the 4-point LSS is the best model for accurate and precise estimation of MPA-AUC_{0-12 h}.

Introduction

Mycophenolate mofetil (MMF), an ester prodrug of mycophenolic acid (MPA), is the first-line immunosuppressant drug used with other drugs in prevention and treatment of graft rejection in solid organ transplantation. Chemical structures of MMF and MPA are shown in Figure 1. Dispersible tablet a new formulation of MMF is widely used in renal transplant patients in China. As MPA has a narrow therapeutic window and large inter-individual variability (Staatz and Tett, 2007; Kuypers et al., 2010; Mathew et al., 2010; Saint-Marcoux et al., 2011), therapeutic drug monitoring (TDM) of the area under the concentration–time curve for the 12-h (h) of exposure (AUC_{0-12 h}) of MPA (MPA-AUC_{0-12 h}) is essential to improve clinical outcomes (van Gelder and Shaw, 2005; de Winter et al., 2007; West-Thielke and Kaplan, 2007). Adequate MPA exposure is especially important to prevent acute rejection in the early renal transplant phase, especially the 1st week after transplantation (Tett et al., 2011; Barraclough et al., 2012b). The MPA-AUC_{0-12 h} has been recommended as the best marker for the MMF dose adjustment to achieve optimal efficacy and keep the toxicity to minimum in clinical application. However, increasing studies reported that the trough level (C₀) of MPA showed a poor correlation with the MPA-AUC_{0-12 h} (Mathew et al., 2010; Barraclough et al., 2012a; Chaabane et al., 2013; Cai et al., 2015). Owing to the requirement for frequent blood sampling, full MPA-AUC_{0-12 h} monitoring is laborious. The limited sampling strategy (LSS), using a limited number of blood samples to estimate the MPA-AUC_{0-12 h} (Baraldo et al., 2009; Barraclough et al., 2012a), solves this problem. The two formulations of MMF (dispersible tablet and capsule) were bioequivalent in healthy Chinese volunteers (Zhang et al., 2010). Due to the differences in the formulations of capsule and dispersible tablet, the PK of MPA may be different. The dynamic profile of the new MMF dispersible tablet in Chinese kidney transplant recipients is unclear. With the increasing use of MMF in renal transplant patients in China, the need for concentration-controlled dosing, based on MPA PK characteristics of MMF dispersible tablet, is becoming important in order to promote optimum safety and efficacy. It is of clinical importance to obtain a better understanding of the dynamic characters of MMF dispersible tablet. Until now, PK data of MPA, based on the new formulation of MMF was lacking in Chinese renal transplant patients. Meanwhile, there was no LSS of MMF dispersible tablet has been suggested for Chinese kidney transplant recipients. So, the accurate and convenient LSS of MMF dispersible tablet for predicting MPA-AUC_{0-12 h} in renal transplant recipients is needed to be developed based on the dynamic characters of the new formulation.

FIGURE 1

FIGURE 1. Chemical structures of MMF (A) and MPA (B).

The aims of the study were to investigate the PK characteristics of MPA in Chinese kidney transplant patients taking MMF dispersible tablet in association with steroids and tacrolimus (Tac) and to develop an accuracy and precision LSS for MPA-AUC_{0-12 h} calculation in large number of samples in such patients.

Materials and Methods

Patients Selection

Consecutive, prospective 60 patients were included in the study. All the patients were over 18 years and underwent the first kidney transplant in this single-center study. Patients were recruited after passing a physical examination and laboratory tests, which included blood biochemistry, hematology, and urine analysis. We excluded patients with cancer, patients with hematologic abnormality, hepatic abnormality, gastrointestinal disorder, pregnancy, systemic infection, or any acute disease, and patients who allergy to MMF dispersible tablet, MPA, steroids or Tac. All patients underwent the same therapeutic scheme of MMF dispersible tablet with Tac and steroids as concomitant immunosuppressive therapy. This research was conducted according to the principles of the Declaration of Helsinki (1964) and its later amendments. The clinical protocol and the informed-consent form were approved by the Ethics Committees of The First Affiliated Hospital of Zhengzhou University. All participants were informed of the details and procedures of the research before they signed a written informed consent.

Study Design

The study was carried out in two steps. We enrolled 40 patients in the first step as the model-building group. The model-building group was used to establish optimal LSS for predicting AUC of MPA.The 20 patients who were enrolled later, were used as the validation group. The validation group was used to assess the accurate and practical applicability of the optimal LSS. The study was carried out in the Clinical Nephrotransplantation Center of the First Affiliated Hospital of Zhengzhou University.

The renal transplant recipients received 500 or 750 mg of MMF dispersible tablet (250 mg/Tab) twice daily with water. The starting Tac dose was 0.1–0.15 mg/kg/day. Targeted Tac whole-blood trough concentration was 8–15 ng/mL and the Tac concentration was assayed by enzyme-multiplied immunotechnique. All patients received two intraoperative corticosteroid of 500 mg of methylprednisolone. Maintenance methylprednisolone dose was tapered to 500 mg on day 2, followed by a stepwise reduction to 375 mg on day 3, 250 mg on day 4, and 120 mg on day 5. Then methylprednisolone (16 mg/d) was administered. Patients fasted overnight before dosing and given standardized meals until 1.5 h after MMF dispersible tablets administration, serial blood samples (2 mL each) were drawn in vials containing EDTA at 0, 0.5, 1, 1.5, 2, 3, 4, 6, 8, 10, and 12 h (C₀, C_0.5, C₁, C_1.5, C₂, C₃, C₄, C₆, C₈, C₁₀, and C₁₂) after the oral morning dose of MMF dispersible tablet on the 5th day after transplantation. The samples were sent to lab on ice. The blood samples were centrifuged at 3000 r/min and plasma samples were separated. MPA concentrations were analyzed using the validated ultra-performance liquid chromatography ultraviolet (UPLC-UV) method (Zhang et al., 2017).

Non-compartmental PK Parameters Analysis

The PK analysis was performed using WinNonlin version 7.0 (Pharsight Corporation, Mountain View, CA, United States). Non-compartmental PK analysis was used to determine the data obtained from individual patients. The area under the concentration-time curve 0 to time (AUC_0-t) calculated by linear trapezoid method. Predose concentration (C₀), post-dose minimum and peak plasma concentration (C_min and C_max), time to reach C_max (T_max) were obtained directly from the concentration vs. time profiles. Total body clearance (CL) was calculated by dose/AUC_0-τ (τ was administration interval). Given the difference in doses, we normalized the AUC_{0-12 h} of 500–750 mg. Dose-normalized estimate of AUC_{0-12 h,norm} was derived by the corresponding value dividing the result that 750 mg divided 500 mg. Inter-individual variability was assessed using coefficient of variation (CV).

Development of LSSs

The study had sufficient sampling points to capture the full characteristics of MPA’s PK. Since, to date, model for estimation of MPA-AUC_{0-12 h} was lacking for MMF dispersible tablet in Chinese renal transplant patients. The time control plasma concentrations of MPA would be the most sensitive to the differences in the drug absorption rate, distribution volume and elimination rate. So, the LSS was performed to estimate MPA-AUC_{0-12 h}. The data were divided into a model-building group (n = 40) and a model validation group (n = 20). All results were presented as mean ± SD. The data were performed using SPSS version 18.0 (SPSS Inc., Chicago, IL, United States). Step 1: the fAUC_{0-12 h} was calculated by linear trapezoid method. Prediction equation models were derived by multiple regression analysis, using MPA concentrations at different time points. The fAUC_{0-12 h} as a dependent variable and the concentration of MPA at each time point were used as independent variable. The correlation coefficient (R²) was used to evaluate the regression level of the equation. The fAUC_{0-12 h} was compared with eAUC_{0-12 h} calculated using the selected equation of Step 1 in the validation group. The prediction bias were quantified as the percentage of mean prediction error (MPE%) and mean absolute prediction error (MAPE%). The two error parameters were calculated by the following equations:

$MPE\%=\frac{1}{n}{\displaystyle \sum _{i=1}^n\{({\text{eAUC}}_{0-12\text{h}}-{\text{fAUC}}_{0-12\text{h}})/{\text{fAUC}}_{0-12\text{h}}\}\times 100}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}(1)$

$MAPE\%=\frac{1}{n}{\displaystyle \sum _{i=1}^n\left\{\right|({\text{eAUC}}_{0-12\text{h}}-{\text{fAUC}}_{0-12\text{h}})/{\text{fAUC}}_{0-12\text{h}}\left|\right\}\times 100}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}(2)$

The Bland–Altman test (Bland and Altman, 1986) was performed to assess the consistency (bias and precision) between the fAUC_{0-12 h} and eAUC_{0-12 h} in the validation group. Bias and precision were visualized by plotting the average AUC_{0-12 h} that results from the estimated and full profile. The eAUC_{0-12 h} was calculated using the selected equation. The correlation between fAUC_{0-12 h} and eAUC_{0-12 h} was expressed by scatter diagram. The best model was selected based on account the values of r², predictive bias, and the Bland–Altman analysis.

Tolerability Evaluation

The investigators including clinicians and clinical pharmacists who monitored the patients for adverse events (AEs). The primary outcomes were 1-week incidence of acute rejection, graft function and other AEs during hospitalization. Delayed graft function was defined as dialysis requirement in the 1st week post-transplantation. Secondary outcomes included 6-month incidence of acute rejection, graft survival and other related AEs. Those events were also determined using patient interview (subjective complaints by patients, objective observation by clinicians) regarding the potential presence of AEs, such as diarrhea, nausea and so on, at any time during hospitalization and by comparing the results of laboratory tests (hematology, blood biochemistry, hepatic function and urinalysis), during the study and 6 months after completion of the study (follow-up visit).

Statistical Analyses

Differences in efficacy and safety end points were evaluated by the Chi-square or Fisher’s exact test. The correlations between AEs and covariates (age, body weight, MPA-AUC_{0-12 h}) were evaluated by Spearman test. All statistical analyses were performed with SPSS version 18.0. P < 0.05 was considered statistically significant.

Results

Patients

A total of 60 patients were enrolled in the study. Table 1 summarizes the baseline characteristics of the patients.

TABLE 1

TABLE 1. The demographic data and clinical parameters.

Non-compartmental Analysis

A total of 60 patients participated in the PK study. The mean plasma concentration–time curves of MPA after administration of MMF dispersible tablets in all patients, model-building group and validation group are shown in Figure 2. The MPA plasma PK parameters are summarized in Table 2. Inter-individual variability (CV) in AUC_{0-12 h}, C_max, C_min, C₀, CL, and T_max, respectively, were 15.19, 27.97, 89.72, 66.53, 23.56, and 80.78%. Due to the enterohepatic recirculation (EHR) most of the patients revealed a second small peak in the 4–12 h after taking MMF dispersible tablets.

FIGURE 2

FIGURE 2. Mean plasma concentration–time curves of MPA after administration of MMF dispersible tablets in all patients (A), model-building group and validation group (B).

TABLE 2

TABLE 2. Pharmacokinetic (PK) parameters of MPA after multiple dose of MMF dispersible tablet in 60 kidney transplant recipients.

Results of LSSs

The results of regression equations obtained in the study group between each sampling time point and the fAUC_{0-12 h} are shown in Table 3. The equation containing only the trough level showed a weak correlation (R² = 0.501) with fAUC_{0-12 h}.

TABLE 3

TABLE 3. Univariate correlation between the MPA plasma concentration at each time point and the AUC_{0-12 h}.

Equations for two-, three-, and four-point estimations of the MPA-AUC_{0-12 h} were obtained by multiple stepwise regression analysis. As seen in Table 4, with the number of sampling time points increasing, the linear correlations between the eAUC_{0-12 h} and fAUC_{0-12 h} also improved. Using the three selected equations obtained in the model building group, we compared the differences between eAUCs and fAUCs in the model validation group. The predictive results are shown in Figures 3, 4. All the three models were found to have good mean residuals and symmetries.

TABLE 4

TABLE 4. Multiple linear regression analysis of the AUC_{0-12 h}.

The differences between the eAUCs obtained with 2, 3, or 4 points and the fAUCs in the 20 validation patients were 3.78 ± 6.17, 3.41 ± 5.51, and 2.11 ± 3.50 μ⋅h/mL, respectively; 95% confidence intervals were (-8.31 to 15.88), (-7.39 to 14.22), and (-4.74 to 8.97), respectively.

For the two-point equation, the MPE% and MAPE% values were 10.28 and 12.99%, respectively. The mean values of fAUC and eAUC were not significantly different (42.13 ± 14.99 μg⋅h/mL and 45.92 ± 15.86 μg⋅h/mL, respectively). The agreements between fAUCs and eAUCs were tested by Bland–Altman Analysis, 10% (2/20) of patients were outside the 95% confidence interval (Figure 3A).

FIGURE 3

FIGURE 3. The Bland–Altman analysis agreement between the fAUC and eAUC using model 2 (A), model 3 (B), and model 4 (C).

For the three-point model, the MPE% and MAPE% values were 8.64 and 12.39%, respectively. The mean values of fAUC and eAUC were not significantly different (42.13 ± 14.99 μg⋅h/mL and 45.55 ± 16.54 μg⋅h/mL, respectively). The Bland–Altman Analysis confirmed the agreement between eAUCs and fAUCs, and 5% (1/20) of patients were outside the 95% confidence interval (Figure 3B).

The four-point model, the MPE% and MAPE% values were 5.26 and 8.35%, respectively. The mean values of fAUC and eAUC were not significantly different (42.13 ± 14.99 μg⋅h/mL and 44.24 ± 15.54 μg⋅h/mL, respectively). The Bland–Altman Analysis showed the best agreements between the eAUCs and fAUCs, and 5% (1/20) of patients were outside the 95% confidence interval (Figure 3C). The abilities of the three models to predict the MPA–AUC_{0-12 h} are depicted in Figure 4. Based on the above analysis, the 4-point method (C₈, C₂, C_4, and C₁) was the best LSS.

FIGURE 4

FIGURE 4. Correlations between eAUC and fAUC calculated using model 2 (A), model 3 (B), and model 4 (C).

Tolerability

No serious AEs occurred in the study. There were no patients who suffered from acute rejection, suspected acute rejection, delayed graft function or systemic infection during hospitalization and 6 months post-transplantation. Results of vital signs, and physical examinations were within normal limits for all participants during hospitalization. The Spearman correlations between AEs and age, body weight or MPA-AUC_{0-12 h} were 0.29, 0.24, and 0.12, respectively. Factor that significantly influenced the incidence of AEs in the Spearman test was MPA-AUC at day 5 (P = 0.023). The incidences of AEs in different MPA-AUC_{0-12 h} ranges were listed in Table 5. There were no significant differences in the incidences of AEs between the patients with MPA-AUC_{0-12 h} levels no more than 60 μg⋅h/mL and more than 60 μg⋅h/mL at the 5 days after transplantation. This may related to the small sample size. Though 60 μg⋅h/mL was not a useful threshold, a relationship was still seen between AUC and adverse effects. Mean MMF daily doses at discharge were similar with hospitalization. Two patients suffered constipation at 2 months post-transplantation. All these AEs were mild. The formulation was well tolerated in transplant patients in the study.

TABLE 5

TABLE 5. Comparison of incidences of AEs in different MPA-AUC_{0-12 h} ranges.

Discussion

The MMF dose management for optimization of post-transplant treatment especially the early postoperative phase has been well recognized. PK studies in the renal transplant population are essential to determine the correct dosage strategies for optimal efficacy. There is an increasing interest in use of MMF dispersible tablet in Chinese renal transplant patients. To date, most studies investigating the PK of MPA in this population applied a sparse sampling protocol, which may have limited the overall precision of the final model. In the present study, the full PK character was available for all the patients and can be used to provide limited sampling strategies information for the early transplant phase patients. The PK characteristics of MMF and equation for estimation of MPA-AUC_{0-12 h} were lacking for the new formulation of MMF in Chinese renal transplant patients. The present study developed LSSs of MMF dispersible tablet in renal transplant patients for the first time.

As previous studies reported, the PK of MPA showed a large interpatient variability (van Hest et al., 2005; Staatz and Tett, 2007; Kiang and Ensom, 2017; Yu et al., 2017). In our study after oral MMF dispersible tablet, the t_max occurred at 0.5 to 6 h after MMF dispersible tablet dosing (mean peak time point: 1.5 h post-dose) and C_max ranged from 4.26 to 38.4 μg/mL. The interindividual variability in AUC_{0-12 h}, C_max, C_min, C₀, CL, and T_max were large. Previous studies (Le Guellec et al., 2004; Yu et al., 2017) found that BW was positively correlated with CL of MPA in adult transplant patients. Kagaya et al. (2017) reported similar results that a lower BW was independently predictive for a higher MPA-AUC_{0-12 h}. CL of MPA was found significant time-dependency in previous study (van Hest et al., 2007). The CL of this study was 18.32 L/h, which was similar to those previously reported in adult kidney transplant recipients (de Winter et al., 2011; Yu et al., 2017). When given 750 mg twice daily, most patients can reach the target exposure in 1st week after renal transplantation.

Mycophenolic acid is primarily metabolized to the inactive metabolite 7-O-MPA glucuronide, which undergoes EHR (Armstrong et al., 2005). Due to EHR a second peak was identified in 6–12 h following MMF oral administration (Bullingham et al., 1998). However, the secondary concentration peaks were relatively small in most of our patients. The reason for little EHR is likely due to the weak continuous secretion of bile into the small intestine. Less 7-O-MPA glucuronide was deglucuronidated in the gut and re-entered the circulation as MPA. The enzyme-multiplied immunoassay technique (EMIT) and UPLC are two principal methods used to monitor MPA concentration. The UPLC-UV method is highly specific for determination of MPA and has low imprecision across the analytical range. In comparison with UPLC, EMIT contributed to a larger measurement errors due to its cross-reaction with acyl-MPA glucuronide (Yu et al., 2017). UPLC method used in this study is an ideal method for rapid detection with high accuracy, and good reproducibility.

The present study shows that the interindividual variability is large, which highlight the TDM of MPA in renal transplant patients. An appropriate AUC_{0-12 h} after drug administration, between 30–60 μg⋅h/mL, was associated with significant decrease in acute graft rejection in kidney transplant patients (Staatz and Tett, 2007). TDM of MPA plays an important role in performing personalized prescription for the prognosis of recipients. The LSS technique has been applied in estimating the MPA exposure of MMF capsule in the past decades. For MMF capsule, substantial research has shown that an LSS can predict MPA-AUC_{0-12 h} in different populations after transplantation (Kuypers et al., 2010; Mathew et al., 2010; Saint-Marcoux et al., 2011; Chaabane et al., 2013). However, those models cannot accurately estimate the AUC values of MMF dispersible tablet in most of our patients. This result may relate to the different PK profiles of MMF dispersible tablet with MMF capsule in Chinese population. MMF dispersible tablet lead a new PK profile which makes it different with MMF capsule in selecting the optimal model for TDM. So, special LSS models are needed for the new formulation of MMF to guide the clinical pharmacists and clinicians to scheme medication regimens for patients. Considering the huge PK variability of MPA, we enrolled 40 patients in the model-building group to establish more practical and precision models. Then we performed external validation with another 20 patients to verify the practicability of the model.

The MPA PK profile is affected by the associated calcineurin inhibitors (Atcheson et al., 2005; Barraclough et al., 2010). The cyclosporine-based studies have been reported to result lower total MPA exposure compared with the Tac-based regimen (Filler et al., 2000). So, in this study, we enrolled the Tac-based regimen patients. The association between the trough level of MPA and clinical outcome such as graft rejection is still controversial. Some researches have shown that the trough level of MPA may be correlated to the patient’s outcomes in kidney transplant recipients (Borrows et al., 2006; Miura et al., 2011). However, many studies have indicated a poor correlation between single-point concentrations and MPA-AUC_{0-12 h} (Barraclough et al., 2012a; Chaabane et al., 2013; Yamaguchi et al., 2013; Cai et al., 2015; Jia et al., 2017). The present study also revealed the same result. The C₀ level was not a reliable parameter to assess the MPA exposure (R² = 0.501). The best correlations between the single-point concentration and fAUC occurred at 3 h (R² = 0.538), 8 h (R² = 0.682), and 10 h (R² = 0.548) treatment with MMF dispersible tablet during the 1st week after kidney transplantation. These results are different with that found by the previous study (R² < 0.4 at any single-point concentration) in Chinese patients treated with MMF capsule in combination with Tac and steroids (Cai et al., 2015). Then, we used the stepwise analysis obtained the equations were listed in Table 4. In the past reports, there were also no models which include single sample concentration could accuracy predict the fAUC of MPA in renal transplant recipients (Pawinski et al., 2002; Baraldo et al., 2009; Miura et al., 2009; Yamaguchi et al., 2013; Cai et al., 2015). As a result, two time points (C₂ and C₈), three time points (C₂, C₄, and C₈) and four time points (C₁, C_2, C₄, and C₈) models were developed in this study. The frequency of sampling times was acceptable in all models.

In this study, the two-sample equation displayed a good agreement with the fAUC. The significant correlation for predicting MPA-AUC_{0-12 h} was C₂ and C₈ with acceptable MPE% and MAPE%. This is the best practicability model with the least blood samples based on the dynamic characters of MMF dispersible tablet in this study. However, in the literature, there were no two-sample models to predict the MPA-AUC_{0-12 h} in early kidney transplant recipients were recommended for clinical practice.

Compared with the two-sample model the 3-time point equation was better correlate with the fAUC. Our equation was different with the published three-point combinations (C₀/C₁/C₄, C₂/C₆/C₁₀, C₂/C₄/C₉, C₀/C₁/C₂, and C₀/C_0.5/C₂) (Pawinski et al., 2002; Muller et al., 2007; Miura et al., 2008; Chaabane et al., 2013; Cai et al., 2015). In the previous published three-point models (C₂/C₆/C₁₀) (Chaabane et al., 2013) the authors choose the patients on 7 days after transplantation. The other models (C₀/C₁/C₄, C₂/C₄/C₉, C₀/C₁/C₂, and C₀/C_0.5/C₂) (Pawinski et al., 2002; Muller et al., 2007; Miura et al., 2008; Chaabane et al., 2013) the patients were on the steady state period. The main reasons that lead to the differences between the published LSSs with our model (C₂/C₄/C₈) possibly were the patient ethnicity, number of patients, pathological conditions and formulation of MMF.

In this study, the best model to assess the fAUC as tested by Bland–Altman analysis was 8.36 + 7.49 ×C₈ + 1.34 ×C₂ + × 1.66 ×C₄ + 0.76 ×C₁ (R² = 0.948). This equation was based on four time points provided the most reliable estimation of the MPA-AUC_{0-12 h}. The model included absorption phase, distribution phase and elimination phase after the MMF dispersible tablet administration.

Considering the MPA kinetic characters that investigated in the present research, three model equations were suggested with strong correlations between the fAUCs and eAUCs. The long time span may restrict the models application for the outpatients, but it was a useful tool to guide the individualization prescription of MMF dispersible tablet dose for the inpatients, especially for early post-transplant phase. The best model was selected by taking into account the values of r², predictive bias, and the Bland–Altman analysis. The 4-point method (C₈, C₂, C₄, and C₁) was the best LSS. This model can be recommended for predicting MPA-AUC_{0-12 h} in early renal transplant recipients in clinical practice.

Conclusion

In this study, the PK characteristics of MMF dispersible tablet was studied in the early renal transplant patients. The result revealed new insights into the PK characters of MPA after multiple dose of MMF dispersible tablet. Large interindividual variability in MPA PK was demonstrated. The accurate and clinically feasible models for evaluating MPA exposure level in renal recipients after administration of MMF dispersible tablet, Tac and steroids were developed. This information can help to apply TDM of MPA more efficient. The 4-time point model was the best LSS. The applicability of this model in predicting MPA-AUC_{0-12 h} had been validated in early Chinese kidney transplant populations. The finding can be recommended in early Chinese adult renal transplant recipients to individualize the dosage regimen in clinical application.

Author Contributions

JZ, ZS, YL, and XZ designed the research. ZZ, JK, and GF performed the experiments. JY, JZ, LZ, and LhZ analyzed the results. All authors approved the final manuscript.

Funding

This work was supported by The First Affiliated Hospital of Zhengzhou University and Medical Science and Technology Project of Henan province (No. 201702044).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer DM and handling Editor declared their shared affiliation.

Acknowledgments

We thank all the patients for their cooperation and participation in this study.

References

Armstrong, V. W., Tenderich, G., Shipkova, M., Parsa, A., Koerfer, R., Schröder, H., et al. (2005). Pharmacokinetics and bioavailability of mycophenolic acid after intravenous administration and oral administration of mycophenolate mofetil to heart transplant recipients. Ther. Drug Monit. 27, 315–321. doi: 10.1097/01.ftd.0000163949.40686.0f

PubMed Abstract | CrossRef Full Text | Google Scholar

Atcheson, B. A., Taylor, P. J., Mudge, D. W., Johnson, D. W., Hawley, C. M., Campbell, S. B., et al. (2005). Mycophenolic acid pharmacokinetics and related outcomes early after renal transplant. Br. J. Clin. Pharmacol. 59, 271–280. doi: 10.1111/j.1365-2125.2004.02235.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Baraldo, M., Cojutti, P. G., Isola, M., Feruglio, M. T., Tursi, V., Livi, U., et al. (2009). Validation of limited sampling strategy for estimation of mycophenolic acid exposure during the first year after heart transplantation. Transplant Proc. 41, 4277–4284. doi: 10.1016/j.transproceed.2009.08.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Barraclough, K. A., Isbel, N. M., Franklin, M. E., Lee, K. J., Taylor, P. J., Campbell, S. B., et al. (2010). Evaluation of limited sampling strategies for mycophenolic acid after mycophenolate mofetil intake in adult kidney transplant recipients. Ther. Drug Monit. 32, 723–733. doi: 10.1097/FTD.0b013e3181fc8fbb

PubMed Abstract | CrossRef Full Text | Google Scholar

Barraclough, K. A., Isbel, N. M., Johnson, D. W., Hawley, C. M., Lee, K. J., McWhinney, B. C., et al. (2012a). A limited sampling strategy for the simultaneous estimation of tacrolimus, mycophenolic acid and unbound prednisolone exposure in adult kidney transplant recipients. Nephrology 17, 294–299. doi: 10.1111/j.1440-1797.2011.01560.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Barraclough, K. A., Staatz, C. E., Johnson, D. W., Lee, K. J., McWhinney, B. C., Ungerer, J. P., et al. (2012b). Kidney transplant outcomes are related to tacrolimus, mycophenolic acid and prednisolone exposure in the first week. Transpl. Int. 25, 1182–1193. doi: 10.1111/j.1432-2277.2012.01553.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bland, J. M., and Altman, D. G. (1986). Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1, 307–310. doi: 10.1016/s0140-6736(86)90837-8

CrossRef Full Text | Google Scholar

Borrows, R., Chusney, G., Loucaidou, M., James, A., Lee, J., Tromp, J. V., et al. (2006). Mycophenolic acid 12-h trough level monitoring in renal transplantation: association with acute rejection and toxicity. Am. J. Transplant. 6, 121–128. doi: 10.1111/j.1600-6143.2005.01151.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bullingham, R. E., Nicholls, A. J., and Kamm, B. R. (1998). Clinical pharmacokinetics of mycophenolate mofetil. Clin. Pharmacokinet. 34, 429–455. doi: 10.2165/00003088-199834060-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, W., Ye, C., Sun, X., Qin, K., Qin, Y., Zhao, D., et al. (2015). Limited sampling strategy for predicting area under the concentration-time curve for mycophenolic acid in Chinese adults receiving mycophenolate mofetil and tacrolimus early after renal transplantation. Ther. Drug Monit. 37, 304–310. doi: 10.1097/FTD.0000000000000165

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaabane, A., Aouam, K., Ben Fredj, N., Hammouda, M., Chadly, Z., El May, M., et al. (2013). Limited sampling strategy of mycophenolic acid in adult kidney transplant recipients: influence of the post-transplant period and the pharmacokinetic profile. J. Clin. Pharmacol. 53, 925–933. doi: 10.1002/jcph.125

PubMed Abstract | CrossRef Full Text | Google Scholar

de Winter, B. C., Mathot, R. A., Sombogaard, F., Vulto, A. G., and van Gelder, T. (2011). Nonlinear relationship between mycophenolate mofetil dose and mycophenolic acid exposure: implications for therapeutic drug monitoring. Clin. J. Am. Soc. Nephrol. 6, 656–663. doi: 10.2215/CJN.05440610

PubMed Abstract | CrossRef Full Text | Google Scholar

de Winter, B. C., Mathot, R. A., van Hest, R. M., and van Gelder, T. (2007). Therapeutic drug monitoring of mycophenolic acid: does it improve patient outcome? Expert Opin. Drug Metab. Toxicol. 3, 251–261. doi: 10.2165/00003088-200746010-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Filler, G., Zimmering, M., and Mai, I. (2000). Pharmacokinetics of mycophenolate mofetil are influenced by concomitant immunosuppression. Pediatr. Nephrol. 14, 100–104. doi: 10.1007/s004670050021

CrossRef Full Text | Google Scholar

Jia, Y., Peng, B., Li, L., Wang, J., Wang, X., Qi, G., et al. (2017). estimation of mycophenolic acid area under the curve with limited-sampling strategy in chinese renal transplant recipients receiving enteric-coated mycophenolate sodium. Ther. Drug Monit. 39, 29–36. doi: 10.1097/FTD.0000000000000360

PubMed Abstract | CrossRef Full Text | Google Scholar

Kagaya, H., Niioka, T., Saito, M., Inoue, T., Numakura, K., Habuchi, T., et al. (2017). Effect of hepatic drug transporter polymorphisms on the pharmacokinetics of mycophenolic acid in patients with severe renal dysfunction before renal transplantation. Xenobiotica 47, 916–922. doi: 10.1080/00498254.2016.1235742

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiang, T. K. L., and Ensom, M. H. H. (2017). Population pharmacokinetics of mycophenolic acid: an update. Clin. Pharmacokinet. 8, 1–12. doi: 10.1007/s40262-017-0593-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuypers, D. R., Le Meur, Y., Cantarovich, M., Tredger, M. J., Tett, S. E., Cattaneo, D., et al. (2010). Consensus report on therapeutic drug monitoring of mycophenolic acid in solid organ transplantation. Clin. J. Am. Soc. Nephrol. 5, 341–358. doi: 10.2215/CJN.07111009

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Guellec, C., Bourgoin, H., Buchler, M., Le Meur, Y., Lebranchu, Y., Marquet, P., et al. (2004). Population pharmacokinetics and Bayesian estimation of mycophenolic acid concentrations in stable renal transplant patients. Clin. Pharmacokinet. 43, 253–266. doi: 10.2165/00003088-200443040-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathew, B. S., Fleming, D. H., Annapandian, V. M., Prasanna, S., Prasanna, C. G., Chandy, S. J., et al. (2010). A reliable limited sampling strategy for the estimation of mycophenolic acid area under the concentration time curve in adult renal transplant patients in the stable posttransplant period. Ther. Drug Monit. 32, 136–140. doi: 10.1097/FTD.0b013e3181cd550f

PubMed Abstract | CrossRef Full Text | Google Scholar

Miura, M., Niioka, T., Kato, S., Kagaya, H., Saito, M., Habuchi, T., et al. (2011). Monitoring of mycophenolic acid predose concentrations in the maintenance phase more than one year after renal transplantation. Ther. Drug Monit. 33, 295–302. doi: 10.1097/FTD.0b013e3182197e38

PubMed Abstract | CrossRef Full Text | Google Scholar

Miura, M., Satoh, S., Niioka, T., Kagaya, H., Saito, M., Hayakari, M., et al. (2008). 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. 30, 52–59. doi: 10.1097/FTD.0b013e31815f5416

PubMed Abstract | CrossRef Full Text | Google Scholar

Miura, M., Satoh, S., Niioka, T., Kagaya, H., Saito, M., Hayakari, M., et al. (2009). Early phase limited sampling strategy characterizing tacrolimus and mycophenolic acid pharmacokinetics adapted to the maintenance phase of renal transplant patients. Ther. Drug Monit. 31, 467–474. doi: 10.1097/FTD.0b013e3181ae44b9

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller, H., Solari, S., Zuniga, C., Guerra, I., Troncoso, J., Ovalle, R., et al. (2007). Therapeutic drug monitoring of mycophenolic acid in kidney transplant patients: a abbreviated sampling strategy. Transplant. Proc. 39, 596–599. doi: 10.1016/j.transproceed.2006.12.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Pawinski, T., Hale, M., Korecka, M., Frrzsimmons, W. E., and Shaw, L. M. (2002). Limited sampling strategy for the estimation of mycophenolic acid area under the curve in adult renal transplant patients treated with concomitant tacrolimus. Clin. Chem. 48, 1497–1504.

Google Scholar

Saint-Marcoux, F., Vandierdonck, S., Premaud, A., Debord, J., Rousseau, A., and Marquet, P. (2011). Large scale analysis of routine dose adjustments of mycophenolate mofetil based on global exposure in renal transplant patients. Ther. Drug Monit. 33, 285–294. doi: 10.1097/FTD.0b013e31821633a6

PubMed Abstract | CrossRef Full Text | Google Scholar

Staatz, C. E., and Tett, S. E. (2007). Clinical pharmacokinetics and pharmacodynamics of mycophenolate in solid organ transplant recipients. Clin. Pharmacokinet. 46, 13–58. doi: 10.2165/00003088-200746010-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Tett, S. E., Saint-Marcoux, F., Staatz, C. E., Brunet, M., Vinks, A. A., Miura, M., et al. (2011). Mycophenolate, clinical pharmacokinetics, formulations, and methods for assessing drug exposure. Transplant. Rev. 25, 47–57. doi: 10.1016/j.trre.2010.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

van Gelder, T., and Shaw, L. M. (2005). The rationale for and limitations of therapeutic drug monitoring for mycophenolate mofetil in transplantation. Transplantation 80, S244–S253. doi: 10.1097/01.tp.0000186380.61251.fc

PubMed Abstract | CrossRef Full Text | Google Scholar

van Hest, R. M., van Gelder, T., Bouw, R., Goggin, T., Gordon, R., Mamelok, R. D., et al. (2007). Time-dependent clearance of mycophenolic acid in renal transplant recipients. Br. J. Clin. Pharmacol. 63, 741–752. doi: 10.1111/j.1365-2125.2006.02841.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van Hest, R. M., van Gelder, T., Vulto, A. G., and Mathot, R. A. (2005). Population pharmacokinetics of mycophenolic acid in renal transplant recipients. Clin. Pharmacokinet. 44, 1083–1096. doi: 10.2165/00003088-200544100-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

West-Thielke, P., and Kaplan, B. (2007). Therapeutic monitoring of mycophenolic acid: is there clinical utility? Am J Transplant. 7, 2441–2442. doi: 10.1111/j.1600-6143.2007.01963.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, K., Fukuoka, N., Kimura, S., Watanabe, M., Tani, K., Tanaka, H., et al. (2013). 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. 36, 1036–1039. doi: 10.1248/bpb.b12-01095

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Z. C., Zhou, P. J., Wang, X. H., Françoise, B., Xu, D., Zhang, W. X., et al. (2017). Population pharmacokinetics and Bayesian estimation of mycophenolic acid concentrations in Chinese adult renal transplant recipients. Acta Pharmacol. Sin. 38, 1566–1579. doi: 10.1038/aps.2017.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Jia, M. M., Zuo, L. H., Li, N., Luo, Y. G., Sun, Z., et al. (2017). Nonlinear relationship between enteric-coated mycophenolate sodium dose and mycophenolic acid exposure in han kidney transplantation recipients. Acta Pharm. Sin. B 7, 347–352. doi: 10.1016/j.apsb.2016.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Tao, Y., Zhu, Y., and Zhu, D. (2010). Bioequivalence and pharmacokinetic comparison of two mycophenolate mofetil formulations in healthy chinese male volunteers: an open-label, randomized-sequence, single-dose, two-way crossover study. Clin. Ther. 32, 171–178. doi: 10.1016/j.clinthera.2010.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar