key: cord-0012888-2guscr3l
authors: Li, Tai-feng; Hu, Lei; Ma, Xiao-lu; Huang, Lin; Liu, Xue-mei; Luo, Xing-xian; Feng, Wan-yu; Wu, Chun-fu
title: Population pharmacokinetics of cyclosporine in Chinese children receiving hematopoietic stem cell transplantation
date: 2019-07-24
journal: Acta Pharmacol Sin
DOI: 10.1038/s41401-019-0277-x
sha: 3e141d5735efd7f8a1281e4f7b9363954c899b27
doc_id: 12888
cord_uid: 2guscr3l

Cyclosporine (CsA) is characterized by a narrow therapeutic window and high interindividual pharmacokinetic variability, particularly in juvenile patients. The aims of this study were to build a population pharmacokinetic model of CsA in Chinese children with hematopathy who received allogeneic hematopoietic stem cell transplantation (allo-HSCT) and to identify covariates affecting CsA pharmacokinetics. A total of 86 Chinese children aged 8.4 ± 3.8 years (range 1.1–16.8 years) who received allo-HSCT were enrolled. Whole blood samples were collected before allo-HSCT. Genotyping was performed using an Agena MassARRAY system. A total of 1010 trough plasma concentration values of CsA and clinical data were collected. The population pharmacokinetic model of CsA was constructed using nonlinear mixed-effects modeling (NONMEM) software. The stability and performance of the final model were validated using bootstrapping and normalized prediction distribution errors. We showed that a one-compartment model with first-order elimination adequately described the pharmacokinetics of CsA. The typical values for clearance (CL) and volume of distribution (V) were 42.3 L/h and 3100 L, respectively. Body weight, postoperative days, CYP3A4*1 G genotype, estimated glomerular filtration rate and coadministration of triazole antifungal drugs were identified as significant covariates for CL. Weight and postoperative days were significant covariates for the V of CsA. Our model can be adopted to optimize the CsA dosing regimen for Chinese children with hematopathy receiving allo-HSCT.

Allogeneic hematopoietic stem cell transplantation (HSCT) has become an effective treatment option for refractory or recurrent hematological malignancies such as leukemia, lymphoma and myelodysplastic syndrome [1, 2] . However, graft-versus-host disease (GVHD) is one of the main reasons for the failure of transplantation. Cyclosporine (CsA), a calcineurin inhibitor, has been widely used to prevent GVHD disease in patients who have undergone allo-HSCT [3] . Indeed, CsA is characterized by a narrow therapeutic window and high interindividual pharmacokinetic variability. Thus, therapeutic drug monitoring (TDM) is advised to adjust doses and reduce toxicity, especially in juvenile patients, who exhibit a wider range of pharmacokinetic variability than adults [4] . However, owing to the hysteresis of TDM, identification of the factors affecting the pharmacokinetics of CsA is essential for improving safety and efficacy.

Population pharmacokinetics (Pop-PK), which can identify and quantify the influences of different factors on drug exposure, is a powerful and potential tool for individualized dosing [5] . Although Pop-PK studies of CsA have been extensively performed in recent years, only two have been undertaken in juvenile allo-HSCT patients [6, 7] . According to reports, genetic polymorphisms have important implications for CsA pharmacokinetics [8] [9] [10] . CsA is mainly metabolized by the cytochrome P450 (CYP) isoenzymes CYP3A4 and CYP3A5 (mainly in the liver, and to a lesser extent in the kidney and intestine) [11, 12] and is also a substrate of the Pglycoprotein (encoded by the ABCB1 gene), which is involved in the absorption and elimination of CsA [13, 14] . Polymorphisms in the genes of those proteins could influence the pharmacokinetics of CsA [8, [15] [16] [17] . Because cytochrome P450 oxidoreductase (POR) is the only obligate electron donor for all hepatic microsomal CYP enzymes [18] , POR genetic polymorphisms could also affect CsA pharmacokinetics [10, 19] . In addition, the pregnane X receptor (PXR) has been identified as the key nuclear receptor regulating the expression of CYP3A4, CYP3A5, and ABCB1 [20, 21] . Furthermore, several studies have also examined the effect of other genetic polymorphisms in ABCC2 [20, 22, 23] , PPARA [19] , RELA [24] , NFKBIA [24] , TNFAIP3 [25] , FOXP3 [26] and CYP2C19 [22, 27] on CsA pharmacokinetics. However, gene polymorphism factors were not included in the existing Pop-PK studies in juvenile allo-HSCT patients conducted by Sarem et al. [6] or Ni et al. [7] . Therefore, a thorough understanding of the effects of genetic polymorphisms is essential for individualized cyclosporine therapy.

The aims of the present study were to establish a population pharmacokinetic model of CsA in children receiving allo-HSCT considering gene polymorphism and to identify the covariates influencing CsA pharmacokinetics.

Patients and data collection A total of 1010 drug monitoring data points were collected retrospectively from clinical reports of 86 child recipients who were first administered CsA and then underwent allo-HSCT because of malignant hematologic disorders at Peking University People's Hospital, Beijing, China.

Data on the following characteristics were collected from electronic medical records: age, sex, body weight (WT), height, body mass index (BMI), hemoglobin (HB), hematocrit (HCT), albumin (Alb), serum creatinine (Scr), total bilirubin (T-BIL), alanine transaminase (ALT), aspartate transaminase (AST), estimated glomerular filtration rate (eGFR), postoperative day (POD), type of disease (TOD), and concomitant medications. This study was approved by the Health Authority Ethics Committee of the Peking University People's Hospital and conducted in accordance with the Declaration of Helsinki. All patients signed an informed consent approved by Ethics Committee.

CsA administration Patients received intravenous CsA (Sandimmun, liquid solution, Neoral, Novartis Pharma, Basel, Switzerland) beginning on Day −7 or −10 of transplantation as an intermittent intravenous infusion over 2 h every 12 h with an initial dose of 2-3 mg/kg. The daily dosage was then adjusted according to trough concentration (C 0 ) to achieve a target concentration of 150-250 ng/mL. The administration route was changed to oral when patients were able to tolerate oral intake after transplantation.

CsA sampling and concentration assay All blood samples were drawn before intravenous CsA infusion in the morning, and the concentrations were trough concentrations. CsA concentrations were detected using a chemiluminescent microparticle immunoassay (CMIA) that was performed on the Architect R i2000SR platform (Abbott, Chicago, IL, USA). This detection method has a lower limit of quantitation of 30 ng/mL and a linearity ranging from 30 to 1500 ng/mL, with a coefficient variation of less than 4%.

DNA extraction and genotyping Genomic DNA was extracted from whole blood samples with a genomic DNA purification kit (Promega, Madison, WI, USA). The genomic DNA was stored at −20°C until use. Single nucleotide polymorphism (SNP) genotyping analysis for CYP3A5*3 (rs776746), CYP3A5 (rs15524), CYP3A4*1 G (rs2242480), CYP3A4 (rs4646437), CYP3A4*18B (rs28371759), ABCB1 c.1236 C > T (rs1128503), ABCB1 c.2677 G > T (rs2032582), ABCB1 c.3435 C > T (rs1045642), POR*28 (rs1057868), POR (rs17685), POR (rs3823884), POR (rs1135612), PXR C-25385T (rs3814055), PXRG24113A (rs2276706), PXR (rs6785049), PXR (rs3842689), PXR (rs1523127), ABCC2 c.3972 C > T (rs3740066), PPARA (rs4823613 A > G), RELA (rs7119750), RELA (rs11820062), NFKBIA G2758A (rs696), TNFAIP3 (rs2230926), FOXP3 (rs3761549), CYP2C19*17 (rs12248560), CYP2C19*2 (rs4244285) and CYP2C19*3 (rs4986893) was carried out using an Agena MassARRAY system (Agena Bioscience, USA). Duplicate samples and negative controls were included to ensure accuracy in genotyping.

Pharmacokinetic analysis Pharmacokinetic analysis was performed using the nonlinear mixed-effects modeling program NONMEM VII (V2.0; Icon Development Solutions, USA). Pharmacokinetic parameters and associated variability were estimated with the first-order conditional estimation (FOCE) method with the interaction.

The one-compartment model with first-order elimination was tested. Interindividual variability of the pharmacokinetic parameters was estimated using an exponential model expressed as follows: (Eq. 1), while the residual error using a mixed model expressed as follows: (Eq. 2).

where P i is the predicted value of the parameter for the i-th patient, P pop is the typical population value of the pharmacokinetic parameter, and η i is a random variable, identically distributed with a mean of 0 and variance of ω 2 , that was used to describe the variation for the individual i.

where C obs is the observed concentration and C pred is the predicted concentration. ε 1 and ε 2 are presumed to obey the distribution of N (0, σ 1 2 ) and N (0, σ 2 2 ), respectively.

Covariate analysis Based on the nonlinear relationship between body weight and drug elimination capacity in children, pharmacokinetic parameters were modified by the allometric scaling theory, which is expressed as follows:

where P i and WT i are the pharmacokinetic parameter and body weight of the i-th individual, respectively. P std is the parameter of an individual with a standard weight of 70 kg. PWR is the allometric coefficient fixed at a value of 0.75 for clearance and a value of 1 for distribution volume [7, [28] [29] [30] .

The covariates screened in this study included the following variables: age, WT, height, BMI, HB, HCT, ALT, AST, Alb, Scr, T-BIL, eGFR, POD, sex, concomitant medications, TOD, and genetic polymorphism.

A forward inclusion and backward elimination process was used to select covariates. The objective function value (OFV) provided by NONMEM was used as a diagnostic criterion. A decrease in OFV greater than 3.84 after introducing a certain covariate indicated that the covariate had a significant effect on the fitting of the model (P < 0.05). The contribution of each covariate was verified by a stepwise backward elimination (P < 0.01, the OFV value increased by more than 6.63).

The quality of the Pop-PK model was extensively evaluated using goodness-of-fit (GOF) plots including observations vs. individual predictions, observations vs. population predictions, conditional weighted residuals vs. population predictions and conditional weighted residuals vs. time.

The bootstrap method was used to evaluate the stability of the final model. During this process, partial data from the original datasets were resampled 500 times. The parameters of each subgroup dataset were calculated. Stability was assessed based on the frequency distribution of these parameters.

The normalized prediction distribution error (NPDE) method was also used to validate the final model [31] . NPDE is appropriate for the design of multiple-dose drug administration with limited sampling points. In this study, a total of 1000 simulations were constructed to generate the corresponding NPDE values for validation of the final model.

Eighty-six subjects (56 boys and 30 girls) were recruited for this population pharmacokinetic study. The mean ± SD age was 8.38 ± 3.78 years (range 1.1-16.8), the mean ± SD body weight was 31.93 ± 16.75 kg (range 6.5-69.0), and the mean ± SD height Table 2 . As shown in Table 2 , the minimum allele frequencies of these SNPs in the Chinese population were very close to those recorded in the National Center for Biotechnology Information.

Population pharmacokinetic model analysis A total of 1010 trough blood CsA concentration values were included in this population pharmacokinetic study (a mean of 12 samples per child). The sampling time of TDM was mostly A one-compartment model with mixed error adequately described the PK profile. The final model was identified: 5 covariates, including WT, POD, CYP3A4*1 G, eGFR and coadministration of triazole antifungals (TAF), were observed to significantly affect the CL, and 2 covariates, including WT and POD, significantly affected the V of CsA. Detailed information for the remaining parameter estimations is shown in Eq. 4, Eq. 5 and Table 3 . The proportional error was 22.9%, and the additive error was 30.3 ng/mL, which was relatively small. The relative standard error (RSE) of all the parameters was lower than 30%, suggesting that the parameters were estimated with good precision and that they gave reliable results.

(If the gene is 1, CL i = CL i *0.984; If the gene is 2 or 3, CL i = CL i *1.22 in Eq. 4, patients were coadministered of TAF, the TAF = 1, otherwise TAF = 0.)

The weight-normalized clearance was 0.75 L·h −1 ·kg −1 . Before and after allometric weight normalization, the comparison of the effect of age on the distribution of the model parameter clearance rate is illustrated in Fig. 1 . After incorporation of body weight with allometric scaling, the effect of age on the eta of clearance rate disappeared.

The GOF plots for the final model are shown in Fig. 2 . The success rate for the sampling of 500 bootstraps was 97.8%, and the estimated values of the parameters were similar to the median values resulting from the bootstrap method. In addition, the 95% CI contained the typical estimates from the final model and had a narrow range, indicating that the results of the model parameters were reliable and that the model was relatively stable (Table 3) . Based on the NPDE results, most of the points in the Q-Q diagram are distributed around the Y = X line (Fig. 3 ). In addition, the average value on the bar distribution graph was close to 0, and the standard deviation was close to 1, showing a near-normal distribution. These results indicate that the model accurately described the dataset.

This study is the first to describe the population pharmacokinetics of cyclosporine in Chinese children with malignant hematologic disorders who received allo-HSCT. The factors that have been reported to have a potential impact on cyclosporine metabolism were evaluated comprehensively and included developmental, biological and clinical factors as well as gene polymorphisms. Our results showed that a one-compartment model incorporating allometric scaling with first-order elimination was adequate for data modeling. Moreover, body weight, POD, CYP3A4*1 G, eGFR and coadministration of TAF had significant impacts on CL, while body weight and POD had significant impacts on V.

Currently, a weight-based dosing regimen is the most common method of administering cyclosporine in patients receiving hematopoietic stem cell transplants, which was confirmed by our study. In this study, allometrically scaled body weight was shown to be significantly associated with both CL and V, which is consistent with many previous studies [7, 32, 33] . However, the effect of age on CL disappeared after incorporating body weight into the model. This result was consistent with previous results in pediatric renal transplant patients [7, 33] . The probable reason is that body weight and age are closely related to children's growth and development.

Although genetic polymorphisms may have an effect on CsA pharmacokinetics, few studies have evaluated the effects of genetic polymorphisms on the population pharmacokinetics of CsA in allo-HSCT patients [22, 34] . Our Pop-PK research comprehensively evaluated genetic polymorphisms as covariates in juvenile allo-HSCT patients, which showed that the clearance rate of CsA in CYP3A4*1 G T allele carriers increased by 24.5% compared with that in CYP3A4*1 G CC carriers. However, we found that CYP3A5 and ABCB1 polymorphisms were not significant covariates in the Pop-PK model of CsA. Dai et al. [35] reported that the in vitro intrinsic metabolic clearance estimated from total 

Population pharmacokinetics of CsA in children with HSCT TF Li et al.

metabolite formation is approximately 2.3-fold greater for CYP3A4 than for CYP3A5. Both Dai's and our results suggested that CYP3A4 could play a more dominant role than CYP3A5 in the metabolism of CsA. Thus far, the effects of CYP3A5 and ABCB1 polymorphisms on the pharmacokinetics of CsA are still contradictory, although many studies have been performed [9, 16, 36, 37] . Differences in study population, sample size, ethnicities, concomitant medication and especially the postoperative day might be possible confounding factors. Children receiving immunosuppressive therapy often take TAF, including voriconazole, itraconazole, fluconazole, and posaconazole, to treat and prevent fungal infections. Consistent with many previous studies [34, 38, 39] , our results showed that coadministration of TAF, a class of potent inhibitors of cytochrome P450 3A enzymes [40, 41] , was significantly associated with lower CL of CsA in Chinese children who received allo-HSCT. However, none of these medicines were identified as having effects on CsA pharmacokinetics by Zhou et al. [42] . An insufficient number of patients for each drug combination may be the cause of this result.

The known CsA PK changed with the course of postoperative day. There were several conflicting studies about the influence of POD on CL or V. Both Irtan et al. [33] and Okada et al. [43] reported a statistically significant increase in CL/F with POD after renal transplantation. In contrast, Wu et al. suggested a decrease in CL/F via an increase in bioavailability with POD in kidney transplantation patients [44] . Jacobson et al. also reported similar findings after hematopoietic stem cell transplantation [45] . However, these studies all involved oral CsA administration, and no final consensus has yet been reached. Our results showed that the CL of CsA in Chinese juvenile allo-HSCT patients during the first 9 d after transplantation increased and then began to decrease significantly. The underlying cause may be that a high dose of steroids, CYP3A inhibitors, were used during the early stages before transplantation, and the dose was tapered down slowly. Thus, the increase in CL in the first few days after transplantation may be the result of a reduction in the steroid dose during that period. Then, the CL significantly decreased with POD, which was presumably related to reduced bile formation/excretion or hepatic cytochrome P450 metabolism due to transplant-or chemotherapy-related toxicity. Therefore, to avoid cyclosporinerelated toxicity, adjustment of the intravenously administered dose is necessary during the early allo-HSCT period. In addition, our study found that POD also has a significant impact on V, which was also proven by Joohan's research [16] . This tendency may be linked to the rebuilding of the recipients' hematopoietic system and immune system after transplantation. At present, the effect of renal function on cyclosporine CL is still controversial. No significant effect was discovered in children after renal and stem cell transplantation [22, 46] . Ni et al. [7] and Fanta et al. [32] reported a significant effect of serum creatinine on CL/F in both renal pretransplant and stem cell transplantation juvenile patients. Similarly, our study showed that the eGFR, which is a more accurate indicator than serum creatinine to reflect renal function, significantly affected cyclosporine clearance in Chinese juvenile allo-HSCT patients. A study in patients with kidney disease reported that renal dysfunction could decrease the activities of hepatic drugmetabolism enzymes, including CYP3A4 [47] . Thus, we speculate that renal dysfunction might decrease the CL of CsA by reducing the activity of CYP3A.

In fact, we also tried stricter statistical standards: P < 0.01 for forward inclusion and P < 0.001 for backward elimination in the modeling process. As a result, only three covariates were retained that had significant impacts on CL, body weight, POD, and coadministration of TAF, while only body weight had significant impacts on V. However, to maximize future research, some possible influencing factors in the new research design could be concerning, so we finally adopted P < 0.05 and P < 0.01 as the criteria for covariate inclusion and elimination, respectively.

Finally, the current study has some limitations. First, only C 0 , the current common practice for CsA drug monitoring, was used in this population pharmacokinetic analysis. However, the CsA concentration measured at 2 h post dose (C 2 ) has been reported to be better related to the area under the curve (AUC) in transplanted patients, which is considered to be the best marker of drug exposure for predicting the efficacy of drug therapy [48, 49] . Second, CsA concentrations were detected using immunoassays, which may overestimate CsA levels due to cross-reactivity compared with high-performance liquid chromatography-tandem mass spectrometry or highperformance liquid chromatography [50] . Third, not all of the mutation sites studied in the present research were genotyped successfully in each sample because of the quality of DNA or detection method, which may reduce the chance of finding genotype-related covariants. Fourth, the number of samples was not sufficient, and studies with larger sample sizes are needed in the future.

In conclusion, a population pharmacokinetic model of cyclosporine was developed in Chinese children receiving allo-HSCT. Body weight, POD, the CYP3A4*1 G genotype, eGFR and coadministration of TAF significantly influenced the cyclosporine CL, and weight and POD significantly influenced the cyclosporine V. These findings might be useful for physicians in optimizing the dosage of intravenously administered cyclosporine in children. 

Allogeneic haematopoietic stem cell transplantation-an overview

Immunotherapy following hematopoietic stem cell transplantation: potential for synergistic effects

Current and emerging strategies for the prevention of graftversus-host disease

Evolution of the therapeutic drug monitoring of cyclosporine

Pharmacometrics: a quantitative tool of pharmacological research

Bayesian approach for the estimation of cyclosporine area under the curve using limited sampling strategies in pediatric hematopoietic stem cell transplantation

Population pharmacokinetics of ciclosporin in Chinese children with aplastic anemia: effects of weight, renal function and stanozolol administration

CYP3A4 *18B and CYP3A5 *3 polymorphisms contribute to pharmacokinetic variability of cyclosporine among healthy Chinese subjects

Associations of ABCB1, NFKB1, CYP3A, and NR1I2 polymorphisms with cyclosporine trough concentrations in Chinese renal transplant recipients

Effect of age and allele variants of CYP3A5, CYP3A4, and POR genes on the pharmacokinetics of cyclosporin A in pediatric renal transplant recipients from Serbia

Drug interactions and interindividual variability of ciclosporin metabolism in the small intestine

CYP3A7, CYP3A5, CYP3A4, and ABCB1 genetic polymorphisms, cyclosporine concentration, and dose requirement in transplant recipients

Biochemical, cellular, and pharmacological aspects of the multidrug transporter

Pharmacogenetics of MDR1 and its impact on the pharmacokinetics and pharmacodynamics of drugs

Association of MDR1, CYP3A4*18B, and CYP3A5*3 polymorphisms with cyclosporine pharmacokinetics in Chinese renal transplant recipients

CYP3A5 polymorphism effect on cyclosporine pharmacokinetics in living donor renal transplant recipients: analysis by population pharmacokinetics

Effect of CYP3A4*1G on the fentanyl consumption for intravenous patient-controlled analgesia after total abdominal hysterectomy in Chinese Han population

Pharmacogenetics of P450 oxidoreductase: implications in drug metabolism and therapy

The influence of CYP3A, PPARA, and POR genetic variants on the pharmacokinetics of tacrolimus and cyclosporine in renal transplant recipients

A comprehensive study of polymorphisms in the ABCB1, ABCC2, ABCG2, NR1I2 genes and lymphoma risk

Role of orphan nuclear receptors in the regulation of drugmetabolising enzymes

Population pharmacokinetics of cyclosporine in hematopoietic stem cell transplant patients: consideration of genetic polymorphisms

Association of four DNA polymorphisms with acute rejection after kidney transplantation

Common genetic polymorphisms within NFkappaB-related genes and the risk of developing invasive

Pharmacogenomics and personalized medicine: a review focused on their application in the Chinese population

FOXP3rs3761548 polymorphism is associated with tacrolimus-induced acute nephrotoxicity in renal transplant patients

Genetic markers in CYP2C19 and CYP2B6 for prediction of cyclophosphamide's 4-hydroxylation, efficacy and side effects in Chinese patients with systemic lupus erythematosus

Population pharmacokinetics and pharmacogenetics of tacrolimus in de novo pediatric kidney transplant recipients

Population pharmacokinetics of micafungin in neonates and young infants

Mechanism-based concepts of size and maturity in pharmacokinetics

Metrics for external model evaluation with an application to the population pharmacokinetics of gliclazide

Developmental pharmacokinetics of ciclosporine population pharmacokinetic study in paediatric renal transplant candidates

Population pharmacokinetics and bayesian estimator of cyclosporine in pediatric renal transplant patients

Population pharmacokinetics and individualized dosage prediction of cyclosporine in allogeneic hematopoietic stem cell transplant patients

In vitro metabolism of cyclosporine A by human kidney CYP3A5

Population pharmacokinetic study of cyclosporine in Chinese renal transplant recipients

Ciclosporin population pharmacokinetics and Bayesian estimation in thoracic transplant recipients

Population pharmacokinetics of cyclosporine in cardiopulmonary transplant recipients

Population pharmacokinetics of cyclosporine A based on NONMEM in Chinese allogeneic hematopoietic stem cell transplantation recipients

Drug-interactions of azole antifungals with selected immunosuppressants in transplant patients: strategies for optimal management in clinical practice

Effect of itraconazole on the concentrations of tacrolimus and cyclosporine in the blood of patients receiving allogeneic hematopoietic stem cell transplants

Population pharmacokinetic study of cyclosporine in the hematopoietic stem cell transplant recipients

Population pharmacokinetics of cyclosporine A in Japanese renal transplant patients: comprehensive analysis in a single center

Population pharmacokinetics of cyclosporine in clinical renal transplant patients

Posttransplant day significantly influences pharmacokinetics of cyclosporine after hematopoietic stem cell transplantation

Ciclosporin kinetics in children after stem cell transplantation

Emerging evidence of the impact of kidney disease on drug metabolism and transport

Cyclosporine monitoring in stable, long-term, pediatric kidney transplant recipients: the value of C2 determination

Developmental pharmacogenetics of immunosuppressants in pediatric organ transplantation

A high-performance liquid chromatography-mass spectrometry method using a novel atmospheric pressure chemical ionization approach for the rapid simultaneous measurement of tacrolimus and cyclosporin in whole blood

The authors gratefully acknowledge financial support from the National Natural Science Foundation of China (Grant No. 8177131171) and the Beijing Municipal Natural Science Foundation (Grant No. 7192218).

TFL, LH, WYF, and CFW designed the research study. TFL, LH, XML, XXL performed the research study. TFL, XLM and LH participated in patient enrollment and genotypic and data collection. TFL, LH, and XLM analyzed the data. TFL wrote the manuscript, which was critically revised by the other author.

Conflict of interest: The authors declare that they have no conflict of interest.