The Population Pharmacokinetic Models of Tacrolimus in Chinese Adult Liver Transplantation Patients

External validation of population pharmacokinetic models of tacrolimus in Thai adult liver transplant recipients

OBJECTIVE

Several population pharmacokinetic models of tacrolimus in liver transplant patients were built, and their predictability was evaluated in their settings. However, the extrapolation in the prediction was unclear. This study aimed to evaluate the predictive performance of published tacrolimus models in adult liver transplant recipients using data from the Thai population as an external dataset.

METHODS

The selected published models were systematically searched and evaluated for their quality. The external dataset of patients who underwent the first liver transplant and received immediate-release tacrolimus was used to assess the predictive performance of each selected model. Trough concentrations between 3 and 6 months were retrospectively collected to evaluate the predictability of each model using prediction-based diagnostics, simulation-based diagnostics, and Bayesian forecasting.

RESULTS

Sixty-seven patients with 360 trough concentrations and eight selected published models were included in this study. None of the models met the predictive precision criteria in prediction-based diagnostics. Meanwhile, four published population pharmacokinetic models showed a normal distribution in NPDE testing. Regarding Bayesian forecasting, all models improved their forecasts with at least one prior information data point.

CONCLUSION

Bayesian forecasting is more accurate and precise than other testing methods for predicting drug concentrations. However, none of the evaluated models provides satisfactory predictive performance for generalization to Thai liver transplant patients. This underscores the need for future research to develop population PK models tailored to the Thai population. Such efforts should consider the inclusion of nonlinear pharmacokinetics and region-specific factors, including genetic variability, to improve model accuracy and applicability.

Development and validation of individualized tacrolimus dosing software for Chinese pediatric liver transplantation patients: a population pharmacokinetic approach

OBJECTIVE

We aim to describe the population pharmacokinetics (PPK) of tacrolimus in Chinese pediatric patients under 4 years old after liver transplantation and to develop individualized tacrolimus dosing software.

METHODS

A total of 663 blood concentrations from 85 patients aged 4.57 months to 3.97 years were collected in this study. PPK analysis was performed using a nonlinear mixed effects modeling approach with the software, Phoenix. Using C#, an individualized tacrolimus dosing software was created. The software was then used to predict the concentrations of another ten pediatric liver transplantation patients to verify the accuracy of said software. The predictive error (PE) and the absolute predictive error (APE) for each predicted time point were computed.

RESULTS

A one-compartment model with first-order elimination best fitted the data. The apparent volume of distribution (V/F) and apparent clearance (CL/F) were 198.65 L and 2.41 L/h. Postoperative days (POD), total bilirubin (TBIL), and the use of voriconazole significantly influenced tacrolimus apparent clearance. The incorporation of an increasing number of actual blood drug concentrations into the prediction resulted in a decrease in both PE (72%, 17%, 7%) and APE (87%, 53%, 26%).

CONCLUSIONS

A qualified PPK model of tacrolimus was developed in Chinese pediatric patients. The individualized tacrolimus dosing software could be used as a suitable tool for the personalization of tacrolimus dosing for pediatric patients after liver transplantation.

Improving prediction of tacrolimus concentration using a combination of population pharmacokinetic modeling and machine learning in chinese renal transplant recipients

Aims

The population pharmacokinetic (PPK) model-based machine learning (ML) approach offers a novel perspective on individual concentration prediction. This study aimed to establish a PPK-based ML model for predicting tacrolimus (TAC) concentrations in Chinese renal transplant recipients.

Methods

Conventional TAC monitoring data from 127 Chinese renal transplant patients were divided into training (80%) and testing (20%) datasets. A PPK model was developed using the training group data. ML models were then established based on individual pharmacokinetic data derived from the PPK basic model. The prediction performances of the PPK-based ML model and Bayesian forecasting approach were compared using data from the test group.

Results

The final PPK model, incorporating hematocrit and CYP3A5 genotypes as covariates, was successfully established. Individual predictions of TAC using the PPK basic model, postoperative date, CYP3A5 genotype, and hematocrit showed improved rankings in ML model construction. XGBoost, based on the TAC PPK, exhibited the best prediction performance.

Conclusion

The PPK-based machine learning approach emerges as a superior option for predicting TAC concentrations in Chinese renal transplant recipients.

Artificial Neural Network Analysis of Determinants of Tacrolimus Pharmacokinetics in Liver Transplant Recipients

BACKGROUND

The efficacy and toxicity of tacrolimus are closely related to its trough blood concentrations. Identifying the influencing factors of pharmacokinetics of tacrolimus in the early postoperative period is conducive to the optimization of the individualized tacrolimus administration protocol and to help liver transplant (LT) recipients achieve the target blood concentrations.

OBJECTIVE

This study aimed to develop an artificial neural network (ANN) for predicting the blood concentration of tacrolimus soon after liver transplantation and for identifying determinants of the concentration based on Shapley additive explanation (SHAP).

METHODS

In this retrospective study, we enrolled 31 recipients who were first treated with liver transplantation from the Department of Liver Transplantation and Hepatic Surgery, the First Affiliated Hospital of Shandong First Medical University (Shandong Provincial Qianfoshan Hospital) from November 2020 to May 2021. The basic information, biochemical indexes, use of concomitant drugs, and genetic factors of organ donors and recipients were used for the ANN model inputs, and the output was the steady-state trough concentration (C0) of tacrolimus after oral administration in LT recipients. The ANN model was established to predict C0 of tacrolimus, SHAP was applied to the trained model, and the SHAP value of each input was calculated to analyze quantitatively the influencing factors for the output C0.

RESULTS

A back-propagation ANN model with 3 hidden layers was established using deep learning. The mean prediction error was 0.27 ± 0.75 ng/mL; mean absolute error, 0.60 ± 0.52 ng/mL; correlation coefficient between predicted and actual C0 values, 0.9677; and absolute prediction error of all blood concentrations obtained by the ANN model, ≤3.0 ng/mL. The results indicated that the following factors had the most significant effect on C0: age, daily drug dose, genotype at CYP3A5 polymorphism rs776746 in both recipient and donor, and concomitant use of caspofungin. The predicted C0 value of tacrolimus in LT recipients increased in a dose-dependent manner when the daily dose exceeded 3 mg, whereas it decreased with age when LT recipients were older than 48 years. The predicted C0 was higher when recipients and donors had the genotype CYP3A5*3*3 than when they had the genotype CYP3A5*1. The predicted C0 value also increased with the use of caspofungin or Wuzhi capsule.

CONCLUSION AND RELEVANCE

The established ANN model can be used to predict the C0 value of tacrolimus in LT recipients with high accuracy and good predictive ability, serving as a reference for personalized treatment in the early stage after liver transplantation.

Modelling changes in the pharmacokinetics of tacrolimus during pregnancy after kidney transplantation: A retrospective cohort study

AIMS

Pregnancy after kidney transplantation is realistic but immunosuppressants should be continued to prevent rejection. Tacrolimus is safe during pregnancy and is routinely dosed based on whole-blood predose concentrations. However, maintaining these concentrations is complicated as physiological changes during pregnancy affect tacrolimus pharmacokinetics. The aim of this study was to describe tacrolimus pharmacokinetics throughout pregnancy and explain the changes by investigating covariates in a population pharmacokinetic model.

METHODS

Data of pregnant women using a twice-daily tacrolimus formulation following kidney transplantation were retrospectively collected from 6 months before conception, throughout gestation and up to 6 months postpartum. Pharmacokinetic analysis was performed using nonlinear mixed effects modelling. Demographic, clinical and genetic parameters were evaluated as covariates. The final model was evaluated using goodness-of-fit plots, visual predictive checks and a bootstrap analysis.

RESULTS

A total of 260 whole-blood tacrolimus predose concentrations from 14 pregnant kidney transplant recipients were included. Clearance increased during pregnancy from 34.5 to 41.7 L/h, by 15, 19 and 21% in the first, second and third trimester, respectively, compared to prior to pregnancy. This indicates a required increase in the tacrolimus dose by the same percentage to maintain the prepregnancy concentration. Haematocrit and gestational age were negatively correlated with tacrolimus clearance (P ≤ 0.01), explaining 18% of interindividual and 85% of interoccasion variability in oral clearance.

CONCLUSIONS

Tacrolimus clearance increases during pregnancy, resulting in decreased exposure to tacrolimus, which is explained by gestational age and haematocrit. To maintain prepregnancy target whole-blood tacrolimus predose concentrations during pregnancy, increasing the dose is required.

Population PK/PD model of tacrolimus for exploring the relationship between accumulated exposure and quantitative scores in myasthenia gravis patients

Tacrolimus is an important immunosuppressant used in the treatment of myasthenia gravis (MG). However, the population pharmacokinetic (PK) characteristics together with the exposure-response of tacrolimus in the treatment of MG remain largely unknown. In this study, we aimed to develop a population PK/pharmacodynamic (PK/PD) model of tacrolimus in patients with MG, in order to explore the relationships among tacrolimus dose, exposure, and its therapeutic efficacy. The genotype of CYP3A5, Osserman's classification, and status of thymus, as well as demographic characteristics and other biomarkers from laboratory testing were tested as covariate, and simulations were performed based on the final model. The population PK model was described using a one-compartment model with first-order elimination and fixed absorption parameters. CYP3A5 genotype significantly influenced the apparent clearance, and total protein (TP) influenced the apparent volume of distribution as covariates. The quantitative MG scores were characterized by the cumulated area under curve of tacrolimus in a maximum effect function. Osserman's classification was a significant covariate on the initial score of patients with MG. The simulations demonstrated that tacrolimus showed an unsatisfying effect possibly due to insufficient exposure in some patients with MG. A starting dose of 2 mg/d and even higher dose for patients with CYP3A5 *1/*1 and *1/*3 and lower TP level were required for the rapid action of tacrolimus. The population PK/PD model quantitatively described the relationships among tacrolimus dose, exposure, and therapeutic efficacy in patients with MG, which could provide reference for the optimization of tacrolimus dosing regimen at the individual patient level.

Population pharmacokinetic analyses of tacrolimus in non-transplant patients: a systematic review

BACKGROUND AND OBJECTIVES

Tacrolimus (TAC) has been increasingly used in patients with non-transplant settings. Because of its large between-subject variability, several population pharmacokinetic (PPK) studies have been performed to facilitate individualized therapy. This review summarized published PPK models of TAC in non-transplant patients, aiming to clarify factors affecting PKs of TAC and identify the knowledge gap that may require further research.

METHODS

The PubMed, Embase databases, and Cochrane Library, as well as related references, were searched from the time of inception of the databases to February 2023, to identify TAC population pharmacokinetic studies modeled in non-transplant patients using a non-linear mixed-effects modeling approach.

RESULTS

Sixteen studies, all from Asian countries (China and Korea), were included in this study. Of these studies, eleven and four were carried out in pediatric and adult patients, respectively. One-compartment models were the commonly used structural models for TAC. The apparent clearance (CL/F) of TAC ranged from 2.05 to 30.9 L·h^-1 (median of 14.9 L·h^-1). Coadministered medication, genetic factors, and weight were the most common covariates affecting TAC-CL/F, and variability in the apparent volume of distribution (V/F) was largely explained by weight. Coadministration with Wuzhi capsules reduced CL/F by about 19 to 43%. For patients with CYP3A5*1*1 and *1*3 genotypes, the CL/F was 39-149% higher CL/F than patients with CYP3A5*1*1.

CONCLUSION

The optimal TAC dosage should be adjusted based on the patient's co-administration, body weight, and genetic information (especially CYP3A5 genotype). Further studies are needed to assess the generalizability of the published models to other ethnic groups. Moreover, external validation should be frequently performed to improve the clinical practicality of the models.

Population pharmacokinetics of tacrolimus in Chinese adult liver transplant patients

Tacrolimus is widely used in organ transplantation to prevent rejection. However, the narrow therapeutic window and the large inter-and intra-individual variability in the pharmacokinetics (PK) of tacrolimus make it difficult for individualization of dosing. This study aimed at developing a population pharmacokinetic model for estimating the oral clearance of tacrolimus in Chinese liver transplant patients, and identifying factors that contribute to the PK variability of tacrolimus. Data of 151 liver transplant patients who received tacrolimus were analyzed in this study. The population PK model was analyzed and the covariates including population demographic and biochemical characteristics, drug combination, and genetic polymorphism were explored using non-linear mixed-effects modeling approach. A single-compartment population PK model was developed, and the final model was CL/F = (14.6-2.38 × cytochrome P450 (CYP) 3A5-3.72 × WZC+1.04 × (POD/9)+2.48 × COR) × Exp(η_i ), where CYP3A5 was 1 for CYP3A5*3/*3, Wuzhi Capsule (WZC) was 1 when patients took tacrolimus combined with WZC, otherwise it was 0, corticosteroids (COR) was 1 when patients take tacrolimus combined with COR, otherwise, it was 0, POD was the post-operative day. Visual inspection and bootstrap indicated that the final model was stable and robust. In this study, we developed the first tacrolimus population PK model in Chinese adult liver transplant patients. We first determined the influence of WZC on tacrolimus in these people, which could provide useful PK information for the drug combination of tacrolimus and WZC. We also revealed the influence of genetic polymorphism of CYP3A5, POD, and a combination of COR on tacrolimus PK. Therefore, these significant factors should be taken into consideration in optimizing dosage regimens.

Evaluation of the performance of a prior tacrolimus population pharmacokinetic kidney transplant model among adult allogeneic hematopoietic stem cell transplant patients

Abstract

Tacrolimus is a calcineurin inhibitor used to prevent acute graft versus host disease in adult patients receiving allogeneic hematopoietic stem cell transplantation (HCT). Previous population pharmacokinetic (PK) models have been developed in solid organ transplant, yet none exists for patients receiving HCT. The primary objectives of this study were to (1) use a previously published population PK model in adult patients who underwent kidney transplant and apply it to allogeneic HCT; (2) evaluate model‐predicted tacrolimus steady‐state trough concentrations and simulations in patients receiving HCT; and (3) evaluate covariates that affect tacrolimus PK in allogeneic HCT. A total of 252 adult patients receiving allogeneic HCT were included in the study. They received oral tacrolimus twice daily (0.03 mg/kg) starting 3 days prior to transplant. Data for these analyses included baseline clinical and demographic data, genotype data for single nucleotide polymorphisms in CYP3A4/5 and ABCB1, and the first tacrolimus steady‐state trough concentration. A dosing simulation strategy based on observed trough concentrations (rather than model‐based predictions) resulted in 12% more patients successfully achieving tacrolimus trough concentrations within the institutional target range (5–10 ng/ml). Stepwise covariate analyses identified HLA match and conditioning regimen (myeloablative vs. reduced intensity) as significant covariates. Ultimately, a previously published tacrolimus population PK model in kidney transplant provided a platform to help establish a model‐based dose adjustment strategy in patients receiving allogenic HCT, and identified HCT‐specific covariates to be considered for future prospective studies.

Impact of Donor and Recipient CYP3A5*3 Genotype on Tacrolimus Population Pharmacokinetics in Chinese Adult Liver Transplant Recipients

Background: Tacrolimus (TAC) is widely used after liver transplantation, but the therapeutic window is narrow. Objective: The purpose was to study both donor and recipient CYP3A5*3 genotypes affecting TAC apparent clearance rate (CL/F) and investigate a TAC population pharmacokinetic (PPK) model in Chinese liver transplant recipients for potential starting-dose individualized medication. Methods: A data set of 721 TAC concentrations was obtained from 43 adult liver transplant recipients. The TAC PPK model was analyzed using nonlinear mixed-effects modeling. Potential covariates, including demographic characteristics, physiological and pathological data, concomitant medications, and CYP3A5*3 genotype, were evaluated. The final model was validated using normalized prediction distribution errors and bootstrapping. Results: A 2-compartment model with first-order absorption and elimination was used to describe TAC disposition. Population estimates of TAC, CL/F, apparent central distribution volume (V2/F), rate of absorption (Ka), and apparent peripheral distribution volume (V3/F) were 18.1 L/h (12%), 72.7 L (34%), 0.163 h−1 (17%), and 412 L (21%), respectively. The model and estimated parameters were found to be stable. Other covariates did not influence TAC CL/F. Both donor and recipient CYP3A5*1 genotypes were significantly correlated with TAC clearance, and CL/F was 1.70-fold higher in both donor and recipient CYP3A5*1 carriers than in noncarriers among Chinese liver transplant recipients. Conclusion and Relevance: A PPK model of TAC was established in Chinese adult liver transplantation recipients for starting-dose individualized medication, which can be expanded to optimize clinical efficacy and minimize toxicity with therapeutic drug monitoring.

Toward a robust tool for pharmacokinetic-based personalization of treatment with tacrolimus in solid organ transplantation: A model-based meta-analysis approach

AIMS

The objective of this study is to develop a generic model for tacrolimus pharmacokinetics modelling using a meta-analysis approach, that could serve as a first step towards a prediction tool to inform pharmacokinetics-based optimal dosing of tacrolimus in different populations and indications.

METHODS

A systematic literature review was performed and a meta-model developed with NONMEM software using a top-down approach. Historical (previously published) data were used for model development and qualification. In-house individual rich and sparse tacrolimus blood concentration profiles from adult and paediatric kidney, liver, lung and heart transplant patients were used for model validation. Model validation was based on successful numerical convergence, adequate precision in parameter estimation, acceptable goodness of fit with respect to measured blood concentrations with no indication of bias, and acceptable performance of visual predictive checks. External validation was performed by fitting the model to independent data from 3 external cohorts and remaining previously published studies.

RESULTS

A total of 76 models were found relevant for meta-model building from the literature and the related parameters recorded. The meta-model developed using patient level data was structurally a 2-compartment model with first-order absorption, absorption lag time and first-time varying elimination. Population values for clearance, intercompartmental clearance, central and peripheral volume were 22.5 L/h, 24.2 L/h, 246.2 L and 109.9 L, respectively. The absorption first-order rate and the lag time were fixed to 3.37/h and 0.33 hours, respectively. Transplanted organ and time after transplantation were found to influence drug apparent clearance whereas body weight influenced both the apparent volume of distribution and the apparent clearance. The model displayed good results as regards the internal and external validation.

CONCLUSION

A meta-model was successfully developed for tacrolimus in solid organ transplantation that can be used as a basis for the prediction of concentrations in different groups of patients, and eventually for effective dose individualization in different subgroups of the population.

A population pharmacokinetic study of tacrolimus in healthy Chinese volunteers and liver transplant patients

AIM

To develop a population pharmacokinetic (PopPK) model of tacrolimus in healthy Chinese volunteers and liver transplant recipients for investigating the difference between the populations, and for potential individualized medication.

METHODS

A set of 1100 sparse trough concentration data points from 112 orthotopic liver transplant recipients, as well as 851 dense data points from 40 healthy volunteers receiving a single dose of tacrolimus (2 mg, p.o.) were collected. PopPK model of tacrolimus was constructed using the program NONMEM. Related covariates such as age, hepatic and renal functions that were potentially associated with tacrolimus disposition were evaluated. The final model was validated using bootstrapping and a visual predictive check.

RESULTS

A two-compartment model of tacrolimus could best describe the data from the two populations. The final model including two covariates, population (liver transplant recipients or volunteers) and serum ALT (alanine aminotransferase) level, was verified and adequately described the pharmacokinetic characteristics of tacrolimus. The estimates of V2/F, Q/F and V3/F were 22.7 L, 76.3 L/h and 916 L, respectively. The estimated CL/F in the volunteers and liver transplant recipients was 32.8 and 18.4 L/h, respectively. Serum ALT level was inversely related to CL/F, whereas age did not influence CL/F. Thus, the elderly (≥65 years) and adult (<65 years) groups in the liver transplant recipients showed no significant difference in the clearance of tacrolimus.

CONCLUSION

Compared with using the sparse data only, the integrating modeling technique combining sparse data from the patients and dense data from the healthy volunteers improved the PopPK analysis of tacrolimus.