Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 03 February 2026

Sec. Alloimmunity and Transplantation

Volume 17 - 2026 | https://doi.org/10.3389/fimmu.2026.1710261

Establishment of accurate estimation equations for area under the concentration curves using simultaneous limited sampling for extended-release tacrolimus and mycophenolic acid in kidney transplant recipients

Kyoko MinamisonoKyoko Minamisono1Toshimasa NakaoToshimasa Nakao1Takehiro OhyamaTakehiro Ohyama1Sho NishidaSho Nishida1Hajime SasakiHajime Sasaki2Takayuki HiroseTakayuki Hirose3Kiyohiko HottaKiyohiko Hotta3Makiko MienoMakiko Mieno4Daiki Iwami*Daiki Iwami1*
  • 1Division of Renal Surgery and Transplantation, Department of Urology, Jichi Medical University, Shimotsuke, Japan
  • 2Department of Kidney Transplant Surgery, Sapporo City General Hospital, Sapporo, Japan
  • 3Department of Renal and Genitourinary surgery, Faculty of Medicine, Hokkaido University, Sapporo, Japan
  • 4Division of Center for Information, Jichi Medical University, Shimotsuke, Japan

Introduction: This study aimed to develop a limited sampling strategy (LSS) and predictive equations to accurately estimate the areas under the concentration-time curves (AUC) of extended-release tacrolimus (TAC-ER) and mycophenolic acid (MPA).

Methods: A retrospective analysis of Japanese kidney transplant recipients yielded 90 TAC-ER AUC0-24 (23 patients) and 80 MPA AUC0-12 (29 patients) datasets, which were randomly split into learning and validation datasets. Training datasets were used to generate the LSS model equations based on multiple linear regression analysis, and the coefficient of determination (R2) was used to assess the goodness of fit of regression models. Validation datasets applied the selected training equations to compute error indices, Passing-Bablok’s Kendall’s τ, and Bland–Altman limits of agreement, thereby assessing predictive bias, accuracy, and precision.

Results and discussion: Four equations (C0-C1-C6, C0-C1-C2-C6, C0-C1-C3-C6, C0-C1-C4-C6) showed strong correlations with the actual AUC (R² > 0.95), with the validation identifying C0-C1-C3-C6 as the most reliable for both TAC-ER and MPA. This study demonstrated that LSS using C0-C1-C3-C6 reliably and accurately estimated both the actual TAC-ER AUC0-24 and MPA AUC0-12 simultaneously in kidney transplant recipients. These equations can be feasibly implemented in outpatient clinical settings to reduce time and cost.

Introduction

The optimization of drug dosing through therapeutic drug monitoring (TDM) is essential for improving graft and patient outcomes after kidney transplantation. Because immunosuppressive drugs used after kidney transplantation have a narrow therapeutic range, patients are prone to developing infections or renal toxicity due to excessive immunosuppression or rejection due to insufficient immunosuppression (13).

Trough concentration (C0), which is mostly used in clinical practice for TDM, may not accurately represent overall drug exposure in all cases, potentially leading to suboptimal dosing, especially for mycophenolate mofetil (MMF) (47). Area under the concentration-time curve (AUC) was informative utilized as one of the key pharmacokinetic parameters (8). Moreover, accurately determining AUC typically requires frequent blood sampling (often 8–10 times), which is impractical in clinical settings because of the burden on patients and associated costs. Limited sampling strategies (LSS) involving 3–5 blood samples have been proposed as valuable tools for accurately estimating the actual AUC while reducing the burden on patients.

Tacrolimus (TAC) and MMF are commonly co-administered in >90% of kidney transplant recipients (913). While many centers administer MMF at a fixed dose, numerous reports have highlighted the importance of monitoring drug exposure using AUC to achieve individualized therapy. Accordingly, MPA exposure is assessed via MPA-AUC in the United States, Europe, China, and Japan, following regional guidelines and recommendations (1418). Similarly, therapeutic drug monitoring of TAC, particularly once-daily extended-release TAC (TAC-ER), remains essential to ensure adequate and personalized immunosuppressive therapy (9, 19).

However, to date, a common set of limited sampling points for simultaneously predicting the AUC of once-daily extended-release TAC (TAC-ER) and mycophenolic acid (MPA), the pharmacologically active metabolites of twice-daily MMF in kidney transplant recipients, has not been established. If the AUCs of both MPA and TAC-ER can be accurately estimated using a limited set of common sampling points, this approach could serve as an effective tool for precise and efficient monitoring of drug exposure in patients.

This study aimed to develop an LSS capable of simultaneously predicting the AUC of TAC-ER and MPA in kidney transplant recipients. Predictive equations for the actual AUC of TAC-ER and MPA, using common blood sampling points and maintaining high precision and accuracy, have the potential to contribute to the prevention of post-transplant complications such as infections and graft rejection.

Materials and method

Patients and protocols

This study included 108 Japanese kidney transplant recipients who underwent transplantation at the Hokkaido University or Sapporo City General Hospital. Between 2006 and 2010, these patients underwent a series of nine blood sample collections before and after drug administration to determine the actual AUC for TAC-ER or MPA.

In this study, data were retrospectively extracted from past clinical records, with a particular focus on nine previously conducted blood sampling points to determine the actual AUC values. Of the 108 included individuals, retrospective data collection yielded 90 actual TAC-ER AUC data points from 57 patients after renal transplantation. Similarly, for the remaining 51 individuals, retrospective analysis provided 80 MPA AUC data points (Figure 1).

Figure 1
Flowchart illustrating the division of TAC-ER and MPA patient data into learning and validation datasets. Learning datasets include 48 TAC-ER and 37 MPA data sets, focusing on equation generation, correlation evaluation, and determination of limited sampling equations. Validation datasets comprise 42 TAC-ER and 43 MPA data sets. The predicted AUC is calculated and evaluated using the learning cohort's equation, with evaluation by MSE, RMSE, MAE, MAPE, Passing-Bablok regression, and Bland-Altman analysis.

Figure 1. Flowchart for dividing the pharmacokinetic data of TAC-ER and MPA in kidney transplant recipients into learning and validation datasets. Pharmacokinetic data were obtained from 52 renal transplant recipients for TAC-ER and 29 recipients for MPA. Patients’ blood concentration data from nine sampling points were randomly divided into learning and validation datasets. For TAC-ER, 90 data points were assigned, while for MPA, 80 data points were included. The learning dataset was used to derive predictive equations for AUC based on 2–4 sampling points using linear regression. The validation dataset evaluated prediction precision using MSE, RMSE, MAE, and MAPE, Passing-Bablok regression, and the Bland–Altman analysis. TAC-ER, extended-release tacrolimus; MPA, mycophenolic acid; MSE, mean squared error; RMSE, root mean squared error; MAE, mean absolute error; MAPE, mean absolute percentage error; AUC, area under the curve;.

Using data from nine-point blood sampling (TAC-ER: pre-dose and post-dose at 0.5, 1, 2, 3, 4, 6, 12, and 24 h; MPA: pre-dose, and post-dose at 0.5, 1, 2, 3, 4, 6, 8, and12 hours), which has already been used in routine clinical practice for measuring TAC-ER or MPA blood concentrations to calculate actual AUC by trapezoidal method, we analyzed the pharmacokinetics of these drugs. The Jichi Medical University Hospital Bioethics Committee for Clinical Research approved the study protocol (approval number: CU23-172).

Analysis

All the datasets were randomly split into two distinct datasets: learning and validation. Initially, multiple regression analysis was applied to the learning datasets to derive predictive equations for the actual AUC of TAC and MPA based on all possible combinations of blood concentrations at various sampling points. The actual AUC was calculated using the linear trapezoidal method. Subsequently, prediction equations were generated through linear regression analysis using all possible non-repetitive combinations of two, three, and four sampling points selected from the nine sampling time points (predicted AUC = α + β0 C0 + β1C1 + … +βnCn).

The coefficient of determination (R2) of each prediction equation was calculated and compared to identify the most accurate model. The R2 value was used to identify the optimal regression equation by selecting those with R2 values >0.95. The validity and reliability of these equations were subsequently evaluated using validation datasets. For validation, error metrics, including the mean squared error (MSE), root mean squared error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE) were calculated and compared. MAPE values <10% were considered highly reliable, whereas those above 20% indicated the need for model improvement. Passing-Bablok regression and Bland–Altman analyses were employed to assess agreement and correlation between predicted and actual values (20, 21). Passing-Bablok regression analysis was used to confirm the monotonic relationship between the two variables. The Bland–Altman analysis is a method for evaluating the agreement between two measurement methods. It is used to compare a new measurement method with a standard measurement method and is suitable for visually evaluating bias (average difference) and measurement variability. Limits of agreement (LOA) were used to evaluate the agreement between the measurement methods, the narrower the limit, the higher the agreement between the measurement methods.

The standard deviation (SD) of differences is an index that indicates the degree of variability between measurement methods, the smaller the limit, the higher the agreement.

Results

Baseline characteristics

Baseline characteristics of the TAC-ER learning (n=48), TAC-ER validation (n=42), MPA learning (n=37), and MPA validation (n=43) datasets included in this study are summarized in Tables 1 and 2. For both TAC-ER and MPA, there were no significant differences in sex ratio, body mass index, dose per body weight, trough levels, or measured AUC between the learning and validation datasets.

Table 1
www.frontiersin.org

Table 1. Clinical characteristics of renal transplant recipients (TAC-ER).

Table 2
www.frontiersin.org

Table 2. Clinical characteristics of renal transplant recipients (MPA).

Overview of the pharmacokinetic profiles

Pharmacokinetic profiles of TAC-ER and MPA were divided into learning and validation datasets, respectively. In most patients treated with TAC-ER, blood concentrations stabilized within 4 h post-administration, with >60% reaching peak levels within 3 h (Figures 2A, C). MPA typically exhibits a biphasic pattern in blood concentration owing to enterohepatic circulation. However, in this study, most MPA blood concentrations reached their peak levels within 2 h, and the average (SD) plasma concentration did not provide sufficient data to clearly demonstrate a biphasic pattern (Figures 2B, D). The correlation between the single-point drug concentration and actual AUC was evaluated for TAC-ER and MPA in the learning datasets. TAC-ER showed the best correlation value at C0 (R2=0.759), followed by C6 (R2=0.672) and C3 (R2=0.618). For MPA, C0 was lower than that of TAC-ER (R2=0.497), whereas C4 demonstrated the best correlation (R2=0.649). Thus, the C0 single-point equation for TAC-ER demonstrated a relatively good correlation, whereas the predictive ability of MPA showed an inferior correlation (Supplementary Table S1).

Figure 2
Graphs labeled A, B, C, and D depict TAC and MPA concentrations over time post-ingestion. Each graph shows multiple line plots with an initial peak followed by a decline. Graphs A and C measure TAC concentrations (ng/dL) over twenty-four hours. Graphs B and D measure MPA concentrations (μg/mL) over twelve hours. Both sets include error bars indicating variability.

Figure 2. Pharmacokinetic profiles of TAC-ER and MPA. Overview of the pharmacokinetic profiles of the learning (A; TAC-ER and B; MPA) and validation (C; TAC-ER and D; MPA) datasets in patients undergoing renal transplantation. The upper graphs represent the pharmacokinetics of individual data, while the lower graphs show their mean ± standard deviation. TAC-ER, extended-release tacrolimus; MPA, mycophenolic acid.

Selection of the best equations for TAC-ER AUC0–24 and MPA AUC0-12 in the learning datasets using limited sampling strategy

The optimal linear regression models for predicting the actual TAC-ER AUC0–24 and MPA AUC0–12 were selected based on the highest R2 values among all the models (Supplementary Table S1). Table 3 summarizes the predictive performance of the equations that have a good correlation with a predictive AUC of R2 > 0.85.

Table 3
www.frontiersin.org

Table 3. Correlation between predicted AUC calculated by various sampling time combination and actual AUC using linear regression analysis.

C0 is easily obtainable for evaluating intra-individual drug fluctuations in outpatient settings. C0 has also been reported to have a strong correlation with the actual AUC, particularly for TAC-ER; therefore, it was included in all selected equations (4). Clinically feasible sampling times for outpatient use were limited to within 6 h after administration, with a maximum of four sampling points, allowing drug concentration measurements in an outpatient setting.

The equations that included C0 did not show inferior R2 values compared with those without C0. Predictive equations incorporating multiple sampling points outperformed the single-point equations. Among the three-point models, C0-C1-C6 (TAC-ER: R2 = 0.958; MPA: R2 = 0.962) exhibited the highest R2 values. For four-point models, C0-C1-C2-C6 (TAC-ER: R2 = 0.959, MPA: R2 = 0.962), C0-C1-C3-C6 (TAC-ER: R2 = 0.962, MPA: R2 = 0.967), and C0-C1-C4-C6 (TAC-ER: R2 = 0.969, MPA: R2 = 0.968) demonstrated superior results to the other four-points model. Therefore, these four prediction models were used to validate the datasets.

Limited sampling with C0-C1-C3-C6 showed the best predictive performance for TAC-ER AUC0–24 and MPA AUC0–12 in the validation datasets

We applied the equations to the validation datasets and evaluated them using MSE, RMSE, MAE, MAPE, Kendall’s τ derived from Passing-Bablok regression analysis, and LOA and SD of differences derived from the Bland–Altman analysis (Table 4). In the Passing-Bablok regression plot, all the equations showed good agreement with a regression line slope close to 1 and an intercept near 0 (Figure 3). All the four models performed similarly for TAC-ER, whereas for MPA, the C0-C1-C6 model yielded a relatively lower value (Kendall’s τ = 0.780). In the Bland–Altman plot, no significant differences were observed between the predicted and actual AUC values with the LOA, including a 95% difference (Figure 4).

Table 4
www.frontiersin.org

Table 4. Prediction performance evaluated on validation data sets with limited sampling formulas for TAC-ER and MPA AUCs elected in learning data sets.

Figure 3
Scatter plots comparing predicted and actual AUC values using different measurement sets labeled C0 to C6. The top row is for TAC-ER and the bottom for MPA. Each plot displays a fitted line and Kendall's Tau values ranging from 0.827 to 0.889, indicating correlation strength.

Figure 3. Passing-Bablok regression analysis comparing predicted AUC with actual AUC. The solid line represents the regression line and the dashed line indicates the line of identity (slope = 1). This analysis assessed the agreement between predicted and observed values. AUC, area under the curve.

Figure 4
Scatter plot grid showing the difference between predicted and measured area under the curve (AUC) against the mean of predicted and actual AUC. There are two rows: the top labeled TAC-ER and the bottom labeled MPA, with varying combinations of C0, C1, C2, C3, C4, and C6 at the top of each column. Each plot includes a horizontal line at zero and two additional parallel lines above and below.

Figure 4. Bland–Altman plots comparing predicted AUC with actual AUC. The dotted lines represent 95% upper and lower limits. The solid line represents an average of the arithmetic differences between predicted AUC and measured AUC. AUC, area under the curve.

For TAC-ER, all the four models demonstrated comparable predictive performance (Table 4).

For MPA, C0-C1-C3-C6 demonstrated the best performance values in all the parameters (R2 = 0.967, MSE = 24.147, RMSE = 4.914, MAE = 3.778, MAPE = 10.284%, and Kendall’s τ = 0.844).

Therefore, there was no clear superiority among the four combinations in TAC-ER; however, the C0-C1-C3-C6 combination was considered the best in MPA. Based on the above, it was thought that the C0-C1-C3-C6 combination was the ideal combination of blood collection points in terms of balancing the burden on the patient and accuracy in the simultaneous LSS of TAC-ER and MPA.

Discussion

Appropriate therapeutic monitoring of immunosuppressive agents is crucial in preventing serious complications such as opportunistic infections due to over-immunosuppression and acute rejection episodes resulting from inadequate drug exposure. While determination of the actual AUC through complete pharmacokinetic profiling would provide the most accurate assessment of drug exposure, obtaining 9 blood samples at multiple time points over 12–24 hours from all patients is impractical in routine clinical practice. Such intensive sampling protocols impose substantial logistical burdens on both patients and healthcare staff, making it unfeasible to implement in everyday transplant management. Therefore, developing a limited sampling strategy (LSS) that maintains adequate predictive accuracy while requiring only 3–4 blood samples obtainable during a single outpatient visit is of considerable clinical importance for optimizing personalized immunosuppressive therapy in the ambulatory care setting.

While many centers administer MMF at a fixed dose, many reports have highlighted the importance of monitoring drug exposure using AUC to achieve individualized therapy. Accordingly, MPA exposure is assessed via MPA-AUC in the United States, Europe, China, and Japan, following regional guidelines and recommendations (1418). Similarly, therapeutic drug monitoring of TAC, particularly once-daily extended-release TAC (TAC-ER), remains essential to ensure adequate and personalized immunosuppressive therapy (9, 19).

We developed and compared predictive equations for TAC-ER and MPA AUCs using various combinations of available sampling time points. Among them, we selected four combinations of sampling points in which the prediction equations using a combination of C0-C1-C6, C0-C1-C2-C6, C0-C1-C3-C6, and C0-C1-C4-C6, demonstrated sufficient predictive accuracy. Although the predictive accuracy improves with additional sampling points, our findings indicate that a clinically feasible number and timing of samples still achieve a sufficiently high predictive performance. There was no clear superiority among the four combinations in TAC-ER; however, the C0-C1-C3-C6 combination was considered the best in MPA. In the validation datasets, predictive assessments demonstrated that for both TAC-ER and MPA, the C0-C1-C3-C6 combination consistently achieved the highest accuracy across all metrics. The combination of blood sampling times can simultaneously measure the AUC0−24 for once-daily TAC-ER and AUC0−12 for twice-daily MMF, which have different drug profiles. While it is scientifically reasonable to construct an AUC estimation formula including C12, in actual clinical practice, TAC and MPA concentrations measurement at night (7pm or 8pm) is required not only for inpatients but also for outpatients. Thus, we prioritized a combination of blood sampling points applicable in both inpatient and outpatient settings. Note that the formulas including C12 are listed in Supplementary Table S1 along with their correlation coefficients. Any of the correlation coefficient (R2) for the formulae including C12 showed no significant improvement compared to formulae using C0, C1, C3, and C6. To ensure clinical feasibility in an outpatient setting, we determined that 4 blood concentrations over a 6-hour period represent the optimal approach. As stated in the manuscript, these constraints are intended to minimize the burden on patients and medical staffs while maintaining practicality. A more extensive (longer) sampling period could potentially improve accuracy of estimation of AUC but is impractical for routine clinical care. The selected model demonstrated high predictive performance (R2 > 0.95) even within this limited sampling framework. Conversely, inclusion of C6 is essential to maintain model accuracy. As shown in Supplementary Table S1, omitting C6 and constructing the predictive formula using only sampling points within a shorter 4-hour window results in reduced predictive performance. As MPA undergoes enterohepatic circulation, it is recommended that at least one blood sampling point after six hours be included for accurate AUC estimation (22). Taking this into account, our formula is also reasonable.

To the best of our knowledge, this is the first study to report the development and validation of an LSS using regression equations to predict the AUC0–24 of TAC-ER and AUC0–12 of MPA. Furthermore, this study included patients with varying postoperative periods, sex, and age, ensuring that the developed equation is applicable to diverse patient backgrounds.

TAC-ER, an immunosuppressive drug containing TAC as its active ingredient, was developed by Astellas Pharma in Japan (23, 24). Compared with PrografR, which requires twice-daily dosing, TAC-ER produces a more gradual peak in blood concentration after post-administration (23). This reduces the risk of side effects associated with fluctuations in blood levels, thereby enhancing patient safety and tolerance (2527). MMF is converted to MPA in the body and metabolized as a glucuronide conjugate in the liver (28, 29). MPA is characterized by a complex metabolism and is significantly affected by enterohepatic recirculation, leading to a biphasic blood concentration profile and considerable interindividual variability (30). As for sample size in the current study, testing Kendall’s correlation coefficient revealed that even with a correlation coefficient of 0.8 and 20 cases, the power exceeded 98%. Although the retrospective nature of the study inherently imposed limitations on sample size, the results obtained demonstrate sufficient power, and the sample size is considered adequate.

This study has some limitations. All cases were taking TAC and MMF orally. However, blood sampling did not measure two drug concentrations simultaneously from the same patients; instead, the datasets were collected from the clinical records where only TAC was measured or where only MPA was measured. This study utilized these datasets. TAC and MMF/MPA have totally different metabolism pathways dependent on CYP3A4/5 (1) and UGT2B7 (28), respectively. Although the TAC and MMF were administered independently, it is considered that measuring the concentrations of the two drugs independently should pose little issue. However, potential pharmacokinetic interactions between the two drugs might not be excluded. Further validation in patients receiving both drugs is warranted. Secondly, this study included multiple AUC measurements obtained from the same subjects at different time points. Consequently, it cannot be ruled out that disparities in contribution may have arisen during the creation of prediction models across individual patients. On the other hand, it may also have enabled the development of estimation models capable of addressing diverse post-transplant periods. Therefore, further validation within a more controlled prospective cohort is desirable. Future studies should evaluate clinical outcomes and assess how post-transplantation periods and patient conditions affect the accuracy of the prediction equations in more-controlled cohort prospectively.

In conclusion, this study demonstrated that LSS using C0-C1-C3-C6 reliably and accurately estimated both the actual TAC-ER AUC0–24 and MPA AUC0–12 in kidney transplant recipients. LSS using C0-C1-C3-C6 can be feasibly implemented in outpatient clinical settings to reduce time and cost. The LSS algorithms validated in this study have proven to be highly useful tools for the TDM of TAC-ER and MPA exposure in clinical settings.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by The Jichi Medical University Hospital Bioethics Committee for Clinical Research. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because To simplify consent procedures in retrospective studies.

Author contributions

KM: Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing. TN: Writing – review & editing. TO: Writing – review & editing. SN: Writing – review & editing. HS: Data curation, Writing – review & editing. TH: Data curation, Writing – review & editing. KH: Data curation, Writing – review & editing. MM: Formal analysis, Methodology, Writing – review & editing. DI: Conceptualization, Data curation, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing.

Funding

The author(s) declared that financial support was not received for this work and/or its publication.

Acknowledgments

The authors would like to express their gratitude to the patients who participated in this study and clinical research nurses at the Hokkaido University and Sapporo City General Hospital for collecting blood samples. We also extend our appreciation to the clinical research teams at these institutions for providing blood samples and their valuable cooperation in conducting this study. We also thank Editage for English language editing.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2026.1710261/full#supplementary-material

Abbreviations

AUC, Area Under the Curve; LOA, Limits of Agreement; LSS, Limited Sampling Strategy; MAE, Mean Absolute Error; MAPE, Mean Absolute Percentage Error; MMF, Mycophenolate mofetil; MPA, Mycophenolic Acid; MSE, Mean Squared Error; RMSE, Root Mean Square Prediction Error; TAC, Tacrolimus; TAC-ER, Extended-Release Tacrolimus; TDM, Therapeutic Drug Monitoring.

References

1. Yu M, Liu M, Zhang W, and Ming Y. Pharmacokinetics, Pharmacodynamics and Pharmacogenetics of Tacrolimus in Kidney Transplantation. Curr Drug Metab. (2018) 19:513–22. doi: 10.2174/1389200219666180129151948

PubMed Abstract | Crossref Full Text | Google Scholar

2. Bauer AC, Franco RF, and Manfro RC. Immunosuppression in Kidney Transplantation: State of the Art and Current Protocols. Curr Pharm Des. (2020) 26:3440–50. doi: 10.2174/1381612826666200521142448

PubMed Abstract | Crossref Full Text | Google Scholar

3. Meneghini M, Bestard O, and Grinyo JM. Immunosuppressive drugs modes of action. Best Pract Res Clin Gastroenterol. (2021) 54-55:101757. doi: 10.1016/j.bpg.2021.101757

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wong KM, Shek CC, Chau KF, and Li CS. Abbreviated tacrolimus area-under-the-curve monitoring for renal transplant recipients. Am J Kidney Dis. (2000) 35:660–6. doi: 10.1016/S0272-6386(00)70013-8

PubMed Abstract | Crossref Full Text | Google Scholar

5. Pisitkun T, Eiam-Ong S, Chusil S, Praditpornsilpa K, Pansin P, and Tungsanga K. The roles of C4 and AUC0–4 in monitoring of tacrolimus in stable kidney transplant patients. Transplant Proc. (2002) 34:3173–5. doi: 10.1016/S0041-1345(02)03684-9

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kuypers DR. Immunosuppressive drug monitoring - what to use in clinical practice today to improve renal graft outcome. Transpl Int. (2005) 18:140–50. doi: 10.1111/j.1432-2277.2004.00041.x

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

8. Wallemacq P, Armstrong VW, Brunet M, Haufroid V, Holt DW, Johnston A, et al. Opportunities to optimize tacrolimus therapy in solid organ transplantation: report of the European consensus conference. Ther Drug Monit. (2009) 31:139–52. doi: 10.1097/FTD.0b013e318198d092

PubMed Abstract | Crossref Full Text | Google Scholar

9. Kidney Disease: Improving Global Outcomes Transplant Work G. KDIGO clinical practice guideline for the care of kidney transplant recipients. Am J Transplant. (2009) 9:S1–155. doi: 10.1111/j.1600-6143.2009.02834.x

PubMed Abstract | Crossref Full Text | Google Scholar

10. Lebranchu Y, Bridoux F, Buchler M, Le Meur Y, Etienne I, Toupance O, et al. Immunoprophylaxis with basiliximab compared with antithymocyte globulin in renal transplant patients receiving MMF-containing triple therapy. Am J Transpl. (2002) 2:48–56. doi: 10.1034/j.1600-6143.2002.020109.x

PubMed Abstract | Crossref Full Text | Google Scholar

11. Webster AC, Woodroffe RC, Taylor RS, Chapman JR, and Craig JC. Tacrolimus versus ciclosporin as primary immunosuppression for kidney transplant recipients: meta-analysis and meta-regression of randomized trial data. BMJ. (2005) 331:810. doi: 10.1136/bmj.38569.471007.AE

PubMed Abstract | Crossref Full Text | Google Scholar

12. Ekberg H, Tedesco-Silva H, Demirbas A, Vitko S, Nashan B, Gurkan A, et al. Reduced exposure to calcineurin inhibitors in renal transplantation. N Engl J Med. (2007) 357:2562–75. doi: 10.1056/NEJMoa067411

PubMed Abstract | Crossref Full Text | Google Scholar

13. Johnson C, Ahsan N, Gonwa T, Halloran P, Stegall M, Hardy M, et al. Randomized trial of tacrolimus (Prograf) in combination with azathioprine or mycophenolate mofetil versus cyclosporine (Neoral) with mycophenolate mofetil after cadaveric kidney transplantation. Transplantation. (2000) 69:834–41. doi: 10.1097/00007890-200003150-00028

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

15. 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. doi: 10.1097/FTD.0b013e31815f5416

PubMed Abstract | Crossref Full Text | Google Scholar

16. Weber LT, Schutz E, Lamersdorf T, Shipkova M, Niedmann PD, Oellerich M, et al. Therapeutic drug monitoring of total and free mycophenolic acid (MPA) and limited sampling strategy for determination of MPA-AUC in pediatric renal transplant recipients. The German Study Group on Mycophenolate Mofetil (MMF) Therapy. Nephrol Dial Transpl. (1999) 14:34–5. doi: 10.1093/ndt/14.suppl_4.34

PubMed Abstract | Crossref Full Text | Google Scholar

17. Xiang H, Zhou H, Zhang J, Sun Y, Wang Y, Han Y, et al. Limited Sampling Strategy for Estimation of Mycophenolic Acid Exposure in Adult Chinese Heart Transplant Recipients. Front Pharmacol. (2021) 12:652333. doi: 10.3389/fphar.2021.652333

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hulin A, Blanchet B, Audard V, Barau C, Furlan V, Durrbach A, et al. Comparison of 3 estimation methods of mycophenolic acid AUC based on a limited sampling strategy in renal transplant patients. Ther Drug Monit. (2009) 31:224–32. doi: 10.1097/FTD.0b013e31819c077c

PubMed Abstract | Crossref Full Text | Google Scholar

19. Durr M, Lachmann N, Zukunft B, Schmidt D, Budde K, and Brakemeier S. Late Conversion to Belatacept After Kidney Transplantation: Outcome and Prognostic Factors. Transplant Proc. (2017) 49:1747–56.e1. doi: 10.1016/j.transproceed.2017.05.013

PubMed Abstract | Crossref Full Text | Google Scholar

20. Bland JM and Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. (1986) 1:307–10.

Google Scholar

21. Passing H and Bablok. A new biometrical procedure for testing the equality of measurements from two different analytical methods. Application of linear regression procedures for method comparison studies in clinical chemistry, Part I. J Clin Chem Clin Biochem. (1983) 21:709–20. doi: 10.1515/cclm.1983.21.11.709

PubMed Abstract | Crossref Full Text | Google Scholar

22. Enokiya T, Nishikawa K, Muraki Y, Iwamoto T, Kanda H, Sugimura Y, et al. 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. J Pharm Health Care Sci. (2017) 3:17. doi: 10.1186/s40780-017-0086-7

PubMed Abstract | Crossref Full Text | Google Scholar

23. Alloway R, Steinberg S, Khalil K, Gourishankar S, Miller J, Norman D, et al. Conversion of stable kidney transplant recipients from a twice daily Prograf-based regimen to a once daily modified release tacrolimus-based regimen. Transplant Proc. (2005) 37:867–70. doi: 10.1016/j.transproceed.2004.12.222

PubMed Abstract | Crossref Full Text | Google Scholar

24. Kino T, Inamura N, Sakai F, Nakahara K, Goto T, Okuhara M, et al. Effect of FK-506 on human mixed lymphocyte reaction in vitro. Transplant Proc. (1987) 19:36–9.

Google Scholar

25. Langone A, Steinberg SM, Gedaly R, Chan LK, Shah T, Sethi KD, et al. Switching STudy of Kidney TRansplant PAtients with Tremor to LCP-TacrO (STRATO): an open-label, multicenter, prospective phase 3b study. Clin Transpl. (2015) 29:796–805. doi: 10.1111/ctr.12581

PubMed Abstract | Crossref Full Text | Google Scholar

26. Barraclough KA, Isbel NM, Johnson DW, Campbell SB, and Staatz CE. Once- versus twice-daily tacrolimus: are the formulations truly equivalent? Drugs. (2011) 71:1561–77. doi: 10.2165/11593890-000000000-00000

PubMed Abstract | Crossref Full Text | Google Scholar

27. van Hest R, Mathot R, Vulto A, Weimar W, and van Gelder T. Predicting the usefulness of therapeutic drug monitoring of mycophenolic acid: a computer simulation. Ther Drug Monit. (2005) 27:163–7. doi: 10.1097/01.ftd.0000158083.45954.97

PubMed Abstract | Crossref Full Text | Google Scholar

28. Staatz CE and Tett SE. Pharmacology and toxicology of mycophenolate in organ transplant recipients: an update. Arch Toxicol. (2014) 88:1351–89. doi: 10.1007/s00204-014-1247-1

PubMed Abstract | Crossref Full Text | Google Scholar

29. Lipsky JJ. Mycophenolate mofetil. Lancet. (1996) 348:1357–9. doi: 10.1016/S0140-6736(96)10310-X

PubMed Abstract | Crossref Full Text | Google Scholar

30. Wagner M, Earley AK, Webster AC, Schmid CH, Balk EM, and Uhlig K. Mycophenolic acid versus azathioprine as primary immunosuppression for kidney transplant recipients. Cochrane Database Syst Rev. (2015) 2015:CD007746. doi: 10.1002/14651858.CD007746.pub2

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: area under the curve, immunosuppressive drug, limited sampling strategy, prediction equations, renal transplantation

Citation: Minamisono K, Nakao T, Ohyama T, Nishida S, Sasaki H, Hirose T, Hotta K, Mieno M and Iwami D (2026) Establishment of accurate estimation equations for area under the concentration curves using simultaneous limited sampling for extended-release tacrolimus and mycophenolic acid in kidney transplant recipients. Front. Immunol. 17:1710261. doi: 10.3389/fimmu.2026.1710261

Received: 22 September 2025; Accepted: 12 January 2026; Revised: 10 December 2025;
Published: 03 February 2026.

Edited by:

Stanislaw Stepkowski, University of Toledo, United States

Reviewed by:

Nikola Stefanović, University of Niš, Serbia
Carla Galvez, Universidad de Chile, Chile

Copyright © 2026 Minamisono, Nakao, Ohyama, Nishida, Sasaki, Hirose, Hotta, Mieno and Iwami. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Daiki Iwami, aXdhbWlkYWlraUB5YmIubmUuanA=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.