ABCB1 C3435T genetic polymorphism on population pharmacokinetics of methotrexate after hematopoietic stem cell transplantation in Korean patients: a prospective analysis

Clinical effects of tacrolimus blood concentrations early after allogeneic hematopoietic stem cell transplantation

BACKGROUND AIMS

Tacrolimus (TAC) plus short-term methotrexate (stMTX) is used for graft-versus-host disease (GVHD) prophylaxis after allogeneic hematopoietic stem cell transplantation (allo-HSCT). TAC blood concentrations are frequently adjusted to enhance the graft-versus-leukemia/lymphoma effect or attenuate severe GVHD. Limited information is available on the clinical impact of these adjustments and the optimal time to perform them in order to achieve good clinical outcomes.

METHODS

We retrospectively analyzed 211 patients who underwent allo-HSCT at our institutes.

RESULTS

Higher TAC concentrations in week 3 correlated with a significantly higher cumulative incidence of relapse (CIR) (P = 0.03) and lower nonrelapse mortality (P = 0.04). The clinical impact of high TAC concentrations in week 3 on CIR was detected in the refined disease risk index: low/intermediate (P = 0.04) and high (P < 0.01), and conditioning regimens other than cyclophosphamide/total body irradiation and busulfan/cyclophosphamide (P = 0.07). Higher TAC concentrations in week 1 correlated with a lower grade 2-4 acute GVHD rate (P = 0.01). Higher TAC concentrations in weeks 2 and 3 correlated with slightly lower (P = 0.05) and significantly lower (P = 0.02) grade 3-4 acute GVHD rates, respectively. Higher TAC concentrations in weeks 1 and 3 were beneficial for severe acute GVHD in patients with a human leukocyte antigen-matched donor (P = 0.03 and P < 0.01, respectively), not treated with anti-thymocyte globulin (P = 0.02 and P = 0.02, respectively), and receiving three stMTX doses (P = 0.03 and P = 0.02, respectively).

CONCLUSIONS

The clinical impact of TAC concentrations varied according to patient characteristics, including disease malignancy, conditioning regimens, donor sources, and GVHD prophylaxis. These results suggest that TAC management needs to be based on patient profiles.

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.

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.

Population Pharmacokinetic Model of Methotrexate in Brazilian Pediatric Patients with Acute Lymphoblastic Leukemia

OBJECTIVES

Methotrexate (MTX) is subject to therapeutic drug monitoring because of its high pharmacokinetic variability and safety risk outside the therapeutic window. This study aimed to develop a population pharmacokinetic model (popPK) of MTX for Brazilian pediatric acute lymphoblastic leukemia (ALL) patients who attended the Hospital de Clínicas de Porto Alegre, Brazil.

METHODS

The model was developed using NONMEM 7.4 (Icon®), ADVAN3 TRANS4, and FOCE-I. To explain inter-individual variability, we evaluated covariates from demographic, biochemical, and genetic data (single nucleotide polymorphisms [SNPs] related to the transport and metabolism of drugs).

RESULTS

A two-compartment model was built using 483 data points from 45 patients (0.33-17.83 years of age) treated with MTX (0.25-5 g/m2) in different cycles. Serum creatinine (SCR), height (HT), blood urea nitrogen (BUN) and a low BMI stratification (according to the z-score defined by the World Health Organization [LowBMI]) were added as clearance covariates. The final model described MTX clearance as [Formula: see text]. In the two-compartment structural model, the central and peripheral compartment volumes were 26.8 L and 8.47 L, respectively, and the inter-compartmental clearance was 0.218 L/h. External validation of the model was performed through a visual predictive test and metrics using data from 15 other pediatric ALL patients.

CONCLUSION

The first popPK model of MTX was developed for Brazilian pediatric ALL patients, which showed that inter-individual variability was explained by renal function and factors related to body size.

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.

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.

Population Pharmacokinetics of High-Dose Methotrexate in Chinese Pediatric Patients With Acute Lymphoblastic Leukemia

High-dose methotrexate (HD-MTX) is widely used in pediatric acute lymphoblastic leukemia (ALL) treatment regimens. In this study, we aimed to develop a population pharmacokinetic (PK) model of HD-MTX in Chinese pediatric patients with ALL for designing personalized dosage regimens. In total, 4,517 MTX serum concentration data for 311 pediatric patients with ALL, aged 0.75–15.2 years and under HD-MTX treatment, were retrospectively collected at a tertiary Children’s Hospital in China. The non-linear mixed-effect model was used to establish the population PK model, using NONMEM software. The potential covariate effects of age, body weight, and biochemical measurements (renal and liver function) on MTX PK disposition were investigated. The model was then evaluated using goodness-of-fit, visual predictive check. MTX PK disposition was described using a three-compartment model reasonable well. Body weight, implemented as a fixed allometric function on all clearance and volume of distribution parameters, showed a substantial improvement in model fit. The final population model demonstrated that the MTX clearance estimate in a typical child with body weight of 19 kg was 6.9 L/h and the central distribution of volume estimate was 20.7 L. The serum creatinine significantly affected the MTX clearance, with a 0.97% decrease in clearance per 1 μmol/L of serum creatinine. Other covariates (e.g., age, sex, bilirubin, albumin, aspartate transaminase, concomitant medication) did not significantly affect PK properties of MTX. The proposed population PK model could describe the MTX concentration data in Chinese pediatric patients with ALL. This population PK model combined with a maximum a posteriori Bayesian approach could be used to estimate individual PK parameters, and optimize personalized MTX therapy in target patients, thus aiming to reduce toxicity and improve treatment outcomes.

Effects of ABCB1, UGT1A1, and UGT1A9 Genetic Polymorphisms on the Pharmacokinetics of Sitafloxacin Granules in Healthy Subjects

Sitafloxacin, a new fluoroquinolone, has strong antibacterial activity. We evaluated the effects of sitafloxacin granules in single‐dose and multidose cohorts and the effects of ABCB1, UGT1A1, and UGT1A9 genetic polymorphisms on the pharmacokinetics (PK) of sitafloxacin in healthy subjects. The single‐dose study included 3 fasted cohorts receiving 50, 100, and 200 mg of sitafloxacin granules and 1 cohort receiving 50 mg of sitafloxacin granules with a high‐fat meal. The multidose study included 1 cohort receiving 100 mg of sitafloxacin granules once daily for 5 days. PK parameters were calculated using noncompartmental parameters based on concentration‐time data. The genotypes for ABCB1, UGT1A1, and UGT1A9 single‐nucleotide polymorphisms were determined using Sanger sequencing. Subsequently, the association between sitafloxacin PK parameters and target single‐nucleotide polymorphisms was analyzed. Sitafloxacin granules were well tolerated up to 200 and 100 mg in the single‐dose and multidose studies, respectively. Sitafloxacin AUC and C_max increased linearly within the detection range, and a steady state was reached within 3 days after the administration of multiple oral doses. Our findings showed that C_max was lower in the ABCB1 (rs1045642) mutation group, whereas t_1/2 was longer in the UGT1A1 (rs2741049) and UGT1A9 (rs3832043) mutation groups. In conclusion, sitafloxacin granules were safe at single doses and multiple doses up to 200 and 100 mg/day, respectively, with a linear plasma PK profile. However, ABCB1 (rs1045642), UGT1A1 (rs2741049), and UGT1A9 (rs3832043) genetic polymorphisms are likely to influence the C_max or t_1/2 and thereby merit further clinical evaluation.

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.

Ethnogeographic and inter-individual variability of human ABC transporters

ATP-binding cassette (ABC) transporters constitute a superfamily of 48 structurally similar membrane transporters that mediate the ATP-dependent cellular export of a plethora of endogenous and xenobiotic substances. Importantly, genetic variants in ABC genes that affect gene function have clinically important effects on drug disposition and can be predictors of the risk of adverse drug reactions and efficacy of chemotherapeutics, calcium channel blockers, and protease inhibitors. Furthermore, loss-of-function of ABC transporters is associated with a variety of congenital disorders. Despite their clinical importance, information about the frequencies and global distribution of functionally relevant ABC variants is limited and little is known about the overall genetic complexity of this important gene family. Here, we systematically mapped the genetic landscape of the entire human ABC superfamily using Next-Generation Sequencing data from 138,632 individuals across seven major populations. Overall, we identified 62,793 exonic variants, 98.5% of which were rare. By integrating five computational prediction algorithms with structural mapping approaches using experimentally determined crystal structures, we found that the functional ABC variability is extensive and highly population-specific. Every individual harbored between 9.3 and 13.9 deleterious ABC variants, 76% of which were found only in a single population. Carrier rates of pathogenic variants in ABC transporter genes associated with autosomal recessive congenital diseases, such as cystic fibrosis or pseudoxanthoma elasticum, closely mirrored the corresponding population-specific disease prevalence, thus providing a novel resource for rare disease epidemiology. Combined, we provide the most comprehensive, systematic, and consolidated overview of ethnogeographic ABC transporter variability with important implications for personalized medicine, clinical genetics, and precision public health.

Germline variant burden in multidrug resistance transporters is a therapy-specific predictor of survival in breast cancer patients

Multidrug resistance due to facilitated drug efflux mediated by ATP-binding cassette (ABC) transporters is a main cause for failure of cancer therapy. Genetic polymorphisms in ABC genes affect the disposition of chemotherapeutics and constitute important biomarkers for therapeutic response and toxicity. Here we correlated germline variability in ABC transporters with disease-specific survival (DSS) in 960 breast cancer (BRCA), 314 clear cell renal cell carcinoma and 325 hepatocellular carcinoma patients. We find that variant burden in ABCC1 is a strong predictor of DSS in BRCA patients, whereas candidate polymorphisms are not associated with DSS. This association is highly drug-specific for subgroups treated with the MRP1 substrates cyclophosphamide (log-rank p = 0.0011) and doxorubicin (log-rank p = 0.0088) independent of age and tumor stage, whereas no association was found in individuals treated with tamoxifen (log-rank p = 0.13). Structural mapping of significant variants revealed multiple variants at residues involved in protein stability, cofactor stabilization or substrate binding. Our results demonstrate that BRCA patients with high variant burden in ABCC1 are less prone to respond appropriately to pharmacological therapy with MRP1 substrates, thus incentivizing the consideration of genomic germline data for precision cancer medicine.

Population Pharmacokinetics of Voriconazole in Chinese Patients with Hematopoietic Stem Cell Transplantation

BACKGROUND AND OBJECTIVE

Voriconazole is widely recommended for the prevention and treatment of invasive fungal infections in hematopoietic stem cell transplantation patients. However, its use is limited by a narrow therapeutic range and large inter-individual variability. This study aimed to characterize the pharmacokinetics of voriconazole in Chinese hematopoietic stem cell transplantation patients, to explore factors affecting its pharmacokinetic parameters, and to provide recommendations for its optimal dosing regimens.

METHODS

A total of 121 serum concentration samples from 23 patients were retrospectively included. Voriconazole concentrations were detected, and patient clinical data were recorded. Population pharmacokinetic analysis was performed by a non-linear, mixed-effect modeling approach. Goodness-of-fit plots, bootstrap method, prediction-corrected visual predictive check and external validation by an independent group of seven patients were performed to evaluate the final model.

RESULTS

A one-compartment model with first-order elimination successfully described the data. The absorption rate constant was fixed at 1.1 h^-1 and bioavailability was fixed at 0.895. The typical values for voriconazole clearance and distribution volume were 9.52 L/h and 155 L, respectively. CYP2C19*2 genotype and mycophenolate mofetil combination presented a significant impact on the clearance. Compared with CYP2C19*2 carriers, voriconazole clearance was proven to be higher in CYP2C19*1/*1 patients.

CONCLUSIONS

A population pharmacokinetic model of voriconazole was successfully established in Chinese hematopoietic stem cell transplantation patients. Based on the final model, CYP2C19*2 genotyping coupled with therapeutic drug monitoring seems to be useful to guide voriconazole dosing and to explain subtherapeutic concentrations in clinical practice.

Pharmacokinetic basis for dosing high-dose methotrexate in infants and young children with malignant brain tumours

AIMS

No population pharmacokinetic studies of high-dose methotrexate (HDMTX) have been conducted in infants with brain tumours, which are a vulnerable population. The aim of this study was to evaluate HDMTX disposition in these children to provide a rational basis for MTX dosing.

METHODS

Patients received 4 monthly courses of HDMTX (5 g/m² or 2.5 g/m² for infants aged ≤31 days) as a 24-h infusion. Serial samples were analysed for MTX by an enzyme immunoassay method. Pharmacokinetic parameters were estimated using nonlinear mixed effects population modelling. Demographics, concomitant medications and genetic polymorphisms were considered as pharmacokinetic covariates while MTX exposure and patient age were considered as covariates for Grade 3 and 4 toxicities.

RESULTS

The population pharmacokinetics of HDMTX were estimated in 178 patients (age range 0.02-4.7 years) in 648 courses. The population clearance and volume were 90 mL/min/m² and 14.4 L/m² , respectively. Significant covariates on body surface area adjusted MTX clearance included estimated glomerular filtration rate and co-treatment with dexamethasone or vancomycin. No significant association was observed between MTX toxicity and MTX exposure, patient age, leucovorin dosage or duration. MTX clearance in infants ≤31 days at enrolment was 44% lower than in older infants, but their incidence of toxicity was not higher since they also received a lower MTX dosage.

CONCLUSIONS

By aggressively following institutional clinical guidelines, HDMTX-related toxicities were low, and using covariates from the population pharmacokinetic model enabled the calculation of a rational dosage for this patient population for future clinical trials.

Variants of FasL and ABCC5 are predictive of outcome after chemotherapy-based treatment in osteosarcoma

Objectives

Previous pharmacogenetics studies showed that genetic variants could be indicative of the response to chemotherapy. We aimed to investigate whether variants of FasL, MSH2, ABCC5, CASP3 and CYP3A4 are associated with the outcome after chemotherapy-based treatment in osteosarcoma.

Methods

132 osteosarcoma patients who had completed the neoadjuvant chemotherapy in our center were included. 5-year progression-free survival (PFS) was assessed from the initial treatment to the earliest sign of disease progression or death from any cause. 5 SNPs were genotyped using TaqMan SNP Genotyping Assay, including rs763110 of FasL, rs4638843 of MSH2, rs939338 of ABCC5, rs2720376 of CASP3 and rs4646437 of CYP3A4. Patients were classified into two groups according to the 5-year PFS (event/no event). The chi-square test was used to analyze difference of genotype frequency. Logistic regression analysis was used to determine the independent predictors of the PFS rate.

Results

The overall 5-year PFS was 61.4% (81/132). Genotype TT/CT of rs763110 and genotype GG/AG of rs939338 were significantly associated with the event of 5-year PFS (p = 0.028 for rs763110; p = 0.039 for rs939338). Patients with no risk allele showed a 5-year PFS of 73.7% (42/57), which was significantly higher than a PFS of 54.2% (26/48) for patients with one risk allele and 48.1% (13/27) for patients with two different risk alleles (p = 0.03). Logistic regression analysis showed that allele T of FasL rs763110 and allele G of ABCC5 rs939338 were independent risk factors of the 5-year PFS. The ORs were 2.14 (95%CI = 1.13–3.35, p = 0.01) for rs763110 and 1.73 (95%CI = 1.05–2.52, p = 0.03) for rs939338, respectively.

Conclusions

The association of variants of FASL and ABCC5 with survival outcome after chemotherapy was validated in patients with osteosarcoma. Our findings may provide a new insight into a more personalized treatment for patients with osteosarcoma.

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

The intra- and inter-individual variances of methotrexate (MTX) pharmacokinetics are extremely large, and the pharmacokinetic property of MTX in patients with primary central nervous system lymphoma (PCNSL) is unestablished. A total of 701 MTX plasma concentrations from 98 patients with PCNSL under high-dose MTX therapy were used to develop the population pharmacokinetic (popPK) model of MTX by using the nonlinear mixed-effects modeling method. A 2-compartment model was employed to describe the pharmacokinetic property of MTX. In the final popPK model, inclusion of serum creatinine and body surface area significantly reduced objective function value for clearance over the base model (p <0.001), and inclusion of age significantly reduced objective function value for distribution volume of central compartment (V_c) over the base model (p <0.001). In the final popPK model, the inter-individual clearance = 6.67 × (SCR/68.1)^-0.48 × (BSA/1.75)^1.17; V_c = 24.46 × (age/57.16)^0.81. The precision of all parameters was acceptable (relative standard error <28.61%). Bootstrap and visual predictive check results indicated that the final popPK model was stable with acceptable predictive ability. The popPK model may be useful for personalized medication in PCNSL patients under high-dose MTX therapy. Further studies are warranted to confirm the results.

Pharmacokinetics, Pharmacodynamics and Pharmacogenomics of Immunosuppressants in Allogeneic Haematopoietic Cell Transplantation: Part I

Although immunosuppressive treatments and target concentration intervention (TCI) have significantly contributed to the success of allogeneic haematopoietic cell transplantation (alloHCT), there is currently no consensus on the best immunosuppressive strategies. Compared with solid organ transplantation, alloHCT is unique because of the potential for bidirectional reactions (i.e. host-versus-graft and graft-versus-host). Postgraft immunosuppression typically includes a calcineurin inhibitor (cyclosporine or tacrolimus) and a short course of methotrexate after high-dose myeloablative conditioning, or a calcineurin inhibitor and mycophenolate mofetil after reduced-intensity conditioning. There are evolving roles for the antithymyocyte globulins (ATGs) and sirolimus as postgraft immunosuppression. A review of the pharmacokinetics and TCI of the main postgraft immunosuppressants is presented in this two-part review. All immunosuppressants are characterized by large intra- and interindividual pharmacokinetic variability and by narrow therapeutic indices. It is essential to understand immunosuppressants' pharmacokinetic properties and how to use them for individualized treatment incorporating TCI to improve outcomes. TCI, which is mandatory for the calcineurin inhibitors and sirolimus, has become an integral part of postgraft immunosuppression. TCI is usually based on trough concentration monitoring, but other approaches include measurement of the area under the concentration-time curve (AUC) over the dosing interval or limited sampling schedules with maximum a posteriori Bayesian personalization approaches. Interpretation of pharmacodynamic results is hindered by the prevalence of studies enrolling only a small number of patients, variability in the allogeneic graft source and variability in postgraft immunosuppression. Given the curative potential of alloHCT, the pharmacodynamics of these immunosuppressants deserves to be explored in depth. Development of sophisticated systems pharmacology models and improved TCI tools are needed to accurately evaluate patients' exposure to drugs in general and to immunosuppressants in particular. Sequential studies, first without and then with TCI, should be conducted to validate the clinical benefit of TCI in homogenous populations; randomized trials are not feasible, because there are higher-priority research questions in alloHCT. In Part I of this article, we review the alloHCT process to facilitate optimal design of pharmacokinetic and pharmacodynamics studies. We also review the pharmacokinetics and TCI of calcineurin inhibitors and methotrexate.

Role of Pharmacogenetics in Hematopoietic Stem Cell Transplantation Outcome in Children

Hematopoietic stem cell transplantation (HSCT) is an established therapeutic procedure for several congenital and acquired disorders, both malignant and nonmalignant. Despite the great improvements in HSCT clinical practices over the last few decades, complications, such as graft vs. host disease (GVHD) and sinusoidal obstructive syndrome (SOS), are still largely unpredictable and remain the major causes of morbidity and mortality. Both donor and patient genetic background might influence the success of bone marrow transplantation and could at least partially explain the inter-individual variability in HSCT outcome. This review summarizes some of the recent studies on candidate gene polymorphisms in HSCT, with particular reference to pediatric cohorts. The interest is especially focused on pharmacogenetic variants affecting myeloablative and immunosuppressive drugs, although genetic traits involved in SOS susceptibility and transplant-related mortality are also reviewed.

An analysis of the therapeutic benefits of genotyping in pediatric hematopoietic stem cell transplantation

Hematopoietic stem cell transplantation is a high-risk procedure that is offered, with curative intent, to patients with malignant and nonmalignant disease. The clinical benefits of personalization of therapy by genotyping have been demonstrated by the reduction in transplant related mortality from donor-recipient HLA matching. However, defining the relationship between genotype and transplant conditioning agents is yet to be translated into clinical practice. A number of the therapeutic agents used in stem cell transplant preparative regimens have pharmacokinetic parameters that predict benefit of incorporating pharmacogenomic data into dosing strategies. Busulfan, cyclophosphamide, thio-TEPA and etoposide have well-described drug metabolism pathways, however candidate gene studies have identified there is a gap in the identification of pharmacogenomic data that can be used to improve transplant outcomes. Incorporating pharmacogenomics into pharmacokinetic modeling may demonstrate the therapeutic benefits of genotyping in transplant preparative regimen agents.

Genotype and allele frequencies of drug-metabolizing enzymes and drug transporter genes affecting immunosuppressants in the Spanish white population

Interpatient variability in drug response can be widely explained by genetically determined differences in metabolizing enzymes, drug transporters, and drug targets, leading to different pharmacokinetic and/or pharmacodynamic behaviors of drugs. Genetic variations affect or do not affect drug responses depending on their influence on protein activity and the relevance of such proteins in the pathway of the drug. Also, the frequency of such genetic variations differs among populations, so the clinical relevance of a specific variation is not the same in all of them. In this study, a panel of 33 single nucleotide polymorphisms in 14 different genes (ABCB1, ABCC2, ABCG2, CYP2B6, CYP2C19, CYP2C9, CYP3A4, CYP3A5, MTHFR, NOD2/CARD15, SLCO1A2, SLCO1B1, TPMT, and UGT1A9), encoding for the most relevant metabolizing enzymes and drug transporters relating to immunosuppressant agents, was analyzed to determine the genotype profile and allele frequencies in comparison with HapMap data. A total of 570 Spanish white recipients and donors of solid organ transplants were included. In 24 single nucleotide polymorphisms, statistically significant differences in allele frequency were observed. The largest differences (>100%) occurred in ABCB1 rs2229109, ABCG2 rs2231137, CYP3A5 rs776746, NOD2/CARD15 rs2066844, TPMT rs1800462, and UGT1A9 rs72551330. In conclusion, differences were recorded between the Spanish and other white populations in terms of allele frequency and genotypic distribution. Such differences may have implications in relation to dose requirements and drug-induced toxicity. These data are important for further research to help explain interindividual pharmacokinetic and pharmacodynamic variability in response to drug therapy.

Cyclosporine and methotrexate-related pharmacogenomic predictors of acute graft-versus-host disease

Effective immunosuppression is mandatory to prevent graft-versus-host disease and to achieve a successful clinical outcome of hematopoietic stem cell transplantation. Here we tested whether germline single nucleotide polymorphisms in 20 candidate genes related to methotrexate and cyclosporine metabolism and activity influence the incidence of graft-versus-host disease in patients who undergo stem cell transplantation for hematologic disorders. Recipient genetic status of the adenosine triphosphate-binding cassette sub-family C1 and adenosine triphosphate-binding cassette sub-family C2 transporters, 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/ inosine monophosphate cyclohydrolase within the methotrexate pathway, and nuclear factor of activated T cells (cytoplasmic 1) loci exhibit a remarkable influence on severe acute graft-versus-host disease prevalence. Indeed, an increased risk of acute graft-versus-host disease was observed in association with single nucleotide polymorphisms located in 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/inosine monophosphate cyclohydrolase (hazard ratio=3.04; P=0.002), nuclear factor of activated T cells (cytoplasmic 1) (hazard ratio=2.69; P=0.004), adenosine triphosphate-binding cassette sub-family C2 (hazard ratio=3.53; P=0.0018) and adenosine triphosphate-binding cassette sub-family C1 (hazard ratio=3.67; P=0.0005). While donor single nucleotide polymorphisms of dihydrofolate reductase and solute carrier family 19 (member 1) genes are associated with a reduced risk of acute graft-versus-host disease (hazard ratio=0.32-0.41; P=0.0009-0.008), those of nuclear factor of activated T cells (cytoplasmic 2) are found to increase such risk (hazard ratio=3.85; P=0.0004). None of the tested single nucleotide polymorphisms was associated with the occurrence of chronic graft-versus-host disease. In conclusion, by targeting drug-related biologically relevant genes, this work emphasizes the potential role of germline biomarkers in predicting acute graft-versus-host disease. Further investigations are warranted to improve our understanding of these relationships to personalize immunosuppressive therapy and optimize outcomes.

A population pharmacokinetic analysis of the influence of nutritional status of digoxin in hospitalized Korean patients

BACKGROUND

Safe and effective use of digoxin in hospitalized populations requires information about the drug's pharmacokinetics and the influence of various factors on drug disposition. However, no attempts have been made to link an individual's digoxin requirements with nutritional status.

OBJECTIVES

The main goal of this study was to estimate the population pharmacokinetics of digoxin and to identify the nutritional status that explains pharmacokinetic variability in hospitalized Korean patients.

METHODS

Routine therapeutic drug-monitoring data from 106 patients who received oral digoxin at Seoul National University Bundang Hospital were retrospectively collected. The pharmacokinetics of digoxin were analyzed with a 1-compartment, open-label pharmacokinetic model by using a nonlinear mixed-effects modeling tool (NONMEM) and a multiple trough screening approach.

RESULTS

The effect of demographic characteristics and biochemical and nutritional indices were explored. Estimates generated by using NONMEM indicated that the CL/F of digoxin was influenced by renal function, serum potassium, age, and percentage of ideal body weight (PIBW). These influences could be modeled by following the equation CL/F (L/h) = 1.36 × (creatinine clearance/50)(1.580) × K(0.835) × 0.055 × (age/65) × (PIBW/100)(0.403). The interindividual %CV for CL/F was 34.3%, and the residual variability (SD) between observed and predicted concentrations was 0.225 μg/L. The median estimates from a bootstrap procedure were comparable and within 5% of the estimates from NONMEM. Correlation analysis with the validation group showed a linear correlation between observed and predicted values.

CONCLUSIONS

The use of this model in routine therapeutic drug monitoring requires that certain conditions be met which are consistent with the conditions of the subpopulations in the present study. Therefore, further studies are needed to clarify the effects of nutritional status on digoxin pharmacokinetics. The present study established important sources of variability in digoxin pharmacokinetics and highlighted the relationship between pharmacokinetic parameters and nutritional status in hospitalized Korean patients.