Population Pharmacokinetics of High-Dose Methotrexate in Patients With Primary Central Nervous System Lymphoma

Population Pharmacokinetic and Toxicity Analysis of High-Dose Methotrexate in Patients with Central Nervous System Lymphoma

BACKGROUND

High-dose methotrexate (HD-MTX)-based polychemotherapy is widely used for patients with central nervous system (CNS) lymphoma. The pharmacokinetic (PK) variability and unpredictable occurrence of toxicity remain major concerns in HD-MTX treatment.

OBJECTIVES

This study aimed to characterize the population PK of HD-MTX in patients with CNS lymphoma and to identify baseline predictors and exposure thresholds that predict a high risk of nephro- and hepatotoxicity.

METHODS

Routinely monitored serum MTX concentrations after intravenous infusion of HD-MTX and MTX dosing information were collected retrospectively. Acute event of toxicity (≥ grade 1) was defined according to the Common Terminology Criteria for Adverse Events (CTCAE) version 5.0 on the basis of serum creatinine and alanine aminotransferase. A population PK model was developed in NONMEM. Toxicity data were analyzed using a logistic regression model, and potential baseline and exposure-related predictors were investigated.

RESULTS

In total, 1584 MTX concentrations from 110 patients were available for analysis. A two-compartment population PK model adequately described the data. Estimated glomerular filtration rate (eGFR), treatment regimen, albumin, alkaline phosphatase, and body weight were identified as significant covariates that explain the PK variability of HD-MTX. Baseline eGFR and sex were identified as significant predictors for renal toxicity, and MTX dose (mg/m2) was the strongest predictor for hepatotoxicity. The MTX area under the concentration-time curve (AUC24-∞) and concentration at 24 h (C24h) were shown to correlate with renal toxicity only, and 49,800 μg/L × h (109.6 μmol/L × h) and C24h > 3930 μg/L (8.65 μmol/L) were potential exposure thresholds predicting high risk (proportion > 60%).

CONCLUSIONS

A population PK model was developed for HD-MTX in patients with CNS lymphoma. The toxicity analysis showed that lower baseline eGFR and male sex, and higher MTX dose are associated with increased risk of acute nephro- and hepatotoxicity, respectively. The proposed exposure thresholds that predict high risk of renal toxicity and the developed models hold the potential to guide HD-MTX dosage individualization and better prevent acute toxicity.

Impact of Methotrexate and 7-Hydroxymethotrexate Exposure on Renal Toxicity in Pediatric Non-Hodgkin Lymphoma

BACKGROUND

7-Hydroxymethotrexate (7-OHMTX) is the main metabolite in plasma following high-dose MTX (HD-MTX), which may result in activity and toxicity of the MTX. Moreover, 7-OHMTX could produce crystalline-like deposits within the renal tubules under acidic conditions or induce renal inflammation, oxidative stress, and cell apoptosis through various signaling pathways, ultimately leading to kidney damage. The objectives of this study were thus to explore the exposure-safety relationship of two compounds and search the most reliable marker for predicting HDMTX nephrotoxicity.

METHOD

A total of 280 plasma concentration data (140 for MTX and 140 for 7-OHMTX) for 60 pediatric patients with non-Hodgkin lymphoma (NHL) were prospectively collected. Plasma MTX and 7-OHMTX concentrations were determined using a high-performance liquid chromatography tandem mass spectrometry (HPLC-MS/MS) method. A nonlinear mixed effect model approach was used to build a joint population pharmacokinetic (PopPK) model. After validation, the model estimated the peak concentration (Cmax) and area under the curve within the initial 48 h (AUC0-48h) of the patients after drug administration by Bayesian feedback. The receiver operating characteristic (ROC) curves were generated to identify an exposure threshold associated with nephrotoxicity.

RESULTS

A three-compartment chain model (central and peripheral compartments for MTX and central compartment 7-OHMTX) with the first-order elimination adequately characterized the in vivo process of MTX and 7-OHMTX. The covariate analysis identified that the aspartate aminotransferase (AST) was strongly associated with the peripheral volume of distribution of MTX. Moreover, the Cmax of MTX and 7-OHMTX showed significant differences (p < 0.0001, p = 0.0472, respectively) among patients with or without nephrotoxicity. Similarly, individuals with nephrotoxicity also exhibited substantially higher ratio of 7-OHMTX to MTX peak concentration and the sum of MTX + 2.25 times the concentration of 7-OHMTX (p < 0.0001, p = 0.0426, respectively). By ROC analysis, the Cmax of MTX and 7-OHMTX had the greatest area under the curve (AUC) values (0.769 and 0.771, respectively). A Cmax threshold of 9.26 μmol/L for MTX or a Cmax threshold of 0.66 μmol/L for 7-OHMTX was associated with the best sensitivity/specificity for toxicity events (MTX: sensitivity = 0.886; specificity = 0.70; 7-OHMTX: sensitivity = 0.886; specificity = 0.70).

CONCLUSIONS

We demonstrated that the Cmax of MTX and 7-OHMTX were the most reliable markers associated with nephrotoxicity and proposed a Cmax threshold of 9.26 μmol/L for MTX and 0.66 μmol/L for 7-OHMTX as the point with a high risk of nephrotoxicity. Altogether, this study may contribute to crucial insights for ensuring the safe administration of drugs in pediatric clinical practice.

Molecular mechanisms underlying methotrexate-induced intestinal injury and protective strategies

Methotrexate (MTX) is a folic acid reductase inhibitor that manages various malignancies as well as immune-mediated inflammatory chronic diseases. Despite being frequently prescribed, MTX's severe multiple toxicities can occasionally limit its therapeutic potential. Intestinal toxicity is a severe adverse effect associated with the administration of MTX, and patients are significantly burdened by MTX-provoked intestinal mucositis. However, the mechanism of such intestinal toxicity is not entirely understood, mechanistic studies demonstrated oxidative stress and inflammatory reactions as key factors that lead to the development of MTX-induced intestinal injury. Besides, MTX causes intestinal cells to express pro-inflammatory cytokines like interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-α), which activate nuclear factor-kappa B (NF-κB). This is followed by the activation of the Janus kinase/signal transducer and activator of the transcription3 (JAK/STAT3) signaling pathway. Moreover, because of its dual anti-inflammatory and antioxidative properties, nuclear factor erythroid-2-related factor 2/heme oxygenase-1 (Nrf2/HO-1) has been considered a critical signaling pathway that counteracts oxidative stress in MTX-induced intestinal injury. Several agents have potential protective effects in counteracting MTX-provoked intestinal injury such as omega-3 polyunsaturated fatty acids, taurine, umbelliferone, vinpocetine, perindopril, rutin, hesperidin, lycopene, quercetin, apocynin, lactobacillus, berberine, zinc, and nifuroxazide. This review aims to summarize the potential redox molecular mechanisms of MTX-induced intestinal injury and how they can be alleviated. In conclusion, studying these molecular pathways might open the way for early alleviation of the intestinal damage and the development of various agent plans to attenuate MTX-mediated intestinal injury.

An interactive dose optimizer based on population pharmacokinetic study to guide dosing of methotrexate in Chinese patients with osteosarcoma

PURPOSE

Osteosarcoma is a rare tumor with an incidence of 4.4 cases per million per year in adolescent. High-dose methotrexate (HD-MTX) is the standard first-line chemotherapeutic agent for osteosarcoma. However, its efficacy can vary significantly among individuals due to wide pharmacokinetic variability. Despite this, only a few population pharmacokinetics (popPK) models based on Chinese patients with osteosarcoma have been reported. Thus, this study aimed to develop a HD-MTX popPK model and an individual model-based dose optimizer for osteosarcoma therapy.

METHOD

A total of 680 MTX serum concentrations from 57 patients with osteosarcoma were measured at the end of MTX infusion and 10 h, 24 h, 48 h, and 72 h after the start of infusion. Using the first-order conditional estimation method with NONMEM, a popPK model was estimated. Goodness-of-fit plots, visual predictive checks, and bootstrap analysis were generated to evaluate the final model. A dose optimizer tool was developed based on the validated models using R Shiny. Additionally, clinical data from 12 patients with newly diagnosed osteosarcoma were collected and used as the validation set to preliminarily verify the predictive ability of the popPK model and the dose optimizer tool.

RESULTS

Body surface area (BSA) was the most significant covariate for compartment distribution. Creatinine clearance (CrCL) and co-administration of NSAIDs were introduced as predictors for central compartmental and peripheral compartmental clearance, respectively. Co-administration of NSAIDs was associated with significantly higher MTX concentrations at 72 h (p = 0.019). The dose optimizer tool exhibited a high consistency in predicting MTX AUC compared to the actual AUC (r = 0.821, p < 0.001) in the validation set.

CONCLUSION

The dose optimizer tool could be used to estimate individual PK parameters, and optimize personalized MTX therapy in particular patients.

Understanding hemoglobin contribution to high-dose methotrexate disposition-population pharmacokinetics in pediatric patients with hematological malignancies

PURPOSE

The aim of the present study was to develop a population pharmacokinetic model for methotrexate (MTX) during high-dose treatment (HDMTX) in pediatric patients with acute lymphoblastic leukemia (ALL) and non-Hodgkin's lymphoma (NHL) and to describe the influence of variability factors.

METHODS

The study included 50 patients of both sexes (aged 1-18 years) who received 3 or 5 g/m2 of HDMTX. A nonlinear mixed effect modeling approach was applied for data analysis. Parameter estimation was performed by first-order conditional estimation method with interaction (FOCEI), whereas stepwise covariate modeling was used to assess variability factors.

RESULTS

The final model is a two-compartment model that incorporates the effect of body surface area and the influence of hemoglobin and serum creatinine on MTX clearance (CL). Population pharmacokinetic values for a typical subject were estimated at 5.75 L/h/m2 for clearance (CL), 21.3 L/m2 for volume of the central compartment (V1), 8.2 L/m2 for volume of the peripheral compartment (V2), and 0.087 L/h/m2 for intercompartmental clearance (Q). According to the final model, MTX CL decreases with increasing serum creatinine, whereas a positive effect was captured for hemoglobin. A difference of almost 32% in MTX CL was observed among patients' hemoglobin values reported in the study.

CONCLUSION

The developed population pharmacokinetic model can contribute to the therapy optimization during HDMTX in pediatric patients with ALL and NHL. In addition to renal function and body weight, it describes the influence of hemoglobin on CL, allowing better understanding of its contribution to the disposition of HDMTX.

External Evaluation of Population Pharmacokinetic Models for High-Dose Methotrexate in Adult Patients with Hematological Tumors

Currently, numerous population pharmacokinetic (popPK) models for methotrexate (MTX) have been published for estimating PK parameters and variability. However, it is unclear whether the accuracy of these models is sufficient for clinical application. The aim of this study is to evaluate published models and assess their predictive performance according to the standards of scientific research. A total of 237 samples from 74 adult patients who underwent high-dose MTX (HDMTX) treatment at Shanghai Changzheng Hospital were collected. The software package NONMEM was used to perform an external evaluation for each model, including prediction-based diagnosis, simulation-based diagnosis, and Bayesian forecasting. The simulation-based diagnosis includes normalized prediction distribution error (NPDE) and visual predictive check (VPC). Following screening, 7 candidate models suitable for external validation were identified for comparison. However, none of these models exhibited excellent predictive performance. Bayesian simulation results indicated that the prediction precision and accuracy of all models significantly improved when incorporating prior concentration information. The published popPK models for MTX exhibit significant differences in their predictive performance, and none of the models were able to accurately predict MTX concentrations in our data set. Therefore, before adopting any model in clinical practice, extensive evaluation should be conducted.

Risk factors associated with high-dose methotrexate induced toxicities

INTRODUCTION

High-dose methotrexate (HDMTX) therapy poses challenges in various neoplasms due to individualized pharmacokinetics and associated adverse effects. Our purpose is to identify early risk factors associated with HDMTX-induced toxicities, paving the way for personalized treatment.

AREAS COVERED

A systematic review of PubMed and Cochrane databases was conducted for articles from inception to July 2023. Eligible studies included reviews, clinical trials, and real-world analyses. Irrelevant studies were excluded, and manual searches and citation reviews were performed. Factors such as MTX exposure, drug interactions, demographics, serum albumin, urine pH, serum calcium, and genetic polymorphisms affecting MTX transport (e.g. SLCO1B1), intracellular folate metabolism (MTHFR), cell development (ARID5B), metabolic pathways (UGT1A1, PNPLA3), as well as epigenetics were identified.

EXPERT OPINION

This comprehensive review aids researchers and clinicians in early identification of HDMTX toxicity risk factors. By understanding the multifaceted risk factors associated with hematologic malignancies, personalized treatment approaches can be tailored to optimize therapeutic outcomes.

An ultrahigh-performance liquid chromatography-tandem mass spectrometry method for quantification of methotrexate and 7-hydroxy-methotrexate and application for therapeutic drug monitoring in patients with central nervous system lymphoma

A method using ultrahigh-performance liquid chromatography-tandem mass spectrometry was developed, validated, and applied to simultaneously determine plasma methotrexate (MTX) and 7-hydroxy-methotrexate (7-OH-MTX) in 117 patients with central nervous system (CNS) lymphoma. The ion transitions utilized were m/z 455.2 > 308.2 for MTX and m/z 471.2 > 324.1 for 7-OH-MTX. Samples were prepared through protein precipitation using methanol. Chromatographic separation was achieved within 3.0 min on a CMS9030 column (Ruixi, 2.1 × 50 mm, 3 μm) through a gradient elution of methanol and a 10% ammonium acetate solution at a flow rate of 0.4 mL/min. The method demonstrated linearity in the concentration range of 0.05-10 μM for MTX and 0.25-50 μM for 7-OH-MTX. The intra- and inter-day inaccuracy ranged from -7.38% to 7.83%, and the imprecision was less than 6.00% for both analytes. The recovery and matrix effect normalized by the internal standard (MTX-D3 ) remained consistent. Both analytes remained stable under nine different storage conditions. In patients with CNS lymphoma, MTX levels at 12 h and 7-OH-MTX levels at 12, 36, and 60 h after dosing in individuals with impaired renal function were significantly higher compared with those with normal renal function. 7-OH-MTX could potentially serve as a superior indicator for nephrotoxicity compared with MTX.

Factors influencing methotrexate pharmacokinetics highlight the need for individualized dose adjustment: a systematic review

PURPOSE

To develop a population pharmacokinetic (PPK) model for methotrexate (MTX) dosage for all ages, assess the association between concentration and clearance, and determine covariates affecting MTX disposition.

METHODS

We compared MTX PK profiles among neonates, children, and adults by performing a systematic literature search for published population MTX models and conducted a Monte Carlo-based meta-analysis. Subsequently, we evaluated study quality and covariates significantly affecting dosage regimens and compared LDMTX and HDMTX PK profiles.

RESULTS

Of the total 40 studies included, 34 were HDMTX, and six were LDMTX studies. For HDMTX, three studies involving neonates reported estimated apparent clearances (median, range) of 0.53 (0.27-0.77) L/kg/h; for 14 studies involving children, 0.23 (0.07-0.23) L/kg/h; and for 13 involving adults, 0.11 (0.03-0.22) L/kg/h. Neonates had a higher volume of distribution than children and adults. For LDMTX studies, apparent clearance was 0.085 (0.05-1.68) L/kg/h, and volume of distribution was 0.25 (0.018-0.47) L/kg, lower than those of HDMTX studies, with large between-subject variability. Bodyweight significantly influenced apparent clearance and volume of distribution, whereas renal function mainly influenced clearance. Mutations in certain genes reduced MTX clearance by 8-35.3%, whereas those in others increased it by 15-48%. Body surface area (BSA) significantly influenced apparent clearance with a median reduction of 51% when BSA increased in pediatric patients.

CONCLUSIONS

Methotrexate dosage regimens were primarily based on body surface area and renal function. Further studies are needed to evaluate MTX pharmacokinetics and pharmacodynamics in both children (especially infants) and adults.

A semimechanistic pharmacokinetic/pharmacodynamic model for alanine aminotransferase-based hepatotoxicity of methotrexate in paediatric patients with acute lymphoid leukaemia

AIMS

Methotrexate (MTX) is recognized for its potential to induce hepatotoxicity, commonly manifested by elevated alanine aminotransferase (ALT) levels. However, the quantitative relationship between the pharmacokinetics (PK) of MTX and ALT-based hepatotoxicity remains unclear. This study aimed to develop a semimechanistic PK/pharmacodynamic (PD) model to characterize the MTX-induced hepatotoxicity based on ALT in paediatric patients with acute lymphoid leukaemia.

METHODS

A retrospective study was conducted on paediatric patients who received high-dose (3-5 g/m2 ) MTX treatment. MTX concentrations were assessed at 24-h intervals until the concentration dropped below 0.1 μmol/L. ALT concentrations were measured both before and after MTX administration. A population PK model was initially developed, which was later connected to a semimechanistic hepatotoxicity model.

RESULTS

The PK model was developed using 354 MTX concentrations obtained from 51 patients, while the PD model was constructed using 379 ALT concentrations collected from 48 patients. The optimal PK model for MTX consisted of a 2-compartment structure, where body surface area served as a covariate for clearance and central volume of distribution. An indirect response model coupled to a liver injury signal transduction model was developed to describe the dynamics of ALT after MTX administration. The drug effect was adequately described by a linear model, exhibiting considerable interoccasion variability for each treatment session. No significant covariates were identified to have an impact on the PD parameters.

CONCLUSION

A semimechanistic model was developed to describe ALT-based hepatotoxicity of MTX, and it has the potential to serve as a valuable tool for characterizing drug-induced hepatotoxicity.

Early clinical indicators of acute kidney injury caused by administering high-dose methotrexate therapy to juvenile pigs

Introduction

Early identification of compromised renal clearance caused by high-dose methotrexate (HDMTX) is essential for initiating timely interventions that can reduce acute kidney injury and MTX-induced systemic toxicity.

Methods

We induced acute kidney injury (AKI) by infusing 42 juvenile pigs with 4 g/kg (80 g/m2) of MTX over 4 hours without high-volume alkalinizing hydration therapy. Concentrations of serum creatinine and MTX were measured at 15 time points up to 148 hours, with 10 samples collected during the first 24 hours after the start of the HDMTX infusion.

Results

During the first 28 hours, 81% of the pigs had increases in the concentrations of serum creatinine in one or more samples indicative of AKI (i.e., > 0.3g/dL increase). A rate of plasma MTX clearance of less than 90% during the initial 4 hours after the HDMTX infusion and a total serum creatinine increase at 6 and 8 hours after starting the infusion greater than 0.3 g/dL were predictive of AKI at 28 hours (p < 0.05 and p < 0.001, respectively). At conclusion of the infusion, pigs with a creatinine concentration more than 0.3 g/dL higher than baseline or serum MTX greater than 5,000 μmol/L had an increased risk of severe AKI.

Conclusions

Our findings suggest that serum samples collected at conclusion and shortly after HDMTX infusion can be used to predict impending AKI. The pig model can be used to identify biological, environmental, and iatrogenic risk factors for HDMTX-induced AKI and to evaluate interventions to preserve renal functions, minimize acute kidney injury, and reduce systemic toxicity.

External evaluation and systematic review of population pharmacokinetic models for high-dose methotrexate in cancer patients

Several population pharmacokinetic (PPK) models have been established to optimize the therapeutic regimen and reduce the toxicity of high-dose methotrexate (HDMTX) in patients with cancer. However, their predictive performance when extrapolated to different clinical centers was unknown. In this study, we aimed to externally evaluate the predictive ability of HDMTX PPK models and determine the potential influencing factors. We searched the literature and determined the predictive performance of the selected models using methotrexate concentrations in 721 samples from 60 patients in the First Affiliated Hospital of the Navy Medical University. Prediction-based diagnostics and simulation-based normalized prediction distribution errors (NPDE) were used to evaluate the predictive performance of the models. The influence of prior information was also assessed using Bayesian forecasting, and the potential factors affecting model predictability were investigated. Thirty models extracted from published PPK studies were assessed. Prediction-based diagnostics showed that the number of compartments potentially influenced model transferability, and simulation-based NPDE indicated model misspecification. Bayesian forecasting significantly improved the predictive performance of the models. Various factors, including bioassays, covariates, and population diagnosis, influence model extrapolation. The published models were unsatisfactory for all prediction-based diagnostics, except for the 24 h methotrexate concentration monitoring and simulation-based diagnostics, making them inappropriate for direct extrapolation. Moreover, Bayesian forecasting combined therapeutic drug monitoring could improve the predictive performance of the models.

Insights from a pharmacometric analysis of HDMTX in adults with cancer: Clinically relevant covariates for application in precision dosing

AIMS

High-dose methotrexate (HDMTX) is an essential part of the treatment of several adult and paediatric malignancies. Despite meticulous supportive care during HDMTX administration, severe toxicities, including acute kidney injury (AKI), may occur contributing to patient morbidity. Population pharmacokinetics provide a powerful tool to predict time to clear HDMTX and adjust subsequent doses. We sought to develop and validate pharmacokinetic models for HDMTX in adults with diverse malignancies and to relate systemic exposure with the occurrence of severe toxicity.

METHODS

Anonymized, de-identified data were provided from 101 US oncology practices that participate in the Guardian Research Network, a non-profit clinical research consortium. Modelled variables included clinical, laboratory, demographic and pharmacological data. Population pharmacokinetic analysis was performed by means of nonlinear mixed effects modelling using MonolixSuite.

RESULTS

A total of 693 HDMTX courses from 243 adults were analysed, of which 62 courses (8.8%) were associated with stage 2/3 acute kidney injury (43 stage 2, 19 stage 3). A three-compartment model adequately fitted the data. Time-dependent serum creatinine, baseline serum albumin and allometrically scaled bodyweight were clinically significant covariates related to methotrexate clearance. External evaluation confirmed a satisfactory predictive performance of the model in adults receiving HDMTX. Dose-normalized methotrexate concentration at 24 and 48 hours correlated with AKI incidence.

CONCLUSION

We developed a population pharmacometric model that considers weight, albumin and time-dependent creatinine that can be used to guide supportive care in adult patients with delayed HDMTX elimination.

Population pharmacokinetics of methotrexate in pediatric patients with acute lymphoblastic leukemia and malignant lymphoma

1. This study aimed to identify physiological and pharmacogenomic covariates and develop a population pharmacokinetic model of high-dose methotrexate (HD-MTX) in Chinese pediatric patients with acute lymphoblastic leukemia (ALL) and malignant lymphoma.2. A total 731 MTX courses and 1658 MTX plasm concentration from 205 pediatric patients with ALL and malignant lymphoma were analyzing using a nonlinear mixed-effects model technique. 47 SNPs in 16 MTX-related gene were genotyped and screened as covariates. A PPK model was established to determine the influence of covariates such as body surface area (BSA), age, laboratory test value, and SNPs on the pharmacokinetic process of HD-MTX.3. Two-compartmental model with allometric scaling using BSA could nicely characterize the in vivo behavior of HD-MTX. After accounting for body size, rs17004785 and rs4148416 were the covariates that influence MTX clearance (CL). The PPK model obtained was: CL =9.33 * (BSA/1.73) ^0.75 * e^{0.13*rs17004785} * e^{0.39*rs4148416} * e^ηCL, Vc =24.98 * (BSA/1.73) * e^ηvc, Q = 0.18 * (BSA/1.73) ^0.75 * e^ηQ and Vp =4.70 * (BSA/1.73) * e^ηvp.4. The established model combined with Bayesian approach could estimate individual pharmacokinetic parameters and optimize personalized HD-MTX therapy for pediatric patients with ALL and malignant lymphoma.

Evaluation and Application of Population Pharmacokinetic Models for Identifying Delayed Methotrexate Elimination in Patients With Primary Central Nervous System Lymphoma

Objective: Several population pharmacokinetic (popPK) models have been developed to determine the sources of methotrexate (MTX) PK variability. It remains unknown if these published models are precise enough for use or if a new model needs to be built. The aims of this study were to 1) assess the predictability of published models and 2) analyze the potential risk factors for delayed MTX elimination.

Methods: A total of 1458 MTX plasma concentrations, including 377 courses (1–17 per patient), were collected from 77 patients who were receiving high-dose MTX for the treatment of primary central nervous system lymphoma in Huashan Hospital. PopPK analysis was performed using the NONMEM® software package. Previously published popPK models were selected and rebuilt. A new popPK model was then constructed to screen potential covariates using a stepwise approach. The covariates were included based on the combination of theoretical mechanisms and data properties. Goodness-of-fit plots, bootstrap, and prediction- and simulation-based diagnostics were used to determine the stability and predictive performance of both the published and newly built models. Monte Carlo simulations were conducted to qualify the influence of risk factors on the incidence of delayed elimination.

Results: Among the eight evaluated published models, none presented acceptable values of bias or inaccuracy. A two-compartment model was employed in the newly built model to describe the PK of MTX. The estimated mean clearance (CL/F) was 4.91 L h⁻¹ (relative standard error: 3.7%). Creatinine clearance, albumin, and age were identified as covariates of MTX CL/F. The median and median absolute prediction errors of the final model were -10.2 and 36.4%, respectively. Results of goodness-of-fit plots, bootstrap, and prediction-corrected visual predictive checks indicated the high predictability of the final model.

Conclusions: Current published models are not sufficiently reliable for cross-center use. The elderly patients and those with renal dysfunction, hypoalbuminemia are at higher risk of delayed elimination.

Extracorporeal Treatment for Methotrexate Poisoning: Systematic Review and Recommendations from the EXTRIP Workgroup

Methotrexate is used in the treatment of many malignancies, rheumatological diseases, and inflammatory bowel disease. Toxicity from use is associated with severe morbidity and mortality. Rescue treatments include intravenous hydration, folinic acid, and, in some centers, glucarpidase. We conducted systematic reviews of the literature following published EXtracorporeal TReatments In Poisoning (EXTRIP) methods to determine the utility of extracorporeal treatments in the management of methotrexate toxicity. The quality of the evidence and the strength of recommendations (either "strong" or "weak/conditional") were graded according to the GRADE approach. A formal voting process using a modified Delphi method assessed the level of agreement between panelists on the final recommendations. A total of 92 articles met inclusion criteria. Toxicokinetic data were available on 90 patients (89 with impaired kidney function). Methotrexate was considered to be moderately dialyzable by intermittent hemodialysis. Data were available for clinical analysis on 109 patients (high-dose methotrexate [>0.5 g/m2]: 91 patients; low-dose [≤0.5 g/m2]: 18). Overall mortality in these publications was 19.5% and 26.7% in those with high-dose and low-dose methotrexate-related toxicity, respectively. Although one observational study reported lower mortality in patients treated with glucarpidase compared with those treated with hemodialysis, there were important limitations in the study. For patients with severe methotrexate toxicity receiving standard care, the EXTRIP workgroup: (1) suggested against extracorporeal treatments when glucarpidase is not administered; (2) recommended against extracorporeal treatments when glucarpidase is administered; and (3) recommended against extracorporeal treatments instead of administering glucarpidase. The quality of evidence for these recommendations was very low. Rationales for these recommendations included: (1) extracorporeal treatments mainly remove drugs in the intravascular compartment, whereas methotrexate rapidly distributes into cells; (2) extracorporeal treatments remove folinic acid; (3) in rare cases where fast removal of methotrexate is required, glucarpidase will outperform any extracorporeal treatment; and (4) extracorporeal treatments do not appear to reduce the incidence and magnitude of methotrexate toxicity.

A Systematic Review of Population Pharmacokinetic Models of Methotrexate

BACKGROUND AND OBJECTIVES

Methotrexate (MTX) is widely used for the treatment of a variety of neoplastic and autoimmune diseases. However, its toxicity and efficacy varied greatly among individuals, and they could be predicted by its pharmacokinetics. Many population pharmacokinetic models have been published to describe MTX pharmacokinetics. The objective of this systematic review was to summarize and discuss covariates with significant influence on MTX pharmacokinetics.

METHODS

We searched PubMed and EMBASE databases from their inception to April 2021 for population pharmacokinetic of MTX. The articles were screened by inclusion and exclusion criteria. The characteristics of studies and information for model construction and validation were extracted, summarized and discussed.

RESULTS

Thirty-five articles were included. The two-compartment model well described the pharmacokinetic behavior of MTX. For inter-individual variability, an exponential distribution error model was usually used for high-dose MTX population pharmacokinetic models, while a proportional distribution error model was used for low-dose MTX population pharmacokinetic models. Proportional and combined proportional and additive error models were used to describe residual error. Renal function was an independent indicator of MTX clearance. Body weight, age, gene polymorphisms (SLCO1B1, ABCC2, ABCB1, ABCG2 and MTHFR) and co-medications (proton pump inhibitors, non-steroidal anti-inflammatory drug, dexamethasone, vancomycin, penicillin and salicylic acid) could influence MTX clearance. Body weight, body surface area, age and dosage regimen have significant influence on MTX central compartment volume. Internal bootstrap test, external validation and visual predictive check were used to evaluate model predictive ability.

CONCLUSIONS

Various covariates could affect MTX pharmacokinetics, and their relationships have been summarized and discussed. This review will be helpful for researchers to develop their own population pharmacokinetic models and select appropriate models for individualized therapy of MTX.

Urine volume to hydration volume ratio is associated with pharmacokinetics of high-dose methotrexate in patients with primary central nervous system lymphoma

High-dose methotrexate (HD-MTX)-based chemotherapy is the first-line treatment for primary central nervous system lymphoma (PCNSL), but is associated with severe adverse effects, including myelosuppression and renal impairment. MTX is primarily excreted by the kidneys. Renal function calculated using serum creatinine (Scr) derived from muscle may be overestimated in elderly PCNSL patients. Therefore, we aimed to construct a population pharmacokinetic model in PCNSL patients and explore the factors associated with MTX clearance. Sixteen PCNSL patients (median age, 66 years) treated with HD-MTX were included, and serum MTX concentrations were measured at 193 points in 49 courses. A population pharmacokinetic analysis was performed using NONMEM. A Monte Carlo simulation was conducted, in which serum MTX concentrations were stratified into three groups of creatine clearance (Ccr) (50, 75, and 100 ml/min) with three groups of the urine volume to hydration volume (UV/HV) ratio (<1, 1-2, and >2). The final model was constructed as follows: MTX clearance = 4.90·(Ccr/94.5)0.456 ·(UV/HV)0.458 . In the Monte Carlo simulation, serum MTX concentrations were below the standard values (10, 1, and 0.1 µM at 24, 48, and 72 h, respectively, after the start of the MTX administration) in most patients with UV/HV >2, even with Ccr of 50 ml/min. Conversely, half of the patients with UV/HV <1 and Ccr of 50 ml/min failed to achieve the standard values. The present results demonstrated that the UV/HV ratio was useful for describing the pharmacokinetics of MTX in PCNSL patients.

Evaluation of body-surface-area adjusted dosing of high-dose methotrexate by population pharmacokinetics in a large cohort of cancer patients

Background

The aim of this study was to identify sources of variability including patient gender and body surface area (BSA) in pharmacokinetic (PK) exposure for high-dose methotrexate (MTX) continuous infusion in a large cohort of patients with hematological and solid malignancies.

Methods

We conducted a retrospective PK analysis of MTX plasma concentration data from hematological/oncological patients treated at the University Hospital of Cologne between 2005 and 2018. Nonlinear mixed effects modeling was performed. Covariate data on patient demographics and clinical chemistry parameters was incorporated to assess relationships with PK parameters. Simulations were conducted to compare exposure and probability of target attainment (PTA) under BSA adjusted, flat and stratified dosing regimens.

Results

Plasma concentration over time data (2182 measurements) from therapeutic drug monitoring from 229 patients was available. PK of MTX were best described by a three-compartment model. Values for clearance (CL) of 4.33 [2.95–5.92] L h^− 1 and central volume of distribution of 4.29 [1.81–7.33] L were estimated. An inter-occasion variability of 23.1% (coefficient of variation) and an inter-individual variability of 29.7% were associated to CL, which was 16 [7–25] % lower in women. Serum creatinine, patient age, sex and BSA were significantly related to CL of MTX. Simulations suggested that differences in PTA between flat and BSA-based dosing were marginal, with stratified dosing performing best overall.

Conclusion

A dosing scheme with doses stratified across BSA quartiles is suggested to optimize target exposure attainment. Influence of patient sex on CL of MTX is present but small in magnitude.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08443-x.

A minimal physiologically based pharmacokinetic model for high-dose methotrexate

Purpose

High-dose methotrexate (HDMTX) is administered for the treatment of a variety of malignant tumors. Wide intra- and inter-individual variabilities characterize the pharmacokinetics of MTX, which is mostly excreted renally. HDMTX dosages are prescribed as a function of body surface area whereas dose adjustments depending on renal function are not well defined. We develop a population pharmacokinetic model with a physiological description of renal excretion as the basis for clinical tools able to suggest model-informed dosages and support therapeutic monitoring.

Methods

This article presents a minimal physiologically based pharmacokinetic (PBPK) model for HDMTX, which specifically accounts for individual characteristics such as body weight, height, gender, age, hematocrit, and serum creatinine to provide individualized predictions. The model supplies a detailed and mechanistic description of capillary and cellular exchanges between plasma, interstitial fluid, and intracellular fluid compartments, and focuses on an individualized description of renal excretion.

Results

The minimal PBPK model is identified and validated with a literature dataset based on Chinese patients suffering from primary central nervous system lymphoma. A comparison with a pharmacokinetic model from the literature suggests that the proposed model provides improved predictions. Remarkably, the model does not present any significant bias in a wide range of degrees of renal function.

Conclusion

Results show that model predictions can capture the wide intra- and inter-individual variability of HDMTX, and highlight the role played by the individual degree of renal function. The proposed model can be the basis for the development of clinical decision-support systems for individualized dosages and therapeutic monitoring.

Comparison of four immunoassays to an HPLC method for the therapeutic drug monitoring of methotrexate: Influence of the hydroxylated metabolite levels and impact on clinical threshold

OBJECTIVES

Methotrexate requires therapeutic drug monitoring in oncology because of narrow therapeutic index, especially the metabolite 7-hydroxymethotrexate exhibits nephrotoxicity. The goal of this study was to evaluate different assays and their impact on clinical decisions.

METHODS

Following routine measurement with Abbott TDxFLx® assay (MTX-TDX), 62 samples were analysed on Architect®i1000 (MTX-ARCHI), Xpand® (ARK/XPND), Indiko® (ARK/INDI), and HPLC (MTX-HPLC) as the reference method. The influence of 7-hydroxymethotrexate was explored on ARK reagent to document the cause of the observed bias. ROC curves were built to study the impact of the method on the discharge thresholds for the patients at three levels.

RESULTS

Total imprecision was below 2.60% for the methotrexate-ARCHI and close to 10% for both ARK assays for plasma pools. The correlation coefficients were 0.93, 0.93, 0.89 and 0.95, the Bland-Altman difference plot revealed a bias of 0.075, 0.037, 0.049 and -0.002, and the number of results exceeding the TE criteria of 0.1 µM was 17 (27%), 13 (21%), 15 (24%) and 15 (24%) for MTX-TDX, ARK/INDI, ARK/XPND and MTX-ARCHI, respectively. Cross reactivity with 7-hydroxymethotrexate was between 1 and 9%. Overestimation of methotrexate concentration was between -4% and +32%. The most robust clinical level was found to be the highest level (0.2 µM) with ROC curves.

CONCLUSIONS

The authors found the best results for imprecision with chemiluminescent microparticle immunoassay method on methotrexate-ARCHI, with bias below to the RICOS recommendations and best correlation to the reference method. Impact on the threshold values for clinical decision need to be clearly exposed to clinicians.

Simultaneous Determination of Urine Methotrexate, 7-Hydroxy Methotrexate, Deoxyaminopteroic Acid, and 7-Hydroxy Deoxyaminopteroic Acid by UHPLC-MS/MS in Patients Receiving High-dose Methotrexate Therapy

Nephrotoxicity, the most important toxicity in high-dose methotrexate (MTX) therapy, is partly caused by the formation of crystal deposits in the kidney due to poor water solubility of MTX and its metabolites 7-hydroxy methotrexate (7-OH MTX), deoxyaminopteroic acid (DAMPA) and 7-hydroxy deoxyaminopteroic acid (7-OH DAMPA). Plasma MTX level-guided urine alkalinization, leucovorin rescue and glucarpidase detoxification are common strategies to overcome MTX-related nephrotoxicity. However, overestimation is a problem for MTX analysis by immunoassays due to the cross-reactivity of MTX metabolites (7-OH MTX and DAMPA). An UHPLC-MS/MS method for the simultaneous determination of MTX, 7-OH MTX, DAMPA and 7-OH DAMPA in human urine was developed, validated and applied in clinical practice. Samples were treated by one-step protein precipitation and analyzed within 3 min. The calibration range was 0.02 to 4 μmol/L for MTX and DAMPA, and 0.1 to 20 μmol/L for 7-OH MTX and 7-OH DAMPA. For all analytes, the intra-day and inter-day bias and imprecision were -8.0 to 7.6 and <9.0%, the internal standard normalized recovery and matrix factor were 92.34 to 109.49 and <20.68%. The plasma MTX and 7-OH MTX levels increased with the urine drug levels, age, serum creatinine and alanine transaminase, but urine could not replace blood for MTX monitoring due to their poor correlation (R², 0.16 to 0.51). Dose-normalized urine and plasma MTX and 7-OH MTX levels were similar between different patient groups (urine pH <7 or ≥7). Due to the large inter-individual variance of the analytes levels in both plasma and urine, these findings should be treated with caution.

Pharmacokinetics and pharmacogenetics of high-dose methotrexate in Chinese adult patients with non-Hodgkin lymphoma: a population analysis

PURPOSE

High-dose methotrexate (HD-MTX) is widely used in the treatment of non-Hodgkin lymphoma (NHL), but the pharmacokinetic properties of HD-MTX in Chinese adult patients with NHL have not yet been established through an approach that integrates genetic covariates. The purposes of this study were to identify both physiological and pharmacogenomic covariates that can explain the inter- and intraindividual pharmacokinetic variability of MTX in Chinese adult patients with NHL and to explore a new sampling strategy for predicting delayed MTX elimination.

METHODS

A total of 852 MTX concentrations from 91 adult patients with NHL were analyzed using the nonlinear mixed-effects modeling method. FPGS, GGH, SLCO1B1, ABCB1 and MTHFR were genotyped using the Sequenom MassARRAY technology platform and were screened as covariates. The ability of different sampling strategies to predict the MTX concentration at 72 h was assessed through maximum a posteriori Bayesian forecasting using a validation dataset (18 patients).

RESULTS

A two-compartment model adequately described the data, and the estimated mean MTX clearance (CL) was 6.03 L/h (9%). Creatinine clearance (CrCL) was identified as a covariate for CL, whereas the intercompartmental clearance (Q) was significantly affected by the body surface area (BSA). However, none of the genotypes exerted a significant effect on the pharmacokinetic properties of MTX. The percentage of patients with concentrations below 0.2 µmol/L at 72 h decreased from 65.6 to 42.6% when the CrCL decreased from 90 to 60 ml/min/1.73 m² with a scheduled dosing of 3 g/m², and the same trend was observed with dose regimens of 1 g/m² and 2 g/m². Bayesian forecasting using the MTX concentrations at 24 and 42 h provided the best predictive performance for estimating the MTX concentration at 72 h after dosing.

CONCLUSIONS

The MTX population pharmacokinetic model developed in this study might provide useful information for establishing personalized therapy involving MTX for the treatment of adult patients with NHL.

Population pharmacokinetics of high-dose methotrexate in Chinese pediatric patients with medulloblastoma

Methotrexate (MTX) pharmacokinetics has substantial inter-individual variability and toxicity. In children with medulloblastoma treated with high-dose methotrexate (HD-MTX), the pharmacokinetic properties of methotrexate have not been established. A total of 660 serum samples from 105 pediatric patients with medulloblastoma were included in a population pharmacokinetic (PPK) analysis of methotrexate by using the nonlinear mixed-effects modeling method. The basic one-compartment population pharmacokinetic model was established by NONMEM software and the first-order conditional estimation (FOCE) method, and the final covariate model was obtained by the stepwise regression method. Weight (WT), creatinine clearance (CrCL), and whether the treatment was combined with dexamethasone (DEX) were covariates that had significant effects on the clearance rate (CL) of the model. The pharmacokinetic equation of CL in the final covariate model was as follows: CL_i = 9.23× (1 + 0.0005× (θ_CrCL -105.78)) × (1 + 0.0017× (θ_WT -16)) × e^ηcl,i (L/h), IF (θ_DEX ) CL_i = 1.19× CL_i (L/h). The estimation accuracy of all pharmacokinetic parameters were acceptable (relative standard error < 14.74%). The goodness-of-fit diagram and bootstrap tests indicated that the final PPK model was stable with acceptable predictive ability. The PPK model may be useful for determining personalized medication levels in pediatric medulloblastoma patients undergoing HD-MTX therapy.

Population Pharmacokinetics of Glasdegib in Patients With Advanced Hematologic Malignancies and Solid Tumors

Glasdegib is an inhibitor of the Hedgehog pathway recently approved in the United States for the treatment of acute myeloid leukemia. A population pharmacokinetic analysis was conducted to characterize the kinetic behavior of glasdegib and its sources of variability (covariates) by utilizing data from 269 patients with cancer treated with oral glasdegib doses ranging from 5 to 640 mg/d. Nonlinear mixed-effects modeling was conducted using NONMEM (v.7.3) and Perl-speaks NONMEM (v.4.2.0). The estimated apparent total clearance, apparent central volume of distribution, and apparent peripheral volume of distribution were 6.27 L/h, 3.32 L, and 279.2 L, respectively. Age, sex, race, and hepatic function were not significant covariates on glasdegib pharmacokinetic parameters. Baseline body weight, percentage bone marrow blasts, creatinine clearance, and use of moderate or strong cytochrome P450 3A inhibitors were statistically significant covariates on apparent total clearance; however, the magnitude of the effects was not considered clinically meaningful.

Population pharmacokinetic analysis of high-dose methotrexate in pediatric and adult oncology patients

PURPOSE

High-dose methotrexate (HD-MTX) is widely used in pediatric and adult oncology treatment regimens. This study aimed to develop a population pharmacokinetic model to characterize pediatric and adult MTX exposure across various disease types and dosing regimens, and to evaluate exposure-toxicity relationships.

METHODS

MTX pharmacokinetic data from pediatric and adult patients were collected. A population pharmacokinetic model was developed to determine the effects of age, liver function, renal function, and demographics on MTX disposition. The final model was used in Monte Carlo simulations to generate expected exposures for different dosing regimens. The association of toxicity, determined through chart review, and MTX area under the curve (AUC) was modeled using logistic regression.

RESULTS

The analysis included 5116 MTX concentrations from 320 patients (135 adult, age 19-79 years; 185 pediatric, age 0.6-19 years). Estimated glomerular filtration rate (eGFR) and treatment cycle number were independent predictors of clearance (CL). CL varied 2.1-fold over the range of study eGFR values and increased 14% for treatment cycle numbers greater than 7. Higher MTX AUC was associated with higher risk of nephrotoxicity in adults, and neurotoxicity and hepatotoxicity in pediatrics.

CONCLUSIONS

This study represents one of the most comprehensive evaluations of HD-MTX PK across a wide range of ages and disease types. After accounting for differences in renal function, age did not impact CL, although toxicity patterns differed by age. The model allows for early identification of patients with slowed MTX clearance and at higher risk of toxicity.