Pharmacokinetic models to assist the prescriber in choosing the best tacrolimus dose

Tacrolimus-Induced Neurotoxicity After Transplant: A Literature Review

Tacrolimus, a calcineurin inhibitor, is an immunosuppressant used globally to prevent rejection after organ transplantation. Although it significantly improves outcomes for solid organ transplant patients, it is associated with various side effects such as nephrotoxicity and neurotoxicity. Tacrolimus-induced neurotoxicity is frequently encountered in clinical practice and can present with a variety of symptoms that may occur even at therapeutic levels. Although tacrolimus-induced neurotoxicity is well documented, there is limited literature available on pharmacologic management. Twenty-eight case reports of tacrolimus-induced neurotoxicity were identified and analyzed in addition to other literature including reviews, retrospective studies, and animal model studies. The severity of cases of tacrolimus-induced neurotoxicity reported ranged from mild symptoms that could be managed with symptomatic treatment to conditions such as posterior reversible encephalopathy syndrome and chronic inflammatory demyelinating polyradiculoneuropathy that may require more immediate intervention. This information was utilized in addition to clinical experience to compile potential management options for prevention and treatment of neurotoxic adverse events. This review is limited by the utilization of primarily retrospective studies and case reports. The available literature on the subject is largely narrative and there are no guidelines on treatment of tacrolimus-induced neurotoxicity at the time of this research. This comprehensive review may guide further studies to investigate the pathophysiology of tacrolimus-induced neurotoxicity and to define patient-specific strategies for mitigation or minimization of neurotoxicity. This is especially important given that management of tacrolimus-induced neurotoxicity can include changes to immunosuppression that can result in an increased risk of rejection.

Population pharmacokinetic models of tacrolimus in paediatric solid organ transplant recipients: A systematic review

AIMS

This study aimed to provide up-to-date information on paediatric population pharmacokinetic models of tacrolimus and to identify factors influencing tacrolimus pharmacokinetic variability.

METHODS

Systematic searches in the Web of Science, PubMed, Scopus, Science Direct, Cochrane, EMBASE databases and reference lists of articles were conducted from inception to March 2023. All population pharmacokinetic studies of tacrolimus using nonlinear mixed-effect modelling in paediatric solid organ transplant patients were included.

RESULTS

Of the 21 studies reviewed, 62% developed from liver transplant recipients and 33% from kidney transplant recipients. Most studies used a 1-compartment model to describe tacrolimus pharmacokinetics. Body weight was a significant predictor for tacrolimus volume of distribution (Vd/F). The estimated Vd/F for 1-compartment models ranged from 20 to 1890 L, whereas the peripheral volume of distribution (Vp/F) for 2-compartment models was between 290 and 1520 L. Body weight, days post-transplant, CYP3A5 genotype or haematocrit were frequently reported as significant predictors of tacrolimus clearance. The estimated apparent clearance values range between 0.12 and 2.18 L/h/kg, with inter-individual variability from 13.5 to 110.0%. Only 29% of the studies assessed the generalizability of the models with external validation.

CONCLUSION

This review highlights the potential factors, modelling approaches and validation methods that impact tacrolimus pharmacokinetics in a paediatric population. The clinician could predict tacrolimus clearance based on body weight, CYP3A5 genotype, days post-transplant or haematocrit. Further research is required to determine the relationship between pharmacogenetics and tacrolimus pharmacodynamics in paediatric patients and confirm the applicability of nonlinear kinetics in this population.

CYP3A5*3 and CYP3A4*22 Cluster Polymorphism Effects on LCP-Tac Tacrolimus Exposure: Population Pharmacokinetic Approach

The aim of the study is to develop a population pharmacokinetic (PopPK) model and to investigate the influence of CYP3A5/CYP3A4 and ABCB1 single nucleotide polymorphisms (SNPs) on the Tacrolimus PK parameters after LCP-Tac formulation in stable adult renal transplant patients. The model was developed, using NONMEM v7.5, from full PK profiles from a clinical study (n = 30) and trough concentrations (C0) from patient follow-up (n = 68). The PK profile of the LCP-Tac formulation was best described by a two-compartment model with linear elimination, parameterized in elimination (CL/F) and distributional (CLD/F) clearances and central compartment (Vc/F) and peripheral compartment (Vp/F) distribution volumes. A time-lagged first-order absorption process was characterized using transit compartment models. According to the structural part of the base model, the LCP-Tac showed an absorption profile characterized by two transit compartments and a mean transit time of 3.02 h. Inter-individual variability was associated with CL/F, Vc/F, and Vp/F. Adding inter-occasion variability (IOV) on CL/F caused a statistically significant reduction in the model minimum objective function MOFV (p < 0.001). Genetic polymorphism of CYP3A5 and a cluster of CYP3A4/A5 SNPs statistically significantly influenced Tac CL/F. In conclusion, a PopPK model was successfully developed for LCP-Tac formulation in stable renal transplant patients. CYP3A4/A5 SNPs as a combined cluster including three different phenotypes (high, intermediate, and poor metabolizers) was the most powerful covariate to describe part of the inter-individual variability associated with apparent elimination clearance. Considering this covariate in the initial dose estimation and during the therapeutic drug monitoring (TDM) would probably optimize Tac exposure attainments.

Bridging the Worlds of Pharmacometrics and Machine Learning

Precision medicine requires individualized modeling of disease and drug dynamics, with machine learning-based computational techniques gaining increasing popularity. The complexity of either field, however, makes current pharmacological problems opaque to machine learning practitioners, and state-of-the-art machine learning methods inaccessible to pharmacometricians. To help bridge the two worlds, we provide an introduction to current problems and techniques in pharmacometrics that ranges from pharmacokinetic and pharmacodynamic modeling to pharmacometric simulations, model-informed precision dosing, and systems pharmacology, and review some of the machine learning approaches to address them. We hope this would facilitate collaboration between experts, with complementary strengths of principled pharmacometric modeling and flexibility of machine learning leading to synergistic effects in pharmacological applications.

Comparing the Performance of Multiple Small-Data Personalized Tacrolimus Dosing Models for Pediatric Liver Transplant: A Retrospective Study

Tacrolimus is a potent immunosuppressant used after pediatric liver transplant. However, tacrolimus's narrow therapeutic window, reliance on physicians' experience for the dose titration, and intra- and inter-patient variability result in liver transplant patients falling out of the target tacrolimus trough levels frequently. Existing personalized dosing models based on the area-under-the-concentration over time curves require a higher frequency of blood draws than the current standard of care and may not be practically feasible. We present a small-data artificial intelligence-derived platform, CURATE.AI, that uses data from individual patients obtained once daily to model the dose and response relationship and identify suitable doses dynamically. Retrospective optimization using 6 models of CURATE.AI and data from 16 patients demonstrated good predictive performance and identified a suitable model for further investigations.Clinical Relevance- This study established and compared the predictive performance of 6 personalized tacrolimus dosing models for pediatric liver transplant patients and identified a suitable model with consistently good predictive performance based on data from pediatric liver transplant patients.

Does the Tacrolimus Trough Level Adequately Predict Drug Exposure in Patients Requiring a High Tacrolimus Dose?

Tacrolimus (Tac) has a narrow therapeutic range. Dosing is generally targeted at Tac trough levels (C ₀), notwithstanding conflicting reports on the correlation between Tac C ₀ and systemic exposure measured by the area-under-the-concentration-over-time curve (AUC). The Tac dose required to meet the target C ₀ varies highly among patients. We hypothesized that patients requiring a relatively high Tac dose for a certain C ₀ may show a higher AUC.

Methods

We retrospectively analyzed data from 53 patients in which a 24-h Tac AUC₂₄ estimation was performed at our center. Patients were divided into those taking a low (≤0.15 mg/kg) or high (>0.15 mg/kg) once-daily Tac dose. Multiple linear regression models were used to investigate if the association between C ₀ and AUC₂₄ changes according to dose level.

Results

Despite the large difference in mean Tac dose between the low- and high-dose group (7 versus 17 mg/d), C ₀ levels were similar. However, the mean AUC₂₄ was substantially higher in the high-dose group (320 ± 96 h·μg/L versus 255 ± 81 h·μg/L, P < 0.001). This difference remained significant after adjusting for age and race. For a same C ₀, every 0.01 mg/kg increase in Tac dose resulted in an AUC₂₄ increase of 3.59 h·μg/L.

Conclusions

This study challenges the general belief that C ₀ levels are sufficiently reliable to estimate systemic drug exposure. We demonstrated that patients requiring a relatively high Tac dose to attain therapeutic C ₀ levels have higher drug exposure and could therefore potentially be overdosed.

A Hybrid Algorithm Combining Population Pharmacokinetic and Machine Learning for Isavuconazole Exposure Prediction

OBJECTIVES

Maximum a posteriori Bayesian estimation (MAP-BE) based on a limited sampling strategy and a population pharmacokinetic (POPPK) model is used to estimate individual pharmacokinetic parameters. Recently, we proposed a methodology that combined population pharmacokinetic and machine learning (ML) to decrease the bias and imprecision in individual iohexol clearance prediction. The aim of this study was to confirm the previous results by developing a hybrid algorithm combining POPPK, MAP-BE and ML that accurately predicts isavuconazole clearance.

METHODS

A total of 1727 isavuconazole rich PK profiles were simulated using a POPPK model from the literature, and MAP-BE was used to estimate the clearance based on: (i) the full PK profiles (refCL); and (ii) C24h only (C24h-CL). Xgboost was trained to correct the error between refCL and C24h-CL in the training dataset (75%). C24h-CL as well as ML-corrected C24h-CL were evaluated in a testing dataset (25%) and then in a set of PK profiles simulated using another published POPPK model.

RESULTS

A strong decrease in mean predictive error (MPE%), imprecision (RMSE%) and the number of profiles outside ± 20% MPE% (n-out20%) was observed with the hybrid algorithm (decreased in MPE% by 95.8% and 85.6%; RMSE% by 69.5% and 69.0%; n-out20% by 97.4% and 100% in the training and testing sets, respectively. In the external validation set, the hybrid algorithm decreased MPE% by 96%, RMSE% by 68% and n-out20% by 100%.

CONCLUSION

The hybrid model proposed significantly improved isavuconazole AUC estimation over MAP-BE based on the sole C24h and may improve dose adjustment.

Application of machine learning to predict tacrolimus exposure in liver and kidney transplant patients given the MeltDose formulation

PURPOSE

Machine Learning (ML) algorithms represent an interesting alternative to maximum a posteriori Bayesian estimators (MAP-BE) for tacrolimus AUC estimation, but it is not known if training an ML model using a lower number of full pharmacokinetic (PK) profiles (= "true" reference AUC) provides better performances than using a larger dataset of less accurate AUC estimates. The objectives of this study were: to develop and benchmark ML algorithms trained using full PK profiles to estimate MeltDose^®-tacrolimus individual AUCs using 2 or 3 blood concentrations; and to compare their performance to MAP-BE.

METHODS

Data from liver (n = 113) and kidney (n = 97) transplant recipients involved in MeltDose-tacrolimus PK studies were used for the training and evaluation of ML algorithms. "True" AUC0-24 h was calculated for each patient using the trapezoidal rule on the full PK profile. ML algorithms were trained to estimate tacrolimus true AUC using 2 or 3 blood concentrations. Performances were evaluated in 2 external sets of 16 (renal) and 48 (liver) transplant patients.

RESULTS

Best estimation performances were obtained with the MARS algorithm and the following limited sampling strategies (LSS): predose (0), 8, and 12 h post-dose (rMPE = - 1.28%, rRMSE = 7.57%), or 0 and 12 h (rMPE = - 1.9%, rRMSE = 10.06%). In the external dataset, the performances of the final ML algorithms based on two samples in kidney (rMPE = - 3.1%, rRMSE = 11.1%) or liver transplant recipients (rMPE = - 3.4%, rRMSE = 9.86%) were as good as or better than those of MAP-BEs based on three time points.

CONCLUSION

The MARS ML models developed using "true" MeltDose^®-tacrolimus AUCs yielded accurate individual estimations using only two blood concentrations.

Dried blood microsampling-assisted therapeutic drug monitoring of immunosuppressants: An overview

In the field of solid organ transplantation, chemotherapy and autoimmune disorders, treatment with immunosuppressant drugs requires intensive follow-up of the blood concentrations via therapeutic drug monitoring (TDM) because of their narrow therapeutic window and high intra- and inter-subject variability. This requires frequent hospital visits and venepunctures to allow the determination of these analytes, putting a high burden on the patients. In the context of patient-centric thinking, it is becoming increasingly established that at least part of these conventional blood draws could be replaced by microsampling, allowing home-sampling and increasing the quality of life for these patients. In this review we discuss the published methods - mostly using liquid chromatography coupled to tandem mass spectrometry - that have utilized (volumetric) dried blood samples as an alternative for conventional liquid whole blood for the TDM of immunosuppressant drugs. Furthermore, some pre-analytical considerations using DBS or volumetric alternatives are considered, as well as the applicability on clinical samples. The implementation status in clinical practice is also discussed, including (1) the cost-effectiveness of this approach compared to venepuncture, (2) the availability of multiplexed methods, (3) the status of harmonization and (4) patient perception. A brief perspective on potential future developments for the dried blood-based TDM of immunosuppressant drugs is provided, by considering how obstacles for the implementation of these strategies into clinical practice might be overcome.

Can the Area Under the Curve/Trough Level Ratio Be Used to Optimize Tacrolimus Individual Dose Adjustment?

BACKGROUND

The aim of this work was to evaluate, in a large data set of renal transplant recipients, the intraindividual variability of the area under the curve (AUC)/predose concentration (C0) ratio in comparison with that of AUC, C0, AUC/dose, and C0/dose.

METHODS

Patients with at least 2 tacrolimus AUC estimation requests were extracted from the Immunosuppressant Bayesian dose Adjustment website, and relative variations between 2 consecutive visits for the different metrics were calculated and compared.

RESULTS

Data from 1325 patients on tacrolimus (3827 measured C0 and estimated AUC) showed that the lowest mean relative variation between 2 consecutives visits was for the AUC/C0 ratio (95% confidence interval [CI] relative fold change = -43% to 44% for AUC/C0; 95% CI, -77% to 72% for AUC; 95% CI, -82% to 98% for AUC/dose; 95% CI, -81% to 80% for C0 and 95% CI, -94% to 117% for C0/dose. The correlation between 2 consecutive requests, whether close or far apart, was also best for the AUC/C0 ratio ( r  = 0.33 and r  = 0.34, respectively) in comparison with C0 ( r  = 0.21 and r  = 0.22, respectively) and AUC ( r  = 0.19 and 0.28, respectively). Regression analysis between AUC0-24 and C0 showed that for some patients, the usual C0 targets translated into some very unusual AUC values. As the AUC/C0 ratio is quite stable during large periods, individualized C0 targets can be derived from the AUC targets, and an algorithm that estimates the individualized C0 was developed for situations in which prior AUC estimates are available or not.

CONCLUSIONS

In this study, we confirmed in a large data set that the AUC/C0 ratio yields low intraindividual variability, whereas C0 shows the largest, and we propose to calculate individualized C0 targets based on this ratio.

Plasma metabolomic profiling reveals factors associated with dose-adjusted trough concentration of tacrolimus in liver transplant recipients

Inter- and intrapatient variability of tacrolimus exposure is a vital prognostic risk factor for the clinical outcome of liver transplantation. New factors or biomarkers characterizing tacrolimus disposition is essential for optimal dose prediction in recipients of liver transplant. The aim of the study was to identify potential plasma metabolites associated with the dose-adjusted trough concentration of tacrolimus in liver transplant recipients by using a global metabolomic approach. A total of 693 plasma samples were collected from 137 liver transplant recipients receiving tacrolimus and regular therapeutic drug monitoring. Untargeted metabolomic analysis was performed by ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry. Univariate and multivariate analyses with a mixed linear model were conducted, and the results showed that the dose-adjusted tacrolimus trough concentration was associated with 31 endogenous metabolites, including medium- and long-chain acylcarnitines such as stearoylcarnitine (β = 0.222, p = 0.001), microbiota-derived uremic retention solutes such as indolelactic acid (β = 0.194, p = 0.007), bile acids such as taurohyodeoxycholic acid (β = -0.056, p = 0.002), and steroid hormones such as testosterone (β = 0.099, p = 0.001). A multiple linear mixed model including 11 metabolites and clinical information was established with a suitable predictive performance (correlation coefficient based on fixed effects = 0.64 and correlation coefficient based on fixed and random effects = 0.78). These data demonstrated that microbiota-derived uremic retention solutes, bile acids, steroid hormones, and medium- and long-chain acylcarnitines were the main metabolites associated with the dose-adjusted trough concentration of tacrolimus in liver transplant recipients.

A Prediction Model for Tacrolimus Daily Dose in Kidney Transplant Recipients With Machine Learning and Deep Learning Techniques

Tacrolimus is a major immunosuppressor against post-transplant rejection in kidney transplant recipients. However, the narrow therapeutic index of tacrolimus and considerable variability among individuals are challenges for therapeutic outcomes. The aim of this study was to compare different machine learning and deep learning algorithms and establish individualized dose prediction models by using the best performing algorithm. Therefore, among the 10 commonly used algorithms we compared, the TabNet algorithm outperformed other algorithms with the highest R² (0.824), the lowest prediction error [mean absolute error (MAE) 0.468, mean square error (MSE) 0.558, and root mean square error (RMSE) 0.745], and good performance of overestimated (5.29%) or underestimated dose percentage (8.52%). In the final prediction model, the last tacrolimus daily dose, the last tacrolimus therapeutic drug monitoring value, time after transplantation, hematocrit, serum creatinine, aspartate aminotransferase, weight, CYP3A5, body mass index, and uric acid were the most influential variables on tacrolimus daily dose. Our study provides a reference for the application of deep learning technique in tacrolimus dose estimation, and the TabNet model with desirable predictive performance is expected to be expanded and applied in future clinical practice.

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.

Point-of-Care Therapeutic Drug Monitoring for Precision Dosing of Immunosuppressive Drugs

BACKGROUND

Immunosuppressive drugs (ISD) are an essential tool in the treatment of transplant rejection and immune-mediated diseases. Therapeutic drug monitoring (TDM) for determination of ISD concentrations in biological samples is an important instrument for dose personalization for improving efficacy while reducing side effects. While currently ISD concentration measurements are performed at specialized, centralized facilities, making the process complex and laborious for the patient, various innovative technical solutions have recently been proposed for bringing TDM to the point-of-care (POC).

CONTENT

In this review, we evaluate current ISD-TDM and its value, limitations, and proposed implementations. Then, we discuss the potential of POC-TDM in the era of personalized medicine, and provide an updated review on the unmet needs and available technological solutions for the development of POC-TDM devices for ISD monitoring. Finally, we provide concrete suggestions for the generation of a meaningful and more patient-centric process for ISD monitoring.

SUMMARY

POC-based ISD monitoring may improve clinical care by reducing turnaround time, by enabling more frequent measurements in order to obtain meaningful pharmacokinetic data (i.e., area under the curve) faster reaction in case of problems and by increasing patient convenience and compliance. The analysis of the ISD-TDM field prompts the evolution of POC testing toward the development of fully integrated platforms able to support clinical decision-making. We identify 4 major areas requiring careful combined implementation: patient usability, data meaningfulness, clinicians' acceptance, and cost-effectiveness.

Evaluation of Bayesian Forecasting Methods for Prediction of Tacrolimus Exposure Using Samples Taken on Two Occasions in Adult Kidney Transplant Recipients

BACKGROUND

Bayesian forecasting-based limited sampling strategies (LSSs) for tacrolimus have not been evaluated for the prediction of subsequent tacrolimus exposure. This study examined the predictive performance of Bayesian forecasting programs/services for the estimation of future tacrolimus area under the curve (AUC) from 0 to 12 hours (AUC0-12) in kidney transplant recipients.

METHODS

Tacrolimus concentrations were measured in 20 adult kidney transplant recipients, 1 month post-transplant, on 2 occasions one week apart. Twelve samples were taken predose and 13 samples were taken postdose at the specified times on the first and second sampling occasions, respectively. The predicted AUC0-12 (AUCpredicted) was estimated using Bayesian forecasting programs/services and data from both sampling occasions for each patient and compared with the fully measured AUC0-12 (AUCmeasured) calculated using the linear trapezoidal rule on the second sampling occasion. The bias (median percentage prediction error [MPPE]) and imprecision (median absolute prediction error [MAPE]) were determined.

RESULTS

Three programs/services were evaluated using different LSSs (C0; C0, C1, C3; C0, C1, C2, C4; and all available concentrations). MPPE and MAPE for the prediction of fully measured AUC0-12 were <15% for each program/service (with the exclusion of when only C0 was used), when using estimated AUC from data on the same (second) occasion. The MPPE and MAPE for the prediction of a future fully measured AUC0-12 were <15% for 2 programs/services (and for the third when participants who had a tacrolimus dose change between sampling days were excluded), when the occasion 1-AUCpredicted, using C0, C1, and C3, was compared with the occasion 2-AUCmeasured.

CONCLUSIONS

All 3 Bayesian forecasting programs/services evaluated had acceptable bias and imprecision for predicting a future AUC0-12, using tacrolimus concentrations at C0, C1, and C3, and could be used for the accurate prediction of tacrolimus exposure in adult kidney transplant recipients.

Advances in personalized medicine and noninvasive diagnostics in solid organ transplantation

Personalized medicine has been a mainstay and in practice in transplant pharmacotherapy since the advent of the field. Decisions pertaining to the diagnosis, selection, and monitoring of transplant pharmacotherapy are aimed toward the individual, the allograft, and the overall immunologic needs of the patient. Recent advances in pharmacogenomics, noninvasive biomarkers, and artificial intelligence (AI) technologies have the promise of transforming the way we individualize treatment and monitor allograft function. Pharmacogenomic testing can provide clinicians with additional data that can minimize toxicity and maximize therapeutic dosing in high-risk patients, leading to more informed decisions that may decrease the risk of rejection and adverse outcomes related to immunosuppressive therapies. Development of noninvasive strategies to monitor allograft function may offer safer and more convenient methods to detect allograft injury. Cell free DNA and gene expression profiling offer the potential to serve as "liquid biopsies" minimizing the risk to patients and providing clinicians with useful molecular data that may help individualize immunosuppression and rejection treatment. Use of big data in transplant and novel AI platforms, such as the iBox, hold tremendous promise in providing clinicians a "glimpse into the future" thereby allowing for a more individualized approach to immunosuppressive therapy that may minimize future adverse outcomes. Advances in diagnostics, laboratory science, and AI have made the application of personalized medicine even more tailored for solid organ transplant recipients. In this perspective, we summarize the current and emerging tools available, literature supporting use, and the horizon for future personalization of transplantation.

Assessment of pharmacokinetic mycophenolic acid clearance models using Monte Carlo numerical analysis

Previously, we performed population pharmacokinetic analysis and indicated age, mycophenolate mofetil (MMF)/mycophenolic acid (MPA) daily dose, and presence of nifedipine in patient therapy as significant predictors of MPA apparent clearance (CL/F) variability. This study aimed to determine the reliability of previously published population pharmacokinetic models derived from similar studies. Furthermore, this study investigated correspondence between chosen population models from the literature.By means of the Monte Carlo simulation method, pharmacokinetic models from different studies are simulated and analysed in the range of standard deviations of measured system parameters as well as the range of observed model parameters taken from the comparison studies.The 1000 numerical simulations were performed for every analysed model in order to calculate the most possible MPA CL/F values according to the expected values from the performed experiment. Fitting our results with other models showed how the presence of nifedipine makes difference in MPA CL/F values.By testing the data from selected studies into our model, a similar range of expected CL/F values was obtained, which may confirm the validity of our model. The results of our population pharmacokinetic study are partially applicable in models by other researchers.

Maintenance Immunosuppression in Solid Organ Transplantation: Integrating Novel Pharmacodynamic Biomarkers to Inform Calcineurin Inhibitor Dose Selection

Calcineurin inhibitors, the primary immunosuppressive therapy used to prevent alloreactivity of transplanted organs, have a narrow therapeutic index. Currently, treatment is individualized based on clinical assessment of the risk of rejection or toxicity guided by trough concentration monitoring. Advances in immune monitoring have identified potential markers that may have value in understanding calcineurin inhibitor pharmacodynamics. Integration of these markers has the potential to complement therapeutic drug monitoring. Existing pharmacokinetic-pharmacodynamic (PK-PD) data is largely limited to correlation between the biomarker and trough concentrations at single time points. Immune related gene expression currently has the most evidence supporting PK-PD integration. Novel biomarker-based approaches to pharmacodynamic monitoring including development of enhanced PK-PD models are proposed to realize the full clinical benefit.

Assessment of the glomerular filtration rate (GFR) in kidney transplant recipients using Bayesian estimation of the iohexol clearance

Background Plasma iohexol clearance (CLiohexol) is a reference technique for glomerular filtration rate (GFR) determination. In routine practice, CLiohexol is calculated using one of several formulas, which have never been evaluated in kidney transplant recipients. We aimed to model iohexol pharmacokinetics in this population, evaluate the predictive performance of three simplified formulas and evaluate whether a Bayesian algorithm improves CLiohexol estimation. Methods After administration of iohexol, six blood samples were drawn from 151 patients at various time points. The dataset was split into two groups, one to develop the population pharmacokinetic (POPPK) model (n = 103) and the other (n = 48) to estimate the predictive performances of the various GFR estimation methods. GFR reference values (GFRref) in the validation dataset were obtained by non-compartmental pharmacokinetic (PK) analysis. Predictive performances of each method were evaluated in terms of bias (ME), imprecision (root mean square error [RMSE]) and number of predictions out of the ±10% or 15% error interval around the GFRref. Results A two-compartment model best fitted the data. The Bayesian estimator with samples drawn at 30, 120 and 270 min allowed accurate prediction of GFRref (ME = 0.47%, RMSE = 3.42%), as did the Brøchner-Mortensen (BM) formula (ME = - 0.0425%, RMSE = 3.40%). With both methods, none of the CL estimates were outside the ±15% interval and only 2.4% were outside the ±10% for the BM formula (and none for the Bayesian estimator). In patients with GFR ≤30 mL/min/1.73 m2, the BM formula performed very well, while the Bayesian method could not be evaluated in depth due to too small a number of patients with adequate sampling times. Conclusions GFR can be estimated with acceptable accuracy in kidney transplant patients using the BM formula, but also using a Bayesian algorithm.

Optimization of tacrolimus in kidney transplantation: New pharmacokinetic perspectives

Tacrolimus is the cornerstone of immunosuppressive therapy after kidney transplantation (KT), but its use is complicated by a narrow therapeutic index and high inter- and intra-patient pharmacokinetic variability. There are three available oral formulations of tacrolimus: immediate-release tacrolimus (IR-Tac), extended-release tacrolimus (ER-Tac) and a MeltDose® (LCPT) formulation, the latter favoring a prolonged drug release and increased bioavailability. The time-concentration curves of these formulations are different. Compared with IR-Tac and ER-Tac, LCPT has a relatively flat pharmacokinetic profile with less fluctuation between trough and peak exposures, and a delayed peak concentration. This translates to a more stable delivery of tacrolimus and may alleviate the risk of underexposure and allograft rejection or overexposure and toxicity. The once-daily formulation of both ER-TAC and LCPT may also offer a potential advantage on patient adherence. Fast metabolizers of tacrolimus, the elderly, and human leukocyte antigen-sensitized patients are at risk of poorer outcomes after KT, possibly associated with a different exhibited pharmacokinetics of tacrolimus or different requirements in terms of exposure. Simple, practical strategies are needed to identify patients at risk of suboptimal KT outcomes and those who would benefit from a more proactively personalized approach to tacrolimus treatment. This review aims to increase awareness of the link between the pharmacokinetics of oral tacrolimus formulations and the clinical needs of patients after KT, particularly among those who have clinically significant pharmacokinetic variation of tacrolimus.

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.

Measured GFR by Utilizing Population Pharmacokinetic Methods to Determine Iohexol Clearance

Introduction

There is an increasing demand for accurately measured glomerular filtration rate (GFR). Iohexol serum clearance has become a new gold standard, but it is challenging when GFR is low and 24-hour sampling is required for accurate results. The primary aim of this study was to develop an iohexol pharmacokinetic population model for accurate determination of individual GFR using limited sampling for up to 5 hours also when renal function is <40 ml/min.

Methods

A nonparametric iohexol population pharmacokinetic model was developed with rich data from 176 patients. In a validation cohort of 43 patients, a model-determined GFR (iohexol clearance) using different limited sampling strategies for up to 5 hours was compared with the strategy currently used in routine care, a log-linear 2-point method. In all, 1526 iohexol concentrations were used, from patients ranging in age from 1 to 82 years and GFR from 14 to 149 ml/min.

Results

The clinical 2-point method showed insufficient agreement compared with reference values; 15% of GFR values had an error of greater than ±10% even when sampling for 24 hours when estimating GFR <40 ml/min per 1.73 m² (standard procedure). Restricted sampling the first 5 hours with the population model required 4 samples to determine GFR accurately. This strategy showed excellent agreement with the reference; <3% of GFR values had an error greater than ±10 %.

Conclusion

Using an iohexol population pharmacokinetic model allows for accurate determination of GFR within 5 hours when applying 4 optimally timed samples, even in patients with GFR <40 ml/min.

A Novel, Dose-Adjusted Tacrolimus Trough-Concentration Model for Predicting and Estimating Variance After Kidney Transplantation

Background and Objective

Given that a high intrapatient variability (IPV) of tacrolimus whole blood concentration increases the risk for a poor kidney transplant outcome, some experts advocate routine IPV monitoring for detection of high-risk patients. However, attempts to estimate the variance of tacrolimus trough concentrations (TTC) are limited by the need for patients to receive a fixed dose over time and/or the use of linear statistical models. A goal of this study is to overcome the current limitations through the novel application of statistical methodology generalizing the relationship between TTC and dose through the use of nonparametric functional regression modeling.

Methods

With TTC as a response and dose as a covariate, the model employs an unknown bivariate function, allowing for the potentially complex, nonlinear relationship between the two parameters. A dose-adjusted variance of TTC is then derived based on standard functional principal component analysis (FPCA). To assess the model, it was compared against an FPCA-based model and linear mixed-effects models using prediction error, bias, and coverage probabilities for simulated data as well as phase III data from the Astellas new drug application studies for extended-release tacrolimus.

Results

Our numerical investigation indicates that the new model better predicts dose-adjusted TTCs compared with the prediction of linear mixed effects models. Estimated coverage probabilities also indicate that the new model accurately accounts for the variance of TTC during the periods of large fluctuation in dose, whereas the linear mixed effects model consistently underestimates the coverage probabilities because of the inaccurate characterization of TTC fluctuation.

Conclusion

This is the first known application of a functional regression model to assess complex relationships between TTC and dose in a real clinical setting. This new method has applicability in future clinical trials including real-world data sets due to flexibility of the nonparametric modeling approach.

Electronic supplementary material

The online version of this article (10.1007/s40268-019-0271-2) contains supplementary material, which is available to authorized users.

Identification of Factors Affecting Tacrolimus Trough Levels in Latin American Pediatric Liver Transplant Patients

Tacrolimus is the cornerstone in pediatric liver transplant immunosuppression. Despite close monitoring, fluctuations in tacrolimus blood levels affect safety and efficacy of immunosuppressive treatments. Identifying the factors related to the variability in tacrolimus exposure may be helpful in tailoring the dose. The aim of the present study was to characterize the clinical, pharmacological, and genetic variables associated with systemic tacrolimus exposure in pediatric liver transplant patients. De novo transplant patients with a survival of more than 1 month were considered for inclusion and were genotyped for cytochrome P450 3A5 (CYP3A5). Peritransplant clinical factors and laboratory covariates were recorded retrospectively between 1 month and 2 years after transplant, including alanine aminotransferase (ALT), aspartate aminotransferase, hematocrit, and tacrolimus predose steady-state blood concentrations collected 12 hours after tacrolimus dosing. A linear mixed effect (LME) model was used to assess the association of these factors and the log-transformed tacrolimus dose-normalized trough concentration (logC0/D) levels. Bootstrapping was used to internally validate the final model. External validation was performed in an independent group of patients who matched the original population. The developed LME model described that logC0/D increases with increases in time after transplant (β = 0.019, 95% confidence interval [CI], 0.010-0.028) and ALT values (β = 0.00030, 95% CI, 0.00002-0.00056), whereas logC0/D is significantly lower in graft CYP3A5 expressers compared with nonexpressers (β = -0.349, 95% CI, -0.631 to -0.062). In conclusion, donor CYP3A5 genotype, time after transplant, and ALT values are associated with tacrolimus disposition between 1 month and 2 years after transplant. A better understanding of tacrolimus exposure is essential to minimize the occurrence of an out-of-range therapeutic window that may lead to adverse drug reactions or acute rejection.

Tacrolimus Levels in the Prophylaxis of Acute Graft-Versus-Host Disease in the Chinese Early After Hematopoietic Stem Cell Transplantation

BACKGROUND

Tacrolimus has been widely accepted as the backbone of acute graft-versus-host disease (aGVHD) prophylaxis in allogeneic hematopoietic stem cell transplantation (alloHSCT). The present work evaluated whether tacrolimus concentrations early after transplant correlate with the incidence of aGVHD in Chinese alloHSCT recipients.

METHODS

One hundred four Chinese alloHSCT recipients were included in this retrospective study. All patients received standard prophylaxis with tacrolimus and short-term methotrexate. Blood samples were taken at steady-state for those on i.v. tacrolimus (Cv) or predose (C0) and 2 hours after the last oral dose (C2).

RESULTS

In the first 8 weeks after alloHSCT, significant variability in Cv, C0, and C2 of Chinese patients was observed. It was found that higher tacrolimus C0 and C2 values tended to be associated with a reduced risk of aGVHD, although this was a nonsignificant trend due to the small sample size involved. Receiver operating characteristic curve analysis indicated that Cv levels of ≥16.52 ng/mL, C0 levels of ≥5.56 ng/mL, and C2 levels of ≥7.83 ng/mL minimized the incidence of treatment failure during weeks 3-4 with intravenous administration and weeks 5-6 with oral administration. There was no statistically significant association of the patient liver and kidney function with the blood concentration of tacrolimus in the desired range of 5-20 ng/mL.

CONCLUSIONS

Tacrolimus therapeutic drug monitoring improved treatment outcomes of Chinese alloHSCT recipients. Cv measurements during weeks 3-4 and C0 or C2 measurements during weeks 5-6 better predicted aGVHD (I-IV) than the concentrations measured at other time points during the first 6 weeks after alloHSCT.

Clinical Evaluation of the Tolerability, Pharmacokinetics, and Inhibition of Platelet Aggregation of Eptifibatide in Healthy Chinese Subjects

The present study aimed to evaluate the pharmacokinetic properties and antiplatelet aggregation activity of eptifibatide in healthy Chinese subjects. Eptifibatide (180 μg/kg) was administrated by 2 bolus injections 10 minutes apart, followed by a 2.0 μg/kg/min infusion for 24 hours. The eptifibatide pharmacokinetic and antiplatelet aggregation activities were evaluated using nonlinear mixed-effects modeling and noncompartmental analysis. Safety assessments included adverse events, hematology, and biochemistry tests. Twelve Chinese healthy subjects were enrolled and completed the study. Steady-state concentrations were achieved at 0.5 to 24 hours after dosing. The median time to maximum concentration was 13 minutes, and the mean terminal elimination half-life was 148.19 minutes. The effective inhibition of platelet aggregation (<20% platelet aggregation) occurred by 3 minutes after starting dosing to 4 hours after termination of the infusion. Eptifibatide concentrations were fitted with a 3-compartment model, and the typical value of clearance was 0.11 L/min, with no significant covariates found. Three mild adverse events were detected in the study. Eptifibatide displays high sensitivity and excellent tolerability in healthy Chinese subjects. The dosage of eptifibatide recommended on the label for whites can effectively inhibit platelet aggregation.