Pharmacokinetic modelling and development of Bayesian estimators for therapeutic drug monitoring of mycophenolate mofetil in reduced-intensity haematopoietic stem cell transplantation

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.

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.

Impact of a Bayesian Individualization of Cyclosporine Dosage Regimen for Children Undergoing Allogeneic Hematopoietic Cell Transplantation: A Cost-Effectiveness Analysis

BACKGROUND

Cyclosporine A (CsA) is the main drug used to prevent graft-versus-host disease in patients undergoing allogeneic hematopoietic stem cell transplantation (HSCT). CsA therapeutic drug monitoring (TDM) has been performed for ages, with studies revealing clinical benefits, but failing to examine its economic impact. In this article, the main objective was to evaluate the economic impact of the CsA TDM strategy, based on a Bayesian approach, by assessing costs related to its clinical impact. Furthermore, TDM effectiveness was analyzed for pharmacokinetics and clinical outcomes.

METHODS

A cost-effective, nonrandomized, retrospective, single-center study compared 2 CsA monitoring and dose adaptation strategies in pediatric patients undergoing HSCT. From 2014 to 2016, CsA TDM was performed using a population pharmacokinetics model-coupled Bayesian approach by a pharmacist ["pharmacist-assisted individualization" (PAI)]. From 2017 to 2018, CsA TDM was performed by the clinician without a Bayesian approach (non-PAI group). HSCT costs were evaluated from the French National Insurance perspective. Economic and clinical outcomes were assessed by measuring incremental cost-effectiveness ratios.

RESULTS

The study included 144 patients: 90 and 54 patients in PAI and non-PAI groups, respectively. Both groups were comparable for sociodemographic and clinical characteristics. The mean total cost per patient was significantly lower (P < 0.01) in the PAI group (€85,947) than in the non-PAI group (€100,435). Multivariate analysis revealed that TDM based on the Bayesian approach was a protective factor (odds ratio = 0.86) for severe acute graft-versus-host disease. We noted that pharmacist-based TDM was the dominant strategy. Bayesian method-based TDM allowed an increase in the percentage of target attainment at any period post-HSCT.

CONCLUSIONS

CsA TDM with a Bayesian approach is a cost-effective procedure, and highlighted clinical benefits encourage the development of new TDM strategies for HSCT.

Personalized Therapy for Mycophenolate: Consensus Report by the International Association of Therapeutic Drug Monitoring and Clinical Toxicology

ABSTRACT

When mycophenolic acid (MPA) was originally marketed for immunosuppressive therapy, fixed doses were recommended by the manufacturer. Awareness of the potential for a more personalized dosing has led to development of methods to estimate MPA area under the curve based on the measurement of drug concentrations in only a few samples. This approach is feasible in the clinical routine and has proven successful in terms of correlation with outcome. However, the search for superior correlates has continued, and numerous studies in search of biomarkers that could better predict the perfect dosage for the individual patient have been published. As it was considered timely for an updated and comprehensive presentation of consensus on the status for personalized treatment with MPA, this report was prepared following an initiative from members of the International Association of Therapeutic Drug Monitoring and Clinical Toxicology (IATDMCT). Topics included are the criteria for analytics, methods to estimate exposure including pharmacometrics, the potential influence of pharmacogenetics, development of biomarkers, and the practical aspects of implementation of target concentration intervention. For selected topics with sufficient evidence, such as the application of limited sampling strategies for MPA area under the curve, graded recommendations on target ranges are presented. To provide a comprehensive review, this report also includes updates on the status of potential biomarkers including those which may be promising but with a low level of evidence. In view of the fact that there are very few new immunosuppressive drugs under development for the transplant field, it is likely that MPA will continue to be prescribed on a large scale in the upcoming years. Discontinuation of therapy due to adverse effects is relatively common, increasing the risk for late rejections, which may contribute to graft loss. Therefore, the continued search for innovative methods to better personalize MPA dosage is warranted.

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.

Population pharmacokinetics and Bayesian estimators for intravenous mycophenolate mofetil in haematopoietic stem cell transplant patients

AIMS

Intravenous mycophenolate mofetil (IV MMF), a prodrug of mycophenolic acid (MPA), is used during nonmyeloablative and reduced-intensity conditioning haematopoetic stem cell transplantation (HCT) to improve engraftment and reduce graft-versus-host disease. The aims of this study were to develop population pharmacokinetic models and Bayesian estimators based on limited sampling strategies to allow for individual dose adjustment of intravenous mycophenolate mofetil administered by infusion in haematopoietic stem cell transplant patients.

METHODS

Sixty-three MPA concentration-time profiles (median [min-max] = 6 [4-7] samples) were collected from 34 HCT recipients transplanted for 14 (1-45) days and administered IV MMF every 8 hours, concomitantly with cyclosporine. The database was split into development (75%) and validation (25%) datasets. Pharmacokinetic models characterized by a single compartment with first-order elimination, combined with two gamma distributions to describe the transformation of MMF into mycophenolic acid, were developed using in parallel nonparametric (Pmetrics) and parametric (ITSIM) approaches. The performances of the models and the derived Bayesian estimators were evaluated in the validation set.

RESULTS

The best limited sampling strategy led to a bias (min, max), root mean square error between observed and modeled interdose areas under the curve in the validation dataset of -11.72% (-31.08%, 5.00%), 14.9% for ITSIM and -2.21% (-23.40%, 30.01%), 12.4% for Pmetrics with three samples collected at 0.33, 2 and 3 hours post dosing.

CONCLUSION

Population pharmacokinetic models and Bayesian estimators for IV MMF in HCT have been developed and are now available online (https://pharmaco.chu-limoges.fr) for individual dose adjustment based on the interdose area under the curve.

Antiviral Properties of Chemical Inhibitors of Cellular Anti-Apoptotic Bcl-2 Proteins

Viral diseases remain serious threats to public health because of the shortage of effective means of control. To combat the surge of viral diseases, new treatments are urgently needed. Here we show that small-molecules, which inhibit cellular anti-apoptotic Bcl-2 proteins (Bcl-2i), induced the premature death of cells infected with different RNA or DNA viruses, whereas, at the same concentrations, no toxicity was observed in mock-infected cells. Moreover, these compounds limited viral replication and spread. Surprisingly, Bcl-2i also induced the premature apoptosis of cells transfected with viral RNA or plasmid DNA but not of mock-transfected cells. These results suggest that Bcl-2i sensitizes cells containing foreign RNA or DNA to apoptosis. A comparison of the toxicity, antiviral activity, and side effects of six Bcl-2i allowed us to select A-1155463 as an antiviral lead candidate. Thus, our results pave the way for the further development of Bcl-2i for the prevention and treatment of viral diseases.

Pharmacokinetics, Pharmacodynamics, and Pharmacogenomics of Immunosuppressants in Allogeneic Hematopoietic Cell Transplantation: Part II

Part I of this article included a pertinent review of allogeneic hematopoietic cell transplantation (alloHCT), the role of postgraft immunosuppression in alloHCT, and the pharmacokinetics, pharmacodynamics, and pharmacogenomics of the calcineurin inhibitors and methotrexate. In this article (Part II), we review the pharmacokinetics, pharmacodynamics, and pharmacogenomics of mycophenolic acid (MPA), sirolimus, and the antithymocyte globulins (ATG). We then discuss target concentration intervention (TCI) of these postgraft immunosuppressants in alloHCT patients, with a focus on current evidence for TCI and on how TCI may improve clinical management in these patients. Currently, TCI using trough concentrations is conducted for sirolimus in alloHCT patients. Several studies demonstrate that MPA plasma exposure is associated with clinical outcomes, with an increasing number of alloHCT patients needing TCI of MPA. Compared with MPA, there are fewer pharmacokinetic/dynamic studies of rabbit ATG and horse ATG in alloHCT patients. Future pharmacokinetic/dynamic research of postgraft immunosuppressants should include '-omics'-based tools: pharmacogenomics may be used to gain an improved understanding of the covariates influencing pharmacokinetics as well as proteomics and metabolomics as novel methods to elucidate pharmacodynamic responses.

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.

Mycophenolic mofetil optimized pharmacokinetic modelling, and exposure-effect associations in adult heart transplant recipients

Mycophenolic acid (MPA) area under the curve (AUC) has been associated with graft outcome.

THE AIMS OF OUR STUDY WERE

(1) to develop pharmacokinetic tools to optimize MPA inter-dose AUC estimation in heart transplant patients; and (2) to investigate the relationships between acute allograft rejection and MPA AUC, trough level (C0) or mycophenolate mofetil (MMF) dose. Two independent modeling approaches (parametric and non parametric) were used to fit 56 rich MPA pharmacokinetic (PK) profiles collected from 40 adult heart transplant recipients enrolled in the PIGREC study, receiving MMF and a calcineurin inhibitor (CNI), in the first year post-transplantation. In addition, associations between drug exposure (MPA C0, AUC and MMF dose) and acute rejection or MMF adverse events were investigated using time-dependent Cox models with stratification on the type of calcineurin inhibitor. Exposure threshold values were investigated using ROC curve analysis. The 2 models developed fit adequately the data and the use of their combination yielded 100% consistency with the measured AUC in terms of strategy of dose adjustment (maintain, increase or decrease). MPA measured AUC adjusted on CNI exposure was significantly associated with rejection (per unit increase: HR [95% CI]=0.97 [0.95-0.99], p=0.0122), while no effect was shown for adverse events attributable to MMF. An AUC threshold of 50 mg×h/L was proposed (sensitivity=77%, specificity=25%) beyond which the risk of rejection was significantly increased (low vs. high: HR=3.48 [1.21-10.0], p=0.0204). The tools developed have already been made available to the heart transplant community on our ISBA website (https://pharmaco.chu-limoges.fr).

Population pharmacokinetics of nefopam in elderly, with or without renal impairment, and its link to treatment response

AIMS

Nefopam is a nonmorphinic central analgesic, for which no recommendation exists concerning adaptation of regimen in aged patients with or without renal impairment. The objective was to describe the pharmacology of nefopam in aged patients to obtain guidelines for practical use.

METHODS

Elderly patients (n = 48), 65-99 years old, with severe or moderate renal impairment or with normal renal function, were recruited. Nefopam (20 mg) was administered as a 30 min infusion postoperatively. Simultaneously, a 1 min intravenous infusion of iohexol was performed, in order to calculate the glomerular filtration rate. Blood samples were drawn to determine nefopam, desmethyl-nefopam and iohexol plasma concentrations. Nefopam and desmethyl-nefopam concentrations were analysed using a nonlinear mixed-effects modelling approach with Monolix version 4.1.3. The association between pharmacokinetic parameters and treatment response was assessed using logistic regression.

RESULTS

A two-compartment open model was selected to describe the pharmacokinetics of nefopam. The typical population estimates (between-subject variability) for clearance, volume of distribution, intercompartmental clearance and peripheral volume were, respectively, 17.3 l h(-1) (53.2%), 114 l (121%), 80.7 l h(-1) (79%) and 208 l (63.6%). Morphine requirement was related to exposure of nefopam. Tachycardia and postoperative nausea and vomiting were best associated with maximal concentration and the rate of increase in nefopam plasma concentration.

CONCLUSIONS

We identified the nefopam pharmacokinetic predictors for morphine requirement and side-effects, such as tachycardia and postoperative nausea and vomiting. In order to maintain morphine sparing and decrease side-effects following a single dose of nefopam (20 mg), simulations suggest an infusion time of >45 min in elderly patients with or without renal impairment.

Clinical mycophenolic acid monitoring in liver transplant recipients

In liver transplantation, the efficacy of mycophenolate mofetil (MMF) has been confirmed in clinical trials and studies. However, therapeutic drug monitoring for mycophenolic acid (MPA) has not been fully accepted in liver transplantation as no long-term prospective study of concentration controlled vs fixed-dose prescribing of MMF has been done. This review addressed MPA measurement, pharmacokinetic variability and reasons of this variation, exposure related to acute rejection and MMF-associated side effects in liver transplant recipients. Limited sampling strategies to predict MPA area under the concentration-time curve have also been described, and the value of clinical use needs to be investigated in future. The published data suggested that a fixed-dosage MMF regimen might not be suitable and monitoring of MPA exposure seems helpful in various clinical settings of liver transplantation.

Nonrelapse mortality and mycophenolic acid exposure in nonmyeloablative hematopoietic cell transplantation

We evaluated the pharmacodynamic relationships between mycophenolic acid (MPA), the active metabolite of mycophenolate mofetil (MMF), and outcomes in 308 patients after nonmyeloablative hematopoietic cell transplantation. Patients were conditioned with total body irradiation ± fludarabine, received grafts from HLA-matched related (n = 132) or unrelated (n = 176) donors, and received postgrafting immunosuppression with MMF and a calcineurin inhibitor. Total and unbound MPA pharmacokinetics were determined to day 25; maximum a posteriori Bayesian estimators were used to estimate total MPA concentration at steady state (Css). Rejection occurred in 9 patients, 8 of whom had a total MPA Css less than 3 μg/mL. In patients receiving a related donor graft, MPA Css was not associated with clinical outcomes. In patients receiving an unrelated donor graft, low total MPA Css was associated with increased grades III to IV acute graft-versus-host disease and increased nonrelapse mortality but not with day 28 T cell chimerism, disease relapse, cytomegalovirus reactivation, or overall survival. We conclude that higher initial oral MMF doses and subsequent targeting of total MPA Css to greater than 2.96 μg/mL could lower grades III to IV acute graft-versus-host disease and nonrelapse mortality in patients receiving an unrelated donor graft.

Population pharmacokinetics and dose optimization of mycophenolic acid in HCT recipients receiving oral mycophenolate mofetil

We sought to create a population pharmacokinetic model for total mycophenolic acid (MPA), to study the effects of different covariates on MPA pharmacokinetics, to create a limited sampling schedule (LSS) to characterize MPA exposure (i.e., area under the curve or AUC) with maximum a posteriori Bayesian estimation, and to simulate an optimized dosing scheme for allogeneic hematopoietic cell transplantation (HCT) recipients. Four thousand four hundred ninety-six MPA concentration-time points from 408 HCT recipients were analyzed retrospectively using a nonlinear mixed effects modeling approach. MPA pharmacokinetics was characterized with a two-compartment model with first-order elimination and a time-lagged first-order absorption process. Concomitant cyclosporine and serum albumin were significant covariates. The median MPA clearance (CL) and volume of the central compartment were 24.2 L/hour and 36.4 L, respectively, for a 70 kg patient receiving tacrolimus with a serum albumin of 3.4 g/dL. Dosing simulations indicated that higher oral MMF doses are needed with concomitant cyclosporine, which increases MPA CL by 33.8%. The optimal LSS was immediately before and at 0.25 hours, 1.25 hours, 2 hours, and 4 hours after oral mycophenolate mofetil administration. MPA AUC in an individual HCT recipient can be accurately estimated using a five-sample LSS and maximum a posteriori Bayesian estimation.

Fundamentals of population pharmacokinetic modelling: validation methods

Population pharmacokinetic modelling is widely used within the field of clinical pharmacology as it helps to define the sources and correlates of pharmacokinetic variability in target patient populations and their impact upon drug disposition; and population pharmacokinetic modelling provides an estimation of drug pharmacokinetic parameters. This method's defined outcome aims to understand how participants in population pharmacokinetic studies are representative of the population as opposed to the healthy volunteers or highly selected patients in traditional pharmacokinetic studies. This review focuses on the fundamentals of population pharmacokinetic modelling and how the results are evaluated and validated. This review defines the common aspects of population pharmacokinetic modelling through a discussion of the literature describing the techniques and placing them in the appropriate context. The concept of validation, as applied to population pharmacokinetic models, is explored focusing on the lack of consensus regarding both terminology and the concept of validation itself. Population pharmacokinetic modelling is a powerful approach where pharmacokinetic variability can be identified in a target patient population receiving a pharmacological agent. Given the lack of consensus on the best approaches in model building and validation, sound fundamentals are required to ensure the selected methodology is suitable for the particular data type and/or patient population. There is a need to further standardize and establish the best approaches in modelling so that any model created can be systematically evaluated and the results relied upon.

A limited sampling schedule to estimate mycophenolic Acid area under the concentration-time curve in hematopoietic cell transplantation recipients

Mycophenolate mofetil (MMF) is a key component of postgrafting immunosuppression in hematopoietic cell transplant (HCT) recipients. The plasma area under the curve (AUC) of its active metabolite, mycophenolic acid (MPA), is associated with MMF efficacy and toxicity. This study developed a population pharmacokinetic model of MPA in HCT recipients and created limited sampling schedules (LSSs) to enable individualized pharmacotherapy. A retrospective evaluation of MPA concentration-time data following a 2-hour MMF intravenous (IV) infusion was conducted in 77 HCT recipients. The final model consisted of 1 and 2 compartments for MMF and MPA pharmacokinetics, respectively. The mean estimated values (coefficient of variation, %) for total systemic clearance, distributional clearance, and central and peripheral compartment volumes of MPA were 36.9 L/h (34.5%), 15.3 L/h (80.4%), 11.9 L (71.7%), and 182 L (127%), respectively. No covariates significantly explained variability among individuals. Optimal LSSs were derived using a simulation approach based on the scaled mean squared error. A 5-sample schedule of 2, 2.5, 3, 5, and 6 hours from the start of the infusion precisely estimated MPA AUC(0-12 h) for Q12-hour IV MMF. A comparable schedule (2, 2.5, 3, 4, and 6 hours) similarly estimated MPA AUC(0-8) (h) for Q8-hour dosing.

Maximum a posteriori Bayesian estimation of mycophenolic Acid area under the concentration-time curve: is this clinically useful for dosage prediction yet?

This review seeks to summarize the available data about Bayesian estimation of area under the plasma concentration-time curve (AUC) and dosage prediction for mycophenolic acid (MPA) and evaluate whether sufficient evidence is available for routine use of Bayesian dosage prediction in clinical practice. A literature search identified 14 studies that assessed the predictive performance of maximum a posteriori Bayesian estimation of MPA AUC and one report that retrospectively evaluated how closely dosage recommendations based on Bayesian forecasting achieved targeted MPA exposure. Studies to date have mostly been undertaken in renal transplant recipients, with limited investigation in patients treated with MPA for autoimmune disease or haematopoietic stem cell transplantation. All of these studies have involved use of the mycophenolate mofetil (MMF) formulation of MPA, rather than the enteric-coated mycophenolate sodium (EC-MPS) formulation. Bias associated with estimation of MPA AUC using Bayesian forecasting was generally less than 10%. However some difficulties with imprecision was evident, with values ranging from 4% to 34% (based on estimation involving two or more concentration measurements). Evaluation of whether MPA dosing decisions based on Bayesian forecasting (by the free website service https://pharmaco.chu-limoges.fr) achieved target drug exposure has only been undertaken once. When MMF dosage recommendations were applied by clinicians, a higher proportion (72-80%) of subsequent estimated MPA AUC values were within the 30-60 mg · h/L target range, compared with when dosage recommendations were not followed (only 39-57% within target range). Such findings provide evidence that Bayesian dosage prediction is clinically useful for achieving target MPA AUC. This study, however, was retrospective and focussed only on adult renal transplant recipients. Furthermore, in this study, Bayesian-generated AUC estimations and dosage predictions were not compared with a later full measured AUC but rather with a further AUC estimate based on a second Bayesian analysis. This study also provided some evidence that a useful monitoring schedule for MPA AUC following adult renal transplant would be every 2 weeks during the first month post-transplant, every 1-3 months between months 1 and 12, and each year thereafter. It will be interesting to see further validations in different patient groups using the free website service. In summary, the predictive performance of Bayesian estimation of MPA, comparing estimated with measured AUC values, has been reported in several studies. However, the next step of predicting dosages based on these Bayesian-estimated AUCs, and prospectively determining how closely these predicted dosages give drug exposure matching targeted AUCs, remains largely unaddressed. Further prospective studies are required, particularly in non-renal transplant patients and with the EC-MPS formulation. Other important questions remain to be answered, such as: do Bayesian forecasting methods devised to date use the best population pharmacokinetic models or most accurate algorithms; are the methods simple to use for routine clinical practice; do the algorithms actually improve dosage estimations beyond empirical recommendations in all groups that receive MPA therapy; and, importantly, do the dosage predictions, when followed, improve patient health outcomes?

Optimal sampling strategy development methodology using maximum a posteriori Bayesian estimation

Maximum a posteriori Bayesian (MAPB) pharmacokinetic parameter estimation is an accurate and flexible method of estimating individual pharmacokinetic parameters using individual blood concentrations and prior information. In the past decade, many studies have developed optimal sampling strategies to estimate pharmacokinetic parameters as accurately as possible using either multiple regression analysis or MAPB estimation. This has been done for many drugs, especially immunosuppressants and anticancer agents. Methods of development for optimal sampling strategies (OSS) are diverse and heterogeneous. This review provides a comprehensive overview of OSS development methodology using MAPB pharmacokinetic parameter estimation, determines the transferability of published OSSs, and compares sampling strategies determined by MAPB estimation and multiple regression analysis. OSS development has the following components: 1) prior distributions; 2) reference value determination; 3) optimal sampling time identification; and 4) validation of the OSS. Published OSSs often lack all data necessary for the OSS to be clinically transferable. MAPB estimation is similar to multiple regression analysis in terms of predictive performance but superior in flexibility.

The evolution of population pharmacokinetic models to describe the enterohepatic recycling of mycophenolic acid in solid organ transplantation and autoimmune disease

With the increasing use of mycophenolic acid (MPA) as an immunosuppressant in solid organ transplantation and in treating autoimmune diseases such as systemic lupus erythematosus, the need for strategies to optimize therapy with this agent has become increasingly apparent. This need is largely based on MPA's significant between-subject and between-occasion (within-subject) pharmacokinetic variability. While there is a strong relationship between MPA exposure and effect, the relationship between drug dose, plasma concentration and exposure (area under the concentration-time curve [AUC]) is very complex and remains to be completely defined. Population pharmacokinetic models using various approaches have been proposed over the past 10 years to further evaluate the pharmacokinetic and pharmacodynamic behaviour of MPA. These models have evolved from simple one-compartment linear iterations to complex multi-compartment versions that try to include various factors, which may influence MPA's pharmacokinetic variability, such as enterohepatic recycling and pharmacogenetic polymorphisms. There have been major advances in the understanding of the roles transport mechanisms, metabolizing and other enzymes, drug-drug interactions and pharmacogenetic polymorphisms play in MPA's pharmacokinetic variability. Given these advances, the usefulness of empirical-based models and the limitations of nonlinear mixed-effects modelling in developing mechanism-based models need to be considered and discussed. If the goal is to individualize MPA dosing, it needs to be determined whether factors which may contribute significantly to variability can be utilized in the population pharmacokinetic models. Some pharmacokinetic models developed to date show promise in being able to describe the impact of physiological processes such as enterohepatic recycling. Most studies have historically been based on retrospective data or poorly designed studies which do not take these factors into consideration. Modelling typically has been undertaken using non-controlled therapeutic drug monitoring data, which do not have the information content to support the development of complex mechanistic models. Only a few recent modelling approaches have moved away from empiricism and have included mechanisms considered important, such as enterohepatic recycling. It is recognized that well thought-out sampling schedules allow for better evaluation of the pharmacokinetic data. It is not possible to undertake complex absorption modelling with very few samples being obtained during the absorption phase (which has often been the case). It is important to utilize robust AUC monitoring which is now being propagated in the latest consensus guideline on MPA therapeutic drug monitoring. This review aims to explore the biological factors that contribute to the clinical pharmacokinetics of MPA and how these have been introduced in the development of population pharmacokinetic models. An overview of the processes involved in the enterohepatic recycling of MPA will be provided. This will summarize the components that complicate absorption and recycling to influence MPA exposure such as biotransformation, transport, bile physiology and gut flora. Already published population pharmacokinetic models will be examined, and the evolution of these models away from empirical approaches to more mechanism-based models will be discussed.

Development of a Bayesian estimator for the therapeutic drug monitoring of mycophenolate mofetil in children with idiopathic nephrotic syndrome

The use of mycophenolate mofetil (MMF) in children with idiopathic nephrotic syndrome (INS) 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 study and model the pharmacokinetics of mycophenolic acid (MPA; the active moiety of MMF) in paediatric patients with INS given MMF, at all stages of the disease; (ii) to develop a Bayesian estimator (MAP-BE) for individual inter-dose area under the concentration-time curve (AUC) prediction in this population, using a limited blood sampling strategy (LSS). Full-pharmacokinetic (PK) profiles of MPA collected in paediatric inpatients with INS already treated with a maintenance immunosuppressive therapy based on MMF (with no calcineurin inhibitors; CNI) were studied. A classical iterative two-stage (ITS) method was applied to model the data and develop MAP-BEs using a one-compartment open model where the absorption is described by a double gamma law allowing the description of a potential enterohepatic recirculation. The performance of the MAP-BE developed for individual exposure assessment was evaluated by the bias and precision of predicted AUCs with respect to measured, trapezoidal AUCs (reference value), and by the proportion of predicted AUCs with absolute error >20%. These PK tools were tested in an independent group of patients. Sixty PK profiles of MPA from children receiving MMF in association to corticosteroids or given alone were included in the study. Forty-five of these PK profiles were used to develop a PK model and a MAP-BE, and 15 for their validation. In the building group, the PK model fitted accurately the PK profiles of MPA: mean residual error of modelled vs. reference AUC was m±SD=-0.015±0.092 (range: -0.153 to 0.204). The MAP-BE which allowed the estimation of MPA AUC on the basis of a 20 min-60 min-180 min LSS was then developed. In the independent group of patients, its mean residual error vs. reference AUCs was m±SD=-0.036±0.145 (range: -0.205 to 0.189). Thus, a PK model and its derived MAP-BE for MMF (without any associated CNI) when given to children with INS have been developed. Clinical trials using these PK tools could test the potential impact of the therapeutic drug monitoring of MMF based on the AUC on the clinical evolution of INS.