Population pharmacokinetics of mycophenolic acid in paediatric patients

Model-informed dose optimization of mycophenolic acid in pediatric kidney transplant patients

PURPOSE

We aimed to develop and evaluate a population PK model of mycophenolic acid (MPA) in pediatric kidney transplant patients to aid MPA dose optimization.

METHODS

Data were collected from pediatric kidney transplant recipients from a Dutch academic hospital (Radboudumc, the Netherlands). Pharmacokinetic model-building and model-validation analyses were performed using NONMEM. Subsequently, we externally evaluated the final model using data from another academic hospital. The final model was used to develop an optimized dosing regimen.

RESULTS

Thirty pediatric patients were included of whom 266 measured MPA plasma concentrations, including 20 full pharmacokinetic (PK) curves and 24 limited sampling curves, were available. A two-compartment model with a transition compartment for Erlang-type absorption best described the data. The final population PK parameter estimates were Ktr (1.48 h-1; 95% CI, 1.15-1.84), CL/F (16.0 L h-1; 95% CI, 10.3-20.4), Vc/F (24.9 L; 95% CI, 93.0-6.71E25), Vp/F (1590 L; 95% CI, 651-2994), and Q/F (36.2 L h-1; 95% CI, 9.63-74.7). The performance of the PK model in the external population was adequate. An optimized initial dose scheme based on bodyweight was developed. With the licensed initial dose, 35% of patients were predicted to achieve the target AUC, compared to 42% using the optimized scheme.

CONCLUSION

We have successfully developed a pharmacokinetic model for MPA in pediatric renal transplant patients. The optimized dosing regimen is expected to result in better target attainment early in treatment. It can be used in combination with model-informed follow-up dosing to further individualize the dose when PK samples become available.

Limited sampling strategy and population pharmacokinetic model of mycophenolic acid in pediatric patients with systemic lupus erythematosus: application of a double gamma absorption model with SAEM algorithm

INTRODUCTION

Mycophenolic acid (MPA), the active metabolite of mycophenolate mofetil (MMF), is widely used in the treatment of systemic lupus erythematosus (SLE). It has been shown that its therapeutic drug monitoring based on the area under the curve (AUC) improves treatment efficacy. MPA exhibits a complex bimodal absorption, and a double gamma distribution model has been already proposed in the past to accurately describe this phenomenon. These previous population pharmacokinetics models (POPPK) have been developed using iterative two stage Bayesian (IT2B) or non-parametric adaptive grid (NPAG) methods. However, non-linear mixed effect (NLME) approaches based on stochastic approximation expectation-maximization (SAEM) algorithms have never been published so far for this particular model. The objectives of this study were (i) to implement the double absorption gamma model in Monolix, (ii) to compare different absorption models to describe the pharmacokinetics of MMF, and (iii) to develop a limited sampling strategy (LSS) to estimate AUC in pediatric SLE patients.

MATERIAL AND METHODS

A data splitting of full pharmacokinetic profiles sampled in 67 children extracted either from the expert system ISBA (n = 34) or the hospital Saint Louis (n = 33) was performed into train (75%) and test (25%) sets. A POPPK was developed for MPA in the train set using a NLME and the SAEM algorithm and different absorption models were implemented and compared (first order, transit, or simple and double gamma). The best limited sampling strategy was then determined in the test set using a maximum-a-posteriori Bayesian method to estimate individual PK parameters and AUC based on three blood samples compared to the reference AUC calculated using the trapezoidal rule applied on all samples and performances were assessed in the test set.

RESULTS

Mean patient age and dose was 13 years old (5-18) and 18.1 mg/kg (7.9-47.6), respectively. MPA concentrations (764) from 107 occasions were included in the analysis. A double gamma absorption with a first-order elimination from the central compartment best fitted the data. The optimal LSS with samples at 30 min, 2 h, and 3 h post-dose exhibited good performances in the test set (mean bias - 0.32% and RMSE 21.0%).

CONCLUSION

The POPPK developed in this study adequately estimated the MPA AUC in pediatric patients with SLE based on three samples. The double absorption gamma model developed with the SAEM algorithm showed very accurate fit and reduced computation time.

Body Surface Area-Based Dosing of Mycophenolate Mofetil in Pediatric Hematopoietic Stem Cell Transplant Recipients: A Prospective Population Pharmacokinetic Study

Mycophenolate mofetil (MMF) is commonly used for acute graft-versus-host disease (aGVHD) after allogeneic hematopoietic stem cell transplantation (HSCT). However, limited population pharmacokinetic (PPK) data are available for pediatric HSCT patients. This study aimed to develop a PPK model and recommend optimal oral MMF dosage in pediatric HSCT patients. This prospective study involved pediatric HSCT patients at a tertiary academic institution. Patients received oral MMF 15-20 mg/kg twice daily for aGVHD prophylaxis and treatment. The PPK analysis was conducted using a nonlinear mixed-effects modeling method. Simulation was performed considering different body surface areas (BSAs) (0.5 m2, 1.0 m2, 1.5 m2) and dosing (400 mg/m2, 600 mg/m2, 900 mg/m2 twice daily). Based on the simulation, an optimal dosage of oral MMF was suggested. A total of 20 patients and 80 samples were included in the PPK model development. A one-compartment model with first-order absorption adequately described the pharmacokinetics of mycophenolic acid (MPA). BSA was a statistically significant covariate on Vd/F. Simulation suggested the optimal dosage of oral MMF as 900 mg/m2 twice daily, respectively. A reliable PPK model was developed with good predictive performance. This model-informed optimal MMF dosage in pediatric HSCT patients can provide valuable dosing guidance in real-world clinical practice.

Significant Effects of Renal Function on Mycophenolic Acid Total Clearance in Pediatric Kidney Transplant Recipients with Population Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is an immunosuppressant commonly prescribed in pediatric kidney transplantation to prevent graft rejection. Large variabilities in MPA plasma exposures have been observed in this population, which could result in severe adverse effects. The majority of the MPA pharmacokinetic data have been reported in adult populations, whereas information in pediatric patients is still very limited. The objective of this study was to establish a novel, nonlinear mixed-effects model for MPA and investigate the clinical variables affecting MPA population pharmacokinetics in pediatric kidney transplant recipients.

METHODS

Data were collected retrospectively from pediatric kidney transplant patients (≤ 18 years when MPA concentrations were initially collected; on oral administration of mycophenolate mofetil) in Calgary, Alberta, Canada. Nonlinear mixed-effect modeling was conducted using stochastic approximation expectation-maximization in Monolix 2021R2 (Lixoft SAS, France) to determine population pharmacokinetic estimates, interindividual variabilities, and interoccasional variabilities. Covariate models were constructed using the Model Proposal function in Monolix in conjunction with a systematic stepwise inclusion/elimination protocol. The best model was selected based on objective function values, relative standard errors, goodness-of-fit plots, prediction-corrected visual predictive checks, and numerical predictive checks.

RESULTS

A total of 50 pediatric kidney transplant patients (25 female) with 219 MPA plasma concentration-time profiles were included. The average age (± standard deviation) and posttransplant time for the sample population were 12.8 ± 4.8 years and 762 ± 1160 days, respectively. The majority of study subjects (i.e., > 85% based on all occasions) were co-administered tacrolimus. A two-compartment, first-order absorption with lag time and linear elimination structural model with lognormal distributed proportional residual errors best described the MPA concentration-time data. The absorption rate constant (2.52 h-1 or 0.042 min-1), lag time (0.166 h or 9.96 min), volumes of distributions of the central (22.8 L) and peripheral (216 L) compartments, and intercompartment clearance (17.6 L h-1 or 0.293 L min-1) were consistent with literature values; whereas total MPA clearance (0.72 L h-1 or 0.012 L min-1) was relatively reduced, likely due to the general lack of cyclosporine interactions and the stabilized graft functions from significantly longer posttransplant time in our sample population. Of the clinical variables tested, only estimated glomerular filtration rate (eGFR) was identified a significant covariate affecting total MPA clearance with a positive, exponential relationship. The final population pharmacokinetic model was successfully evaluated/validated using a variety of complementary methods.

CONCLUSION

We have successfully constructed and validated a novel population pharmacokinetic model of MPA in pediatric kidney transplant patients. A positive, nonlinear relationship between eGFR and total MPA clearance identified in our model is likely attributed to multiple concurrent mechanisms, which warrant further systematic investigations.

Identification of biomarkers that predict renal flare in childhood-onset lupus nephritis with mycophenolate acid

BACKGROUND

Children-onset lupus nephritis (cLN) occurs > 50 % of patients with systemic lupus erythematosus. Mycophenolic acid (MPA) is the first-line agent for the induction and maintenance treatment of LN. This study was to explore the predictors of renal flare in cLN.

METHODS

Data from 90 patients were included in population pharmacokinetic (PK) models to predict MPA exposure. Cox regression models and restricted cubic spline were performed in 61 patients to identify the risk factors for renal flare, baseline clinical characteristics and MPA exposures as potential covariates.

RESULTS

PK best fitted a two-compartment model of first-order absorption and linear elimination, with delayed absorption. Clearance increased with weight and immunoglobulin G (IgG), but decreased with albumin and serum creatinine. During follow-up of 1040 (658-1359) days, 18 patients experienced a renal flare, after a median time of 932.5 (663.5-1316) days. Each 1-mg·h/L increase of MPA-AUC was associated with a 6 % decreased risk of an event (HR = 0.94; 95 % CI: 0.90-0.98), while IgG significantly increased this risk (HR = 1.17; 95 % CI: 1.08-1.26). ROC analysis showed that MPA-AUC_0-12h < 35 mg·h/L and IgG > 17.6 g/L had a good prediction of renal flare. Of restricted cubic spline, the risk of renal flares decreased with higher MPA exposure but reached a plateau when AUC_0-12h > 55 mg·h/L, while substantially increases when IgG is > 18.2 g/L.

CONCLUSIONS

Monitoring MPA exposure together with IgG could be very useful during clinical practice to identify patients with a potential high risk of renal flare. This early risk assessment would allow for the treat-to-target and tailored medicine.

Limited sampling strategy to predict free mycophenolic acid area under the concentration-time curve in paediatric patients with nephrotic syndrome

In paediatric patients, there is no data on the recommended area under the concentration-time curve from 0 to 12 h (AUC_0-12 ) for free mycophenolic acid (fMPA), which is the active form of the drug, responsible for the pharmacological effect. We decided to establish the limited sampling strategy (LSS) for fMPA for its use in MPA therapeutic monitoring in children with nephrotic syndrome treated with mycophenolate mofetil (MMF). This study included 23 children (aged 11 ± 4 years) from whom eight blood samples were collected within 12 h after MMF administration. The fMPA was determined using the high-performance liquid chromatography with fluorescence detection method. LSSs were estimated with the use of R software and bootstrap procedure. The best model was chosen based on a number of profiles with AUC predicted within ± 20% of AUC_0-12 (good guess), r² , mean prediction error (%MPE) of ±10% and mean absolute error (%MAE) of less than 25%. The fMPA AUC_0-12 was 0.1669 ± 0.0697 μg h/mL and the free fraction was within 0.16%-0.81%. In total, there were 92 equations developed of which five fulfilled the acceptance criteria for %MPE, %MAE, good guess >80% and r²  > 0.900. These equations consisted of three time points: model 1 (C₁ , C₂ , C₆ ), model 2 (C₁ , C₃ , C₆ ), model 3 (C₁ , C₄ , C₆ ), model 5 (C₀ , C₁ , C₂ ), and model 6 (C₁ , C₂ , C₉ ). Although blood sampling up to 9 h after MMF dosing is impractical, it is crucial to include C₆ or C₉ in LSS to assess fMPA AUC_pred correctly. The most practical fMPA LSS, which fulfilled the acceptance criteria in the estimation group, was fMPA AUC_pred  = 0.040 + 2.220 × C₀  + 1.130 × C₁  + 1.742 × C₂ . Further studies should define the recommended fMPA AUC_0-12 value in children with nephrotic syndrome.

Clinical Evidence on the Purported Pharmacokinetic Interactions between Corticosteroids and Mycophenolic Acid

Corticosteroids (steroids) are commonly used concurrently with mycophenolic acid (MPA) as the first-line immunosuppression therapy for the prevention of rejection in solid organ transplantations. Steroids are also commonly administered with MPA in various autoimmune disorders such as systemic lupus erythematosus and idiopathic nephrotic syndrome. Despite various review articles having suggested the presence of pharmacokinetic interactions between MPA and steroids, definitive data have not yet been demonstrated. The aim of this Current Opinion is to critically evaluate the available clinical data and propose the optimal study design for characterising the MPA-steroid pharmacokinetic interactions. The PubMed and Embase databases were searched for relevant clinical articles in English as of September 29, 2022, where a total of 8 papers have been identified as supporting and 22 as non-supporting the purported drug interaction. To objectively evaluate the data, novel assessment criteria to effectively diagnose the interaction based on known MPA pharmacology were formulated, including the availability of independent control groups, prednisolone concentrations, MPA metabolite data, unbound MPA concentrations, and the characterisations of entero-hepatic recirculation and MPA renal clearance. Overall, the majority of the identified corticosteroid data were pertaining to prednisone or prednisolone. Our assessment indicated that no conclusive mechanistic data supporting the interaction are available in the current clinical literature, and further studies are required to quantify the effects/mechanisms of steroid-tapering or withdrawal on MPA pharmacokinetics. This current opinion provides justification for further translational investigations, as this particular drug interaction has the potential to exert significant adverse outcomes in patients prescribed MPA.

Population pharmacokinetics of mycophenolate mofetil in pediatric patients early after liver transplantation

Objective: To investigate the factors influencing the pharmacokinetics of mycophenolate mofetil (MMF) in pediatric patients after liver transplantation, and to establish a population pharmacokinetics model, which can provide a reference for clinical dosage adjustment.

Methods: A prospective study in a single center was performed on pediatric patients who were administrated with mycophenolate mofetil dispersible tablets (MMFdt) for at least 4 days after liver transplantation continuously. Blood samples were collected in ethylene diamine tetraacetic acid anticoagulant tubes before dosing and 0.5, 1, 2, 4, 8, and 12 h after the morning intake of MMFdt. The concentrations of mycophenolic acid (MPA) in plasma were assayed with a validated reverse-phase high-performance liquid chromatography method. UGT1A8 518C > G, UGT1A9 -275T > A, UGT1A9 -2152C > T, UGT2B7 211G > T, SLC O 1B1 521T > C polymorphism were determined by Sanger sequencing. Nonlinear mixed effects modeling was used to establish the population pharmacokinetics (PPK) model. The predictability and stability of the model were internally evaluated by the goodness of fit plots, visual prediction check, normalized prediction errors, and bootstraps.

Results: A two-compartment model with first-order absorption and first-order elimination was established with 115 MPA concentrations from 20 pediatric patients. The final model were: CL/F (L/h) = 14.8×(WT/7.5)^0.75×(DOSE/11.16)^0.452×е^0.06, Ka (h⁻¹) = 2.02×(WT/7.5)^−0.25, Vc/F (L) = 6.01×(WT/7.5), Vp/F (L) = 269 (fixed), Q/F (L/h) = 15.4×(WT/7.5)^0.75×е^1.39. Where CL/F was the apparent clearance rate, Ka was the absorption rate constant, Vc/F was the apparent distribution volume of the central compartment, Vp/F was the apparent distribution volume of the peripheral compartment, Q/F was the atrioventricular clearance rate, WT was the body weight of the subject, and DOSE was the MMFdt administered dose. The model indicated there was large inter-individual variability in CL/F and Q/F after multiple dosing of MMFdt. Internal evaluation results showed that the final model had good stability and prediction performance.

Conclusion: A stable and predictive population pharmacokinetic model of MMFdt in pediatric patients after the early stage of liver transplantation was established. The pediatric patient’s weight and the dose of MMFdt can be a reference to adjust the MMFdt dose.

Population Pharmacokinetics of Mycophenolic Acid in Renal Transplant Patients: A Comparison of the Early and Stable Posttransplant Stages

Mycophenolic acid (MPA) is an antimetabolic immunosuppressive drug widely used in solid organ transplantation and autoimmune diseases. Pharmacokinetics (PK) of MPA demonstrates high inter- and intra-variability. The aim of this study was to compare the population PK properties of MPA in adult renal transplant patients in the early and stable post-transplant stages and to simulate an optimal dosing regimen for patients at different stages. A total of 51 patients in the early post-transplant period (1 week after surgery) and 48 patients in the stable state (5.5–10 years after surgery) were included in the study. In the two-compartment population PK model, CL/F (23.36 L/h vs. 10.25 L/h) and V/F (78.07 vs. 16.24 L) were significantly different between the two stages. The dose-adjusted area under the concentration time curve (AUC_ss,12h/dose) for patients in the early stage were significantly lower than those for patients in the stable state (40.83 ± 22.26 mg h/L vs. 77.86 ± 21.34 mg h/L; p < 0.001). According to Monte Carlo simulations, patients with 1.0–1.5 g of mycophenolate mofetil twice daily in the early phase and 0.50–0.75 g twice daily in the stable phase had a high probability of achieving an AUC_ss,12h of 30–60 mg h/L. In addition, limited sampling strategies showed that two 4-point models (C0-C1-C2-C4 and C1-C2-C3-C6) performed well in predicting MPA exposure by both Bayesian estimate and regression equation and could be applied in clinical practice to assist therapeutic drug monitoring of MPA.

Significant Correlations between p-Cresol Sulfate and Mycophenolic Acid Plasma Concentrations in Adult Kidney Transplant Recipients

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is a commonly prescribed life-long immunosuppressant for kidney transplant recipients. The frequently observed large variations in MPA plasma exposure may lead to severe adverse outcomes; therefore, characterizations of contributing factors can potentially improve the precision dosing of MPA. Our group recently reported the potent inhibitory effects of p-cresol (a protein-bound uremic toxin that can be accumulated in kidney transplant patients) on the hepatic metabolism of MPA in human in vitro models. Based on these data, the hypothesis for this clinical investigation was that a direct correlation between p-cresol and MPA plasma exposure should be evident in adult kidney transplant recipients.

METHODS

Using a prospective and observational approach, adult kidney transplant recipients within the first year after transplant on oral mycophenolate mofetil (with tacrolimus ± prednisone) were screened for recruitment. The exclusion criteria were cold ischemia time > 30 h, malignancy, pregnancy, severe renal dysfunction (i.e., estimated glomerular filtration rate, eGFR, < 10 mL/min/1.73 m²), active graft rejection, or MPA intolerance. Patients' demographic and biochemistry data were collected. Total and free plasma concentrations of MPA, MPA glucuronide (MPAG), and total p-cresol sulfate (the predominant, quantifiable form of p-cresol in the plasma) were quantified using validated assays. Correlational and categorical analyses were performed using GraphPad Prism.

RESULTS

Forty patients (11 females) were included: donor type (living/deceased: 20/20), induction regimen (basiliximab/thymoglobulin/basiliximab followed by thymoglobulin: 35/3/2), post-transplant time (74 ± 60 days, mean ± standard deviation), age (53.7 ± 12.4 years), bodyweight (79.8 ± 18.5 kg), eGFR (51.9 ± 18.0 mL/min/1.73 m²), serum albumin (3.6 ± 0.5 g/dL), prednisone dose (18.5 ± 13.2 mg, n = 33), and tacrolimus trough concentration (9.4 ± 2.4 µg/L). Based on Spearman analysis, significant control correlations supporting the validity of our dataset were observed between total MPA trough concentration (C₀) and total MPAG C₀ (correlation coefficient [R] = 0.39), ratio of total MPAG C₀-to-total MPA C₀ and post-transplant time (R = - 0.56), total MPAG C₀ and eGFR (R = - 0.35), and p-cresol sulfate concentration and eGFR (R = - 0.70). Our primary analysis indicated the novel observation that total MPA C₀ (R = 0.39), daily dose-normalized total MPA C₀ (R = 0.32), and bodyweight-normalized total MPA C₀ (R = 0.32) were significantly correlated with plasma p-cresol sulfate concentrations. Consistently, patients categorized with elevated p-cresol sulfate concentrations (i.e., ≥ median of 3.2 µg/mL) also exhibited increased total MPA C₀ (by 57 % vs those below median), daily dose-normalized total MPA C₀ (by 89 %), and bodyweight-normalized total MPA C₀ (by 62 %). Our secondary analyses with MPA metabolites, unbound concentrations, free fractions, and MPA metabolite ratios supported additional potential interacting mechanisms.

CONCLUSION

We have identified a novel, positive association between p-cresol sulfate exposure and total MPA C₀ in adult kidney transplant recipients, which is supported by published mechanistic in vitro data. Our findings confirm a potential role of p-cresol as a significant clinical variable affecting the pharmacokinetics of MPA. These data also provide the justifications for conducting subsequent full-scale pharmacokinetic-pharmacodynamic studies to further characterize the cause-effect relationships of this interaction, which could also rule out potential confounding variables not adequately controlled in this correlational study.

Recent lessons learned from population pharmacokinetic studies of mycophenolic acid: physiological, genomic, and drug interactions leading to the prediction of drug effects

INTRODUCTION

Mycophenolic acid (MPA) is a widely used immunosuppressant in transplantation and autoimmune disease. Highly variable pharmacokinetics have been observed with MPA, but the exact mechanisms remain largely unknown.

AREAS COVERED

The current review provided a critical, comprehensive update of recently published population pharmacokinetic/dynamic models of MPA (n=16 papers identified from PubMed and Embase, inclusive from January 2017 to August 2021), with specific emphases on the intrinsic and extrinsic factors influencing the pharmacology of MPA. The significance of the identified covariates, potential mechanisms, and comparisons to historical literature have been provided.

EXPERT OPINION

While select covariates affecting the population pharmacokinetics of MPA are consistently observed and mechanistically supported, some variables have not been regularly reported and/or lacked mechanistic explanation. Very few pharmacodynamic models were available, pointing to the need to extrapolate pharmacokinetic findings. Ideal models of MPA should consist of: i) utilizing optimal sampling points to allow the characterizations of absorption, re-absorption, and elimination phases; ii) characterizing unbound/total MPA, MPA metabolites, plasma/urinary concentrations, and genetic polymorphisms to facilitate mechanistic interpretations; and iii) incorporating actual outcomes and pharmacodynamic data to establish clinical relevance. We anticipate the field will continue to expand in the next 5 to 10 years.

Exposure levels of mycophenolic acid are associated with comorbidities in children with systemic lupus erythematosus

INTRODUCTION

Little is known about the relationship between exposure levels of mycophenolic acid (MPA), the active metabolite of mycophenolate mofetil (MMF), and comorbidities of systemic lupus erythematosus (SLE) in children. This study aims to explore this association.

METHODS

Longitudinal data from SLE children, who were taking MMF for immunosuppression and under therapeutic drug monitoring (TDM), were retrospectively collected. Area under the concentration-time curve of mycophenolic acid (MPA) over 24 hours (AUC_0-24h) was estimated with Bayesian methods. Logistic regression and random forest models were used to explore the association between comorbidities and MPA exposure levels.

RESULTS

This study included 107 children with 358 times of follow-up (median age 169.02 months). The incidence of diabetes, acute kidney injury (AKI), or pneumonia was significantly associated with AUC_0-24h (odds ratio [OR] 0.991, 95% confidence interval [CI] 0.982-0.999), SLE duration (OR 1.012, 95% CI 1.002-1.022), lymphocyte percentage (OR 0.959, 95% CI 0.925-0.991), plasma albumin levels (OR 0.891, 95% CI 0.843-0.940), use of aspirin (OR 0.292, 95% CI 0.126-0.633) and hydroxychloroquine (OR 0.407, 95% CI 0.184-0.906). The random forest model showed that albumin and AUC_0-24h were two important predictors. The case group (with the three comorbidities) had a mean AUC_0-24h of 73.63 mg · h/L, while the control group had a mean AUC_0-24h of 100.39 mg · h/L.

CONCLUSIONS

Increased levels of MPA exposure are associated with decreased incidence odds of diabetes, AKI or pneumonia in SLE children. An AUC_0-24h of 100.39 mg · h/L or an AUC_0-12h of 50.20 mg · h/L could be used as the targeted exposure level for clinical practice.

Mycophenolic Acid Exposure Optimization Based on Vitamin D Status in Children with Systemic Lupus Erythematosus: A Single-Center Retrospective Study

Introduction

Systemic lupus erythematosus (SLE) can affect bone metabolism and homeostasis of serum electrolytes that are associated with abnormal levels of vitamin D. Mycophenolate mofetil (MMF) is a commonly used immunosuppressant with the active metabolite mycophenolic acid (MPA). The area under the plasma concentration–time curve (AUC) of MPA is often monitored during the treatment to assess the exposure levels. This study aims to explore the association between exposure levels of MPA and 25-hydroxyvitamin D [25(OH)D] levels in children with SLE.

Methods

Repeated measured data of children with SLE who were treated with MMF and under therapeutic drug monitoring (TDM) were retrospectively collected from the electronic medical records. MPA exposure levels were reflected by the area under the concentration–time curve over 24 h (AUC_0–24h). Univariate and multivariate linear regression models were employed to analyze factors associated with 25(OH)D levels. Hierarchical linear models were developed to analyze the intra- and inter-individual effects of AUC_0–24h on the variance of 25(OH)D levels.

Results

Data from 184 children with SLE (142 female and 42 male) with 518 follow-ups were collected. The median age was 14 years (range 3–18 years) at TDM. Children with normal 25(OH)D levels had significantly higher AUC_0–24h than children with low 25(OH)D levels (98.71 vs. 84.05 mg·h/L, P = 0.004). Intra- and inter-individual effects of AUC_0–24h on 25(OH)D levels were similar (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\gamma }_{10}$$\end{document}γ10 = 0.034 vs. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\gamma }_{01}$$\end{document}γ01 = 0.037) but only the intra-individual effect was significant (P = 0.001) in hierarchical models. Other associated factors include age, sex, season at measurement, glucocorticoid daily dose, and external vitamin D₃ supplements.

Conclusion

25(OH)D levels are associated with MPA exposure levels, and may serve as a potential indicator to optimize the exposure level of MPA during treatment. AUC_0–24h of 98.71 mg·h/L or AUC_0–12h of 49.36 mg·h/L could be the targeted exposure level for children with SLE.