Development of population PK model with enterohepatic circulation for mycophenolic acid in patients with childhood-onset systemic lupus erythematosus

Efficacy and Safety of Pharmacokinetically-Driven Dosing of Mycophenolate Mofetil for the Treatment of Pediatric Proliferative Lupus Nephritis-A Double-Blind Placebo Controlled Clinical Trial (The Pediatric Lupus Nephritis Mycophenolate Mofetil Study)

Background

The safety and efficacy of mycophenolate mofetil (MMF) for lupus nephritis (LN) treatment is established in adults and in some children. MMF is rapidly converted to the biologically active metabolite mycophenolic acid (MPA) whose pharmacokinetics (PK) is characterized by large inter- and intra-individual variability.

Methods/Design

This randomized, double-blind, active comparator, controlled clinical trial of pediatric subjects with proliferative LN compares pharmacokinetically-guided precision-dosing of MMF (MMFPK, i.e. the dose is adjusted to the target area under the concentration-time curve (AUC0-12h) of MPA ≥ 60-70 mg*h/L) and MMF dosed per body surface area (MMFBSA, i.e. MMF dosed 600 mg/m2 body surface area), with MMF dosage taken about 12 hours apart. At baseline, subjects are randomized 1:1 to receive blinded treatment with MMFPK or MMFBSA for up to 53 weeks. The primary outcome is partial clinical remission of LN (partial renal response, PRR) at week 26, and the major secondary outcome is complete renal response (CRR) at week 26. Subjects in the MMFBSA arm with PRR at week 26 will receive MMFPK from week 26 onwards, while subjects with CRR will continue MMFBSA or MMFPK treatment until week 53. Subjects who achieve PRR at week 26 are discontinued from study intervention.

Discussion

The Pediatric Lupus Nephritis Mycophenolate Mofetil (PLUMM) study will provide a thorough evaluation of the PK of MMF in pediatric LN patients, yielding a head-to-head comparison of MMFBSA and MMFPK for both safety and efficacy. This study has the potential to change current treatment recommendations for pediatric LN, thereby significantly impacting childhood-onset SLE (cSLE) disease prognosis and current clinical practice.

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.

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.

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.

Drug monitoring in systemic lupus erythematosus

Therapeutic drug monitoring (TDM) is not yet accepted by systemic lupus erythematosus (SLE) treatment guidelines. Studies in SLE, however, have proven benefit in three areas: identification of non-adherence or poor adherence; targets for clinical benefit; and ranges of toxicity. This review covers the data on three medications commonly used for SLE, drawing on studies from both the SLE and non-SLE literature.

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.

Integration of advanced methods and models to study drug absorption and related processes: An UNGAP perspective

This collection of contributions from the European Network on Understanding Gastrointestinal Absorption-related Processes (UNGAP) community assembly aims to provide information on some of the current and newer methods employed to study the behaviour of medicines. It is the product of interactions in the immediate pre-Covid period when UNGAP members were able to meet and set up workshops and to discuss progress across the disciplines. UNGAP activities are divided into work packages that cover special treatment populations, absorption processes in different regions of the gut, the development of advanced formulations and the integration of food and pharmaceutical scientists in the food-drug interface. This involves both new and established technical approaches in which we have attempted to define best practice and highlight areas where further research is needed. Over the last months we have been able to reflect on some of the key innovative approaches which we were tasked with mapping, including theoretical, in silico, in vitro, in vivo and ex vivo, preclinical and clinical approaches. This is the product of some of us in a snapshot of where UNGAP has travelled and what aspects of innovative technologies are important. It is not a comprehensive review of all methods used in research to study drug dissolution and absorption, but provides an ample panorama of current and advanced methods generally and potentially useful in this area. This collection starts from a consideration of advances in a priori approaches: an understanding of the molecular properties of the compound to predict biological characteristics relevant to absorption. The next four sections discuss a major activity in the UNGAP initiative, the pursuit of more representative conditions to study lumenal dissolution of drug formulations developed independently by academic teams. They are important because they illustrate examples of in vitro simulation systems that have begun to provide a useful understanding of formulation behaviour in the upper GI tract for industry. The Leuven team highlights the importance of the physiology of the digestive tract, as they describe the relevance of gastric and intestinal fluids on the behaviour of drugs along the tract. This provides the introduction to microdosing as an early tool to study drug disposition. Microdosing in oncology is starting to use gamma-emitting tracers, which provides a link through SPECT to the next section on nuclear medicine. The last two papers link the modelling approaches used by the pharmaceutical industry, in silico to Pop-PK linking to Darwich and Aarons, who provide discussion on pharmacometric modelling, completing the loop of molecule to man.

Precision Medicine: Towards Individualized Dosing in Pediatric Rheumatology

Despite an increase in the number of available therapeutics, many children with rheumatic disease continue to experience active inflammatory disease and treatment failure. One reason for treatment failure is the lack of dosing paradigms to account for the wide between-patient variability in drug pharmacokinetics because of developmental changes or genetic polymorphisms that effect drug absorption, distribution, metabolism, and elimination. This review highlights several strategies to optimize dosing for biologic and nonbiologic disease-modifying antirheumatic drugs, including therapeutic drug monitoring, pharmacogenomics, and the use of pharmacokinetic/pharmacodynamic modeling.

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.

Influence of SLCO1B1 521T>C, UGT2B7 802C>T and IMPDH1 -106G>A Genetic Polymorphisms on Mycophenolic Acid Levels and Adverse Reactions in Chinese Autoimmune Disease Patients

Introduction

Mycophenolate mofetil (MMF), a new type of immunosuppressant, has emerged as a frontline agent for treating autoimmune diseases. Mycophenolic acid (MPA) is an active metabolite of MMF. MPA exposure varies greatly among individuals, which may lead to adverse drug reactions such as gastrointestinal side effects, infection, and leukopenia. Genetic factors play an important role in the variation of MPA levels and its side effects. Although many published studies have focused on MMF use in patients after organ transplant, studies that examine the use of MMF in patients with autoimmune diseases are still lacking.

Methods

This study will not only explore the genetic factors affecting MPA levels and adverse reactions but also investigate the relationships between UGT1A9 −118(dT)9/10, UGT1A9 - 1818T>C, UGT2B7 802C>T, SLCO1B1 521T>C, SLCO1B3 334T>G, IMPDH1 −106G>A and MPA trough concentration (MPA C₀), along with adverse reactions among Chinese patients with autoimmune diseases. A total of 120 patients with autoimmune diseases were recruited. The MPA trough concentration was detected using the enzyme multiplied immunoassay technique (EMIT). Genotyping was performed using a real-time polymerase chain reaction (PCR) system and validated allelic discrimination assays. Clinical data were collected for the determination of side effects.

Results

SLCO1B1 521T>C demonstrated a significant association with MPA C₀/d (p=0.003), in which patients with the CC type showed a higher MPA C₀/d than patients with the TT type (p=0.001) or the CT type (p=0.000). No significant differences were found in MPA C₀/d among the other SNPs. IMPDH1 −106G>A was found to be significantly related to infections (p=0.006). Subgroup analysis revealed that UGT2B7 802C>T was significantly related to Pneumocystis carinii pneumonia infection (p=0.036), while SLCO1B1 521T>C was associated with anemia (p=0.029).

Conclusion

For Chinese autoimmune disease patients, SLCO1B1 521T>C was correlated with MPA C₀/d and anemia. IMPDH1 −106G>A was significantly related to infections. UGT2B7 802C>T was significantly related to Pneumocystis carinii pneumonia infection.

Pharmacokinetic Model Analysis of Supralingual, Oral and Intravenous Deliveries of Mycophenolic Acid

Mycophenolic acid (MPA) is commonly used for organ rejection prophylaxis via oral administration in the clinic. Recent studies have shown that MPA also has anticancer activities. To explore new therapeutic options for oral precancerous/cancerous lesions, MPA was designed to release topically on the dorsal tongue surface via a mucoadhesive patch. The objective of this study was to establish the pharmacokinetic (PK) and tongue tissue distribution of mucoadhesive MPA patch formulation after supralingual administration in rats and also compare the PK differences between oral, intravenous, and supralingual administration of MPA. Blood samples were collected from Sprague Dawley rats before and after a single intravenous bolus injection, a single oral dose, or a mucoadhesive patch administration on the dorsal tongue surface for 4 h, all with a dose of 0.5 mg/kg of MPA. Plots of MPA plasma concentration versus time were obtained. As multiple peaks were found in all three curves, the enterohepatic recycling (EHR) model in the Phoenix software was adapted to describe their PK parameters with an individual PK analysis method. The mean half-lives of intravenous and oral administrations were 10.5 h and 7.4 h, respectively. The estimated bioavailability after oral and supralingual administration was 72.4% and 7.6%, respectively. There was a 0.5 h lag-time presented after supralingual administration. The results suggest that the systemic plasma MPA concentrations were much lower in rats receiving supralingual administration compared to those receiving doses from the other two routes, and the amount of MPA accumulated in the tongue after patch application showed a sustained drug release pattern. Studies on the dynamic of drug retention in the tongue after supralingual administration showed that ~3.8% of the dose was accumulated inside of tongue right after the patch removal, ~0.11% of the dose remained after 20 h, and ~20.6% of MPA was not released from the patches 4 h after application. The data demonstrate that supralingual application of an MPA patch can deliver a high amount of drug at the site of administration with little systemic circulation exposure, hence lowering the potential gastrointestinal side effects associated with oral administration. Thus, supralingual administration is a potential alternative route for treating oral lesions.

Population Pharmacokinetic Model of Plasma and Cellular Mycophenolic Acid in Kidney Transplant Patients from the CIMTRE Study

Background and Objective

Mycophenolate mofetil is widely used in kidney transplant recipients. Mycophenolate mofetil is hydrolysed by blood esterases to mycophenolic acid (MPA), the active drug. Although MPA therapeutic drug monitoring has been recommended to optimise the treatment efficacy by the area under the plasma concentration vs time curve, little is known regarding MPA concentrations in peripheral blood mononuclear cells, where MPA inhibits inosine monophosphate dehydrogenase. This study aimed to build a pharmacokinetic model using a population approach to describe MPA total and unbound concentrations in plasma and into peripheral blood mononuclear cells in 78 adult kidney transplant recipients receiving mycophenolate mofetil therapy combined with tacrolimus and prednisone.

Methods

Total and unbound plasma concentrations and peripheral blood mononuclear cell concentrations were assayed. A three-compartment model, two for plasma MPA and one for peripheral blood mononuclear cell MPA, with a zero-order absorption and a first-order elimination was used to describe the data.

Results

Mycophenolic acid average concentrations in peripheral blood mononuclear cells were well above half-maximal effective concentration for inosine monophosphate dehydrogenase and no relationship was found with the occurrence of graft rejection. Three covariates affected unbound and intracellular MPA pharmacokinetics: creatinine clearance, which has an effect on unbound MPA clearance, human serum albumin, which influences fraction unbound MPA and the ABCB1 3435 C>T (rs1045642) genetic polymorphism, which has an effect on MPA efflux transport from peripheral blood mononuclear cells.

Conclusion

This population pharmacokinetic model demonstrated the intracellular accumulation of MPA, the efflux of MPA out of the cells being dependent on P-glycoprotein transporters. Nevertheless, further studies are warranted to investigate the relevance of MPA concentrations in peripheral blood mononuclear cells to dosing regimen optimisation.

Electronic supplementary material

The online version of this article (10.1007/s40268-020-00319-y) contains supplementary material, which is available to authorized users.

Population pharmacokinetics of mycophenolic acid in paediatric patients

Mycophenolic acid (MPA) is widely used in paediatric kidney transplant patients and sometimes prescribed for additional indications. Population pharmacokinetic or pharmacodynamic modelling has been frequently used to characterize the fixed, random and covariate effects of MPA in adult patients. However, MPA population pharmacokinetic data in the paediatric population have not been systematically summarized. The objective of this narrative review was to provide an up-to-date critique of currently available paediatric MPA population pharmacokinetic models, with emphases on modelling techniques, pharmacological findings and clinical relevance. PubMed and EMBASE were searched from inception of database to May 2020, where a total of 11 studies have been identified representing kidney transplant (n = 4), liver transplant (n = 1), haematopoietic stem cell transplant (n = 1), idiopathic nephrotic syndrome (n = 2), systemic lupus erythematosus (n = 2), and a combined population consisted of kidney, liver and haematopoietic stem cell transplant patients (n = 1). Critical analyses were provided in the context of MPA absorption, distribution, metabolism, excretion and bioavailability in this paediatric database. Comparisons to adult patients were also provided. With respect to clinical utility, Bayesian estimation models (n = 6) with acceptable accuracy and precision for MPA exposure determination have also been identified and systematically evaluated. Overall, our analyses have identified unique features of MPA clinical pharmacology in the paediatric population, while recognizing several gaps that still warrant further investigations. This review can be used by pharmacologists and clinicians for improving MPA pharmacokinetic-pharmacodynamic modelling and patient care.

Mycophenolic acid area under the concentration-time curve is associated with therapeutic response in childhood-onset lupus nephritis

BACKGROUND

Mycophenolic acid (MPA), the active compound of mycophenolate mofetil (MMF), is widely used in lupus nephritis treatment. Therapeutic drug monitoring of adults suggests that area under the concentration-time curve (AUC) of MPA (MPA-AUC) is associated with clinical outcomes, but childhood data are scarce.

METHODS

Retrospective study of 27 children with biopsy-proven lupus nephritis treated with MMF between 2008 and 2016. In 25 children, MPA-AUC was performed within 6 months after kidney biopsy and MMF initiation. Treatment response at 6 months was defined as normal or improved GFR by 25% compared with baseline, 50% reduction of proteinuria to < 0.5 g/day or 50 mg/mmol, and no hematuria.

RESULTS

A total of 62 MPA-AUC were analyzed in 27 patients. Overall median was 44 mg h/L (interquartile range [IQR] 33-54). Individual dose adaptation was required in 32 cases (52%) to achieve target AUC of 30-60 mg h/L. At 6 months, 14/25 patients were defined as responders (56%, median MPA-AUC 49 mg h/L (40-59)) and 11/25 as non-responders (44%, 29 mg h/L (24-38)). Patients with MPA-AUC levels > 45, 30-45, and < 30 mg h/L had 6-month response rates of 89% (8/9), 60% (6/10), and 0% (0/6), respectively. In a logistic regression model adjusted for age, sex, lupus nephritis classification, and time since MMF initiation, an MPA-AUC > 45 mg h/L was significantly associated with therapeutic response (OR 3.6, 95% CI 2.4-9.5, p = 0.03).

CONCLUSIONS

Therapeutic drug monitoring leading to individualized dosing may improve efficacy of MMF. MPA-AUC > 45 mg h/L is associated with better response rate and may be considered as a target value in pediatric lupus nephritis.

PK/PD Study of Mycophenolate Mofetil in Children With Systemic Lupus Erythematosus to Inform Model-Based Precision Dosing

Objectives: To evaluate the mycophenolic acid [MPA, the active form of mycophenolate mofetil (MMF)] pharmacokinetic parameters in relation to clinical response to identify target exposure ranges in pediatric patients with systemic lupus erythematosus (SLE).

Methods: This was a retrospective study using pharmacokinetic data collected in 67 pediatric patients aged 4–18 years with SLE. Target MPA exposures for effective inhibition of SLE activity (as measured by SLE disease Activity Index (SLEDAI), active SLE was defined as a SLEDAI score of ≥6, and a controlled disease was defined as a SLEDAI score of ≤4) were assessed by receiver operating characteristic (ROC) curve and logistic regression. Exposure-response models were developed to quantitatively describe the relationship between SLEDAI score and AUC_0–12 or C_trough, respectively.

Results: The MPA AUC_0-12 in patients with active SLE was significantly lower than that in patients with inactive SLE. ROC analysis revealed that an AUC_0–12 threshold of 39 μg h/ml or a C_trough of 1.01 μg/ml was associated with the lowest risk of active SLE. Logistic regression analysis revealed that an AUC_0–12 of less than 34 μg h/ml or a C_trough of less than 1.2 μg/ml probably is associated with active SLE. The results of the exposure-response modeling also indicated that an AUC_0-12 less than 32 μg h/ml or a C_trough less than 1.1 μg/ml was associated with suboptimal clinical outcome. An AUC_0-12 above 50 μg h/ml or a C_trough above 1.7 ug/ml was associated with disease control.

Conclusion: Both AUC_0–12 and C_trough of MPA are predictive of the likelihood of active SLE in pediatric patients receiving MMF. An individualized dosing regimen of MMF, with a target AUC_0–12 or C_trough, should be considered for SLE patients.

Prediction of Free from Total Mycophenolic Acid Concentrations in Stable Renal Transplant Patients: A Population-Based Approach

BACKGROUND

A population pharmacokinetic (PK) protein-binding model was developed to (1) predict free mycophenolic acid (fMPA) based on total MPA (tMPA) concentrations in renal transplant patients, to establish the therapeutic range of fMPA through pharmacokinetic-pharmacodynamic studies; and (2) provide a guideline for dosing mycophenolate mofetil (MMF).

METHODS

Full PK profiles of 56 patients (from five different occasions) during the first year after transplantation who were treated with oral MMF and cyclosporine, or macrolides (either tacrolimus or sirolimus), were analysed. fMPA protein-binding was modelled using nonlinear mixed effects modelling (NONMEM). The influence of physiological factors and coadministered immunosupressant was studied.

RESULTS

A two-compartment model with first-order absorption and elimination, linear protein binding and enterohepatic circulation (EHC) best described the PK of MPA. Different recycling rate constants were considered depending on the coadministered immunosuppressant. The protein-binding rate constant (K_B [relative standard error, RSE%]) increased nonlinearly with renal function according to K _B = 43.1 (3.13)·(CL_CR/59.51)^0.394(10.66) h^-1. Furthermore, fMPA plasma clearance, given by clearance of the free mycophenolic acid (CL_fMPA), CL_fMPA = 410 (RSE%3.00)·(1+CsA·0.594 (22.39)) L/h, was 59.4% greater in cyclosporine-treated patients than in macrolide-treated patients, leading to lower MPA exposures. External evaluation proved acceptable area under the plasma concentration-time curve and trough concentration predictions.

CONCLUSIONS

A reliable protein-binding population PK model was developed for prediction of fMPA or tMPA from each other and for dose guiding in stable renal transplant recipients.

Development of a joint population pharmacokinetic model of ezetimibe and its conjugated metabolite

Ezetimibe (EZE) is an extensively used antihyperlipidemic drug with an important cholesterol lowering activity. It undergoes extensive first-pass metabolism to form its active glucuronide metabolite (EZEG). Both drugs exhibit complex pharmacokinetic profiles attributed mainly to repetitive enterohepatic kinetics. The aim of the present study was the investigation of EZE and EZEG pharmacokinetics (PK), through the development of a joint population pharmacokinetic model able to characterize their kinetic processes and enterohepatic recirculation simultaneously. Concentration-time data derived from a bioequivalence study in 28 healthy subjects were used for the analysis. Population PK modeling was performed on the obtained data using nonlinear mixed effect modeling approach, where different methodologies were applied for the description of the complex metabolism and recirculation processes of the two compounds. EZE and EZEG concentrations were best described by a population PK model incorporating first-pass metabolism and an enterohepatic recirculation loop, accounting for the recycling process of the two moieties. This is the first joint population pharmacokinetic model describing the kinetics of both EZE and EZEG.

Pharmacokinetics of mycophenolate mofetil in juvenile patients with autoimmune diseases

Objectives: This study aimed to determine the association between the dosage and pharmacokinetics of mycophenolate mofetil (MMF) in juvenile patients with autoimmune diseases.Methods: Totally, 29 patients were administered oral MMF. The blood concentrations of mycophenolate acid (MPA) at seven points, the area under the time-concentration curve (MPA-AUC_0-12h), the peak concentration (C_max), and the time to peak concentration (T_max) were measured. To obtain a dose-normalized MPA-AUC_0-12h value, the actual measured MPA-AUC_0-12h value was divided by the dose value of the morning administration corrected for body weight (BW) or body surface area (BSA). The patients were classified into three age groups (group 1, ≤10 years; group 2, >10-≤15 years; and group 3, >15 years), and pharmacokinetic parameters were compared among the groups.Results: In total, we obtained 37 measurements. The actual measured MPA-AUC_0-12h values and the MPA-AUC_0-12h values corrected for dose per BW and T_max were lower in young patients. The MPA-AUC_0-12h values corrected for dose per BSA and C_max were comparable among all the groups.Conclusion: In patients with juvenile autoimmune diseases, determining MMF administration dosage according to BSA may facilitate MPA-AUC_0-12h value prediction.

Innovative Study Designs Optimizing Clinical Pharmacology Research in Infants and Children

Almost half of recent pediatric trials failed to achieve labeling indications, in large part because of inadequate study design. Therefore, innovative study methods are crucial to optimizing trial design while also reducing the potential harms inherent with drug investigation. Several methods exist to optimize the amount of pharmacokinetic data collected from the smallest possible volume and with the fewest number of procedures, including the use of opportunistic and sparse sampling, alternative and noninvasive matrices, and microvolume assays. In addition, large research networks using master protocols promote collaboration, reduce regulatory burden, and increase trial efficiency for both early- and late-phase trials. Large pragmatic trials that leverage electronic health records can capitalize on central management strategies to reduce costs, enroll patients with rare diseases on a large scale, and augment study generalizability. Further, trial efficiency and safety can be optimized through Bayesian adaptive techniques that permit planned protocol changes based on analyses of prior and accumulated data. In addition to these trial design features, advances in modeling and simulation have paved the way for systems-based and physiologically based models that individualize pediatric dosing recommendations and support drug approval. Last, given the low prevalence of many pediatric diseases, collecting deidentified genetic and clinical data on a large scale is a potentially transformative way to augment clinical pharmacology research in children.

Population Pharmacokinetic Modeling With Enterohepatic Circulation for AZD3241 in Healthy Subjects and Patients With Multiple System Atrophy

AZD3241 is a potent and selective myeloperoxidase inhibitor potentially for the treatment of a number of neurodegenerative disorders, including multiple system atrophy (MSA). The objectives of this work were to develop a population pharmacokinetic (PopPK) model for AZD3241 and to investigate the correlation between AZD3241 exposure and myeloperoxidase inhibition. The PopPK model was developed using AZD3241 data from one phase 1 study in healthy subjects and one phase 2 study in patients with MSA. A one-compartment model incorporating a gallbladder compartment for enterohepatic circulation, sequential zero-first order absorption, and first-order elimination adequately described the AZD3241 concentration profiles. The apparent clearance and central volume of distribution were 63.1 L/h (interindividual variability: 34.8%) and 121.9 L (interindividual variability: 44.0%), respectively. The enterohepatic circulation model reasonably captured the second peak of AZD3241, and high-fat food increased the absorption rate by 69%. A linear regression model was applied to describe the relationship between AZD3241 exposure and percentage change from baseline in myeloperoxidase-specific activity. The developed PopPK model was consistent with known pharmacokinetic characteristics of AZD3241. This model can be used to estimate AZD3241 exposure in patients with MSA and could be applied to future pharmacokinetic-pharmacodynamic analyses of AZD3241 in clinical development.

Increased Exposure of Tacrolimus by Co-administered Mycophenolate Mofetil: Population Pharmacokinetic Analysis in Healthy Volunteers

The objective of the study was to investigate the pharmacokinetic drug-drug interactions between tacrolimus (TAC) and mycophenolate mofetil (MMF) in healthy Korean male volunteers. Seventeen volunteers participated in a three-period, single-dose, and fixed sequence study. They sequentially received MMF, TAC, and the combination. Concentrations of TAC, mycophenolic acid (MPA), and its metabolites MPA 7-O-glucuronide and MPA acyl glucuronide were measured. The variants of CYP3A4, CYP3A5, SLCO1B1, SLCO1B3, ABCC2, UGT1A9, and UGT2B7 were genotyped. Drug interaction was evaluated with a non-compartmental analysis and population pharmacokinetic modelling to quantify the interaction effect. A total of 1,082 concentrations of those analytes were analysed. AUC_0-inf of TAC increased by 22.1% (322.4 ± 174.1 to 393.6 ± 121.7 ng·h/mL; P < 0.05) when co-administered with MMF, whereas the pharmacokinetic parameters of MPA and its metabolites were not changed by TAC. Apparent clearance (CL/F) of TAC was 17.8 L/h [relative standard error (RSE) 11%] or 13.8 L/h (RSE 11%) without or with MMF, respectively. Interaction was explained by the exponential model. The CYP3A5 genotype was the only significant covariate. The population estimate of CL/F of TAC was 1.48-fold (RSE 16%) in CYP3A5 expressers when compared to nonexpressers. CL/F of TAC was decreased when co-administered with MMF in these subjects.

Pharmacokinetics of mycophenolate mofetil in children with lupus and clinical findings in favour of therapeutic drug monitoring

AIMS

The use of mycophenolate mofetil (MMF) in children with systemic lupus erythematosus (SLE) is increasing. However, the clinical benefit of its monitoring has been scarcely studied, and little is known about its pharmacokinetics in this context. The objectives of the present study were: (i) to describe mycophenolic acid (MPA, the active moiety of MMF) pharmacokinetics, (ii) to develop a Bayesian estimator (BE) allowing the determination AUC (area under the curve) from a limited number of blood samples and (iii) to explore the relationships between exposure indices to MPA and the clinical status in children with SLE.

METHODS

This was a retrospective study including 36 children with SLE, extracted from the expert system ISBA, for whom full- pharmacokinetic profiles of MPA were collected together with clinical data. A pharmacokinetic model and a BE were developed using an iterative two stage Bayesian approach. ROC curve analyses and logistic regressions were used to investigate the association of exposure and active disease.

RESULTS

A pharmacokinetic model and a BE were developed that allowed good AUC estimation performance (bias ± SD = -0.02 ± 0.15). ROC curve analyses showed that AUC/dose <0.06 and AUC <4 mg l(-1)  h were associated with a good sensitivity and specificity for active disease (78%/94% and 94%/56%, respectively). When introduced in a logistic regression model, AUC <44 mg l(-1)  h and AUC/dose <0.06 were associated with an increased risk of active disease (OR = 21.2, 95% CI 2.3, 196.1, P = 0.007 and OR = 59.5, 95% CI 5.9, 588.2, P = 0.0005 respectively].

CONCLUSIONS

The developed pharmacokinetic BE could be used to test prospectively the interest of MPA monitoring for limiting relapse of the disease or its progression.

Population pharmacokinetic-pharmacodynamic modelling of mycophenolic acid in paediatric renal transplant recipients in the early post-transplant period

AIM

The purpose of this study was to develop a population pharmacokinetic and pharmacodynamic (PK-PD) model for mycophenolic acid (MPA) in paediatric renal transplant recipients in the early post-transplant period.

METHODS

A total of 214 MPA plasma concentrations-time data points from 24 patients were available for PK model development. In 17 out of a total of 24 patients, inosine monophosphate dehydrogenase (IMPDH) enzyme activity measurements (n = 97) in peripheral blood mononuclear cells were available for PK-PD modelling. The PK-PD model was developed using non-linear mixed effects modelling sequentially by 1) developing a population PK model and 2) incorporating IMPDH activity into a PK-PD model using post hoc Bayesian PK parameter estimates. Covariate analysis included patient demographics, co-medication and clinical laboratory data. Non-parametric bootstrapping and prediction-corrected visual predictive checks were performed to evaluate the final models.

RESULTS

A two compartment model with a transit compartment absorption best described MPA PK. A non-linear relationship between dose and MPA exposure was observed and was described by a power function in the model. The final population PK parameter estimates (and their 95% confidence intervals) were CL/F, 22 (14.8, 25.2) l h(-1) 70 kg(-1) ; Vc /F, 45.4 (29.6, 55.6) l; Vp /F, 411 (152.6, 1472.6)l; Q/F, 22.4 (16.0, 32.5) l h(-1) ; Ka , 2.5 (1.45, 4.93) h(-1) . Covariate analysis in the PK study identified body weight to be significantly correlated with CL/F. A simplified inhibitory Emax model adequately described the relationship between MPA concentration and IMPDH activity. The final population PK-PD parameter estimates (and their 95% confidence intervals) were: E0 , 3.45 (2.61, 4.56) nmol h(-1)  mg(-1) protein and EC50 , 1.73 (1.16, 3.01) mg l(-1) . Emax was fixed to 0. There were two African-American patients in our study cohorts and both had low IMPDH baseline activities (E0 ) compared with Caucasian patients (mean value 2.13 mg l(-1) vs. 3.86 mg l(-1) ).

CONCLUSION

An integrated population PK-PD model of MPA has been developed in paediatric renal transplant recipients. The current model provides information that will facilitate future studies and may be implemented in a Bayesian algorithm to allow a PK-PD guided therapeutic drug monitoring strategy.

Pharmacokinetic modeling of therapies for systemic lupus erythematosus

With the increasing use of different types of therapies in treating autoimmune diseases such as systemic lupus erythematosus (SLE), there is a need to utilize pharmacokinetic (PK) strategies to optimize the clinical outcome of these treatments. Various PK analysis approaches, including population PK modeling and physiologically based PK modeling, have been used to evaluate drug PK characteristics and population variability or to predict drug PK profiles in a mechanistic manner. This review outlines the PK modeling of major SLE therapies including immunosuppressants (methotrexate, azathioprine, mycophenolate and cyclophosphamide, among others) and immunomodulators (intravenous immunoglobulin). It summarizes the population PK modeling, physiologically based PK modeling and model-based individualized dosing strategies to improve the therapeutic outcomes in SLE patients.

Population Pharmacokinetic Modeling of the Enterohepatic Recirculation of Fimasartan in Rats, Dogs, and Humans

Enterohepatic recirculation (EHC) can greatly enhance plasma drug exposures and therapeutic effects. This study aimed to develop a population pharmacokinetic model that can simultaneously characterize the extent and time-course of EHC in three species using fimasartan, a novel angiotensin II receptor blocker, as a model drug. All fimasartan plasma concentration profiles in 32 rats (intravenous doses, 0.3-3 mg/kg; oral doses, 1-10 mg/kg), 34 dogs (intravenous doses, 0.3-1 mg/kg; oral doses, 1-10 mg/kg), and 42 healthy volunteers (single or multiple oral doses, 20-480 mg) were determined via liquid chromatography-tandem mass spectrometry (LC-MS/MS) and simultaneously modeled in S-ADAPT. The proposed model quantitatively characterized EHC in three species after oral and intravenous dosing. The median (range) fraction of drug undergoing recirculation was 76.3% (64.9-88.7%) in rats, 33.3% (24.0-45.9%) in dogs, and 65.6% (56.5-72.0%) in humans. In the presence compared with the absence of EHC, the area under the curve in plasma was predicted to be 4.22-fold (2.85-8.85) as high in rats, 1.50-fold (1.32-1.85) in dogs, and 2.91-fold (2.30-3.57) in humans. The modeled oral bioavailability in rats (median (range), 38.7% (20.0-59.8%)) and dogs (median, 7.13% to 15.4%, depending on the formulation) matched the non-compartmental estimates well. In humans, the predicted oral bioavailability was 25.1% (15.1-43.9%) under fasting and 18.2% (12.2-31.0%) under fed conditions. The allometrically scaled area under the curve predicted from rats was 420 ng·h/mL for 60 mg fimasartan compared with 424 ± 63 ng·h/mL observed in humans. The developed population pharmacokinetic model can be utilized to characterize the impact of EHC on plasma drug exposure in animals and humans.

Population pharmacokinetics of mycophenolic acid and its main glucuronide metabolite: a comparison between healthy Chinese and Caucasian subjects receiving mycophenolate mofetil

PURPOSE

Mycophenolate mofetil (MMF), a prodrug of the immunosuppressive agent mycophenolic acid (MPA), is widely used for prophylaxis of solid organ transplant rejection. MPA is primarily metabolized to 7-O-mycophenolic acid glucuronide (MPAG), an inactive metabolite that undergoes enterohepatic recirculation (EHC). This study assessed ethnic differences in the pharmacokinetics (PK) of MPA and MPAG between healthy Chinese and Caucasian subjects using population PK analysis.

METHODS

Data were pooled from 132 healthy subjects (80 Chinese, 52 Caucasians) in eight clinical studies in which MMF was administered in a single oral dose. Population PK analysis was performed using NONMEM®.

RESULTS

The PK of MPA and MPAG were best described by a five-chain compartment model, including a gallbladder compartment for EHC and a transit absorption model. Ethnicity was significantly correlated with the apparent clearance (CL/F) and volume of distribution (V/F) of MPAG but not those of MPA. Weight was identified as a covariate and was correlated with the PK of MPA and MPAG. MPA CL/F was 11.5 L/h for a 70-kg healthy subject, and the MPAG CL/F values were 1.36 and 1.90 L/h for 70-kg Chinese and Caucasian individuals, respectively. Internal and external evaluation indicated model validity.

CONCLUSIONS

This is the first population PK analysis to evaluate ethnic differences in the PK of MPA and MPAG in healthy Chinese and Caucasian subjects. No differences were observed in the PK of MPA between healthy Chinese and Caucasian subjects. Although, the MPAG CL/F was approximately 40 % higher in Caucasians, this finding may not be clinically relevant.