Pharmacometabonomic Prediction of Busulfan Clearance in Hematopoetic Cell Transplant Recipients

The presence of busulfan metabolites and pharmacometabolomics in plasma drawn immediately before allograft infusion in hematopoietic cell transplant recipients

Busulfan is hepatically metabolized through glutathione (GSH) conjugation; in vitro, this process depletes hepatocyte GSH stores and generates the cytotoxic metabolite γ-glutamyldehydroalanylglycine, which is too unstable to be quantitated in vivo. We sought to evaluate if pre-graft (i.e., immediately before allograft infusion) concentrations of busulfan metabolites' and of endogenous metabolomic compounds (EMCs) representing the glutathione pathway were associated with clinical outcomes in hematopoietic cell transplant (HCT) recipients receiving busulfan. The clinical outcomes evaluated were relapse, acute graft versus host disease (GVHD), chronic GVHD, non-relapse mortality, and neutrophil nadir. In pre-graft samples obtained from patients immediately before allograft infusion, our objectives were to evaluate for: (1) the presence of busulfan and its metabolites tetrahydrothiophenium ion (THT+), tetrahydrothiophene 1-oxide, sulfolane, and 3-hydroxysulfolane (N = 124); (2) EMCs using a global metabolomics assay (N = 77); and (3) the association of the busulfan metabolites and the EMCs with clinical outcomes. In the pre-graft samples, busulfan and THT+ could not be detected. THT 1-oxide, sulfolane, and 3-hydroxysulfolane were quantitated in 9.6%, 26%, and 58% of pre-graft samples; their concentrations were not associated with clinical outcomes. Four pre-graft EMCs were statistically significantly associated with the neutrophil nadir. The pre-graft EMCs were not associated with the other clinical outcomes. In conclusion, busulfan's metabolites are present in patients' plasma immediately before allograft infusion; the neutrophil nadir is associated with pre-graft EMCs. Future research should investigate the association of clinical outcomes with the concentrations of busulfan's metabolites and EMCs in the pre-graft plasma from allogeneic HCT recipients.

Pharmacokinetics in Pharmacometabolomics: Towards Personalized Medication

Indiscriminate drug administration may lead to drug therapy results with varying effects on patients, and the proposal of personalized medication can help patients to receive effective drug therapy. Conventional ways of personalized medication, such as pharmacogenomics and therapeutic drug monitoring (TDM), can only be implemented from a single perspective. The development of pharmacometabolomics provides a research method for the realization of precise drug administration, which integrates the environmental and genetic factors, and applies metabolomics technology to study how to predict different drug therapeutic responses of organisms based on baseline metabolic levels. The published research on pharmacometabolomics has achieved satisfactory results in predicting the pharmacokinetics, pharmacodynamics, and the discovery of biomarkers of drugs. Among them, the pharmacokinetics related to pharmacometabolomics are used to explore individual variability in drug metabolism from the level of metabolism of the drugs in vivo and the level of endogenous metabolite changes. By searching for relevant literature with the keyword "pharmacometabolomics" on the two major literature retrieval websites, PubMed and Web of Science, from 2006 to 2023, we reviewed articles in the field of pharmacometabolomics that incorporated pharmacokinetics into their research. This review explains the therapeutic effects of drugs on the body from the perspective of endogenous metabolites and pharmacokinetic principles, and reports the latest advances in pharmacometabolomics related to pharmacokinetics to provide research ideas and methods for advancing the implementation of personalized medication.

Precision sirolimus dosing in children: The potential for model-informed dosing and novel drug monitoring

The mTOR inhibitor sirolimus is prescribed to treat children with varying diseases, ranging from vascular anomalies to sporadic lymphangioleiomyomatosis to transplantation (solid organ or hematopoietic cell). Precision dosing of sirolimus using therapeutic drug monitoring (TDM) of sirolimus concentrations in whole blood drawn at the trough (before the next dose) time-point is the current standard of care. For sirolimus, trough concentrations are only modestly correlated with the area under the curve, with R 2 values ranging from 0.52 to 0.84. Thus, it should not be surprising, even with the use of sirolimus TDM, that patients treated with sirolimus have variable pharmacokinetics, toxicity, and effectiveness. Model-informed precision dosing (MIPD) will be beneficial and should be implemented. The data do not suggest dried blood spots point-of-care sampling of sirolimus concentrations for precision dosing of sirolimus. Future research on precision dosing of sirolimus should focus on pharmacogenomic and pharmacometabolomic tools to predict sirolimus pharmacokinetics and wearables for point-of-care quantitation and MIPD.

Prediction of Busulfan Clearance by Predose Plasma Metabolomic Profiling

Intravenous busulfan doses are often personalized to a target plasma exposure (targeted busulfan) using an individual's busulfan clearance (BuCL). We evaluated whether BuCL could be predicted by a predose plasma panel of 841 endogenous metabolomic compounds (EMCs). In this prospective cohort of 132 hematopoietic cell transplantation (HCT) patients, all had samples collected immediately before busulfan administration (preBU) and 96 had samples collected 2 weeks before busulfan (2-week-preBU). BuCL was significantly associated with 37 EMCs after univariate linear regression analysis and controlling for false discovery (< 0.05) in the 132 preBU samples. In parallel, with preBU samples, we included all 841 EMCs in a least absolute shrinkage and selection operator-penalized regression which selected 13 EMCs as predominantly associated with BuCL. Then, we constructed a prediction model by estimating coefficients for these 13 EMCs, along with sex, using ordinary least-squares. When the resulting linear prediction model was applied to the 2-week-preBU samples, it explained 40% of the variation in BuCL (adjusted R²  = 0.40). Pathway enrichment analysis revealed 18 pathways associated with BuCL. Lysine degradation followed by steroid biosynthesis, which aligned with the univariate analysis, were the top two pathways. BuCL can be predicted before busulfan administration with a linear regression model of 13 EMCs. This pharmacometabolomics method should be prioritized over use of a busulfan test dose or pharmacogenomics to guide busulfan dosing. These results highlight the potential of pharmacometabolomics as a precision medicine tool to improve or replace pharmacokinetics to personalize busulfan doses.

Pharmacometabonomic Association of Cyclophosphamide 4-hydroxylation in Hematopoietic Cell Transplant Recipients

The widely used alkylating agent cyclophosphamide (CY) has substantive interpatient variability in the area under the curve (AUC) of it and its metabolites. Numerous factors may influence the drug‐metabolizing enzymes that metabolize CY to 4‐hydroxycyclophosphamide (4HCY), the principal precursor to CY’s cytotoxic metabolite. We sought to identify endogenous metabolomics compounds (EMCs) associated with 4HCY formation clearance (ratio of 4HCY/CY AUC) using global metabolomics. Patients who undergo hematopoietic cell transplantation receiving post‐transplant CY (PT‐CY) were enrolled, cohort 1 (n = 26) and cohort 2 (n = 25) donating longitudinal blood samples before they started HCT (pre‐HCT), before infusion of the donor allograft (pre‐graft), before the first dose of PT‐CY (pre‐CY), and 24 h after the first dose of PT‐CY (24‐h post‐CY), which is also immediately before the second dose of CY. A total of 512 and 498 EMCs were quantitated in two cohorts, respectively. Both univariate linear regression with false discovery rate (FDR), and pathway enrichment analyses using a global association test were performed. At the pre‐CY time point, no EMCs were associated at FDR less than 0.1. At pre‐HCT, cohort 1 had one EMC (levoglucosan) survive the FDR threshold. At pre‐graft, cohort 1 and cohort 2 had 20 and 13 EMCs, respectively, exhibiting unadjusted p values less than 0.05, with the only EMCs having an FDR less than 0.1 being two unknown EMCs. At 24‐h post‐CY, there were three EMCs, two ketones, and threitol, at FDR less than 0.1 in cohort 2. These results demonstrate the potential of pharmacometabonomics, but future studies in larger samples are needed to optimize CY.

Multi-omics profiling: the way towards precision medicine in metabolic diseases

Metabolic diseases including type 2 diabetes mellitus (T2DM), non-alcoholic fatty liver disease (NAFLD), and metabolic syndrome (MetS) are alarming health burdens around the world, while therapies for these diseases are far from satisfying as their etiologies are not completely clear yet. T2DM, NAFLD, and MetS are all complex and multifactorial metabolic disorders based on the interactions between genetics and environment. Omics studies such as genetics, transcriptomics, epigenetics, proteomics, and metabolomics are all promising approaches in accurately characterizing these diseases. And the most effective treatments for individuals can be achieved via omics pathways, which is the theme of precision medicine. In this review, we summarized the multi-omics studies of T2DM, NAFLD, and MetS in recent years, provided a theoretical basis for their pathogenesis and the effective prevention and treatment, and highlighted the biomarkers and future strategies for precision medicine.

Feasibility of pharmacometabolomics to identify potential predictors of paclitaxel pharmacokinetic variability

PURPOSE

Paclitaxel is a commonly used chemotherapy drug with substantial variability in pharmacokinetics (PK) that affects treatment efficacy and toxicity. Pharmacometabolomic signatures that explain PK variability could be used to individualize dosing to improve therapeutic outcomes. The objective of this study was to identify pretreatment metabolites or metabolomic signatures that explain variability in paclitaxel PK.

METHODS

This analysis was conducted using data previously collected on a prospective observational study of 48 patients with breast cancer receiving weekly 80 mg/m² paclitaxel infusions. Paclitaxel plasma concentrations were measured during the first infusion to estimate paclitaxel time above threshold (T_c>0.05) and maximum concentration (C_max). Metabolites measured in pretreatment whole blood by nuclear magnetic resonance spectrometry were analyzed for an association with T_c>0.05 and C_max using Pearson correlation followed by stepwise linear regression.

RESULTS

Pretreatment creatinine, glucose, and lysine concentrations were positively correlated with T_c>0.05, while pretreatment betaine was negatively correlated and lactate was positively correlated with C_max (all uncorrected p < 0.05). After stepwise elimination, creatinine was associated with T_c>0.05, while betaine and lactate were associated with C_max (all p < 0.05).

CONCLUSION

This study identified pretreatment metabolites that may be associated with paclitaxel PK variability demonstrating feasibility of a pharmacometabolomics approach for understanding paclitaxel PK. However, identification of more robust pharmacometabolomic predictors will be required for broad and routine application for the clinical dosing of paclitaxel.

Short Report: Using Targeted Urine Metabolomics to Distinguish Between Manganese Exposed and Unexposed Workers in a Small Occupational Cohort

Objectives: Despite the widespread use of manganese (Mn) in industrial settings and its association with adverse neurological outcomes, a validated and reliable biomarker for Mn exposure is still elusive. Here, we utilize targeted metabolomics to investigate metabolic differences between Mn-exposed and -unexposed workers, which could inform a putative biomarker for Mn and lead to increased understanding of Mn toxicity. Methods: End of shift spot urine samples collected from Mn exposed (n = 17) and unexposed (n = 15) workers underwent a targeted assay of 362 metabolites using LC-MS/MS; 224 were quantified and retained for analysis. Differences in metabolite abundances between exposed and unexposed workers were tested with a Benjamini-Hochberg adjusted Wilcoxon Rank-Sum test. We explored perturbed pathways related to exposure using a pathway analysis. Results: Seven metabolites were significantly differentially abundant between exposed and unexposed workers (FDR ≤ 0.1), including n-isobutyrylglycine, cholic acid, anserine, beta-alanine, methionine, n-isovalerylglycine, and threonine. Three pathways were significantly perturbed in exposed workers and had an impact score >0.5: beta-alanine metabolism, histidine metabolism, and glycine, serine, and threonine metabolism. Conclusion: This is one of few studies utilizing targeted metabolomics to explore differences between Mn-exposed and -unexposed workers. Metabolite and pathway analysis showed amino acid metabolism was perturbed in these Mn-exposed workers. Amino acids have also been shown to be perturbed in other occupational cohorts exposed to Mn. Additional research is needed to characterize the biological importance of amino acids in the Mn exposure-disease continuum, and to determine how to appropriately utilize and interpret metabolomics data collected from occupational cohorts.

Pharmacometabonomics: data processing and statistical analysis

Individual variations in drug efficacy, side effects and adverse drug reactions are still challenging that cannot be ignored in drug research and development. The aim of pharmacometabonomics is to better understand the pharmacokinetic properties of drugs and monitor the drug effects on specific metabolic pathways. Here, we systematically reviewed the recent technological advances in pharmacometabonomics for better understanding the pathophysiological mechanisms of diseases as well as the metabolic effects of drugs on bodies. First, the advantages and disadvantages of all mainstream analytical techniques were compared. Second, many data processing strategies including filtering, missing value imputation, quality control-based correction, transformation, normalization together with the methods implemented in each step were discussed. Third, various feature selection and feature extraction algorithms commonly applied in pharmacometabonomics were described. Finally, the databases that facilitate current pharmacometabonomics were collected and discussed. All in all, this review provided guidance for researchers engaged in pharmacometabonomics and metabolomics, and it would promote the wide application of metabolomics in drug research and personalized medicine.

Persistent metabolomic alterations characterize chronic critical illness after severe trauma

BACKGROUND

Following trauma, persistent inflammation, immunosuppression, and catabolism may characterize delayed recovery or failure to recover. Understanding the metabolic response associated with these adverse outcomes may facilitate earlier identification and intervention. We characterized the metabolic profiles of trauma victims who died or developed chronic critical illness (CCI) and hypothesized that differences would be evident within 1-week postinjury.

METHODS

Venous blood samples from trauma victims with shock who survived at least 7 days were analyzed using mass spectrometry. Subjects who died or developed CCI (intensive care unit length of stay of ≥14 days with persistent organ dysfunction) were compared with subjects who recovered rapidly (intensive care unit length of stay, ≤7 days) and uninjured controls. We used partial least squares discriminant analysis, t tests, linear mixed effects regression, and pathway enrichment analyses to make broad comparisons and identify differences in metabolite concentrations and pathways.

RESULTS

We identified 27 patients who died or developed CCI and 33 who recovered rapidly. Subjects were predominantly male (65%) with a median age of 53 years and Injury Severity Score of 36. Healthy controls (n = 48) had similar age and sex distributions. Overall, from the 163 metabolites detected in the samples, 56 metabolites and 21 pathways differed between injury outcome groups, and partial least squares discriminant analysis models distinguished injury outcome groups as early as 1-day postinjury. Differences were observed in tryptophan, phenylalanine, and tyrosine metabolism; metabolites associated with oxidative stress via methionine metabolism; inflammatory mediators including kynurenine, arachidonate, and glucuronic acid; and products of the gut microbiome including indole-3-propionate.

CONCLUSIONS

The metabolic profiles in subjects who ultimately die or develop CCI differ from those who have recovered. In particular, we have identified differences in markers of inflammation, oxidative stress, amino acid metabolism, and alterations in the gut microbiome. Targeted metabolomics has the potential to identify important metabolic changes postinjury to improve early diagnosis and targeted intervention.

LEVEL OF EVIDENCE

Prognostic/epidemiologic, level III.

Longitudinal Pharmacometabonomics for Predicting Malignant Tumor Patient Responses to Anlotinib Therapy: Phenotype, Efficacy, and Toxicity

Anlotinib is an oral small molecule inhibitor of multiple receptor tyrosine kinases (RTKs), which was approved by the National Medical Products Administration (NMPA) of China in 2018 for the third-line treatment of non-small cell lung cancer (NSCLC). Here, for the first time, the longitudinal pharmacometabonomics was explored for predicting malignant tumor patient responses to anlotinib, including the metabolic phenotype variation, drug efficacy, and toxicity. A total of 393 plasma samples from 16 subjects collected from a phase I additional study of anlotinib (NCT02752516) were submitted to targeted metabolomics analysis. The orthogonal partial least-squares discriminant analysis (OPLS-DA) models were constructed for the predication of anlotinib efficacy and toxicity based on the longitudinal pharmacometabonomics data. Statistical results showed that 38 metabolites, mainly involved in aminoacyl-tRNA biosynthesis, alanine, aspartate, and glutamate metabolism, and steroid hormone biosynthesis, were all significantly upregulated attributing to anlotinib treatment. The anti-tumor efficacy and occurrence of proteinuria after anlotinib administration can be predicted with 100% accuracy using the established OPLS-DA models. Glycodeoxycholic acid and glycocholic acid possessed the most excellent sensitivity and specificity in predicting the efficacy of anlotinib, with area under the receiver operating characteristic curve (AUC of ROC curve) 0.847 and 0.828, respectively. NG, NG-dimethylarginine was the most promising biomarker for the prediction of proteinuria occurrence after anlotinib administration, with AUC of ROC curve 0.814. In conclusion, this work developed efficient and convenient discriminant models that can accurately predict the efficacy and toxicity of anlotinib based on longitudinal pharmacometabonomics study.

Prediction of Acute Graft versus Host Disease and Relapse by Endogenous Metabolomic Compounds in Patients Receiving Personalized Busulfan-Based Conditioning

Busulfan-based conditioning is the most commonly used high-dose conditioning regimen for allogeneic hematopoietic cell transplant (HCT). The alkylating agent busulfan has a narrow therapeutic index, with busulfan doses personalized to a target plasma exposure (targeted busulfan). Using a global pharmacometabonomics approach, we sought to identify novel biomarkers of relapse or acute graft versus host disease (GVHD) in a cohort of 84 patients receiving targeted busulfan before allogeneic HCT. A total of 763 endogenous metabolomic compounds (EMCs) were quantitated in 230 longitudinal blood samples before, during, and shortly after intravenous busulfan administration. We performed both univariate linear regression and pathway enrichment analyses using global testing. The cysteine/methionine pathway and the glycine, serine, and threonine metabolism pathway were most associated with relapse. The latter be explained by the fact that glutathione S-transferases conjugate both busulfan and glutathione, which contains glycine as a component. The d-arginine and d-ornithine metabolism pathway and arginine and proline metabolism pathway were most associated with acute GVHD. None of these associations were significant after correcting for false discovery rate (FDR) with a strict cutoff of FDR-adjusted p < 0.1. Although larger studies are needed to substantiate these findings, the results show that EMCs may be used as predictive biomarkers in HCT patients.

A semi-mechanistic model based on glutathione depletion to describe intra-individual reduction in busulfan clearance

Aim

To develop a semi‐mechanistic model, based on glutathione depletion and predict a previously identified intra‐individual reduction in busulfan clearance to aid in more precise dosing.

Methods

Busulfan concentration data, measured as part of regular care for allogeneic hematopoietic cell transplantation (HCT) patients, were used to develop a semi‐mechanistic model and compare it to a previously developed empirical model. The latter included an empirically estimated time effect, where the semi‐mechanistic model included theoretical glutathione depletion. As older age has been related to lower glutathione levels, this was tested as a covariate in the semi‐mechanistic model. Lastly, a therapeutic drug monitoring (TDM) simulation was performed comparing the two models in target attainment.

Results

In both models, a similar clearance decrease of 7% (range −82% to 44%), with a proportionality to busulfan metabolism, was found. After 40 years of age, the time effect increased with 4% per year of age (0.6–8%, P = 0.009), causing the effect to increase more than a 2‐fold over the observed age‐range (0–73 years). Compared to the empirical model, the final semi‐mechanistic model increased target attainment from 74% to 76%, mainly through better predictions for adult patients.

Conclusion

These results suggest that the time‐dependent decrease in busulfan clearance may be related to gluthathione depletion. This effect increased with older age (>40 years) and was proportional to busulfan metabolism. The newly constructed semi‐mechanistic model could be used to further improve TDM‐guided exposure target attainment of busulfan in patients undergoing HCT.

Pharmacometabonomics: The Prediction of Drug Effects Using Metabolic Profiling

Metabonomics, also known as metabolomics, is concerned with the study of metabolite profiles in humans, animals, plants and other systems in order to assess their health or other status and their responses to experimental interventions. Metabonomics is thus widely used in disease diagnosis and in understanding responses to therapies such as drug administration. Pharmacometabonomics, also known as pharmacometabolomics, is a related methodology but with a prognostic as opposed to diagnostic thrust. Pharmacometabonomics aims to predict drug effects including efficacy, safety, metabolism and pharmacokinetics, prior to drug administration, via an analysis of pre-dose metabolite profiles. This article will review the development of pharmacometabonomics as a new field of science that has much promise in helping to deliver more effective personalised medicine, a major goal of twenty-first century healthcare.

NMR-based pharmacometabonomics: A new paradigm for personalised or precision medicine

Metabolic profiling by NMR spectroscopy or hyphenated mass spectrometry, known as metabonomics or metabolomics, is an important tool for systems-based approaches in biology and medicine. The experiments are typically done in a diagnostic fashion where changes in metabolite profiles are interpreted as a consequence of an intervention or event; be that a change in diet, the administration of a drug, physical exertion or the onset of a disease. By contrast, pharmacometabonomics takes a prognostic approach to metabolic profiling, in order to predict the effects of drug dosing before it occurs. Differences in pre-dose metabolite profiles between groups of subjects are used to predict post-dose differences in response to drug administration. Thus the paradigm is inverted and pharmacometabonomics is the metabolic equivalent of pharmacogenomics. Although the field is still in its infancy, it is expected that pharmacometabonomics, alongside pharmacogenomics, will assist with the delivery of personalised or precision medicine to patients, which is a critical goal of 21st century healthcare.

Pharmacometabolomics identifies dodecanamide and leukotriene B4 dimethylamide as a predictor of chemosensitivity for patients with acute myeloid leukemia treated with cytarabine and anthracycline

Clinical responses to standard cytarabine plus anthracycline regimen in acute myeloid leukemia (AML) are heterogeneous and there is an unmet need for biological predictors of response to this regimen. Here, we applied a pharmacometabolomics approach to identify potential biomarkers associated with response to this regimen in AML patients. Based on clinical response the enrolled 82 patients were subdivided into two groups: complete remission(CR) responders (n=42) and non-responders (n=40). Metabolic profiles of pre-treatment serum from patients were analyzed by ultra-high performance liquid chromatography coupled with mass spectrometry and the metabolic differences between the two groups were investigated by multivariate statistical analysis. A metabolite panel containing dodecanamide and leukotriene B4 dimethylamide (LTB4-DMA) had the power capacity to differentiate the two groups of patients, yielding an area under the receiver operating characteristic of 0.945 (85.2% sensitivity and 88.9% specificity) in the training set and 0.944(84.6% sensitivity and 80.0% specificity) in the test set. The patients with high levels of LTB4-DMA and low amounts of dodecanamide had good sensitivity to chemotherapeutic agents. The possible reasons were that dodecanamide was produced by leukemic cells as a lipolytic factor to fuel their growth with a potential role in drug resistance and LTB4-DMA was a potent leukotriene B4 antagonist that could be applicable in the treatment of AML. These preliminary results demonstrates the feasibility of relating chemotherapy responses with pre-treatment metabolic profiles of AML patients, and pharmacometabolomics may be a useful tool to select patients that are more likely to benefit from cytarabine plus anthracycline chemotherapy.