Physiologically based pharmacokinetic modelling of Busulfan: a new approach to describe and predict the pharmacokinetics in adults

Pharmacokinetic Modeling and Simulation with Pharmacogenetic Insights Support the Relevance of Therapeutic Drug Monitoring for Myeloablative Busulfan Dosing in Adult HSCT

Therapeutic drug monitoring (TDM) of busulfan (Bu) is well-established in pediatric hematopoietic stem cell transplantation (HSCT), but its use in adults is limited due to a lack of clear recommendations and scarcity of evidence regarding its utility. GSTA1 promoter variants are reported to affect Bu clearance in both adults and pediatric patients. This study aimed to evaluate the value of preemptive genotyping GSTA1 and body composition (obesity) in individualizing Bu dosing in adults, through pharmacokinetic (PK) modeling and simulations. A population pharmacokinetic (PopPK) model was developed and validated with data from 60 adults who underwent HSCT. Simulations assessed different dosing scenarios based on body size metrics and GSTA1 genotypes. Due to the limited number of obese patients in the cohort, the effect of obesity on Bu pharmacokinetics (PK) was evaluated in silico using a physiologically-based pharmacokinetic (PBPK) model and relevant virtual populations from Simcyp software. Patients with at least 1 GSTA1*B haplotype had 17% lower clearance on average. PopPK simulations indicated that adjusting doses based on genotype increased the probability of achieving the target exposure (3.7 to 5.5 mg.h/L) from 53% to 60 % in GSTA1*A homozygous patients, and from 50% to 61% in *B carriers. Still, Approximately 40% of patients would not achieve this therapeutic window without TDM. A 2-sample optimal design was validated for routine model-based Bu first dose AUC0-∞ estimation, and the model was implemented in the Tucuxi user-friendly TDM software. PBPK simulations confirmed body surface area-based doses of 29 to 31 mg/m2/6h as the most appropriate, regardless of obesity status. This study emphasizes the importance of individualized Bu dosing strategies in adults to achieve therapeutic targets. Preemptive genotyping alone may not have a significant clinical impact, and routine TDM may be necessary for optimal transplantation outcomes.

Limited Sampling Strategies for Estimating Busulfan Area Under the Concentration-Time Curve: Based on Peak and Trough Concentrations in Saliva

Therapeutic drug monitoring for busulfan is currently performed by multiple plasma sampling. Saliva is considered a noninvasive therapeutic drug monitoring matrix. This study aimed to investigate intravenous busulfan pharmacokinetics (PK) in plasma and saliva, and establish a limited sampling strategy (LSS) for predicting the area under the concentration-time curve from time zero to infinity in plasma (AUC0-∞,p) by using saliva samples. Therefore, the PK of busulfan was studied in 37 Chinese patients. Pearson correlation analysis was used to evaluate the correlation between the AUC of busulfan in plasma and saliva. LSS models were established by the multiple linear regression analysis. The prediction error, the mean prediction error, and the root mean square error were used to evaluate the predictive accuracy. The agreement between the predicted and observed AUC0-∞ in saliva was investigated by the intraclass correlation coefficient and Bland-Altman analysis. The accuracy and robustness of the models were evaluated by using the bootstrap procedure. The result of PK analysis 62.2% of patients (23/37) was within the target range of AUC0-∞,p . A good correlation between saliva and plasma busulfan AUC0-∞ was observed (r = 0.63, p < .01). The bias and precision of the models 7 and 13 were less than 15%. The intraclass correlation coefficient exceeded 0.9, and the limits of agreement were within ±15%. The 2-point LSS model in saliva is a convenient and desirable approach to predict the AUC0-∞ of 4 times daily intravenous busulfan in plasma, which can be used to design personalized dosing for busulfan.

Evaluating the Pharmacokinetics of Fentanyl in the Brain Extracellular Fluid, Saliva, Urine, and Plasma of Newborns from Transplacental Exposure from Parturient Mothers Dosed with Epidural Fentanyl Utilizing PBPK Modeling

BACKGROUND AND OBJECTIVE

Fentanyl can mitigate the mother and newborn complications resulting from labor pain. However, fentanyl shows a narrow therapeutic index between its respiratory depressive and analgesic effects. Thus, prenatally acquired high fentanyl levels in the newborn brain extracellular fluid (bECF) may induce respiratory depression which requires therapeutic drug monitoring (TDM). TDM using saliva and urine in newborns can reduce the possibility of infections and distress associated with TDM using blood. The objective of this study was to develop  a physiologically based pharmacokinetic (PBPK) model to predict fentanyl concentrations in different newborn tissues due to intrauterine exposure.

METHODS

A fentanyl PBPK model in adults after intravenous and epidural administration was built, validated, and scaled to pregnancy and newborn populations. The dose that the newborn received transplacentally at birth was calculated using the pregnancy model. Then, the newborn bECF, saliva, plasma, and urine concentrations after such a dose were predicted using the newborn PBPK model.

RESULTS

After a maternal epidural dose of fentanyl 245 µg, the predicted newborn plasma and bECF levels were below the toxicity thresholds. Furthermore, the salivary threshold levels in newborns for fentanyl analgesic and respiratory depression effects were estimated to be 0.39 and 14.7-18.2 ng/ml, respectively.

CONCLUSION

The salivary TDM of fentanyl in newborns can be useful in newborns exposed to intrauterine exposure from parturient females dosed with epidural fentanyl. However, newborn-specific values of µ-opioid receptors IC50 for respiratory depression are needed.

The Analysis of Pethidine Pharmacokinetics in Newborn Saliva, Plasma, and Brain Extracellular Fluid After Prenatal Intrauterine Exposure from Pregnant Mothers Receiving Intramuscular Dose Using PBPK Modeling

BACKGROUND AND OBJECTIVE

Pethidine (meperidine) can decrease labor pain-associated mother's hyperventilation and high cortisol-induced newborn complications. However, prenatal transplacentally acquired pethidine can cause side effects in newborns. High pethidine concentrations in the newborn brain extracellular fluid (bECF) can cause a serotonin crisis. Therapeutic drug monitoring (TDM) in newborns' blood distresses them and increases infection incidence, which can be overcome by using salivary TDM. Physiologically based pharmacokinetic (PBPK) modeling can predict drug concentrations in newborn plasma, saliva, and bECF after intrauterine pethidine exposure.

METHODS

A healthy adult PBPK model was constructed, verified, and scaled to newborn and pregnant populations after intravenous and intramuscular pethidine administration. The pregnancy PBPK model was used to predict the newborn dose received transplacentally at birth, which was used as input to the newborn PBPK model to predict newborn plasma, saliva, and bECF pethidine concentrations and set correlation equations between them.

RESULTS

Pethidine can be classified as a Salivary Excretion Classification System class II drug. The developed PBPK model predicted that, after maternal pethidine intramuscular doses of 100 mg and 150 mg, the newborn plasma and bECF concentrations were below the toxicity thresholds. Moreover, it was estimated that newborn saliva concentrations of 4.7 µM, 11.4 µM, and 57.7 µM can be used as salivary threshold concentrations for pethidine analgesic effects, side effects, and the risk for serotonin crisis, respectively, in newborns.

CONCLUSION

It was shown that saliva can be used for pethidine TDM in newborns during the first few days after delivery to mothers receiving pethidine.

Drug Exposure to Establish Pharmacokinetic-Response Relationships in Oncology

In the oncology field, understanding the relationship between the dose administered and the exerted effect is particularly important because of the narrow therapeutic index associated with anti-cancer drugs and the high interpatient variability. Therefore, in this review, we provide a critical perspective of the different methods of characterising treatment exposure in the oncology setting. The increasing number of modelling applications in oncology reflects the applicability and the impact of pharmacometrics on all phases of the drug development process and patient management as well. Pharmacometric modelling is a worthy component within the current paradigm of model-based drug development, but pharmacometric modelling techniques are also accessible for the clinician in the optimisation of current oncology therapies. Consequently, the application of population models in a hospital setting by generating close collaborations between physicians and pharmacometricians is highly recommended, providing a systematic means of developing and assessing model-based metrics as 'drivers' for various responses to treatments, which can then be evaluated as predictors for treatment success. Characterising the key determinants of variability in exposure is of particular importance for anticancer agents, as efficacy and toxicity are associated with exposure. We present the different strategies to describe and predict drug exposure that can be applied depending on the data available, with the objective of obtaining the most useful information in the patients' favour throughout the full drug cycle. Therefore, the objective of the present article is to review the different approaches used to characterise a patient's exposure to oncology drugs, which will result in a better understanding of the time course of the response and the magnitude of interpatient variability.

Physiologically-based pharmacokinetic model for alectinib, ruxolitinib, and panobinostat in the presence of cancer, renal impairment, and hepatic impairment

Renal (RIP) and hepatic (HIP) impairments are prevalent conditions in cancer patients. They can cause changes in gastric emptying time, albumin levels, hematocrit, glomerular filtration rate, hepatic functional volume, blood flow rates, and metabolic activity that can modify drug pharmacokinetics. Performing clinical studies in such populations has ethical and practical issues. Using predictive physiologically-based pharmacokinetic (PBPK) models in the evaluation of the PK of alectinib, ruxolitinib, and panobinostat exposures in the presence of cancer, RIP, and HIP can help in using optimal doses with lower toxicity in these populations. Verified PBPK models were customized under scrutiny to account for the pathophysiological changes induced in these diseases. The PBPK model-predicted plasma exposures in patients with different health conditions within average 2-fold error. The PBPK model predicted an area under the curve ratio (AUCR) of 1, and 1.8, for ruxolitinib and panobinostat, respectively, in the presence of severe RIP. On the other hand, the severe HIP was associated with AUCR of 1.4, 2.9, and 1.8 for alectinib, ruxolitinib, and panobinostat, respectively, in agreement with the observed AUCR. Moreover, the PBPK model predicted that alectinib therapeutic cerebrospinal fluid levels are achieved in patients with non-small cell lung cancer, moderate HIP, and severe HIP at 1-, 1.5-, and 1.8-fold that of healthy subjects. The customized PBPK models showed promising ethical alternatives for simulating clinical studies in patients with cancer, RIP, and HIP. More work is needed to quantify other pathophysiological changes induced by simultaneous affliction by cancer and RIP or HIP.

Use of normalized prediction distribution errors for assessing population physiologically-based pharmacokinetic model adequacy

Currently employed methods for qualifying population physiologically-based pharmacokinetic (Pop-PBPK) model predictions of continuous outcomes (e.g., concentration-time data) fail to account for within-subject correlations and the presence of residual error. In this study, we propose a new method for evaluating Pop-PBPK model predictions that account for such features. The approach focuses on deriving Pop-PBPK-specific normalized prediction distribution errors (NPDE), a metric that is commonly used for population pharmacokinetic model validation. We describe specific methodological steps for computing NPDE for Pop-PBPK models and define three measures for evaluating model performance: mean of NPDE, goodness-of-fit plots, and the magnitude of residual error. Utility of the proposed evaluation approach was demonstrated using two simulation-based study designs (positive and negative control studies) as well as pharmacokinetic data from a real-world clinical trial. For the positive-control simulation study, where observations and model simulations were generated under the same Pop-PBPK model, the NPDE-based approach denoted a congruency between model predictions and observed data (mean of NPDE =  - 0.01). In contrast, for the negative-control simulation study, where model simulations and observed data were generated under different Pop-PBPK models, the NPDE-based method asserted that model simulations and observed data were incongruent (mean of NPDE =  - 0.29). When employed to evaluate a previously developed clindamycin PBPK model against prospectively collected plasma concentration data from 29 children, the NPDE-based method qualified the model predictions as successful (mean of NPDE = 0). However, when pediatric subpopulations (e.g., infants) were evaluated, the approach revealed potential biases that should be explored.

Utility of physiologically based pharmacokinetic (PBPK) modeling in oncology drug development and its accuracy: a systematic review

PURPOSE

Physiologically based pharmacokinetic (PBPK) modeling, a mathematical modeling approach which uses a pharmacokinetic model to mimick human physiology to predict drug concentration-time profiles, has been used for the discover and development of drugs in various fields, including oncology, since 2000. There have been a few general review articles on the utilization of PBPK in the development of oncology drugs, but these do not include an evaluation of model prediction accuracy. We therefore conducted a systematic review to define the accuracy of PBPK model prediction and its utility throughout all the developmental phases of oncology drugs.

METHODS

A systematic search was performed in the PubMed, PubMed Central and Cochrane Library databases from 1980 to February 2017 for articles (1) written in English, (2) focused on oncology or antineoplastic or anticancer drugs, tumor or cancer or anticancer drugs listed in the U.S. National Institutes of Health and (3) involving a PBPK model. The absolute-average-folding-errors (AAFEs) of the area under the curve (AUC) between predicted and observed values in each article were calculated to assess model prediction accuracy.

RESULTS

Of the 2341 articles initially identified by our search of the databases, 40 were included in the review analysis. These articles reported on six types of studies, i.e. in vivo (n = 4), first-in-human (n = 5), phase II/III clinical trials (n = 9), organ impairment (n = 3), pediatrics (n = 4) and drug-drug interactions (n = 15). AAFEs of the predicted AUC for all groups of studies were within 1.3-fold of each other despite variations in experimental methodologies.

CONCLUSION

PBPK modeling is a potential tool which can be effectively applied throughout all phases of oncology drug development. The number of experimental animals and human participants enrolled in the studies can be reduced using PBPK modeling and PBPK-population-PK modeling. The limited number of publications of unsuccessful model application to date may contribute to bias toward the usefulness of modeling.

GSTA1 diplotypes affect busulfan clearance and toxicity in children undergoing allogeneic hematopoietic stem cell transplantation: a multicenter study

Busulfan (BU) dose adjustment following therapeutic drug monitoring contributes to better outcome of hematopoietic stem cell transplantation (HSCT). Further improvement could be achieved through genotype-guided BU dose adjustments. To investigate this aspect, polymorphism within glutathione S transferase genes were assessed. Particularly, promoter haplotypes of the glutathione S transferase A1 (GSTA1) were evaluated in vitro, with reporter gene assays and clinically, in a pediatric multi-center study (N =138) through association with BU pharmacokinetics (PK) and clinical outcomes. Promoter activity significantly differed between the GSTA1 haplotypes (p<0.001) supporting their importance in capturing PK variability. Four GSTA1 diplotype groups that significantly correlated with clearance (p=0.009) were distinguished. Diplotypes underlying fast and slow metabolizing capacity showed higher and lower BU clearance (ml/min/kg), respectively. GSTA1 diplotypes with slow metabolizing capacity were associated with higher incidence of sinusoidal obstruction syndrome, acute graft versus host disease and combined treatment-related toxicity (p<0.0005). Among other GST genes investigated, GSTP1 313GG correlated with acute graft versus host disease grade 1-4 (p=0.01) and GSTM1 non-null genotype was associated with hemorrhagic cystitis (p=0.003). This study further strengthens the hypothesis that GST diplotypes/genotypes could be incorporated into already existing population pharmacokinetic models for improving first BU dose prediction and HSCT outcomes. (N^o Clinicaltrials.gov identifier: NCT01257854. Registered 8 December 2010, retrospectively registered).

The busulfan metabolite EdAG irreversibly glutathionylates glutaredoxins

The DNA alkylating agent busulfan is used to 'precondition' patients with leukemia, lymphomas and other hematological disorders prior to hematopoietic stem cell transplants. Busulfan is metabolized via conjugation with glutathione (GSH) followed by intramolecular rearrangement to the GSH analog γ-glutamyl-dehydroalanyl -glycine (EdAG). EdAG contains the electrophilic dehydroalanine, which is expected to react with protein nucleophiles, particularly proteins with GSH binding sites such as glutaredoxins (Grx's). Incubation of EdAG with human Grx-1 or Grx-2 results in facile adduction of cys-23 and cys-77, respectively, as determined by ESI-MS/MS. The resulting modified proteins are catalytically inactive. In contrast, the glutathione transferase A1-1 includes a GSH binding site with a potentially reactive tyrosinate (Tyr-9) but it does not react with EdAG. Similarly, Cys-112 of GSTA1-1, which lies outside the active site and is known to form disulfides with GSH, does not react with EdAG. The results provide the first demonstration of the reactivity of any busulfan metabolites with intact proteins, and they suggest that GSH-binding sites containing thiolates are most susceptible. The adduction of Grx's by EdAG suggests the possible alteration of proteins that are normally regulated via Grx-dependent reversible glutathionylation or deglutathionylation. Dysregulation of Grx-dependent processes could contribute to cellular toxicity of busulfan.

Predictive performance of a physiologically based pharmacokinetic model of busulfan in children

A physiologically based pharmacokinetic (PBPK) model of the DNA-alkylating agent busulfan was slightly modified and scaled from adults to children in order to predict the systemic busulfan drug exposure in children. Capitalizing on the recent major software release of PK-Sim®, we refined our PBPK model by implementing glutathione S transferase (GST) in 11 organs using the software integrated enzyme expression database. In addition, two irreversible binding processes (i.e., DNA and plasma protein binding) were applied by using Koff and KD values. The model was scaled from adults to children. Simulations were computed and compared to concentration-time data after intravenous (i.v.) busulfan administration to 36 children. Based on the results, an age-dependent enzyme activity and maturation ratio was tailored and evaluated with an external dataset consisting of 23 children. Initial adult to children scaling indicated lower clearance values for children in comparison to adults. Subsequent age-dependent maturation ratio resulted in three different age groups: Activity of busulfan-glutathione conjugate formation was 80%, 61%, and 89% in comparison to adults for children with an age of up to 2 years, > 2-6 years, and > 6-18 years, respectively. Patients of the evaluation dataset were simulated with a mean percentage error (MPE) for all patients of 3.9% with 3/23 children demonstrating a MPE of > ±30%. The PBPK model parameterization sufficiently described the observed concentration-time data of the validation dataset while showing an adequate predictive performance. This PBPK model could be helpful to determine the first dose of busulfan in children.