To be or not to be, when synthetic data meet clinical pharmacology: A focused study on pharmacogenetics.

The use of synthetic data in pharmacology research has gained significant attention due to its potential to address privacy concerns and promote open science. In this study, we implemented and compared three synthetic data generation methods, CT-GAN, TVAE, and a simplified implementation of Avatar, for a previously published pharmacogenetic dataset of 253 patients with one measurement per patient (non-longitudinal). The aim of this study was to evaluate the performance of these methods in terms of data utility and privacy trade off. Our results showed that CT-GAN and Avatar used with k = 10 (number of patients used to create the local model of generation) had the best overall performance in terms of data utility and privacy preservation. However, the TVAE method showed a relatively lower level of performance in these aspects. In terms of Hazard ratio estimation, Avatar with k = 10 produced HR estimates closest to the original data, whereas CT-GAN slightly underestimated the HR and TVAE showed the most significant deviation from the original HR. We also investigated the effect of applying the algorithms multiple times to improve results stability in terms of HR estimation. Our findings suggested that this approach could be beneficial, especially in the case of small datasets, to achieve more reliable and robust results. In conclusion, our study provides valuable insights into the performance of CT-GAN, TVAE, and Avatar methods for synthetic data generation in pharmacogenetic research. The application to other type of data and analyses (data driven) used in pharmacology should be further investigated.

Everolimus Personalized Therapy: Second Consensus Report by the International Association of Therapeutic Drug Monitoring and Clinical Toxicology.

ABSTRACT: The Immunosuppressive Drugs Scientific Committee of the International Association of Therapeutic Drug Monitoring and Clinical Toxicology established the second consensus report to guide Therapeutic Drug Monitoring (TDM) of everolimus (EVR) and its optimal use in clinical practice 7 years after the first version was published in 2016. This version provides information focused on new developments that have arisen in the last 7 years. For the general aspects of the pharmacology and TDM of EVR that have retained their relevance, readers can refer to the 2016 document. This edition includes new evidence from the literature, focusing on the topics updated during the last 7 years, including indirect pharmacological effects of EVR on the mammalian target of rapamycin complex 2 with the major mechanism of direct inhibition of the mammalian target of rapamycin complex 1. In addition, various concepts and technical options to monitor EVR concentrations, improve analytical performance, and increase the number of options available for immunochemical analytical methods have been included. Only limited new pharmacogenetic information regarding EVR has emerged; however, pharmacometrics and model-informed precision dosing have been constructed using physiological parameters as covariates, including pharmacogenetic information. In clinical settings, EVR is combined with a decreased dose of calcineurin inhibitors, such as tacrolimus and cyclosporine, instead of mycophenolic acid. The literature and recommendations for specific organ transplantations, such as that of the kidneys, liver, heart, and lungs, as well as for oncology and pediatrics have been updated. EVR TDM for pancreatic and islet transplantation has been added to this edition. The pharmacodynamic monitoring of EVR in organ transplantation has also been updated. These updates and additions, along with the previous version of this consensus document, will be helpful to clinicians and researchers treating patients receiving EVR.

Daptomycin Dosage Optimization in Renal Impairment Using Model-Informed Precision Dosing.

BACKGROUND: Daptomycin's efficacy and toxicity are closely related to its exposure, which can vary widely among individuals. The patient, a 59-year-old male with an estimated glomerular filtration rate (eGFR) of 12 mL/min/1.73 m² and a weight of 64 kg, was treated with 850 mg of daptomycin every other day for infective endocarditis caused by methicillin-resistant Staphylococcus aureus (MRSA). For patients with an estimated glomerular filtration rate of less than 30 mL/min/1.73 m², the dosing recommendations are not explicitly defined in the endocarditis guidelines. Subsequently, the pharmacology department was contacted to adjust the dosage. METHODS: A population pharmacokinetic model developed by Dvorchik et al. was used for Bayesian estimation of the patient's pharmacokinetic parameters. The 24-hour area under the curve (AUC24) of daptomycin was calculated at steady state using peak and trough plasma samples. RESULTS: The minimum inhibitory concentration (MIC) of the MRSA strain was 0.25 mg/L. An AUC24/MIC ratio below 666 is associated with higher mortality risk, while an AUC24 above 939 h·mg/L correlates with increased risk of muscular toxicity. Initial AUC24 estimation was 1091 h·mg/L. Following a dosage reduction to 700 mg every other day, the AUC24 increased to 1600 h·mg/L. Further reduction to 500 mg every other day brought the AUC24 down to 750 h mg/L, with two subsequent measurements showing consistent AUC24 values of 500 h·mg/L, which is within the target range. CONCLUSIONS: Daptomycin ended 6 weeks after the initial negative blood culture, with no adverse effects or recurrence of MRSA infection. This case underscores the need for therapeutic drug monitoring and a multidisciplinary approach to adjust daptomycin doses in patients with renal impairment.

Population Pharmacokinetic of Vancomycin Administered by Continuous Infusion in Critically Ill Patients.

INTRODUCTION: Administration of vancomycin dose by continuous infusion (CI) according to population pharmacokinetic (Pop Pk) models is highly recommended in critically ill patients who exhibit pathophysiological changes. OBJECTIVE: The objective of this study was to develop and validate a Pop Pk model of vancomycin administered by CI in critically ill patients with normal and impaired renal functions. METHODS: The Pop Pk study was performed using a nonparametric approach (Pmetrics*). The influence of covariates (gender, age, weight, height, and creatinine clearance [Cr-Cl]) was tested on the model's Pk parameters. The performance of the final model was assessed using an external dataset. RESULTS: A one-compartment model (volume of distribution [Vd], elimination from compartment [Ke]) was found to show a good prediction performance. The influence of covariates has shown that age and Cr-Cl affected significantly Vd and Ke, respectively. The distribution of simulated vancomycin clearance (CLv) according to different renal function levels showed a negative correlation between CLv and the severity of the renal impairment. The internal validation of the final model showed that the plot of individual-predicted concentration versus observed concentration resulted in r2 = 0.86 in the final model. The external validation of the final model showed an acceptable predictive performance. CONCLUSION: We developed a Pop Pk model for vancomycin administered by CI in critically ill patients. A significant impact of Cr-Cl and different stages of renal failure on CLv has been demonstrated. The establishment of an individualized proposal dose based on this model may be helpful to achieve the target range which is critical in optimizing the efficacy and safety of this antibiotic.

Application of machine-learning models to predict the ganciclovir and valganciclovir exposure in children using a limited sampling strategy.

Intravenous ganciclovir and oral valganciclovir display significant variability in ganciclovir pharmacokinetics, particularly in children. Therapeutic drug monitoring currently relies on the area under the concentration-time (AUC). Machine-learning (ML) algorithms represent an interesting alternative to Maximum-a-Posteriori Bayesian-estimators for AUC estimation. The goal of our study was to develop and validate an ML-based limited sampling strategy (LSS) approach to determine ganciclovir AUC0-24 after administration of either intravenous ganciclovir or oral valganciclovir in children. Pharmacokinetic parameters from four published population pharmacokinetic models, in addition to the World Health Organization growth curve for children, were used in the mrgsolve R package to simulate 10,800 pharmacokinetic profiles of children. Different ML algorithms were trained to predict AUC0-24 based on different combinations of two or three samples. Performances were evaluated in a simulated test set and in an external data set of real patients. The best estimation performances in the test set were obtained with the Xgboost algorithm using a 2 and 6 hours post dose LSS for oral valganciclovir (relative mean prediction error [rMPE] = 0.4% and relative root mean square error [rRMSE] = 5.7%) and 0 and 2 hours post dose LSS for intravenous ganciclovir (rMPE = 0.9% and rRMSE = 12.4%). In the external data set, the performance based on these two sample LSS was acceptable: rMPE = 0.2% and rRMSE = 16.5% for valganciclovir and rMPE = -9.7% and rRMSE = 17.2% for intravenous ganciclovir. The Xgboost algorithm developed resulted in a clinically relevant individual estimation using only two blood samples. This will improve the implementation of AUC-targeted ganciclovir therapeutic drug monitoring in children.

A Machine Learning Algorithm to Predict the Starting Dose of Daptomycin.

BACKGROUND AND OBJECTIVE: The dosage of daptomycin is usually based on body weight. However, it has been shown that this approach yields too high an exposure in obese patients. Pharmacokinetic and pharmacodynamic indexes (PK/PD) have been proposed for daptomycin's antibacterial effect (AUC/CMI >666) and toxicity (C0 > 24.3 mg/L). We previously developed machine learning (ML) algorithms to predict starting doses based on Monte Carlo simulations. We propose a new way to perform probability of target attainment based on an ML algorithm to predict the daptomycin starting dose. METHODS: The Dvorchik model of daptomycin was implemented in the mrgsolve R package and 4950 pharmacokinetic profiles were simulated with doses ranging from 4 to 12 mg/kg. We trained and benchmarked four machine learning algorithms and selected the best to iteratively search for the optimal dose of daptomycin maximizing the event (AUC/CMI > 666 and C0 < 24.3 mg/L). The ML algorithm was evaluated in simulations and an external database of real patients in comparison with population pharmacokinetics. RESULTS: The performance of the Xgboost algorithms developed to predict the event (ROC AUC) in the training and test set were 0.762 and 0.761, respectively. The most important prediction variables were dose, creatinine clearance, body weight and sex. In the external database of real patients, the starting dose administered based on the ML algorithm significantly improved the target attainment by 7.9% (p-value = 0.02929) in comparison with the dose administered based on body weight. CONCLUSION: The developed algorithm improved the target attainment for daptomycin in comparison with weight-based dosing. We built a Shiny app to calculate the optimal starting dose.

Synergizing pharmacometrics and pharmacovigilance for medication error management: the case of secukinumab.

PURPOSE: To demonstrate the effective integration of pharmacometrics and pharmacovigilance in managing medication errors, highlighted by a case involving secukinumab in a patient with hidradenitis suppurativa. METHODS: We present the case of a 41-year-old male with progressive hidradenitis suppurativa, unresponsive to multiple antibiotic regimens and infliximab treatment. Due to a medication error, the patient received 300 mg of secukinumab daily for 4 days instead of weekly, totaling 1200 mg. The regional pharmacovigilance center assessed potential toxicity, and a pharmacometric analysis using a population pharmacokinetic model was performed to inform dosing adjustments. RESULTS: Clinical data indicated that the received doses were within a non-toxic range. No adverse effects were observed. Pharmacometric simulations revealed a risk of underexposure due to the dosing error. Based on these simulations, it was recommended to restart monthly secukinumab injections on day 35 after the initial dose. Measured plasma concentrations before re-administration confirmed the model's accuracy. CONCLUSION: This case highlights the crucial collaboration between clinical services, pharmacovigilance, and pharmacometrics in managing medication errors. Such interdisciplinary efforts ensure therapeutic efficacy and patient safety by maintaining appropriate drug exposure levels.

Population Pharmacokinetics and Limited Sampling Strategy of Mycophenolate Mofetil for Indian Patients With Lupus Nephritis.

BACKGROUND: Mycophenolic acid is widely used to treat lupus nephritis (LN). However, it exhibits complex pharmacokinetics with large interindividual variability. This study aimed to develop a population pharmacokinetic (popPK) model and a 3-sample limited sampling strategy (LSS) to optimize therapeutic drug monitoring in Indian patients with LN. METHODS: Five blood samples from each LN patient treated with mycophenolic acid were collected at steady-state predose and 1, 2, 4, and 6 hours postdose. Demographic parameters were tested as covariates to explain interindividual variability. PopPK analysis was performed using Monolix and the stochastic approximation expectation-maximization algorithm. An LSS was derived from 500 simulated pharmacokinetic (PK) profiles using maximum a posteriori Bayesian estimation to estimate individual PK parameters and area under the curve (AUC). The LSS-calculated AUC was compared with the AUC calculated using the trapezoidal rule and all the simulated samples. RESULTS: A total of 51 patients were included in this study. Based on the 245 mycophenolic acid concentrations, a 1-compartmental model with double absorption using gamma distributions best fitted the data. None of the covariates improved the model significantly. The model was internally validated using diagnostic plots, prediction-corrected visual predictive checks, and bootstrapping. The best LSS included samples at 1, 2, and 4 hours postdose and exhibited good performances in an external dataset (root mean squared error, 21.9%; mean bias, -4.20%). CONCLUSIONS: The popPK model developed in this study adequately estimated the PK of mycophenolic acid in adult Indian patients with LN. This simple LSS can optimize TDM based on the AUC in routine practice.

A machine learning approach to predict daptomycin exposure from two concentrations based on Monte Carlo simulations.

Daptomycin is a concentration-dependent lipopeptide antibiotic for which exposure/effect relationships have been shown. Machine learning (ML) algorithms, developed to predict the individual exposure to drugs, have shown very good performances in comparison to maximum a posteriori Bayesian estimation (MAP-BE). The aim of this work was to predict the area under the blood concentration curve (AUC) of daptomycin from two samples and a few covariates using XGBoost ML algorithm trained on Monte Carlo simulations. Five thousand one hundred fifty patients were simulated from two literature population pharmacokinetics models. Data from the first model were split into a training set (75%) and a testing set (25%). Four ML algorithms were built to learn AUC based on daptomycin blood concentration samples at pre-dose and 1 h post-dose. The XGBoost model (best ML algorithm) with the lowest root mean square error (RMSE) in a 10-fold cross-validation experiment was evaluated in both the test set and the simulations from the second population pharmacokinetic model (validation). The ML model based on the two concentrations, the differences between these concentrations, and five other covariates (sex, weight, daptomycin dose, creatinine clearance, and body temperature) yielded very good AUC estimation in the test (relative bias/RMSE = 0.43/7.69%) and validation sets (relative bias/RMSE = 4.61/6.63%). The XGBoost ML model developed allowed accurate estimation of daptomycin AUC using C0, C1h, and a few covariates and could be used for exposure estimation and dose adjustment. This ML approach can facilitate the conduct of future therapeutic drug monitoring (TDM) studies.

Optimization of Ganciclovir and Valganciclovir Starting Dose in Children by Machine Learning.

BACKGROUND AND OBJECTIVES: Ganciclovir (GCV) and valganciclovir (VGCV) show large interindividual pharmacokinetic variability, particularly in children. The objectives of this study were (1) to develop machine learning (ML) algorithms trained on simulated pharmacokinetics profiles obtained by Monte Carlo simulations to estimate the best ganciclovir or valganciclovir starting dose in children and (2) to compare its performances on real-world profiles to previously published equation derived from literature population pharmacokinetic (POPPK) models achieving about 20% of profiles within the target. MATERIALS AND METHODS: The pharmacokinetic parameters of four literature POPPK models in addition to the World Health Organization (WHO) growth curve for children were used in the mrgsolve R package to simulate 10,800 pharmacokinetic profiles. ML algorithms were developed and benchmarked to predict the probability to reach the steady-state, area-under-the-curve target (AUC0-24 within 40-60 mg × h/L) based on demographic characteristics only. The best ML algorithm was then used to calculate the starting dose maximizing the target attainment. Performances were evaluated for ML and literature formula in a test set and in an external set of 32 and 31 actual patients (GCV and VGCV, respectively). RESULTS: A combination of Xgboost, neural network, and random forest algorithms yielded the best performances and highest target attainment in the test set (36.8% for GCV and 35.3% for the VGCV). In actual patients, the best GCV ML starting dose yielded the highest target attainment rate (25.8%) and performed equally for VGCV with the Franck model formula (35.3% for both). CONCLUSION: The ML algorithms exhibit good performances in comparison with previously validated models and should be evaluated prospectively.

Limited sampling strategy and population pharmacokinetic model of mycophenolic acid in pediatric patients with systemic lupus erythematosus: application of a double gamma absorption model with SAEM algorithm.

INTRODUCTION: Mycophenolic acid (MPA), the active metabolite of mycophenolate mofetil (MMF), is widely used in the treatment of systemic lupus erythematosus (SLE). It has been shown that its therapeutic drug monitoring based on the area under the curve (AUC) improves treatment efficacy. MPA exhibits a complex bimodal absorption, and a double gamma distribution model has been already proposed in the past to accurately describe this phenomenon. These previous population pharmacokinetics models (POPPK) have been developed using iterative two stage Bayesian (IT2B) or non-parametric adaptive grid (NPAG) methods. However, non-linear mixed effect (NLME) approaches based on stochastic approximation expectation-maximization (SAEM) algorithms have never been published so far for this particular model. The objectives of this study were (i) to implement the double absorption gamma model in Monolix, (ii) to compare different absorption models to describe the pharmacokinetics of MMF, and (iii) to develop a limited sampling strategy (LSS) to estimate AUC in pediatric SLE patients. MATERIAL AND METHODS: A data splitting of full pharmacokinetic profiles sampled in 67 children extracted either from the expert system ISBA (n = 34) or the hospital Saint Louis (n = 33) was performed into train (75%) and test (25%) sets. A POPPK was developed for MPA in the train set using a NLME and the SAEM algorithm and different absorption models were implemented and compared (first order, transit, or simple and double gamma). The best limited sampling strategy was then determined in the test set using a maximum-a-posteriori Bayesian method to estimate individual PK parameters and AUC based on three blood samples compared to the reference AUC calculated using the trapezoidal rule applied on all samples and performances were assessed in the test set. RESULTS: Mean patient age and dose was 13 years old (5-18) and 18.1 mg/kg (7.9-47.6), respectively. MPA concentrations (764) from 107 occasions were included in the analysis. A double gamma absorption with a first-order elimination from the central compartment best fitted the data. The optimal LSS with samples at 30 min, 2 h, and 3 h post-dose exhibited good performances in the test set (mean bias - 0.32% and RMSE 21.0%). CONCLUSION: The POPPK developed in this study adequately estimated the MPA AUC in pediatric patients with SLE based on three samples. The double absorption gamma model developed with the SAEM algorithm showed very accurate fit and reduced computation time.

In Silico Pharmacokinetics Evaluation of Forgiveness for Doravirine and Rilpivirine.

BACKGROUND: This study aimed to evaluate the concentrations of rilpivirine (RLP) and doravirine (DOR) after 3 days-off using simulations from population pharmacokinetics models. METHODS: The authors conducted a series of 500 sets of 10,000 Monte Carlo simulations to examine the steady-state conditions for 2 common dosage levels: 25 mg/d for RLP and 100 mg/d for DOR. These simulations were conducted under 2 scenarios: 1 without drug cessation and another after a 3-day break. The validity of the implementation was established through a comparison of median trough concentrations (C24h) with previously reported data. Subsequently, the proportion of simulated patients with C24h and C72h after 3 days-off (C72h/3do) that exceeded the inhibitory concentration 50 (IC50), 5.2 mcg/L for DOR and 20.5 mcg/L for RLP respectively, was calculated. The inhibitory quotient (IQ) was also computed, which was 6 times IC50 for DOR and 4.5 times IC50 for RLP. Finally, nomograms were constructed to estimate the probability of having C72h/3do > IC50 or > IQ for different ranges of C24h. RESULTS: Simulated C24h median ± SD for RLP were 61.8 ± 0.4 mcg/L and for DOR 397 ± 0 mcg/L. For RLP, 99.3 ± 0.1% exceeded IC50 at C24h, 16.4 ± 0.4% at C72h/3do, and none surpassed the IQ threshold. In contrast, DOR had 100% ± 0% above IC50 at C24h, 93.6 ± 0.2% at C72h/3do, and 58.6 ± 0.5% exceeded the IQ. CONCLUSIONS: These findings suggest that treatment with DOR may offer a more forgiving therapeutic profile than RLP, given the larger proportion of patients achieving effective drug exposure with DOR. However, it is important to acknowledge a significant limitation of this study, namely, the assumption that drug concentration is a perfect surrogate for drug effectiveness.

Diagnostic, prognostic and clinical value of left ventricular radial strain to identify paradoxical septal motion in ventilated patients with the acute respiratory distress syndrome: an observational prospective multicenter study.

BACKGROUND: Acute cor pulmonale (ACP) is prognostic in patients with acute respiratory distress syndrome (ARDS). Identification of paradoxical septal motion (PSM) using two-dimensional echocardiography is highly subjective. We sought to describe feature-engineered metrics derived from LV radial strain changes related to PSM in ARDS patients with ACP of various severity and to illustrate potential diagnostic and prognostic yield. METHODS: This prospective bicentric study included patients under protective ventilation for ARDS related to COVID-19 who were assessed using transesophageal echocardiography (TEE). Transgastric short-axis view at mid-papillary level was used to visually grade septal motion, using two-dimensional imaging, solely and combined with LV radial strain: normal (grade 0), transient end-systolic septal flattening (grade 1), prolonged end-systolic septal flattening or reversed septal curvature (grade 2). Inter-observer variability was calculated. Feature engineering was performed to calculate the time-to-peak and area under the strain curve in 6 LV segments. In the subset of patients with serial TEE examinations, a multivariate Cox model analysis accounting for new-onset of PSM as a time-dependent variable was used to identify parameters associated with ICU mortality. RESULTS: Overall, 310 TEE examinations performed in 182 patients were analyzed (age: 67 [60-72] years; men: 66%; SAPSII: 35 [29-40]). Two-dimensional assessment identified a grade 1 and grade 2 PSM in 100 (32%) and 48 (15%) examinations, respectively. Inter-rater reliability was weak using two-dimensional imaging alone (kappa = 0.49; 95% CI 0.40-0.58; p < 0.001) and increased with associated LV radial strain (kappa = 0.84, 95% CI 0.79-0.90, p < 0.001). The time-to-peak of mid-septal and mid-lateral segments occurred significantly later in systole and increased with the grade of PSM. Similarly, the area under the strain curve of these segments increased significantly with the grade of PSM, compared with mid-anterior or mid-inferior segments. Severe acute cor pulmonale with a grade 2 PSM was significantly associated with mortality. Requalification in an upper PSM grade using LV radial strain allowed to better identify patients at risk of death (HR: 6.27 [95% CI 2.28-17.2] vs. 2.80 [95% CI 1.11-7.09]). CONCLUSIONS: In objectively depicting PSM and quantitatively assessing its severity, TEE LV radial strain appears as a valuable adjunct to conventional two-dimensional imaging.

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.

Within-subject pharmacokinetic variability has a strong influence on individual exposure ratios in bioequivalence studies, hence on drug formulation interchangeability.

INTRODUCTION: Bioequivalence between a reference and a generic drug is based on the hypothesis that a ± 20% change in blood exposure (or ± 10% for drugs with narrow therapeutic index, NTI) following the generic/reference switch will not have any therapeutic consequences. However, the individual exposure ratio between generic and reference can be higher than 1.20 (or 1.10). This study aims to analyse the different parameters influencing the individual exposure ratio, hence the conditions for reference/generic interchangeability. METHODS: Bioequivalence studies with a double cross-over design for a virtual drug were simulated using 100 random sets of 12, 24, 48 or 100 pairs of areas under the curve (AUC), varying the generic/reference AUC geometric mean ratios between 0.80 and 1.25 and the within-subject exposure variance of the reference and the generic formulations. RESULTS: The proportion of subjects with an exposure generic/reference ratio outside the ± 10% or ± 20% acceptance intervals increases when (1) the reference within-subject variance increases; (2) the ratio of the generic within-subject variance on the reference within-subject variance increases; and (3) the generic/reference mean AUC ratio diverges from 1.0. When only considering replicated administrations of the reference, the individual exposure ratio increases with the within-subject variance, yielding values outside the usually accepted individual exposure ratio range of 0.5 to 2 for drugs with narrow therapeutic index as soon as the within-subject variance standard deviation is ≥ 0.25 (equivalent to within-patient CV% > 25%). CONCLUSIONS: Interchangeability between reference and generic formulations, especially for drugs with narrow therapeutic index can only be assumed if, the within-subject variance of generic is less or equal to the within-subject variance of reference or, if this is not the case, if the distribution of the generic/generic individual exposure ratios is included within the therapeutic margins of the reference drug.

Tacrolimus Exposure Before and After a Switch From Twice-Daily Immediate-Release to Once-Daily Prolonged Release Tacrolimus: The ENVARSWITCH Study.

LCP-tacrolimus displays enhanced oral bioavailability compared to immediate-release (IR-) tacrolimus. The ENVARSWITCH study aimed to compare tacrolimus AUC0-24 h in stable kidney (KTR) and liver transplant recipients (LTR) on IR-tacrolimus converted to LCP-tacrolimus, in order to re-evaluate the 1:0.7 dose ratio recommended in the context of a switch and the efficiency of the subsequent dose adjustment. Tacrolimus AUC0-24 h was obtained by Bayesian estimation based on three concentrations measured in dried blood spots before (V2), after the switch (V3), and after LCP-tacrolimus dose adjustment intended to reach the pre-switch AUC0-24 h (V4). AUC0-24 h estimates and distributions were compared using the bioequivalence rule for narrow therapeutic range drugs (Westlake 90% CI within 0.90-1.11). Fifty-three KTR and 48 LTR completed the study with no major deviation. AUC0-24 h bioequivalence was met in the entire population and in KTR between V2 and V4 and between V2 and V3. In LTR, the Westlake 90% CI was close to the acceptance limits between V2 and V4 (90% CI = [0.96-1.14]) and between V2 and V3 (90% CI = [0.96-1.15]). The 1:0.7 dose ratio is convenient for KTR but may be adjusted individually for LTR. The combination of DBS and Bayesian estimation for tacrolimus dose adjustment may help with reaching appropriate exposure to tacrolimus rapidly after a switch.

Tacrolimus population pharmacokinetics in adult heart transplant patients.

INTRODUCTION: Tacrolimus is an immunosuppressant largely used in heart transplantation. However, the calculation of its exposure based on the area under the curve (AUC) requires the use of a population pharmacokinetic (PK) model. The aims of this work were (i) to develop a population PK model for tacrolimus in heart transplant patients, (ii) to derive a maximum a posteriori Bayesian estimator (MAP-BE) based on a limited sampling strategy (LSS) and (iii) to estimate probabilities of target attainment (PTAs) for AUC and trough concentration (C0). MATERIAL AND METHODS: Forty-seven PK profiles (546 concentrations) of 18 heart transplant patients of the Pharmacocinétique des Immunosuppresseurs chez les patients GREffés Cardiaques study receiving tacrolimus (Prograf®) were included. The database was split into a development (80%) and a validation (20%) set. PK parameters were estimated in MONOLIX® and based on this model a Bayesian estimator using an LSS was built. Simulations were performed to calculate the PTA for AUC and C0. RESULTS: The best model to describe the tacrolimus PK was a two-compartment model with a transit absorption and a linear elimination. Only the CYP3A5 covariate was kept in the final model. The derived MAP-BE based on the LSS (0-1-2 h postdose) yielded an AUC bias ± SD = 2.7 ± 10.2% and an imprecision of 9.9% in comparison to the reference AUC calculated using the trapezoidal rule. PTAs based on AUC or C0 allowed new recommendations to be proposed for starting doses (0.11 mg·kg-1 ·12 h-1 for the CYP3A5 nonexpressor and 0.22 mg·kg1 ·12 h-1 for the CYP3A5 expressor). CONCLUSION: The MAP-BE developed should facilitate estimation of tacrolimus AUC in heart transplant patients.

A Method for Evaluating Robustness of Limited Sampling Strategies-Exemplified by Serum Iohexol Clearance for Determination of Measured Glomerular Filtration Rate.

In combination with Bayesian estimates based on a population pharmacokinetic model, limited sampling strategies (LSS) may reduce the number of samples required for individual pharmacokinetic parameter estimations. Such strategies reduce the burden when assessing the area under the concentration versus time curves (AUC) in therapeutic drug monitoring. However, it is not uncommon for the actual sample time to deviate from the optimal one. In this work, we evaluate the robustness of parameter estimations to such deviations in an LSS. A previously developed 4-point LSS for estimation of serum iohexol clearance (i.e., dose/AUC) was used to exemplify the effect of sample time deviations. Two parallel strategies were used: (a) shifting the exact sampling time by an empirical amount of time for each of the four individual sample points, and (b) introducing a random error across all sample points. The investigated iohexol LSS appeared robust to deviations from optimal sample times, both across individual and multiple sample points. The proportion of individuals with a relative error greater than 15% (P15) was 5.3% in the reference run with optimally timed sampling, which increased to a maximum of 8.3% following the introduction of random error in sample time across all four time points. We propose to apply the present method for the validation of LSS developed for clinical use.