Enterohepatic circulation model for population pharmacokinetic analysis

Effect of Meal Timings and Meal Content on the AUC0-12h of Mycophenolic Acid: A Simulation Study

Meal timings and content related to gallbladder emptying in the enterohepatic circulation are important for explaining the high variability in mycophenolic acid exposure. The limited sampling strategy (LSS) was established to estimate the area under the plasma concentration-time curve from time 0 to 12 hours (AUC_0-12h ) of mycophenolic acid in therapeutic drug monitoring. The aim of this study is to investigate the effect of meal timings and content on the AUC_0-12h of mycophenolic acid and to assess the influence of meals on LSS. A mycophenolic acid pharmacokinetic model with a mechanism-based enterohepatic circulation process was employed to perform simulations under various assumed meal scenarios. The simulations were compared to evaluate the effect of meal timings and meal content on mycophenolic acid AUC_0-12h . Monte Carlo simulations were performed using the meal parameter with the greatest impact on mycophenolic acid AUC_0-12h as a variable. The corresponding LSS equations were established, and the predictive performance was assessed. Both the meal timings and meal content affected the mycophenolic acid AUC_0-12h , and the postdose fasting period had the greatest impact. The predictive performance of the LSS is sensitive to the postdose fasting period. Therefore, meal timings may improve the estimation of mycophenolic acid AUC_0-12h and the efficacy of therapeutic drug monitoring.

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.

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.

Pharmacokinetics, bioavailability and withdrawal period of antibiotic oxytetracycline in catfish Pangasianodon hypophthalmus

The antibiotic oxytetracycline (OTC) has been widely used for therapeutic and preventive management of bacterial diseases in finfish and shellfish. In the present study the bioavailability, pharmacokinetics, and withdrawal period of the OTC have been determined following in-feed administration in intensively cultured catfish Pangasianodon hypophthalmus. Furthermore, the pharmacokinetic parameters of oral route were also compared with parenteral route. Drug concentrations were measured in various tissues at different time intervals by LC-MS/MS. The study revealed the drug kinetics best followed the enterohepatic circulation model with very poor bioavailability and low blood concentration after oral administration. In the withdrawal study, after 10-days of in-feed administration at the therapeutic dose the drug reached very high concentrations in the liver and kidneys but did not attain minimum inhibitory concentrations (MICs) in blood or flesh. OTC concentration also did not exceed the recommended MRL value in flesh; however, considering high amounts of the chemical in the liver and kidneys a withdrawal period of 4 days (at 28 ± 1.5 ^oC) is recommended for consumer safety. Poor bioavailability and non-attainment of minimum therapeutic concentration in blood and flesh do not warrant in-feed administration of OTC for control of bacterial diseases in P. hypophthalmus. AVAILABILITY OF DATA AND MATERIALS: All data generated and analyzed during this study are included in this article. Raw data may be shared upon reasonable request.

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.

Effect of Protein Binding on Exposure of Unbound and Total Mycophenolic Acid: A Population Pharmacokinetic Analysis in Chinese Adult Kidney Transplant Recipients

OBJECTIVES

The population pharmacokinetic (popPK) characteristics of total mycophenolic acid (tMPA) have been investigated in various ethnic populations. However, investigations of popPK of unbound MPA (uMPA) are few. Thus, a popPK analysis was performed to: (1) characterize the PK of uMPA and tMPA and its 7-O-mycophenolic acid glucuronide (MPAG) metabolite in kidney transplant patients cotreated with cyclosporine (CsA), and (2) identify the clinically significant covariates that explain variability in the dose-exposure relationship.

METHODS

A total of 740 uMPA, 741 tMPA, and 734 total MPAG (tMPAG) concentration-time data from 58 Chinese kidney transplant patients receiving MPA in combination with CsA were analyzed using NONMEM® software with the stochastic approximation expectation maximization (SAEM) followed by the important sampling (IMP) method. The influence of covariates was tested using a stepwise procedure.

RESULTS

The PK of uMPA and unbound MPAG (uMPAG) were characterized by a two- and one-compartment model with first-order elimination, respectively. A linear protein binding model was used to link uMPA and tMPA. Apparent clearance (CL/F) and central volume of distribution (VC/F) of uMPA (CLuMPA/F and VCuMPA/F, respectively) and protein binding rate constant (k B) were estimated to be 851 L/h [relative standard error (RSE), 7.1%], 718 L (18.5%) and 53.4/h (2.3%), respectively. For uMPAG, the population values (RSE) of CL/F (CLuMPAG) and VC/F (VCuMPAG/F) were 5.71 L/h (4.4%) and 29.9 L (7.7%), respectively. Between-subject variability (BSVs) on CLuMPA/F, VCuMPA/F, CLuMPAG/F, and VCuMPAG/F were 51.0, 80.0, 31.8 and 48.4%, respectively, whereas residual unexplained variability (RUVs) for uMPA, tMPA, and uMPAG were 47.0, 45.9, and 22.0%, respectively. Significant relationships were found between k B and serum albumin (ALB) and between CLuMPAG/F and glomerular filtration rate (GFR). Additionally, model-based simulation showed that changes in ALB concentrations substantially affected tMPA but not uMPA exposure.

CONCLUSIONS

The established model adequately described the popPK characteristics of the uMPA, tMPA, and MPAG. The estimated CLuMPA/F and unbound fraction of MPA (FUMPA) in Chinese kidney transplant recipients cotreated with CsA were comparable to those published previously in Caucasians. We recommend monitoring uMPA instead of tMPA to optimize mycophenolate mofetil (MMF) dosing for patients with lower ALB levels.

A compartmental approach to isosteviol's disposition in Sprague-Dawley rats

Isosteviol has been reported to reverse hypertrophy and related inflammatory responses in in vitro models representative of cardiac muscle cells. The disposition of isosteviol is, however, characterized by secondary peaks and long plasma residence time despite reports of a relatively short half-life in liver fractions. The present study describes a compartmental approach to modelling the secondary peaks characteristic of isosteviol's concentration-time data in Sprague-Dawley rats. Oral (4 mg/kg) and intravenous (4 mg/kg) doses of isosteviol were administered to male and female Sprague-Dawley rats. Plasma samples collected between 0 and 72 h, and total bile secreted in 24 h, were analysed for isosteviol content with LC-MS/MS techniques. The disposition of isosteviol was, thereafter, described with a structural model that accounted for the sampling, liver and biliary secretion compartments, with a gap-time characterizing the accumulation and subsequent emptying of isosteviol for re-absorption. The half-life of isosteviol following oral dosing was about 103% greater in female rats than in the male, and the model-derived area under the concentration-time curve (AUC) in 72 h was about 756% greater in female animals than in males. Following the administration of intravenous doses of isosteviol, half-life and AUC in 24 h were about 332% and 595%, respectively, higher in female rats than in males. Isosteviol equivalent secreted into bile over 24 h accounted for about 94% of orally administered dose in male rats, and about 59% of oral dose in females. These findings show a differential systemic removal of isosteviol in Sprague-Dawley rats, likely explainable by gender-related differences in the glucuronidation-capacity of isosteviol.

Optimization and application of an HPLC method for quantification of inosine-5'-monophosphate dehydrogenase activity as a pharmacodynamic biomarker of mycophenolic acid in Chinese renal transplant patients

BACKGROUND

The immunosuppressive agent mycophenolic acid (MPA) is a non-competitive, reversible inhibitor of inosine-5'-monophosphate dehydrogenase (IMPDH). Thus, IMPDH activity can serve as a potential pharmacodynamic biomarker to optimize dosing of MPA.

METHODS

Peripheral blood mononuclear cells were isolated from 2 mL blood samples and an in vitro enzymatic reaction was subsequently performed for 120 min. To determine IMPDH activity in Chinese healthy volunteers and renal transplant patients, a high performance liquid chromatography assay was established and validated by subtracting adenosine monophosphate (AMP) from blank samples for eliminating exogenous AMP interference.

RESULTS

The accuracy of our method ranged between -0.8% and 12.5%, and the precision ranged between 0.7% and 6.3%. The mean value of IMPDH activity across 11 healthy volunteers was 46.60 ± 14.28 μmol/s/mol AMP. A negative relationship between MPA concentration and IMPDH activity was observed in four renal transplant patients treated with MPA 13 days post-transplantation, while the inhibitory rate of IMPDH activity ranged from 24% to 42%.

CONCLUSION

A bioanalytical assay for IMPDH quantification was optimized and evaluated. The differences in the pharmacodynamics of MPA between Asians and Caucasians may provide some evidence for dosing differences among ethnicities.

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.

On the pharmacokinetics of two inhaled budesonide/formoterol combinations in asthma patients using modeling approaches

Dry powder inhalers containing the budesonide/formoterol combination have currently a well-established position among other inhaled products. Even though their efficacy mainly depends on the local concentrations of the drug they deliver within the lungs, their safety profile is directly related to their total systemic exposure. The aim of the present investigation was to explore the absorption and disposition kinetics of the budesonide/formoterol combination delivered via two different dry powder inhalers in asthma patients. Plasma concentration-time data were obtained from a single-dose, crossover bioequivalence study in asthma patients. Non-compartmental and population compartmental approaches were applied to the available datasets. The non-compartmental analysis allowed for an initial characterization of the primary pharmacokinetic (PK) parameters of the two inhaled drugs and subsequently the bioequivalence assessment of the two different dry powder inhalers. The population pharmacokinetic analysis further explored the complex absorption and disposition characteristics of the two drugs. In case of inhaled FOR, a five-compartment PK model including an enterohepatic re-circulation process was developed. For inhaled BUD, the incorporation of two parallel first-order absorption rate constants (fast and slow) for lung absorption in a two-compartment PK model emphasized the importance of pulmonary anatomical features and underlying physiological processes during model development. The role of potential covariates on the variability of the PK parameters was also investigated.

Role of enterohepatic recirculation in drug disposition: cooperation and complications

Enterohepatic recirculation (EHC) concerns many physiological processes and notably affects pharmacokinetic parameters such as plasma half-life and AUC as well as estimates of bioavailability of drugs. Also, EHC plays a detrimental role as the compounds/drugs are allowed to recycle. An in-depth comprehension of this phenomenon and its consequences on the pharmacological effects of affected drugs is important and decisive in the design and development of new candidate drugs. EHC of a compound/drug occurs by biliary excretion and intestinal reabsorption, sometimes with hepatic conjugation and intestinal deconjugation. EHC leads to prolonged elimination half-life of the drugs, altered pharmacokinetics and pharmacodynamics. Study of the EHC of any drug is complicated due to unavailability of the apposite model, sophisticated procedures and ethical concerns. Different in vitro and in vivo methods for studies in experimental animals and humans have been devised, each having its own merits and demerits. Involvement of the different transporters in biliary excretion, intra- and inter-species, pathological and biochemical variabilities obscure the study of the phenomenon. Modeling of drugs undergoing EHC has always been intricate and exigent models have been exploited to interpret the pharmacokinetic profiles of drugs witnessing multiple peaks due to EHC. Here, we critically appraise the mechanisms of bile formation, factors affecting biliary drug elimination, methods to estimate biliary excretion of drugs, EHC, multiple peak phenomenon and its modeling.

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.

Empirical and semi-mechanistic modelling of double-peaked pharmacokinetic profile phenomenon due to gastric emptying

Models have been developed to explain double-peaked plasma concentration-time profiles using mechanisms such as variable absorption and enterohepatic recirculation. Interruption of gastric emptying has also been shown to produce double-peaks, and this work proposes models for analysis of such data. In the presence of levodopa, gastric emptying is interrupted at times associated with double-peaks in pharmacokinetic profiles. Data from a simultaneous scintigraphy and paracetamol absorption study with levodopa was obtained, and models with compartments for stomach, intestine, central and peripheral tissue were developed to describe levodopa and paracetamol pharmacokinetics, including the double-peak phenomenon. The empirical model uses two gastric emptying parameter rates which are applied over separate time periods to describe the varying gastric emptying rate. The semi-mechanistic model uses a feedback mechanism acting via an effect compartment to link the plasma concentration of levodopa to the rate of gastric emptying, allowing levodopa pharmacokinetics to vary the rate of gastric emptying and give rise to a multiple-peaked plasma pharmacokinetic profile. The models were applied to plasma levodopa and paracetamol pharmacokinetic data with and without simultaneous analysis of scintigraphy data, in both cases giving a good fit and in the absence of scintigraphy data adequately predicting the stomach profile. For the semi-mechanistic model, the first-order constant governing gastric emptying was shown to switch between fast and slow values at a critical levodopa effect compartment concentration. New models have thus been proposed for analysis of plasma concentration profiles that exhibit double-peak phenomenon and applied successfully to levodopa data.

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

AIM

This study aimed to develop a population pharmacokinetic (PK) enterohepatic recycling model for MPA in patients with childhood-onset systemic lupus erythematosus (cSLE).

METHODS

MPA concentration-time data were from outpatients on stable oral mycophenolate mofetil (MMF) and collected under fasting conditions, with standardized meals (1 and 4 h post-dose). Sampling times were pre-dose, 20, 40 min, 1, 1.5, 2, 3, 4, 6 and 9 h, post dose. The population PK analysis simultaneously modelled MPA and 7-O-MPA-β-glucuronide (MPAG) concentrations using nonlinear mixed effect modelling.

RESULTS

PK analysis included 186 MPA and MPAG concentrations (mg l(-1)) from 19 patients. cSLE patients, age range 10-28 years, median 16.5 years were included. Mean ± SD disease duration was 3.8 ± 3.7 years. The final PK model included a gallbladder compartment for enterohepatic recycling and bile release time related to meal times, with first order absorption and single series of transit compartments. The PK estimates for MPA were CL(1) /F 25.3 l h(-1), V(3) /F 20.9 l, V(4) /F 234 l and CL(2) /F 19.8 l h(-1).

CONCLUSION

The final model fitted the complex processes of absorption and enterohepatic circulation (EHC) in those treated with MMF for cSLE and could be applied in Bayesian dose optimization algorithms.

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.

PKSolver: An add-in program for pharmacokinetic and pharmacodynamic data analysis in Microsoft Excel

This study presents PKSolver, a freely available menu-driven add-in program for Microsoft Excel written in Visual Basic for Applications (VBA), for solving basic problems in pharmacokinetic (PK) and pharmacodynamic (PD) data analysis. The program provides a range of modules for PK and PD analysis including noncompartmental analysis (NCA), compartmental analysis (CA), and pharmacodynamic modeling. Two special built-in modules, multiple absorption sites (MAS) and enterohepatic circulation (EHC), were developed for fitting the double-peak concentration-time profile based on the classical one-compartment model. In addition, twenty frequently used pharmacokinetic functions were encoded as a macro and can be directly accessed in an Excel spreadsheet. To evaluate the program, a detailed comparison of modeling PK data using PKSolver and professional PK/PD software package WinNonlin and Scientist was performed. The results showed that the parameters estimated with PKSolver were satisfactory. In conclusion, the PKSolver simplified the PK and PD data analysis process and its output could be generated in Microsoft Word in the form of an integrated report. The program provides pharmacokinetic researchers with a fast and easy-to-use tool for routine and basic PK and PD data analysis with a more user-friendly interface.

The pharmacokinetics of ceftriaxone based on population pharmacokinetics and the prediction of efficacy in Japanese adults

OBJECTIVE

By using a population pharmacokinetic analysis method, we predicted the efficacy of Ceftriaxone (CTRX) based on the pharmacokinetics of CTRX in Japanese adults and the sensitivity of infective organisms to CTRX in 2004. In addition, we clarified the difference in efficacy between once-a-day administration and twice-a-day administration.

METHODS

The population pharmacokinetic analysis was based on the serum concentrations of CTRX already published by NONMEM. The possible effect of body weight and age on the pharmacokinetics of CTRX was examined using a model which incorporated the change of a specific protein-binding ratio of CTRX. A Monte Carlo simulation was conducted based on the population pharmacokinetic parameters obtained by this analysis, and thereby the time above MIC (TAM) was determined from the MIC values of CTRX administered once at 0.5, 1, and 2 g and twice at 1 g. The efficacy ratio was predicted from the TAM thus obtained.

RESULTS

Because the time course of serum concentration of CTRX in adult subjects was fitted to a 2-compartment model and both body weight and age were not incorporated as the covariate, the dosing method by which a fixed amount of CTRX is administered to patients has been thought to be adequate. Based on the efficacy ratio estimated from the MIC of CTRX, we have predicted that the once-a-day administration of CTRX even at 0.5g is effective on various infecting organisms.

Population pharmacokinetics of mycophenolic acid and its 2 glucuronidated metabolites in kidney transplant recipients

The population pharmacokinetics of mycophenolic acid (MPA) and its phenolic (MPAG) and acyl (AcMPAG) glucuronide metabolites were studied in patients taking enteric-coated mycophenolate sodium. Plasma samples (n = 232), obtained from 18 renal transplant recipients, were analyzed for MPA, MPAG, and AcMPAG using a validated high-performance liquid chromatography/ultraviolet assay. Population pharmacokinetic analysis was performed using NONMEM. The pharmacokinetics of MPA were best described by a 2-compartment model, with MPAG and AcMPAG produced from the central compartment and with enterohepatic recirculation of MPA via these 2 metabolites. Population mean estimates for MPA were apparent clearance (CL/F) of 10.6 L/h (interindividual variability [IIV] = 21.4%) and apparent volume of distribution of the central compartment (V(1)/F) of 25.9 L (IIV = 87.8%). Mean elimination rate constants of MPAG and AcMPAG were 0.323 h(-1) (IIV = 29.1%) and 0.206 h(-1) (IIV = 48.8%), respectively. The mean fraction of MPA converted to MPAG and AcMPAG, normalized by their volumes of distribution (FM(AG) and FM(AC), respectively), was also estimated. The elimination rate constant for MPAG and FM(AC) was influenced by glomerular filtration rate in patients with renal impairment. The visual predictive check, based on 100 simulated data sets each for MPA, MPAG, and AcMPAG, found that the final pharmacokinetic model adequately predicts the observed concentrations of all 3 species.

Prediction of the Possibility of the Second Peak of Drug Plasma Concentration Time Curve after iv Bolus Administration from the Standpoint of the Traditional Multi-Compartmental Linear Pharmacokinetics

It is shown that the existence of the second peak on the drug plasma concentration time curve C(p)(t) after iv bolus dosing can be explained by considering the traditional multi-compartmental linear pharmacokinetics. It was found that a direct solution of the general three-compartment model yields the second peak of C(p)(t) for the certain values of the rate constants, and C(p)(t) includes the term with oscillating preexponent, that is, K sin(omegat + phi) exp(-lambdat), in this case. The considered model describes the drug entero-hepatic recirculation in the species which do not have gall bladder (rats). The model fit of the experimental data from rat pharmacokinetic studies where the second peak of C(p)(t) was observed, yields the rate of bile production that is consistent with the physiological value ( approximately 0.7 mL/h).

Population pharmacokinetic modelling for enterohepatic circulation of mycophenolic acid in healthy Chinese and the influence of polymorphisms in UGT1A9

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT

Mycophenolic acid (MPA) undergoes enterohepatic circulation (EHC) in the body and several population models have been proposed to describe this process using sparse data. Recent studies in Whites have found that polymorphism in UGT1A9 could partly explain the large interindividual variability associated with the pharmacokinetics of MPA.

WHAT THIS STUDY ADDS

A new population pharmacokinetic model for EHC combining MPA and its main glucuronide metabolite (MPAG) simultaneously was established based on physiological aspects of biliary excretion using intensive sampling data. Pharmacokinetic profiles of MPA and MPAG with the UGT1A9 polymorphism in healthy Chinese were characterized. AIMS To establish a population pharmacokinetic model that describes enterohepatic circulation (EHC) of mycophenolic acid (MPA) based on physiological considerations and to investigate the influence of polymorphisms of UGT1A9 on the pharmacokinetics of MPA.

METHODS

Pharmacokinetic data were obtained from two comparative bioavailability studies of oral mycophenolic mofetil formulations. Nonlinear mixed effects modelling was employed to develop an EHC model including both MPA and its main glucuronide metabolite (MPAG) simultaneously. Demographic characteristics and UGT1A9 polymorphisms were screened as covariates.

RESULTS

In total, 590 MPA and 589 MPAG concentration-time points from 42 healthy male volunteers were employed in this study. The chain compartment model included an intestinal compartment, a gallbladder compartment, a central and a peripheral compartment for MPA and a central compartment for MPAG. The typical population clearance (CL/F) estimates with its relative standard error for MPA and MPAG were 10.2 l h(-1) (5.7%) and 1.38 l h(-1) (6.9%), respectively. The amount of MPA recycled in the body was estimated to be 29.1% of the total amount absorbed. Covariate analysis showed that body weight was positively correlated with CL/F of MPA, intercompartment CL/F of MPA and distribution volume of MPA peripheral compartment. Polymorphisms of UGT1A9 did not show any effect on the pharmacokinetics of MPA and MPAG. The model evaluation tests indicated that the proposed model can describe the pharmacokinetic profiles of MPA and MPAG in healthy Chinese subjects.

CONCLUSIONS

The proposed model may provide a valuable approach for planning future pharmacokinetic-pharmacodynamic studies and for designing proper dosage regimens of MPA.

Population pharmacokinetic modeling for enterohepatic recirculation in Rhesus monkey

Enterohepatic recirculation (EHR) occurs via biliary excretion and intestinal reabsorption of a drug. Drug recycling through EHR can lead to a change in pharmacokinetic (PK) properties, such as reduced clearance (CL), extended half-life (T(1/2)) and increased plasma exposure (AUC). As a result, EHR may prolong the pharmacological effect of drugs. In the present study, the compound (Cpd A) was found to exhibit EHR in Rhesus monkeys associated with a reduction in CL (from 3.8 to 0.33 Lh(-1), IV; from 2.3 to 0.4 Lh(-1), PO), and an increase in T(1/2) (from 0.9 to 18 h, IV) and in AUC (from 1.5 to 17.4 microg h/mL, IV; from 2.8 to 16.3 microg h/mL, PO), by comparing the PK in the monkeys via the interruption of EHR (bile-duct cannulation) with that in the intact monkeys. A population four-compartment model was constructed based on recirculation loops incorporating all possible inputs (bile secretion, a lag-time model for gall bladder emptying, routes and amounts of a single dose administration) to fully evaluate the EHR of Cpd A. The plasma concentrations versus time profiles predicted from the model had a good fit to the values observed in the subjects and were further simulated with 90% confidence interval to demonstrate its utility. Thus, the model could be applied as a useful tool to evaluate the drugs or compounds that undergo EHR in different species.

Population pharmacokinetics of mycophenolic acid in kidney transplant pediatric and adolescent patients

Current data on mycophenolate mofetil (MMF) suggest that there is a pharmacokinetic/pharmacodynamic relationship between the mycophenolic acid (MPA) area under the curve (AUC) during treatment and both the risk of acute rejection and the occurrence of side effects. The aim of this study was to characterize the population pharmacokinetics of MPA in kidney transplant patients between the ages of 2 and 21 years and to propose a limited sampling strategy to estimate individual MPA AUCs. Forty-one patients received long-term oral MMF continuous therapy as part of a triple immunosuppressive regimen, which also included cyclosporine or tacrolimus (n=3) and corticosteroids. Therapy was initiated at a dose of 600 mg/m twice daily. The population parameters were calculated from an initial group of 32 patients. The data were analyzed by nonlinear mixed-effect modeling using a 2-compartment structural model with first-order absorption and a lag time. The interindividual variability in the initial volume of distribution was partially explained by the fact that this parameter was weight-dependent. Fifteen concentration-time profiles from 13 patients were used to evaluate the predictive performance of the Bayesian approach and to devise a limited sampling strategy. The protocol, involving two sampling times, 1 and 4 hours after oral administration, allows the precise and accurate determination of MPA AUCs (bias -0.9 microg.h/mL; precision 6.02 microg.h/mL). The results of this study combine the relationships between the pharmacokinetic parameters of MPA and patient covariates, which may be useful for dose adjustment, with a convenient sampling procedure that may aid in optimizing pediatric patient care.

Raloxifene pharmacokinetics in males with normal and impaired renal function

AIMS

To compare raloxifene pharmacokinetics between renally impaired and healthy subjects.

METHODS

Raloxifene 120 mg was administered to 10 males with renal impairment (creatinine 2-4 mg dl(-1)) and to 10 healthy males. Data were analysed by two noncompartmental and one compartmental nonlinear regression methods.

RESULTS

The medians (95% confidence interval) of the area under the curves (AUC) were 35.1 (25.8, 74) and 20.5 (16.8, 28.0) h ng ml(-1) per mg kg(-1), P < 0.01, and of the clearances (CL/F) were 28.5 (13.5, 38.8) and 48.8 (35.8, 59.4) l h(-1) kg(-1), P < 0.01, in renally impaired and healthy subjects, respectively. 95% Confidence intervals on the differences for AUC and CL/F were 6.5-44.1 and -35.1 to -7.9, respectively.

CONCLUSION

Exposure to raloxifene was twice as high in males with renal impairment compared with healthy subjects.

A Double Absorption-Phase Model Adequately Describes Mycophenolic Acid Plasma Profiles in De Novo Renal Transplant Recipients Given Oral Mycophenolate Mofetil

BACKGROUND

Mycophenolic acid (MPA) shows complex plasma concentration-time profiles, particularly in the immediate (first month) post-transplantation phase for which no relevant pharmacokinetic model has been proposed thus far.

OBJECTIVE

The aim of this study was to develop a model to accurately describe the time profile of plasma MPA concentrations after oral administration of mycophenolate mofetil in adult kidney transplant patients, in any post-transplantation period.

METHOD

Full interdose pharmacokinetic profiles were collected in 45 adult renal transplant patients who were orally administered mycophenolate mofetil and ciclosporin; 25 patients were de novo transplant patients for whom individual pharmacokinetics were assessed at three post-transplantation periods (days 3, 7 and 30) and 20 patients were stable transplant patients (>3 months post-transplantation). MPA was determined in plasma by liquid chromatography-mass spectrometry. Models combining a single- or double-input (described as single or double gamma distributions) with one- or two-compartments were developed using in-house software and fitted to the individual profiles by nonlinear regression.

RESULTS

Visual inspection of the pharmacokinetic profiles showed highly variable absorption profiles and secondary peaks of various intensity. The pharmacokinetic models including a double gamma distribution best fitted these various profiles in the immediate post-transplantation period (mean bias and precision of -0.92% and 20.19%; -1.5% and 18.02%, on day 7 and day 30, respectively), while in the stable post-grafting phase (beyond 3 months), the single- and double-absorption models performed similarly (mean bias and precision of -3.37% and 17.64%; -3.12% and 18.44%, on day 7 and day 30, respectively).

CONCLUSION

The proposed pharmacokinetic models adequately describe the concentration-time profiles of MPA in renal transplant patients and could be helpful in the development of tools for MPA monitoring.

Using Established Immunosuppressant Therapy Effectively

Getting the most effective use of immunosuppressant medications in transplant patients continues to be a major challenge in clinical practice. This need applies to established immunosuppressants just as it does to new agents. In this review this principle is illustrated for mycophenolic acid, the active metabolite of mycophenolate mofetil, the most commonly used immunosuppressant, in various combinations with other immunosuppressants, in current clinical practice. Defining, as rigorously as possible, the requirements for effective therapeutic monitoring of MPA is an important goal given all of the changes in immunosuppressive drug regimens. This review will focus on the major factors known to influence MPA clearance including: UDP-glucuronyltransferases, enterohepatic circulation, MPA free fraction, the effect of time posttransplantation, concomitant administration of immunosuppressant drugs. The significant variability of MPA pharmacokinetics and the need to deepen our understanding of the influence of these factors on MPA clearance are strong reasons why additional clinical trials are needed to define best practice therapeutic drug monitoring of this drug.

Population pharmacokinetics and Bayesian estimation of mycophenolic acid concentrations in stable renal transplant patients

BACKGROUND

Therapeutic drug monitoring of mycophenolic acid (MPA) may minimise the risk of acute rejection after transplantation. Area under the curve (AUC) rather than trough concentration-based monitoring is recommended and models for AUC estimation are needed.

OBJECTIVES

To develop a population pharmacokinetic model suitable for Bayesian estimation of individual AUC in stable renal transplant patients.

PATIENTS AND METHODS

The population pharmacokinetics of MPA were studied using nonlinear mixed effects modelling (NONMEM) in 60 patients (index group) receiving MPA on a twice-daily basis. Ten blood samples were collected at fixed timepoints from ten patients and four blood samples were collected at sparse timepoints from 50 patients. Bayesian estimation of individual AUC was made on the basis of three blood concentration measurements and covariates. The predictive performances of the Bayesian procedure were evaluated in an independent group of patients (test group) comprising ten subjects in whom ten blood samples were collected at fixed timepoints.

RESULTS

A two-compartment model with zero-order absorption best fitted the data. Covariate analysis showed that bodyweight was positively correlated with oral clearance. However, the weak magnitude of the reduction in variability (from 34.8 to 28.2%) indicates that administration on a per kilogram basis would be of limited value in decreasing interindividual variability in MPA exposure. Bayesian estimation of pharmacokinetic parameters using samples drawn at 20 minutes and 1 and 3 hours enabled estimation of individual AUC with satisfactory accuracy (bias 7.7%, range of prediction errors 0.43-15.1%) and precision (root mean squared error 12.4%) as compared with the reference value obtained using the trapezoidal method.

CONCLUSION

This paper reports for the first time population pharmacokinetic data for MPA in stable renal transplant patients, and shows that Bayesian estimation can allow accurate prediction of AUC with only three samples. This method provides a tool for therapeutic drug monitoring of MPA or for concentration-effect studies. Its application to MPA monitoring in the early period post-transplantation needs to be evaluated.

A recirculatory model with enterohepatic circulation by measuring portal and systemic blood concentration difference

The present study describes a recirculatory model for the evaluation of pharmacokinetic characteristics of drugs possessing enterohepatic circulation (EHC). The advantage of the model is to separately define the extent and rate of absorption for the dosage and EHC after oral administration. Cephradine was used as a model drug and was intravenously or orally administered to rats. Portal and systemic bloods were simultaneously collected in order to estimate various local moments after defining the global moments obtained by non-compartment analysis. For the zero-order moments, bioavailability (BA), the hepatic recovery ratio (Fh), the sum of the local absorption ratio for the dosage and recirculatory local absorption ratio for EHC (F(a)po), and the recirculatory local absorption ratio for EHC (F(a)ehc) after oral administration were estimated to be 95.6, 77.9, 172, and 71.5%, respectively. These data indicate that a complete absorption and substantial EHC contribute high oral exposure of cephradine. For the first-order moments, the sum of the mean local absorption times for the dosage and EHC (t(a)po) and the mean transit time for a single pass of EHC (tc) were 2.50 and 0.117 hr, suggesting a rapid EHC of cephradine compared with the absorption from the dosage. With this model, the absorption rate-time profiles for the dosage and EHC were separately simulated by using a program of nonlinear least squares (MULTI) with fast inverse Laplace transform (FILT). The cumulative biliary excretion ratio (Fbile) calculated by the model was in good agreement with the experimental value obtained in the bile ductcannulated rats. These results suggest that the model proposed in this study would be useful for evaluating the extent and rate of ECH along with absorption from the dosage after oral administration of drugs.

Enterohepatic circulation: physiological, pharmacokinetic and clinical implications

Enterohepatic recycling occurs by biliary excretion and intestinal reabsorption of a solute, sometimes with hepatic conjugation and intestinal deconjugation. Cycling is often associated with multiple peaks and a longer apparent half-life in a plasma concentration-time profile. Factors affecting biliary excretion include drug characteristics (chemical structure, polarity and molecular size), transport across sinusoidal plasma membrane and canniculae membranes, biotransformation and possible reabsorption from intrahepatic bile ductules. Intestinal reabsorption to complete the enterohepatic cycle may depend on hydrolysis of a drug conjugate by gut bacteria. Bioavailability is also affected by the extent of intestinal absorption, gut-wall P-glycoprotein efflux and gut-wall metabolism. Recently, there has been a considerable increase in our understanding of the role of transporters, of gene expression of intestinal and hepatic enzymes, and of hepatic zonation. Drugs, disease and genetics may result in induced or inhibited activity of transporters and metabolising enzymes. Reduced expression of one transporter, for example hepatic canalicular multidrug resistance-associated protein (MRP) 2, is often associated with enhanced expression of others, for example the usually quiescent basolateral efflux MRP3, to limit hepatic toxicity. In addition, physiologically relevant pharmacokinetic models, which describe enterohepatic recirculation in terms of its determinants (such as sporadic gall bladder emptying), have been developed. In general, enterohepatic recirculation may prolong the pharmacological effect of certain drugs and drug metabolites. Of particular importance is the potential amplifying effect of enterohepatic variability in defining differences in the bioavailability, apparent volume of distribution and clearance of a given compound. Genetic abnormalities, disease states, orally administered adsorbents and certain coadministered drugs all affect enterohepatic recycling.

A pharmacokinetic model for analysis of drug disposition profiles undergoing enterohepatic circulation

A new and simple pharmacokinetic model that can explain enterohepatic circulation profiles for both single and repeated dosing was developed, and its applicability and usefulness were assessed by an actual published data set and simulation study. The model is basically a conventional compartment model, and the transfer rate from the bile compartment to the central compartment is assumed to change periodically, with the sine function being used to describe this periodical change. Using this model, the effect of the parameter values on plasma time-course profiles was examined by simulation, and the applicability of the model was tested by curve fitting to obtain the parameter estimates using an actual published data set. These studies confirmed that our model can simulate the periodical increase of the concentration due to re-absorption. By averaging the sine function in the transfer rate from the bile compartment to the central compartment, a smoothed time-course profile in the elimination phase that is independent of the enterohepatic cycle can be obtained. Also, the apparent half-life in the elimination phase can be estimated, whith is useful especially for evaluating drug accumulation during repeated dosing. It was suggested that the present model can be used to evaluate the drug disposition profile with enterohepatic circulation. The effects of sampling points and sampling time on parameter estimation are also discussed.