Mechanism-Based Modeling of Gastric Emptying Rate and Gallbladder Emptying in Response to Caloric Intake

Usefulness of the Beagle Model in the Evaluation of Paracetamol and Ibuprofen Exposure after Oral Administration to Pediatric Populations: An Exploratory Study

The present study aimed to explore the usefulness of beagle dogs in combination with physiologically based pharmacokinetic (PBPK) modeling in the evaluation of drug exposure after oral administration to pediatric populations at an early stage of pharmaceutical product development. An exploratory, single-dose, crossover bioavailability study in six beagles was performed. A paracetamol suspension and an ibuprofen suspension were coadministered in the fasted-state conditions, under reference-meal fed-state conditions, and under infant-formula fed-state conditions. PBPK models developed with GastroPlus v9.7 were used to inform the extrapolation of beagle data to human infants and children. Beagle-based simulation outcomes were compared with published human-adult-based simulations. For paracetamol, fasted-state conditions and reference-meal fed-state conditions in beagles appeared to provide adequate information for the applied scaling approach. Fasted-state and/or reference-meal fed-state conditions in beagles appeared suitable to simulate the performance of ibuprofen suspension in pediatric populations. Contrary to human-adult-based translations, extrapolations based on beagle data collected under infant-formula fed-state conditions appeared less useful for informing simulations of plasma levels in pediatric populations. Beagle data collected under fasted and/or reference-meal fed-state conditions appeared to be useful in the investigation of pediatric product performance of the two investigated highly permeable and highly soluble drugs in the upper small intestine. The suitability of the beagle as a preclinical model to understand pediatric drug product performance under different dosing conditions deserves further evaluation with a broader spectrum of drugs and drug products and comparisons with pediatric in vivo data.

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.

A glucose-insulin-glucagon coupled model of the isoglycemic intravenous glucose infusion experiment

Type 2 diabetes (T2D) is a pathophysiology that is characterized by insulin resistance, beta- and alpha-cell dysfunction. Mathematical models of various glucose challenge experiments have been developed to quantify the contribution of insulin and beta-cell dysfunction to the pathophysiology of T2D. There is a need for effective extended models that also capture the impact of alpha-cell dysregulation on T2D. In this paper a delay differential equation-based model is developed to describe the coupled glucose-insulin-glucagon dynamics in the isoglycemic intravenous glucose infusion (IIGI) experiment. As the glucose profile in IIGI is tailored to match that of a corresponding oral glucose tolerance test (OGTT), it provides a perfect method for studying hormone responses that are in the normal physiological domain and without the confounding effect of incretins and other gut mediated factors. The model was fit to IIGI data from individuals with and without T2D. Parameters related to glucagon action, suppression, and secretion as well as measures of insulin sensitivity, and glucose stimulated response were determined simultaneously. Significant impairment in glucose dependent glucagon suppression was observed in patients with T2D (duration of T2D: 8 (6–36) months) relative to weight matched control subjects (CS) without diabetes (k1 (mM)−1: 0.16 ± 0.015 (T2D, n = 7); 0.26 ± 0.047 (CS, n = 7)). Insulin action was significantly lower in patients with T2D (a1 (10 pM min)−1: 0.000084 ± 0.0000075 (T2D); 0.00052 ± 0.00015 (CS)) and the Hill coefficient in the equation for glucose dependent insulin response was found to be significantly different in T2D patients relative to CS (h: 1.4 ± 0.15; 1.9 ± 0.14). Trends in parameters with respect to fasting plasma glucose, HbA1c and 2-h glucose values are also presented. Significantly, a negative linear relationship is observed between the glucagon suppression parameter, k1, and the three markers for diabetes and is thus indicative of the role of glucagon in exacerbating the pathophysiology of diabetes (Spearman Rank Correlation: (n = 12; (−0.79, 0.002), (−0.73,.007), (−0.86,.0003)) respectively).

A Population Pharmacokinetic-Pharmacodynamic Model of Navtemadlin, its Major Active Metabolite (M1) and Serum Macrophage Inhibitory Cykokine-1 (MIC-1)

Navtemadlin is a potent, selective, orally available inhibitor of murine double minute 2 that restores p53 activity to induce apoptosis in TP53 wild-type malignancies. Using richly sampled pharmacokinetic (PK) and pharmacodynamic (PD) data from healthy volunteers, a population PK/PD model was developed. A population PK (PPK) model described the PK characteristics of navtemadlin and its major metabolite acyl glucuronide (M1) and quantified enterohepatic recirculation (EHR). Post hoc individual PK parameters from this model were coupled with PD data for serum macrophage inhibitory cytokine-1 (MIC-1, GDF15), a cytokine biomarker of p53 activation, to construct a population PK/PD model that described plasma concentration-driven MIC-1 excursions and enabled simulation of the extent and duration of navtemadlin PD effects. The median apparent clearance (CL/F) and apparent central volume (V2/F) of navtemadlin were 36.4 L/hr and 159 L. The typical maximum stimulatory effect (S_max) was close to the median maximum MIC-1 ratio to baseline of 7.29 in observed data. Simulation revealed a dose-dependent increase of MIC-1 with steady state attained in approximately 7 days, in a 7-day-on/21-day-off dose regimen. Elevated MIC-1 concentrations persist through 17-19 days, leaving about 9-11 PD-free days in a 28-day cycle.

The Use of Physiologically Based Pharmacokinetic Analyses-in Biopharmaceutics Applications -Regulatory and Industry Perspectives

The use of physiologically based pharmacokinetic (PBPK) modeling to support the drug product quality attributes, also known as physiologically based biopharmaceutics modeling (PBBM) is an evolving field and the interest in using PBBM is increasing. The US-FDA has emphasized on the use of patient centric quality standards and clinically relevant drug product specifications over the years. Establishing an in vitro in vivo link is an important step towards achieving the goal of patient centric quality standard. Such a link can aid in constructing a bioequivalence safe space and establishing clinically relevant drug product specifications. PBBM is an important tool to construct a safe space which can be used during the drug product development and lifecycle management. There are several advantages of using the PBBM approach, though there are also a few challenges, both with in vitro methods and in vivo understanding of drug absorption and disposition, that preclude using this approach and therefore further improvements are needed. In this review we have provided an overview of experience gained so far and the current perspective from regulatory and industry point of view. Collaboration between scientists from regulatory, industry and academic fields can further help to advance this field and deliver on promises that PBBM can offer towards establishing patient centric quality standards.

Prior information for population pharmacokinetic and pharmacokinetic/pharmacodynamic analysis: overview and guidance with a focus on the NONMEM PRIOR subroutine

Abstract

Population pharmacokinetic analysis is used to estimate pharmacokinetic parameters and their variability from concentration data. Due to data sparseness issues, available datasets often do not allow the estimation of all parameters of the suitable model. The PRIOR subroutine in NONMEM supports the estimation of some or all parameters with values from previous models, as an alternative to fixing them or adding data to the dataset. From a literature review, the best practices were compiled to provide a practical guidance for the use of the PRIOR subroutine in NONMEM. Thirty-three articles reported the use of the PRIOR subroutine in NONMEM, mostly in special populations. This approach allowed fast, stable and satisfying modelling. The guidance provides general advice on how to select the most appropriate reference model when there are several previous models available, and to implement and weight the selected parameter values in the PRIOR function. On the model built with PRIOR, the similarity of estimates with the ones of the reference model and the sensitivity of the model to the PRIOR values should be checked. Covariates could be implemented a priori (from the reference model) or a posteriori, only on parameters estimated without prior (search for new covariates).

Graphic abstract

Electronic supplementary material

The online version of this article (10.1007/s10928-020-09695-z) contains supplementary material, which is available to authorized users.

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.

Prediction of Clinical Transporter-Mediated Drug-Drug Interactions via Comeasurement of Pitavastatin and Eltrombopag in Human Hepatocyte Models

A structurally identifiable micro‐rate constant mechanistic model was used to describe the interaction between pitavastatin and eltrombopag, with improved goodness‐of‐fit values through comeasurement of pitavastatin and eltrombopag. Transporter association and dissociation rate constants and passive rates out of the cell were similar between pitavastatin and eltrombopag. Translocation into the cell through transporter‐mediated uptake was six times greater for pitavastatin, leading to pronounced inhibition of pitavastatin uptake by eltrombopag. The passive rate into the cell was 91 times smaller for pitavastatin compared with eltrombopag. A semimechanistic physiologically‐based pharmacokinetic (PBPK) model was developed to evaluate the potential for clinical drug–drug interactions (DDIs). The PBPK model predicted a twofold increase in the pitavastatin peak blood concentration and area under the concentration‐time curve in the presence of eltrombopag in simulated healthy volunteers. The use of structural identifiability supporting experimental design combined with robust micro‐rate constant parameter estimates and a semimechanistic PBPK model gave more informed predictions of transporter‐mediated DDIs.

A Physiology-Based Model of Bile Acid Distribution and Metabolism Under Healthy and Pathologic Conditions in Human Beings

BACKGROUND & AIMS

Disturbances of the enterohepatic circulation of bile acids (BAs) are seen in a number of clinically important conditions, including metabolic disorders, hepatic impairment, diarrhea, and gallstone disease. To facilitate the exploration of underlying pathogenic mechanisms, we developed a mathematical model built on quantitative physiological observations across different organs.

METHODS

The model consists of a set of kinetic equations describing the syntheses of cholic, chenodeoxycholic, and deoxycholic acids, as well as time-related changes of their respective free and conjugated forms in the systemic circulation, the hepatoportal region, and the gastrointestinal tract. The core structure of the model was adapted from previous modeling research and updated based on recent mechanistic insights, including farnesoid X receptor-mediated autoregulation of BA synthesis and selective transport mechanisms. The model was calibrated against existing data on BA distribution and feedback regulation.

RESULTS

According to model-based predictions, changes in intestinal motility, BA absorption, and biotransformation rates affected BA composition and distribution differently, as follows: (1) inhibition of transintestinal BA flux (eg, in patients with BA malabsorption) or acceleration of intestinal motility, followed by farnesoid X receptor down-regulation, was associated with colonic BA accumulation; (2) in contrast, modulation of the colonic absorption process was predicted to not affect the BA pool significantly; and (3) activation of ileal deconjugation (eg, in patents with small intestinal bacterial overgrowth) was associated with an increase in the BA pool, owing to higher ileal permeability of unconjugated BA species.

CONCLUSIONS

This model will be useful in further studying how BA enterohepatic circulation modulation may be exploited for therapeutic benefits.

Model-Based Prediction of Plasma Concentration and Enterohepatic Circulation of Total Bile Acids in Humans

Bile acids released postprandially can modify the rate and extent of lipophilic compounds' absorption. This study aimed to predict the enterohepatic circulation (EHC) of total bile acids (TBAs) in response to caloric intake from their spillover in plasma. A model for TBA EHC was combined with a previously developed gastric emptying (GE) model. Longitudinal gallbladder volumes and TBA plasma concentration data from 30 subjects studied after ingestion of four different test drinks were supplemented with literature data. Postprandial gallbladder refilling periods were implemented to improve model predictions. The TBA hepatic extraction was reduced with the high-fat drink. Basal and nutrient-induced gallbladder emptying rates were altered by type 2 diabetes (T2D). The model was predictive of the central trend and the variability of gallbladder volume and TBA plasma concentration for all test drinks. Integration of this model within physiological pharmacokinetic modeling frameworks could improve the predictions for lipophilic compounds' absorption considerably.