Total plasma protein effect on tacrolimus elimination in kidney transplant patients – Population pharmacokinetic approach

Customizing Tacrolimus Dosing in Kidney Transplantation: Focus on Pharmacogenetics

ABSTRACT

Different polymorphisms in genes encoding metabolizing enzymes and drug transporters have been associated with tacrolimus pharmacokinetics. In particular, studies on CYP3A4 and CYP3A5, and their combined cluster have demonstrated their significance in adjusting tacrolimus dosing to minimize under- and overexposure thereby increasing the proportion of patients who achieve tacrolimus therapeutic target. Many factors influence the pharmacokinetics of tacrolimus, contributing to inter-patient variability affecting individual dosing requirements. On the other hand, the growing use of population pharmacokinetic models in solid organ transplantation, including different tacrolimus formulations, has facilitated the integration of pharmacogenetic data and other variables into algorithms to easier implement the personalized dose adjustment in transplant centers. The future of personalized medicine in transplantation lies in implementing these models in clinical practice, with pharmacogenetics as a key factor to account for the high inter-patient variability in tacrolimus exposure. To date, three clinical trials have validated the clinical application of these approaches. The aim of this review is to provide an overview of the current studies regarding the different population pharmacokinetic including pharmacogenetics and those translated to the clinical practice for individualizing tacrolimus dose adjustment in kidney transplantation.

Model-Informed Precision Dosing of Tacrolimus: A Systematic Review of Population Pharmacokinetic Models and a Benchmark Study of Software Tools

BACKGROUND AND OBJECTIVE

Tacrolimus is an immunosuppressant commonly administered after solid organ transplantation. It is characterized by a narrow therapeutic window and high variability in exposure, demanding personalized dosing. In recent years, population pharmacokinetic models have been suggested to guide model-informed precision dosing of tacrolimus. We aimed to provide a comprehensive overview of population pharmacokinetic models and model-informed precision dosing software modules of tacrolimus in all solid organ transplant settings, including a simulation-based investigation of the impact of covariates on exposure and target attainment.

METHODS

We performed a systematic literature search to identify population pharmacokinetic models of tacrolimus in solid organ transplant recipients. We integrated selected population pharmacokinetic models into an interactive software tool that allows dosing simulations, Bayesian forecasting, and investigation of the impact of covariates on exposure and target attainment. We conducted a web survey amongst model-informed precision dosing software tool providers and benchmarked publicly available tools in terms of models, target populations, and clinical integration.

RESULTS

We identified 80 population pharmacokinetic models, including 44 one-compartment and 36 two-compartment models. The most frequently retained covariates on clearance and distribution parameters were cytochrome P450 3A5 polymorphisms and body weight, respectively. Our simulation tool, hosted at https://lpmx.shinyapps.io/tacrolimus/ , allows thorough investigation of the impact of covariates on exposure and target attainment. We identified 15 model-informed precision dosing software tool providers, of which ten offer a tacrolimus solution and nine completed the survey.

CONCLUSIONS

Our work provides a comprehensive overview of the landscape of available tacrolimus population pharmacokinetic models and model-informed precision dosing software modules. Our simulation tool allows an interactive thorough exploration of covariates on exposure and target attainment.

A prospective controlled, randomized clinical trial of kidney transplant recipients developed personalized tacrolimus dosing using model-based Bayesian Prediction

For three decades, tacrolimus (Tac) dose adjustment in clinical practice has been calculated empirically according to the manufacturer's labeling based on a patient's body weight. Here, we developed and validated a Population pharmacokinetic (PPK) model including pharmacogenetics (cluster CYP3A4/CYP3A5), age, and hematocrit. Our study aimed to assess the clinical applicability of this PPK model in the achievement of Tac Co (therapeutic trough Tac concentration) compared to the manufacturer's labelling dosage. A prospective two-arm, randomized, clinical trial was conducted to determine Tac starting and subsequent dose adjustments in 90 kidney transplant recipients. Patients were randomized to a control group with Tac adjustment according to the manufacturer's labeling or the PPK group adjusted to reach target Co (6-10 ng/ml) after the first steady state (primary endpoint) using a Bayesian prediction model (NONMEM). A significantly higher percentage of patients from the PPK group (54.8%) compared with the control group (20.8%) achieved the therapeutic target fulfilling 30% of the established superiority margin defined. Patients receiving PPK showed significantly less intra-patient variability compared to the control group, reached the Tac Co target sooner (5 days vs 10 days), and required significantly fewer Tac dose modifications compared to the control group within 90 days following kidney transplant. No statistically significant differences occurred in clinical outcomes. Thus, PPK-based Tac dosing offers significant superiority for starting Tac prescription over classical labeling-based dosing according to the body weight, which may optimize Tac-based therapy in the first days following transplantation.

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.

Adaptation of a population pharmacokinetic model to inform tacrolimus therapy in heart transplant recipients

AIM

Existing tacrolimus population pharmacokinetic models are unsuitable for guiding tacrolimus dosing in heart transplant recipients. This study aimed to develop and evaluate a population pharmacokinetic model for tacrolimus in heart transplant recipients that considers the tacrolimus-azole antifungal interaction.

METHODS

Data from heart transplant recipients (n = 87) administered the oral immediate-release formulation of tacrolimus (Prograf®) were collected. Routine drug monitoring data, principally trough concentrations, were used for model building (n = 1099). A published tacrolimus model was used to inform the estimation of Ka , V2 /F, Q/F and V3 /F. The effect of concomitant azole antifungal use on tacrolimus CL/F was quantified. Fat-free mass was implemented as a covariate on CL/F, V2 /F, V3 /F and Q/F on an allometry scale. Subsequently, stepwise covariate modelling was performed. Significant covariates influencing tacrolimus CL/F were included in the final model. Robustness of the final model was confirmed using prediction-corrected visual predictive check (pcVPC). The final model was externally evaluated for prediction of tacrolimus concentrations of the fourth dosing occasion (n = 87) from one to three prior dosing occasions.

RESULTS

Concomitant azole antifungal therapy reduced tacrolimus CL/F by 80%. Haematocrit (∆OFV = -44, P < .001) was included in the final model. The pcVPC of the final model displayed good model adequacy. One recent drug concentration is sufficient for the model to guide tacrolimus dosing.

CONCLUSION

A population pharmacokinetic model that adequately describes tacrolimus pharmacokinetics in heart transplant recipients, considering the tacrolimus-azole antifungal interaction was developed. Prospective evaluation is required to assess its clinical utility to improve patient outcomes.

Model-Based Tacrolimus Follow-up Dosing in Adult Renal Transplant Recipients: A Simulation Trial

BACKGROUND

Initial algorithm-based dosing appears to be effective in predicting tacrolimus dose requirement. However, achieving and maintaining the target concentrations is challenging. Model-based follow-up dosing, which considers patient characteristics and pharmacological data, may further personalize treatment. This study investigated whether model-based follow-up dosing could lead to more accurate tacrolimus exposure than standard therapeutic drug monitoring (TDM) in kidney transplant recipients after an initial algorithm-based dose.

METHODS

This simulation trial included patients from a prospective trial that received an algorithm-based tacrolimus starting dose followed by TDM. For every measured tacrolimus predose concentration (C 0,obs ), model-based dosing advice was simulated using the InsightRX software. Based on previous tacrolimus doses and C 0 , age, body surface area, CYP3A4 and CYP3A5 genotypes, hematocrit, albumin, and creatinine, the optimal next dose, and corresponding tacrolimus concentration (C 0,pred ) were predicted.

RESULTS

Of 190 tacrolimus C 0 values measured in 59 patients, 121 (63.7%; 95% CI 56.8-70.5) C 0,obs were within the therapeutic range (7.5-12.5 ng/mL) versus 126 (66.3%, 95% CI 59.6-73.0) for C 0,pred ( P = 0.89). The median absolute difference between the tacrolimus C 0 and the target tacrolimus concentration (10.0 ng/mL) was 1.9 ng/mL for C 0,obs versus 1.6 ng/mL for C 0,pred . In a historical cohort of 114 kidney transplant recipients who received a body weight-based starting dose followed by TDM, 172 of 335 tacrolimus C 0 (51.3%) were within the therapeutic range (10.0-15.0 ng/mL).

CONCLUSIONS

The combination of an algorithm-based tacrolimus starting dose with model-based follow-up dosing has the potential to minimize under- and overexposure to tacrolimus in the early posttransplant phase, although the additional effect of model-based follow-up dosing on initial algorithm-based dosing seems small.

Important lack of difference in tacrolimus and mycophenolic acid pharmacokinetics between Aboriginal and Caucasian kidney transplant recipients

AIM

To examine whether differences in tacrolimus and mycophenolic acid (MPA) pharmacokinetics contribute to the poorer kidney transplant outcomes experienced by Aboriginal Australians.

METHODS

Concentration-time profiles for tacrolimus and MPA were prospectively collected from 43 kidney transplant recipients: 27 Aboriginal and 16 Caucasian. Apparent clearance (CL/F) and distribution volume (V/F) for each individual were derived from concentration-time profiles combined with population pharmacokinetic priors, with subsequent assessment for between-group difference in pharmacokinetics. In addition, population pharmacokinetic models were developed using the prospective dataset supplemented by previously developed structural models for tacrolimus and MPA. The change in NONMEM objective function was used to assess improvement in goodness of model fit.

RESULTS

No differences were found between Aboriginal and Caucasian groups or empirical Bayes estimates, for CL/F or V/F of MPA or tacrolimus. However, a higher prevalence of CYP3A5 expressers (26% compared with 0%) and wider between-subject variability in tacrolimus CL/F (SD = 5.00 compared with 3.25 L/h/70 kg) were observed in the Aboriginal group, though these differences failed to reach statistical significance (p = .07 and p = .08).

CONCLUSION

There were no differences in typical tacrolimus or MPA pharmacokinetics between Aboriginal and Caucasian kidney transplant recipients. This means that Bayesian dosing tools developed to optimise tacrolimus and MPA dosing in Caucasian recipients may be applied to Aboriginal recipients. In turn, this may improve drug exposure and thereby transplant outcomes in this group. Aboriginal recipients appeared to have greater between-subject variability in tacrolimus CL/F and a higher prevalence of CYP3A5 expressers, attributes that have been linked with inferior outcomes.

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.

Predicting model‐informed precision dosing: A test‐case in tacrolimus dose adaptation for kidney transplant recipients

Before investing resources into the development of a precision dosing (model‐informed precision dosing [MIPD]) tool for tacrolimus, the performance of the tool was evaluated in silico. A retrospective dataset of 315 de novo kidney transplant recipients was first used to identify a one‐compartment pharmacokinetic (PK) model with time‐dependent clearance. MIPD performance was subsequently evaluated by calculating errors to predict future concentrations, which is directly related to dosing precision and probability of target attainment (PTA). Based on the identified model residual error, the theoretical upper limit was 45% PTA for a target of 13.5 ng/ml and an acceptable range of 12–15 ng/ml. Using empirical Bayesian estimation, this limit was reached on day 5 post‐transplant and beyond. By incorporating correlated within‐patient variability when predicting future individual concentrations, PTA improved beyond the theoretical upper limit. This yielded a Bayesian feedback dosing algorithm accurately predicting future trough concentrations and adapting each dose to reach a target concentration. Simulated concentration‐time profiles were then used to quantify MIPD‐based improvement on three end points: average PTA increased from 28% to 39%, median time to three concentrations in target decreased from 10 to 7 days, and mean log‐squared distance to target decreased from 0.080 to 0.055. A study of 200 patients was predicted to have sufficient power to demonstrate these nuanced PK end points reliably. These simulations supported our decision to develop a precision dosing tool for tacrolimus and test it in a prospective trial.

Predictive Performance of Published Tacrolimus Population Pharmacokinetic Models in Thai Kidney Transplant Patients

BACKGROUND AND OBJECTIVE

Tacrolimus is a narrow therapeutic index drug with high pharmacokinetic variability, and several tacrolimus population pharmacokinetic (PopPK) models were developed to guide individualized drug dosing. These models, however, may not perform well in other clinical settings. Therefore, we aimed to assess the predictive ability of published tacrolimus PopPK models using a dataset of Thai kidney transplant patients.

METHODS

The external dataset was retrospectively collected from medical records of Bhumibol Adulyadej Hospital, Thailand. Published tacrolimus PopPK models were systematically searched from PubMed, Science Direct, CINAHL Complete, and Scopus databases. Models conducted using a nonlinear mixed-effects approach with covariate resemblance to our external dataset were selected. The external dataset consisted of Thai kidney transplant patients receiving oral immediate- or extended-release tacrolimus formulations twice or once daily, respectively. Accuracy and precision of predicted concentrations were evaluated using mean absolute prediction error (MAPE), root mean square error (RMSE), and goodness of fit plots.

RESULTS

Only three models produced acceptable population predictions with the MAPE of < 50%. By using the Bayesian posthoc estimate of individual pharmacokinetic parameters, all models well performed with the MAPE and RMSE of < 30% and 40%, respectively, except two models; one could not successfully converge and the other substantially underpredicted tacrolimus concentrations.

CONCLUSION

We evaluated ten tacrolimus PopPK models, and eight models resulted in satisfactorily individual predicted tacrolimus concentrations in Thai kidney transplant patients and may be used to aid tacrolimus dose adjustment along with a clinical judgment.

Machine learning-based method for tacrolimus dose predictions in Chinese kidney transplant perioperative patients

WHAT IS KNOWN AND OBJECTIVES

Tacrolimus (TAC), a first-line immunosuppressant in solid-organ transplant, has a narrow therapeutic window and large inter-individual variability, which affects its use in clinical practice. Successful predictions using machine learning algorithms have been reported in several fields. However, a comparison of 10 machine learning model-based TAC pharmacogenetic and pharmacokinetic dosing algorithms for kidney transplant perioperative patients of Chinese descent has not been reported. The objective of this study was to screen and establish an appropriate machine learning method to predict the individualized dosages of TAC for perioperative kidney transplant patients.

METHODS

The records of 2551 patients were collected from three transplant centres, 80% of which were randomly selected as a 'derivation cohort' to develop the dose prediction algorithm, while the remaining 20% constituted a 'validation cohort' to validate the final algorithm selected. Important features were screened according to our previously established population pharmacokinetic model of tacrolimus. The performances of the algorithms were evaluated and compared using R-squared and the mean percentage in the remaining 20% of patients.

RESULTS AND DISCUSSION

This study identified several factors influencing TAC dosage, including CYP3A5 rs776746, CYP3A4 rs4646437, haematocrit, Wuzhi capsules, TAC daily dose, age, height, weight, post-operative time, nifedipine and the medication history of the patient. According to our results, among the 10 machine learning models, the extra trees regressor (ETR) algorithm showed the best performance in the training set (R-squared: 1, mean percentage within 20%: 100%) and test set (R-squared: 0.85, mean percentage within 20%: 92.77%) of the derivation cohort. The ETR model successfully predicted the ideal TAC dosage in 97.73% of patients, especially in the intermediate dosage range (>5 mg/day to <8 mg/day), whereby the ideal TAC dosage could be successfully predicted in 99% of the patients.

WHAT IS NEW AND CONCLUSION

The results indicated that the ETR algorithm, which was chosen to establish the dose prediction model, performed better than the other nine machine learning models. This study is the first to establish ETR algorithms to predict TAC dosage. This study will further promote the individualized medication of TAC in kidney transplant patients in the future, which has great significance in ensuring the safety and effectiveness of drug use.

Population Pharmacokinetic Models of Tacrolimus in Adult Transplant Recipients: A Systematic Review

BACKGROUND AND OBJECTIVES

Numerous population pharmacokinetic (PK) models of tacrolimus in adult transplant recipients have been published to characterize tacrolimus PK and facilitate dose individualization. This study aimed to (1) investigate clinical determinants influencing tacrolimus PK, and (2) identify areas requiring additional research to facilitate the use of population PK models to guide tacrolimus dosing decisions.

METHODS

The MEDLINE and EMBASE databases, as well as the reference lists of all articles, were searched to identify population PK models of tacrolimus developed from adult transplant recipients published from the inception of the databases to 29 February 2020.

RESULTS

Of the 69 studies identified, 55% were developed from kidney transplant recipients and 30% from liver transplant recipients. Most studies (91%) investigated the oral immediate-release formulation of tacrolimus. Few studies (17%) explained the effect of drug-drug interactions on tacrolimus PK. Only 35% of the studies performed an external evaluation to assess the generalizability of the models. Studies related variability in tacrolimus whole blood clearance among transplant recipients to either cytochrome P450 (CYP) 3A5 genotype (41%), days post-transplant (30%), or hematocrit (29%). Variability in the central volume of distribution was mainly explained by body weight (20% of studies).

CONCLUSION

The effect of clinically significant drug-drug interactions and different formulations and brands of tacrolimus should be considered for any future tacrolimus population PK model development. Further work is required to assess the generalizability of existing models and identify key factors that influence both initial and maintenance doses of tacrolimus, particularly in heart and lung transplant recipients.

Avoiding Tacrolimus Underexposure and Overexposure with a Dosing Algorithm for Renal Transplant Recipients: A Single Arm Prospective Intervention Trial

Bodyweight‐based tacrolimus dosing followed by therapeutic drug monitoring is standard clinical care after renal transplantation. However, after transplantation, a meager 38% of patients are on target at first steady‐state and it can take up to 3 weeks to reach the target tacrolimus predose concentration (C₀). Tacrolimus underexposure and overexposure is associated with an increased risk of rejection and drug‐related toxicity, respectively. To minimize subtherapeutic and supratherapeutic tacrolimus exposure in the immediate post‐transplant phase, a previously developed dosing algorithm to predict an individual’s tacrolimus starting dose was tested prospectively. In this single‐arm, prospective, therapeutic intervention trial, 60 de novo kidney transplant recipients received a tacrolimus starting dose based on a dosing algorithm instead of a standard, bodyweight‐based dose. The algorithm included cytochrome P450 (CYP)3A4 and CYP3A5 genotype, body surface area, and age as covariates. The target tacrolimus C₀, measured for the first time at day 3, was 7.5–12.5 ng/mL. Between February 23, 2019, and July 7, 2020, 60 patients were included. One patient was excluded because of a protocol violation. On day 3 post‐transplantation, 34 of 59 patients (58%, 90% CI 47–68%) had a tacrolimus C₀ within the therapeutic range. Markedly subtherapeutic (< 5.0 ng/mL) and supratherapeutic (> 20 ng/mL) tacrolimus concentrations were observed in 7% and 3% of the patients, respectively. Biopsy‐proven acute rejection occurred in three patients (5%). In conclusion, algorithm‐based tacrolimus dosing leads to the achievement of the tacrolimus target C₀ in as many as 58% of the patients on day 3 after kidney transplantation.

Concept and utility of population pharmacokinetic and pharmacokinetic/pharmacodynamic models in drug development and clinical practice

Due to frequent clinical trial failures and consequently fewer new drug approvals, the need for improvement in drug development has, to a certain extent, been met using model-based drug development. Pharmacometrics is a part of pharmacology that quantifies drug behaviour, treatment response and disease progression based on different models (pharmacokinetic - PK, pharmacodynamic - PD, PK/PD models, etc.) and simulations. Regulatory bodies (European Medicines Agency, Food and Drug Administration) encourage the use of modelling and simulations to facilitate decision-making throughout all drug development phases. Moreover, the identification of factors that contribute to variability provides a basis for dose individualisation in routine clinical practice. This review summarises current knowledge regarding the application of pharmacometrics in drug development and clinical practice with emphasis on the population modelling approach.

Predictors of tacrolimus pharmacokinetic variability: current evidences and future perspectives

INTRODUCTION

In kidney transplantation, tacrolimus (TAC) is at the cornerstone of current immunosuppressive strategies. Though because of its narrow therapeutic index, it is critical to ensure that TAC levels are maintained within this sharp window through reactive adjustments. This would allow maximizing efficiency while limiting drug-associated toxicity. However, TAC high intra- and inter-patient pharmacokinetic (PK) variability makes it more laborious to accurately predict the appropriate dosage required for a given patient.

AREAS COVERED

This review summarizes the state-of-the-art knowledge regarding drug interactions, demographic and pharmacogenetics factors as predictors of TAC PK. We provide a scoring index for each association to grade its relevance and we present practical recommendations, when possible for clinical practice.

EXPERT OPINION

The management of TAC concentration in transplanted kidney patients is as critical as it is challenging. Recommendations based on rigorous scientific evidences are lacking as knowledge of potential predictors remains limited outside of DDIs. Awareness of these limitations should pave the way for studies looking at demographic and pharmacogenetic factors as well as gut microbiota composition in order to promote tailored treatment plans. Therapeutic approaches considering patients' clinical singularities may help allowing to maintain appropriate concentration of TAC.

Effects of red blood cell concentrate transfusion on blood tacrolimus concentration

Background Elevated blood concentration of tacrolimus is frequently observed following transfusion of red blood cell concentrate in patients after allogeneic hematopoietic stem cell transplantation. Objective The aim of this retrospective study was to clarify the effects of transfusion of red blood cell concentrate on the blood concentration of tacrolimus. Setting Chiba University Hospital in Japan. Method Fifty-two patients (aged 0-65 years) receiving both tacrolimus and transfusion after allogeneic hematopoietic stem cell transplantation were enrolled. The ratio of measurement after transfusion to measurement before transfusion was calculated for hematocrit and blood concentration/dose ratio of tacrolimus (termed the hematocrit ratio and the tacrolimus ratio, respectively). Main outcome measure Change in blood concentration/dose ratio of tacrolimus and variable factors associated with variation in tacrolimus ratio. Results The blood concentration/dose ratio of tacrolimus was increased after transfusion compared with before transfusion (p < 0.001). A statistically significant correlation was seen between the hematocrit ratio and tacrolimus ratio (r = 0.32, p < 0.001). Hematocrit ratio, age or body surface area, and difference in aspartate aminotransferase level before and after transfusion were associated with the variation in tacrolimus ratio. There was no correlation between tacrolimus ratio and change in serum creatinine or potassium level in the short term. Conclusion Change in the blood concentration/dose ratio of tacrolimus was associated with change in the hematocrit ratio after transfusion, and more attention is required for children or patients with small body surface area. Dose adjustment of tacrolimus is required if the blood concentration of tacrolimus is much higher than the target concentration.

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.

Population Pharmacokinetic Analysis of Immediate-Release Oral Tacrolimus Co-administered with Mycophenolate Mofetil in Corticosteroid-Free Adult Kidney Transplant Recipients

BACKGROUND AND OBJECTIVE

Tacrolimus is the mainstay calcineurin inhibitor frequently administered with mycophenolic acid with or without corticosteroids to prevent graft rejection in adult kidney transplant recipients. The primary objective of this study was to develop and evaluate a population pharmacokinetic model characterizing immediate-release oral tacrolimus co-administered with mycophenolate mofetil (a pro-drug of mycophenolic acid) in adult kidney transplant recipients on corticosteroid-free regimens. The secondary objective was to investigate the effects of clinical covariates on the pharmacokinetics of tacrolimus, emphasizing the interacting effects of mycophenolic acid.

METHODS

Population modeling and evaluation were conducted with Monolix (Suite-2018R1) using the stochastic approximation expectation-maximization algorithm in 49 adult subjects (a total of 320 tacrolimus whole-blood concentrations). Effects of clinical variables on tacrolimus pharmacokinetics were determined by population covariate modeling, regression modeling, and categorical analyses.

RESULTS

A two-compartment, first-order absorption with a lag-time, linear elimination, and constant error model best represented the population pharmacokinetics of tacrolimus. The apparent clearance value for tacrolimus was 17.9 l/h (6.95% relative standard error) in our model, which is lower compared with similar subjects on corticosteroid-based therapy. The glomerular filtration rate had significant effects on the apparent clearance and central compartment volume of distribution. Conversely, mycophenolic acid did not affect the apparent clearance of tacrolimus.

CONCLUSION

We have developed and internally evaluated a novel population pharmacokinetic model for tacrolimus co-administered with mycophenolate mofetil in corticosteroid-free adult kidney transplant patients. These findings are clinically important and provide further reasons for conducting therapeutic drug monitoring in this specific population.

Exploring Sirolimus Pharmacokinetic Variability Using Data Available from the Routine Clinical Care of Renal Transplant Patients - Population Pharmacokinetic Approach

Background

Due to wide intra- and inter-individual pharmacokinetic variability and narrow therapeutic index of sirolimus, the therapeutic drug monitoring (TDM) of sirolimus with detailed biochemical and clinical monitoring is necessary for dose individualization in kidney transplant patients. The purpose of the study was to explore and identify factors that contribute to pharmacokinetic variability by developing and validating a population model using routine TDM data and routinely monitored biochemical and clinical parameters.

Methods

The data obtained by routine monitoring of 38 patients over a period of one year from the sirolimus treatment initiation, were collected from patients' records. Population analysis was performed using the software NONMEM^®. The validity of the model was tested by the internal and external validation techniques.

Results

The pharmacokinetic variability was partially explained with patient's age and liver function. CL/F was found to decrease with age. According to the developed model, sirolimus CL/F decreases by, in average, 37% in patients with aspartate aminotransferase (AST) greater than 37 IU/L. The internal and external validation confirmed the satisfactory prediction of the developed model.

Conclusions

The population modeling of routinely monitored data allowed quantification of the age and liver function influence on sirolimus CL/F. According to the final model, patients with compromised liver function expressed via AST values require careful monitoring and dosing adjustments. Proven good predictive performance makes this model a useful tool in everyday clinical practice.

A population pharmacokinetic model to predict the individual starting dose of tacrolimus in adult renal transplant recipients

Aims

The aims of this study were to describe the pharmacokinetics of tacrolimus immediately after kidney transplantation, and to develop a clinical tool for selecting the best starting dose for each patient.

Methods

Data on tacrolimus exposure were collected for the first 3 months following renal transplantation. A population pharmacokinetic analysis was conducted using nonlinear mixed‐effects modelling. Demographic, clinical and genetic parameters were evaluated as covariates.

Results

A total of 4527 tacrolimus blood samples collected from 337 kidney transplant recipients were available. Data were best described using a two‐compartment model. The mean absorption rate was 3.6 h⁻¹, clearance was 23.0 l h^–1 (39% interindividual variability, IIV), central volume of distribution was 692 l (49% IIV) and the peripheral volume of distribution 5340 l (53% IIV). Interoccasion variability was added to clearance (14%). Higher body surface area (BSA), lower serum creatinine, younger age, higher albumin and lower haematocrit levels were identified as covariates enhancing tacrolimus clearance. Cytochrome P450 (CYP) 3A5 expressers had a significantly higher tacrolimus clearance (160%), whereas CYP3A4*22 carriers had a significantly lower clearance (80%). From these significant covariates, age, BSA, CYP3A4 and CYP3A5 genotype were incorporated in a second model to individualize the tacrolimus starting dose: Dosemg=222nghml–1*22.5lh–1*1.0ifCYP3A5*3/*3or1.62ifCYP3A5*1/*3orCYP3A5*1/*1*1.0ifCYP3A4*1or unknownor0.814ifCYP3A4*22*Age56−0.50*BSA1.930.72/1000Both models were successfully internally and externally validated. A clinical trial was simulated to demonstrate the added value of the starting dose model.

Conclusions

For a good prediction of tacrolimus pharmacokinetics, age, BSA, CYP3A4 and CYP3A5 genotype are important covariates. These covariates explained 30% of the variability in CL/F. The model proved effective in calculating the optimal tacrolimus dose based on these parameters and can be used to individualize the tacrolimus dose in the early period after transplantation.

Impact of POR and CYP3A5 Polymorphisms on Trough Concentration to Dose Ratio of Tacrolimus in the Early Post-operative Period Following Kidney Transplantation

BACKGROUND

Tacrolimus, a critical dose drug, is widely used in transplantation. Knowing the contribution of genetic factors, which significantly influence tacrolimus variability, is beneficial in the personalization of its starting dose. The significant impact of CYP3A5*3 polymorphisms on tacrolimus exposure has been reported. Conflicting results of the additional influence of POR*28 polymorphisms on tacrolimus pharmacokinetic interindividual variability have been observed among different populations. The objective of this study was to explore the interaction between POR*28 and CYP3A5*3 polymorphisms and their main effects on tacrolimus trough concentration to dose ratios on day 7 after kidney transplantation.

METHODS

Two hundred sixteen adult kidney transplant recipients participated in this retrospective study. All participants received a twice daily tacrolimus regimen. Blood samples and data were collected on day 7 after transplantation. A 2-way analysis of covariance was performed. Tested covariates were age, hemoglobin, serum albumin, and prednisolone dose.

RESULTS

A 2 × 2 analysis of covariance revealed that the interaction between CYP3A5 polymorphisms (CYP3A5 expresser and CYP3A5 nonexpresser) and POR polymorphisms (POR*28 carrier and POR*28 noncarrier) was not significant (F(1, 209) = 2.473, P = 0.117, (Equation is included in full-text article.)= 0.012). The predicted main effect of CYP3A5 and POR polymorphisms was significant (F(1, 209) = 105.565, P < 0.001, (Equation is included in full-text article.)= 0.336 and F(1, 209) = 4.007, P = 0.047, (Equation is included in full-text article.)= 0.019, respectively). Hemoglobin, age, and steroid dose influenced log C0/dose of tacrolimus (F(1, 209) = 20.612, P < 0.001, (Equation is included in full-text article.)= 0.090; F(1, 209) = 14.360, P < 0.001, (Equation is included in full-text article.)= 0.064; and F(1, 209) = 5.512, P = 0.020, (Equation is included in full-text article.)= 0.026, respectively).

CONCLUSIONS

After adjusting for the influences of hemoglobin, age, and prednisolone dose, significant impacts of the CYP3A5 and POR polymorphisms on tacrolimus exposure were found. The effect of POR*28 and CYP3A5*3 polymorphisms during the very early period after kidney transplantation is independent of each other.

Evaluating tacrolimus pharmacokinetic models in adult renal transplant recipients with different CYP3A5 genotypes

PURPOSE

Numerous studies have been conducted on the population pharmacokinetics of tacrolimus in adult renal transplant recipients. It has been reported that the cytochrome P450 (CYP) 3A5 genotype is an important cause of variability in tacrolimus pharmacokinetics. However, the predictive performance of population pharmacokinetic (PK) models of tacrolimus should be evaluated prior to their implementation in clinical practice. The aim of the study reported here was to test the predictive performance of these published PK models of tacrolimus.

METHODS

A literature search of the PubMed and Web of Science databases ultimately led to the inclusion of eight one-compartment models in our analysis. We collected a total of 1715 trough concentrations from 174 patients. Predictive performance was assessed based on visual and numerical comparison bias and imprecision and by the use of simulation-based diagnostics and Bayesian forecasting.

RESULTS

Of the eight one-compartment models assessed, seven showed better predictive performance in CYP3A5 extensive metabolizers in terms of bias and imprecision. Results of the simulation-based diagnostics also supported the findings. The model based on a Chinese population in 2013 (model 3) showed the best and most stable predictive performance in all the tests and was more informative in CYP3A5 extensive metabolizers. As expected, Bayesian forecasting improved model predictability. Diversity among models and between different CYP3A5 genotypes of the same model was also narrowed by Bayesian forecasting.

CONCLUSIONS

Based on our results, we recommend using model 3 in CYP3A5 extensive metabolizers in clinical practice. All models had a poor predictive performance in CYP3A5 poor metabolizers, and they should be used with caution in this patient population. However, Bayesian forecasting improved the predictability and reduced differences, and thus the models could be applied in this latter patient population for the design of maintenance dose.

Tacrolimus Population Pharmacokinetics and Multiple CYP3A5 Genotypes in Black and White Renal Transplant Recipients

Tacrolimus exhibits inter-patient pharmacokinetic variability attributed to CYP3A5 isoenzymes and the efflux transporter, P-glycoprotein. Most black renal transplant recipients require higher tacrolimus doses compared to whites to achieve similar troughs when race-adjusted recommendations are used. An established guideline provides tacrolimus genotype dosing recommendations based on CYP3A5*1(W/T) and loss of protein function variants: CYP3A5*3 (rs776746), CYP3A5*6 (rs10264272), CYP3A5*7 (rs41303343) and may provide more comprehensive race-adjusted dosing recommendations. Our objective was to develop a tacrolimus population pharmacokinetic model evaluating demographic, clinical, and genomic factors in stable black and white renal transplant recipients. A secondary objective investigated race-based tacrolimus regimens and genotype-specific dosing. Sixty-seven recipients receiving oral tacrolimus and mycophenolic acid ≥6 months completed a 12-hour pharmacokinetic study. CYP3A5*3,*6,*7 and ABCB1 1236C>T, 2677G>T/A, 3435C>T polymorphisms were characterized. Patients were classified as extensive, intermediate, and poor metabolizers using a novel CYP3A5*3*6*7 metabolic composite. Modeling and simulation was performed with computer software (NONMEM 7.3, ICON Development Solutions; Ellicott City, Maryland). A 2-compartment model with first-order elimination and absorption with lag time best described the data. The CYP3A5*3*6*7 metabolic composite was significantly associated with tacrolimus clearance (P value < .05), which was faster in extensive (mean: 45.0 L/hr) and intermediate (29.5 L/hr) metabolizers than poor metabolizers (19.8 L/hr). Simulations support CYP3A5*3*6*7 genotype-based tacrolimus dosing to enhance general race-adjusted regimens, with dose increases of 1.5-fold and 2-fold, respectively, in intermediate and extensive metabolizers for comparable exposures to poor metabolizers. This model offers a novel approach to determine tacrolimus dosing adjustments that maintain comparable therapeutic exposure between black and white recipients with different CYP3A5 genotypes.

External evaluation of published population pharmacokinetic models of tacrolimus in adult renal transplant recipients

AIM

Several tacrolimus population pharmacokinetic models in adult renal transplant recipients have been established to facilitate dose individualization. However, their applicability when extrapolated to other clinical centres is not clear. This study aimed to (1) evaluate model external predictability and (2) analyze potential influencing factors.

METHODS

Published models were screened from the literature and were evaluated using an external dataset with 52 patients (609 trough samples) collected by postoperative day 90 via methods that included (1) prediction-based prediction error (PE%), (2) simulation-based prediction- and variability-corrected visual predictive check (pvcVPC) and normalized prediction distribution error (NPDE) tests and (3) Bayesian forecasting to assess the influence of prior observations on model predictability. The factors influencing model predictability, particularly the impact of structural models, were evaluated.

RESULTS

Sixteen published models were evaluated. In prediction-based diagnostics, the PE% within ±30% was less than 50% in all models, indicating unsatisfactory predictability. In simulation-based diagnostics, both the pvcVPC and the NPDE indicated model misspecification. Bayesian forecasting improved model predictability significantly with prior 2-3 observations. The various factors influencing model extrapolation included bioassays, the covariates involved (CYP3A5*3 polymorphism, postoperative time and haematocrit) and whether non-linear kinetics were used.

CONCLUSIONS

The published models were unsatisfactory in prediction- and simulation-based diagnostics, thus inappropriate for direct extrapolation correspondingly. However Bayesian forecasting could improve the predictability considerably with priors. The incorporation of non-linear pharmacokinetics in modelling might be a promising approach to improving model predictability.