Evolution of a detailed physiological model to simulate the gastrointestinal transit and absorption process in humans, Part 1: Oral solutions

Analyzing Parametric Influences Driving Age-Associated Changes in Absorption using a PBPK-GSA Approach

The advanced age population may be susceptible to an increased risk of adverse effects due to increased drug exposure after oral dosing. Factors such as high-interindividual variability and lack of data has led to poor characterization of absorption's role in pharmacokinetic changes in this population. Physiologically based pharmacokinetic (PBPK) models are increasingly being used during the drug development process, as their unique qualities are advantageous in atypical scenarios such as drug-drug interactions or special populations such as older people. Along with relying on various sources of data, auxiliary tools including parameter estimation and sensitivity analysis techniques are employed to support model development and other applications. However, sensitivity analyses have mostly been limited to localized techniques in the majority of reported PBPK models using them. This is disadvantageous, since local sensitivity analyses are unsuitable for risk analysis, which require assessment of parametric interactions and proper coverage of the input space to better estimate and subsequently mitigate the effects of the phenomenon of interest. For this reason, this study seeks to integrate a global sensitivity analysis screening method with PBPK models based in PK-Sim® to characterize the consequences of potential changes in absorption that are often associated with advanced age. The Elementary Effects (Morris) method and visualization of the results are implemented in R and three model drugs representing Biopharmaceutical Classification System classes I-III that are expected to exhibit some sensitivity to three age-associated hypotheses were successfully tested.

A physiologically-based pharmacokinetic modeling approach for dosing amiodarone in children on ECMO

Extracorporeal membrane oxygenation (ECMO) is a cardiopulmonary bypass device commonly used to treat cardiac arrest in children. The American Heart Association guidelines for cardiopulmonary resuscitation (CPR) and emergency cardiovascular care recommend using amiodarone as a first-line agent to treat ventricular arrhythmias in children with cardiac arrest. However, there are no dosing recommendations for amiodarone to treat ventricular arrhythmias in pediatric patients on ECMO. Amiodarone has a high propensity for adsorption to the ECMO components due to its physicochemical properties leading to altered pharmacokinetics (PK) in ECMO patients. The change in amiodarone PK due to interaction with ECMO components may result in a difference in optimal dosing in patients on ECMO when compared with non-ECMO patients. To address this clinical knowledge gap, a physiologically-based pharmacokinetic model of amiodarone was developed in adults and scaled to children, followed by the addition of an ECMO compartment. The pediatric model included ontogeny functions of cytochrome P450 (CYP450) enzyme maturation across various age groups. The ECMO compartment was parameterized using the adsorption data of amiodarone obtained from ex vivo studies. Model predictions captured observed concentrations of amiodarone in pediatric patients with ECMO well with an average fold error between 0.5 and 2. Model simulations support an amiodarone intravenous (i.v) bolus dose of 22 mg/kg (neonates), 13 mg/kg (infants), 8 mg/kg (children), and 6 mg/kg (adolescents). This PBPK modeling approach can be applied to explore the dosing of other drugs used in children on ECMO.

Is there a fast track ("Darmstrasse") for fluids in the small intestine? Evidence from magnetic resonance imaging

BACKGROUND

The transit and distribution pattern of fluids in the small intestine is a key parameter for the dissolution and absorption of drugs. Although some information is known about the small intestinal water content after administration of fluid volumes and meals, the intestinal transit of orally ingested fluids and solutions has been barely investigated. The aim of this three-arm, cross-over, 9-subject human study was to investigate the transit of orally ingested water in the small intestine under fasting and postprandial conditions using MRI. To identify the ingested water, manganese gluconate, which can be identified with T1-weighted MRI sequences, was added as a marker. Using Horos (DICOM software), quantification of the distribution of Mn2+ ions in the gastrointestinal tract in fasted versus fed state (standard meal by FDA guidance and a light meal) was possible. The distribution and approximate wetted intestinal length was very similar in the fasting and postprandial states, suggesting rapid transport of water ingested after a meal through the chyme-filled small intestine in continuation of the "Magenstrasse" (stomach road). In some subjects, manganese gluconate reached deeper parts of the small intestine even more quickly in the postprandial state than in the fasting arm of the study. A deeper understanding of the behaviour of solutes in the gastrointestinal tract is fundamental to a mechanistic explanation for the kinetic interaction between food and drug intake (food effects).

Physiologically Based Pharmacokinetic Modeling to Unravel the Drug-gene Interactions of Venlafaxine: Based on Activity Score-dependent Metabolism by CYP2D6 and CYP2C19 Polymorphisms

BACKGROUND

Venlafaxine (VEN) is a commonly utilized medication for alleviating depression and anxiety disorders. The presence of genetic polymorphisms gives rise to considerable variations in plasma concentrations across different phenotypes. This divergence in phenotypic responses leads to notable differences in both the efficacy and tolerance of the drug.

PURPOSE

A physiologically based pharmacokinetic (PBPK) model for VEN and its metabolite O-desmethylvenlafaxine (ODV) to predict the impact of CYP2D6 and CYP2C19 gene polymorphisms on VEN pharmacokinetics (PK).

METHODS

The parent-metabolite PBPK models for VEN and ODV were developed using PK-Sim® and MoBi®. Leveraging prior research, derived and implemented CYP2D6 and CYP2C19 activity score (AS)-dependent metabolism to simulate exposure in the drug-gene interactions (DGIs) scenarios. The model's performance was evaluated by comparing predicted and observed values of plasma concentration-time (PCT) curves and PK parameters values.

RESULTS

In the base models, 91.1%, 94.8%, and 94.6% of the predicted plasma concentrations for VEN, ODV, and VEN + ODV, respectively, fell within a twofold error range of the corresponding observed concentrations. For DGI scenarios, these values were 81.4% and 85% for VEN and ODV, respectively. Comparing CYP2D6 AS = 2 (normal metabolizers, NM) populations to AS = 0 (poor metabolizers, PM), 0.25, 0.5, 0.75, 1.0 (intermediate metabolizers, IM), 1.25, 1.5 (NM), and 3.0 (ultrarapid metabolizers, UM) populations in CYP2C19 AS = 2.0 group, the predicted DGI AUC0-96 h ratios for VEN were 3.65, 3.09, 2.60, 2.18, 1.84, 1.56, 1.34, 0.61, and for ODV, they were 0.17, 0.35, 0.51, 0.64, 0.75, 0.83, 0.90, 1.11, and the results were similar in other CYP2C19 groups. It should be noted that PK differences in CYP2C19 phenotypes were not similar across different CYP2D6 groups.

CONCLUSIONS

In clinical practice, the impact of genotyping on the in vivo disposition process of VEN should be considered to ensure the safety and efficacy of treatment.

Physiologically based pharmacokinetic modeling of candesartan to predict the exposure in hepatic and renal impairment and elderly populations

Background

Candesartan cilexetil is a widely used angiotensin II receptor blocker with minimal adverse effects and high tolerability for the treatment of hypertension. Candesartan is administered orally as the prodrug candesartan cilexetil, which is wholly and swiftly converted to the active metabolite candesartan by carboxylesterase during absorption in the intestinal tract. In populations with renal or hepatic impairment, candesartan's pharmacokinetic (PK) behavior may be altered, necessitating dosage adjustments.

Objectives

This study was conducted to examine how the physiologically based PK (PBPK) model characterizes the PKs of candesartan in adult and geriatric populations and to predict the PKs of candesartan in elderly populations with renal and hepatic impairment.

Design

After developing PBPK models using the reported physicochemical properties of candesartan and clinical data, these models were validated using data from clinical investigations involving various dose ranges.

Methods

Comparing predicted and observed blood concentration data and PK parameters was used to assess the fit performance of the models.

Results

Doses should be reduced to approximately 94% of Chinese healthy adults for the Chinese healthy elderly population; approximately 92%, 68%, and 64% of that of the Chinese healthy adult dose in elderly populations with mild, moderate, and severe renal impairment, respectively; and approximately 72%, 71%, and 52% of that of the Chinese healthy adult dose in elderly populations with Child-Pugh-A, Child-Pugh-B, and Child-Pugh-C hepatic impairment, respectively.

Conclusion

The results suggest that the PBPK model of candesartan can be utilized to optimize dosage regimens for special populations.

Physiologically-based pharmacokinetic model to predict doxorubicin and paclitaxel exposure in infants through breast milk

Limited information is available concerning infant exposure and safety when breastfed by mothers receiving chemotherapy. Whereas defining distribution to breast milk is important to infer drug exposure, infant pharmacokinetics also determine to what extent the infant will be exposed to potential toxic effects. We aimed to assess the impact of chemotherapy containing breast milk on infants by predicting systemic and local (intestinal) exposure of paclitaxel and doxorubicin in infants through breast milk using a physiologically-based pharmacokinetic (PBPK) approach. Whole-body PBPK models of i.v. paclitaxel and doxorubicin were extended from the literature, with an oral absorption component to enable predictions in infants receiving paclitaxel or doxorubicin-containing breast milk. For safety considerations, worst-case scenarios were explored. Finally, paclitaxel and doxorubicin exposures in plasma and intestinal tissue of infants following feeding of breast milk from paclitaxel- or doxorubicin-treated mothers were simulated and breast milk discarding strategies were evaluated. The upper 95th percentile of the predicted peak concentrations in peripheral venous blood were 3.48 and 0.74 nM (0.4%-1.7% and 0.1%-1.8% of on-treatment) for paclitaxel and doxorubicin, respectively. Intestinal exposure reached peak concentrations of 1.0 and 140 μM for paclitaxel and doxorubicin, respectively. Discarding breast milk for the first 3 days after maternal chemotherapy administration reduced systemic and tissue exposures even further, to over 90% and 80% for paclitaxel and doxorubicin, respectively. PBPK simulations of chemotherapy exposure in infants after breastfeeding with chemotherapy containing breast milk suggest that particularly local gastrointestinal adverse events should be monitored, whereas systemic adverse events are not expected.

Computational Amendment of Parenteral In Situ Forming Particulates' Characteristics: Design of Experiment and PBPK Physiological Modeling

Lipid and/or polymer-based drug conjugates can potentially minimize side effects by increasing drug accumulation at target sites and thus augment patient compliance. Formulation factors can present a potent influence on the characteristics of the obtained systems. The selection of an appropriate solvent with satisfactory rheological properties, miscibility, and biocompatibility is essential to optimize drug release. This work presents a computational study of the effect of the basic formulation factors on the characteristics of the obtained in situ-forming particulates (IFPs) encapsulating a model drug using a 21.31 full factorial experimental design. The emulsion method was employed for the preparation of lipid and/or polymer-based IFPs. The IFP release profiles and parameters were computed. Additionally, a desirability study was carried out to choose the optimum formulation for further morphological examination, rheological study, and PBPK physiological modeling. Results revealed that the type of particulate forming agent (lipid/polymer) and the incorporation of structure additives like Brij 52 and Eudragit RL can effectively augment the release profile as well as the burst of the drug. The optimized formulation exhibited a pseudoplastic rheological behavior and yielded uniformly spherical-shaped dense particulates with a PS of 573.92 ± 23.5 nm upon injection. Physiological modeling simulation revealed the pioneer pharmacokinetic properties of the optimized formulation compared to the observed data. These results assure the importance of controlling the formulation factors during drug development, the potentiality of the optimized IFPs for the intramuscular delivery of piroxicam, and the reliability of PBPK physiological modeling in predicting the biological performance of new formulations with effective cost management.

A physiologically based model of bile acid metabolism in mice

Bile acid (BA) metabolism is a complex system that includes a wide variety of primary and secondary, as well as conjugated and unconjugated BAs that undergo continuous enterohepatic circulation (EHC). Alterations in both composition and dynamics of BAs have been associated with various diseases. However, a mechanistic understanding of the relationship between altered BA metabolism and related diseases is lacking. Computational modeling may support functional analyses of the physiological processes involved in the EHC of BAs along the gut-liver axis. In this study, we developed a physiologically based model of murine BA metabolism describing synthesis, hepatic and microbial transformations, systemic distribution, excretion, and EHC of BAs at the whole-body level. For model development, BA metabolism of specific pathogen-free (SPF) mice was characterized in vivo by measuring BA levels and composition in various organs, expression of transporters along the gut, and cecal microbiota composition. We found significantly different BA levels between male and female mice that could only be explained by adjusted expression of the hepatic enzymes and transporters in the model. Of note, this finding was in agreement with experimental observations. The model for SPF mice could also describe equivalent experimental data in germ-free mice by specifically switching off microbial activity in the intestine. The here presented model can therefore facilitate and guide functional analyses of BA metabolism in mice, e.g., the effect of pathophysiological alterations on BA metabolism and translation of results from mouse studies to a clinically relevant context through cross-species extrapolation.

Development and Validation of Physiologically Based Pharmacokinetic Model of Levetiracetam to Predict Exposure and Dose Optimization in Pediatrics

Levetiracetam (Lev) is an antiepileptic drug that has been increasingly used in the epilepsy pediatric population in recent years, but its pharmacokinetic behavior in pediatric population needs to be characterized clearly. Clinical trials for the pediatric drug remain difficult to conduct due to ethical and practical factors. The purpose of this study was to use the physiologically based pharmacokinetic (PBPK) model to predict changes in plasma exposure of Lev in pediatric patients and to provide recommendations for dose adjustment. A PBPK model of Lev in adults was developed using PK-Sim® software and extrapolated to the entire age range of the pediatric population. The model was evaluated using clinical pharmacokinetic data. The results showed the good fit between predictions and observations of the adult and pediatric models. The recommended doses for neonates, infants and children are 0.78, 1.67 and 1.22 times that of adults, respectively. Moreover, at the same dose, plasma exposure in adolescents was similar to that of adults. The PBPK models of Lev for adults and pediatrics were successfully developed and validated to provide a reference for the rational administration of drugs in the pediatric population.

Physiologically based pharmacokinetic (PBPK) modelling of oral drug absorption in older adults - an AGePOP review

The older population consisting of persons aged 65 years or older is the fastest-growing population group and also the major consumer of pharmaceutical products. Due to the heterogenous ageing process, this age group shows high interindividual variability in the dose-exposure-response relationship and, thus, a prediction of drug safety and efficacy is challenging. Although physiologically based pharmacokinetic (PBPK) modelling is a well-established tool to inform and confirm drug dosing strategies during drug development for special population groups, age-related changes in absorption are poorly accounted for in current PBPK models. The purpose of this review is to summarise the current state-of-knowledge in terms of physiological changes with increasing age that can influence the oral absorption of dosage forms. The capacity of common PBPK platforms to incorporate these changes and describe the older population is also discussed, as well as the implications of extrinsic factors such as drug-drug interactions associated with polypharmacy on the model development process. The future potential of this field will rely on addressing the gaps identified in this article, which can subsequently supplement in-vitro and in-vivo data for more robust decision-making on the adequacy of the formulation for use in older adults and inform pharmacotherapy.

The In Vitro, In Vivo, and PBPK Evaluation of a Novel Lung-Targeted Cardiac-Safe Hydroxychloroquine Inhalation Aerogel

Hydroxychloroquine (HCQ) was repurposed for COVID-19 treatment. Subtherapeutic HCQ lung levels and cardiac toxicity of oral HCQ were overcome by intratracheal (IT) administration of lower HCQ doses. The crosslinker-free supercritical fluid technology (SFT) produces aerogels and impregnates them with drugs in their amorphous form with efficient controlled release. Mechanistic physiologically based pharmacokinetic (PBPK) modeling can predict the lung's epithelial lining fluid (ELF) drug levels. This study aimed to develop a novel HCQ SFT formulation for IT administration to achieve maximal ELF levels and minimal cardiac toxicity. HCQ SFT formulation was prepared and evaluated for physicochemical, in vitro release, pharmacokinetics, and cardiac toxicity. Finally, the rat HCQ ELF concentrations were predicted using PBPK modeling. HCQ was amorphous after loading into the chitosan-alginate nanoporous microparticles (22.7±7.6 μm). The formulation showed a zero-order release, with only 40% released over 30 min compared to 94% for raw HCQ. The formulation had a tapped density of 0.28 g/cm3 and a loading efficiency of 35.3±1.3%. The IT administration of SFT HCQ at 1 mg/kg resulted in 23.7-fold higher bioavailability, fourfold longer MRT, and eightfold faster absorption but lower CK-MB and LDH levels than oral raw HCQ at 4 mg/kg. The PBPK model predicted 6 h of therapeutic ELF levels for IT SFT HCQ and a 100-fold higher ELF-to-heart concentration ratio than oral HCQ. Our findings support the feasibility of lung-targeted and more effective SFT HCQ IT administration for COVID-19 compared to oral HCQ with less cardiac toxicity. Graphical abstract.

Advanced In Vivo Prediction by Introducing Biphasic Dissolution Data into PBPK Models

Coupling biorelevant in vitro dissolution with in silico physiological-based pharmacokinetic (PBPK) tools represents a promising method to describe and predict the in vivo performance of drug candidates in formulation development including non-passive transport, prodrug activation, and first-pass metabolism. The objective of the present study was to assess the predictability of human pharmacokinetics by using biphasic dissolution results obtained with the previously established BiPHa+ assay and PBPK tools. For six commercial drug products, formulated by different enabling technologies, the respective organic partitioning profiles were processed with two PBPK in silico modeling tools, namely PK-Sim and GastroPlus^®, similar to extended-release dissolution profiles. Thus, a mechanistic dissolution/precipitation model of the assessed drug products was not required. The developed elimination/distribution models were used to simulate the pharmacokinetics of the evaluated drug products and compared with available human data. In essence, an in vitro to in vivo extrapolation (IVIVE) was successfully developed. Organic partitioning profiles obtained from the BiPHa+ dissolution analysis enabled highly accurate predictions of the pharmacokinetic behavior of the investigated drug products. In addition, PBPK models of (pro-)drugs with pronounced first-pass metabolism enabled adjustment of the solely passive diffusion predicting organic partitioning profiles, and increased prediction accuracy further.

Physiologically based pharmacokinetic modeling of levetiracetam to predict the exposure in hepatic and renal impairment and elderly populations

Levetiracetam (LEV) is an anti-epileptic drug approved for use in various populations. The pharmacokinetic (PK) behavior of LEV may be altered in the elderly and patients with renal and hepatic impairment. Thus, dosage adjustment is required. This study was conducted to investigate how the physiologically-based PK (PBPK) model describes the PKs of LEV in adult and elderly populations, as well as to predict the PKs of LEV in patients with renal and hepatic impairment in both populations. The whole-body PBPK models were developed using the reported physicochemical properties of LEV and clinical data. The models were validated using data from clinical studies with different dose ranges and different routes and intervals of administration. The fit performance of the models was assessed by comparing predicted and observed blood concentration data and PK parameters. It is recommended that the doses be reduced to ~70%, 60%, and 45% of the adult dose for the mild, moderate, and severe renal impairment populations and ~95%, 80%, and 57% of the adult dose for the Child Pugh-A (CP-A), Child Pugh-B (CP-B), and Child Pugh-C (CP-C) hepatic impairment populations, respectively. No dose adjustment is required for the healthy elderly population, but dose reduction is required for the elderly with organ dysfunction accordingly, on a scale similar to that of adults. A PBPK model of LEV was successfully developed to optimize dosing regimens for special populations.

Develop adult extrapolation to pediatrics and pediatric dose optimization based on the physiological pharmacokinetic model of azithromycin

Physiologically-based pharmacokinetic (PBPK) models are more frequently used for supporting pediatric dose selection in small-molecule drugs. Through literature research, drug parameters of azithromycin and clinical data from different studies were obtained. Through parameter optimization of the absorption and dissolution process, the adult intravenous model was extended to the adult oral model. The adult intravenous and oral PBPK models are precise to meet the AAFE<2 standard, and the pharmacokinetic parameters of the predicted values of the model are all within the mean standard deviation of the clinical observations. The values of plasma protein unbound fraction, renal clearance, and gastric juice pH between adults and pediatrics were changed by using the age-dependent pediatric organ maturity formula, and the adult model was extrapolated to the pediatric model. The final developed pediatric PBPK model was used to evaluate optimal dosing for children of different developmental ages. The relationship between the frist dose and age was as follows: 8.8 mg/kg/day from 0.5 to 2 years old, 9.2 mg/kg/day from 3 to 6 years old, 9.4 mg/kg/day from 7 to 12 years old, and 8.2 mg/kg/day from 13 to 18 years old, taken in half for 2-5 days. Simultaneously, the simulated exposures achieved with the dosing regimen proposed were comparable to adult plasma exposures for treatment of community-acquired pneumonia. A reasonable azithromycin pharmacokinetic-pharmacodynamic model for adults and pediatrics has been established, which can be demonstrated by the use of literature pediatric data to develop pediatric PBPK models, expanding the scope of this powerful modeling tool.

The Analysis of Pethidine Pharmacokinetics in Newborn Saliva, Plasma, and Brain Extracellular Fluid After Prenatal Intrauterine Exposure from Pregnant Mothers Receiving Intramuscular Dose Using PBPK Modeling

BACKGROUND AND OBJECTIVE

Pethidine (meperidine) can decrease labor pain-associated mother's hyperventilation and high cortisol-induced newborn complications. However, prenatal transplacentally acquired pethidine can cause side effects in newborns. High pethidine concentrations in the newborn brain extracellular fluid (bECF) can cause a serotonin crisis. Therapeutic drug monitoring (TDM) in newborns' blood distresses them and increases infection incidence, which can be overcome by using salivary TDM. Physiologically based pharmacokinetic (PBPK) modeling can predict drug concentrations in newborn plasma, saliva, and bECF after intrauterine pethidine exposure.

METHODS

A healthy adult PBPK model was constructed, verified, and scaled to newborn and pregnant populations after intravenous and intramuscular pethidine administration. The pregnancy PBPK model was used to predict the newborn dose received transplacentally at birth, which was used as input to the newborn PBPK model to predict newborn plasma, saliva, and bECF pethidine concentrations and set correlation equations between them.

RESULTS

Pethidine can be classified as a Salivary Excretion Classification System class II drug. The developed PBPK model predicted that, after maternal pethidine intramuscular doses of 100 mg and 150 mg, the newborn plasma and bECF concentrations were below the toxicity thresholds. Moreover, it was estimated that newborn saliva concentrations of 4.7 µM, 11.4 µM, and 57.7 µM can be used as salivary threshold concentrations for pethidine analgesic effects, side effects, and the risk for serotonin crisis, respectively, in newborns.

CONCLUSION

It was shown that saliva can be used for pethidine TDM in newborns during the first few days after delivery to mothers receiving pethidine.

Physiologically Based Biopharmaceutics Modeling of regional and colon absorption in humans

Colon absorption is a key determinant for successful development of extended release and colon targeted drug products. This is the first systematic evaluation of the ability to predict in vivo regional differences in absorption and the extent of colon absorption in humans using mechanistic physiologically based biopharmaceutics modeling (PBBM). A new dataset, consisting of 19 drugs with a wide range of biopharmaceutics properties and extent of colon absorption in humans, was established. Mechanistic predictions of the extent of absorption and plasma exposure after oral, or jejunal and direct colon administration were performed in GastroPlus and GI-Sim using an a priori approach. Two new colon models developed in GI-Sim, were also evaluated to assess if the prediction performance could be improved. Both GastroPlus and GI-Sim met the pre-defined criteria for accurate predictions of regional and colon absorption for high permeability drugs irrespective of formulation type, while the prediction performance was poor for low permeability drugs. For solutions, the two new GI-Sim colon models improved the colon absorption prediction performance for the low permeability drugs while maintaining the accurate prediction performance for the high permeability drugs. In contrast, the prediction performance decreased for non-solutions using the two new colon models. In conclusion, PBBM can be used with sufficient accuracy to predict regional and colon absorption in humans for high permeability drugs in candidate selection as well as early design and development of extended release or colon targeted drug products. The prediction performance of the current models needs to be improved to allow high accuracy predictions for commercial drug product applications including highly accurate predictions of the entire plasma concentration-time profiles as well as for low permeability drugs.

Physiologically Based Pharmacokinetic Modelling to Predict Pharmacokinetics of Enavogliflozin, a Sodium-Dependent Glucose Transporter 2 Inhibitor, in Humans

Enavogliflozin is a sodium-dependent glucose cotransporter 2 (SGLT2) inhibitor approved for clinical use in South Korea. As SGLT2 inhibitors are a treatment option for patients with diabetes, enavogliflozin is expected to be prescribed in various populations. Physiologically based pharmacokinetic (PBPK) modelling can rationally predict the concentration–time profiles under altered physiological conditions. In previous studies, one of the metabolites (M1) appeared to have a metabolic ratio between 0.20 and 0.25. In this study, PBPK models for enavogliflozin and M1 were developed using published clinical trial data. The PBPK model for enavogliflozin incorporated a non-linear urinary excretion in a mechanistically arranged kidney model and a non-linear formation of M1 in the liver. The PBPK model was evaluated, and the simulated pharmacokinetic characteristics were in a two-fold range from those of the observations. The pharmacokinetic parameters of enavogliflozin were predicted using the PBPK model under pathophysiological conditions. PBPK models for enavogliflozin and M1 were developed and validated, and they seemed useful for logical prediction.

The investigation of the complex population-drug-drug interaction between ritonavir-boosted lopinavir and chloroquine or ivermectin using physiologically-based pharmacokinetic modeling

OBJECTIVES

Therapy failure caused by complex population-drug-drug (PDDI) interactions including CYP3A4 can be predicted using mechanistic physiologically-based pharmacokinetic (PBPK) modeling. A synergy between ritonavir-boosted lopinavir (LPVr), ivermectin, and chloroquine was suggested to improve COVID-19 treatment. This work aimed to study the PDDI of the two CYP3A4 substrates (ivermectin and chloroquine) with LPVr in mild-to-moderate COVID-19 adults, geriatrics, and pregnancy populations.

METHODS

The PDDI of LPVr with ivermectin or chloroquine was investigated. Pearson's correlations between plasma, saliva, and lung interstitial fluid (ISF) levels were evaluated. Target site (lung epithelial lining fluid [ELF]) levels of ivermectin and chloroquine were estimated.

RESULTS

Upon LPVr coadministration, while the chloroquine plasma levels were reduced by 30, 40, and 20%, the ivermectin plasma levels were increased by a minimum of 425, 234, and 453% in adults, geriatrics, and pregnancy populations, respectively. The established correlation equations can be useful in therapeutic drug monitoring (TDM) and dosing regimen optimization.

CONCLUSIONS

Neither chloroquine nor ivermectin reached therapeutic ELF levels in the presence of LPVr despite reaching toxic ivermectin plasma levels. PBPK modeling, guided with TDM in saliva, can be advantageous to evaluate the probability of reaching therapeutic ELF levels in the presence of PDDI, especially in home-treated patients.

Experimental and computational models to investigate intestinal drug permeability and metabolism

Oral administration is the preferred route for drug administration that leads to better therapy compliance. The intestine plays a key role in the absorption and metabolism of oral drugs, therefore, new intestinal models are being continuously proposed, which contribute to the study of intestinal physiology, drug screening, drug side effects, and drug-drug interactions.Advances in pharmaceutical processes have produced more drug formulations, causing challenges for intestinal models. To adapt to the rapid evolution of pharmaceuticals, more intestinal models have been created. However, because of the complexity of the intestine, few models can take all aspects of the intestine into account, and some functions must be sacrificed to investigate other areas. Therefore, investigators need to choose appropriate models according to the experimental stage and other requirements to obtain the desired results.To help researchers achieve this goal, this review summarised the advantages and disadvantages of current commonly used intestinal models and discusses possible future directions, providing a better understanding of intestinal models.

A Physiologically-Based Pharmacokinetic Model of Ruxolitinib and Posaconazole to Predict CYP3A4-Mediated Drug-Drug Interaction Frequently Observed in Graft versus Host Disease Patients

Ruxolitinib (RUX) is approved for the treatment of steroid-refractory acute and chronic graft versus host disease (GvHD). It is predominantly metabolized via cytochrome P450 (CYP) 3A4. As patients with GvHD have an increased risk of invasive fungal infections, RUX is frequently combined with posaconazole (POS), a strong CYP3A4 inhibitor. Knowledge of RUX exposure under concomitant POS treatment is scarce and recommendations on dose modifications are inconsistent. A physiologically based pharmacokinetic (PBPK) model was developed to investigate the drug-drug interaction (DDI) between POS and RUX. The predicted RUX exposure was compared to observed concentrations in patients with GvHD in the clinical routine. PBPK models for RUX and POS were independently set up using PK-Sim^® Version 11. Plasma concentration-time profiles were described successfully and all predicted area under the curve (AUC) values were within 2-fold of the observed values. The increase in RUX exposure was predicted with a DDI ratio of 1.21 (C_max) and 1.59 (AUC). Standard dosing in patients with GvHD led to higher RUX exposure than expected, suggesting further dose reduction if combined with POS. The developed model can serve as a starting point for further simulations of the implemented DDI and can be extended to further perpetrators of CYP-mediated PK-DDIs or disease-specific physiological changes.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Lornoxicam: Exploration of doses for CYP2C9 Genotypes and Patients with Cirrhosis

Lornoxicam physiologically based pharmacokinetic (PBPK) models were developed and validated on the basis of clinical pharmacokinetic results obtained by considering CYP2C9 genetic polymorphisms in healthy adult populations. PBPK models were extended to predict lornoxicam pharmacokinetics for patients with cirrhosis by quantitatively examining the pathophysiological information associated with cirrhosis. The predicted plasma exposure to lornoxicam was approximately 1.12-2.83 times higher in the CYP2C9*1/*3 and *1/*13 groups than in the CYP2C9*1/*1 group of healthy adult populations and patients with cirrhosis. The predicted plasma exposure to lornoxicam was approximately 1.28-3.61 times higher in patients with cirrhosis than in healthy adult populations. If the relationship between lornoxicam exposure in plasma and drug efficacy was proportional, then the proposed adjusted doses of lornoxicam for each group varied from 1.25 mg to 8 mg. As the severity of cirrhosis increased, or when the CYP2C9 genotype was *1 heterozygous, the dose adjustment range of lornoxicam increased. Therefore, the effect of pathophysiological factors (cirrhosis severity) on the pharmacokinetics of lornoxicam might be more important than that of CYP2C9 genetic factors. These results could be useful for broadening the scope of clinical application of lornoxicam by enabling dose selection based on CYP2C9 genotypes and liver cirrhosis degree.

Predictive performance and verification of physiologically based pharmacokinetic model of propylthiouracil

Background: Propylthiouracil (PTU) treats hyperthyroidism and thyroid crisis in all age groups. A variety of serious adverse effects can occur during clinical use and require attention to its pharmacokinetic and pharmacodynamic characteristics in various populations.

Objective: To provide information for individualized dosing and clinical evaluation of PTU in the clinical setting by developing a physiologically based pharmacokinetic (PBPK) model, predicting ADME characteristics, and extrapolating to elderly and pediatric populations.

Methods: Relevant databases and literature were retrieved to collect PTU’s pharmacochemical properties and ADME parameters, etc. A PBPK model for adults was developed using PK-Sim® software to predict tissue distribution and extrapolated to elderly and pediatric populations. The mean fold error (MFE) method was used to compare the differences between predicted and observed values to assess the accuracy of the PBPK model. The model was validated using PTU pharmacokinetic data in healthy adult populations.

Result: The MFE ratios of predicted to observed values of AUC0-t, Cmax, and Tmax were mainly within 0.5 and 2. PTU concentrations in various tissues are lower than venous plasma concentrations. Compared to healthy adults, the pediatric population requires quantitative adjustment to the appropriate dose to achieve the same plasma exposure levels, while the elderly do not require dose adjustments.

Conclusion: The PBPK model of PTU was successfully developed, externally validated, and applied to tissue distribution prediction and special population extrapolation, which provides a reference for clinical individualized drug administration and evaluation.

Does the choice of applied physiologically‐based pharmacokinetics platform matter? A case study on simvastatin disposition and drug–drug interaction

Physiologically‐based pharmacokinetic (PBPK) models have an important role in drug discovery/development and decision making in regulatory submissions. This is facilitated by predefined PBPK platforms with user‐friendly graphical interface, such as Simcyp and PK‐Sim. However, evaluations of platform differences and the potential implications for disposition‐related applications are still lacking. The aim of this study was to assess how PBPK model development, input parameters, and model output are affected by the selection of PBPK platform. This is exemplified via the establishment of simvastatin PBPK models (workflow, final models, and output) in PK‐Sim and Simcyp as representatives of established whole‐body PBPK platforms. The major finding was that the choice of PBPK platform influenced the model development strategy and the final model input parameters, however, the predictive performance of the simvastatin models was still comparable between the platforms. The main differences between the structure and implementation of Simcyp and PK‐Sim were found in the absorption and distribution models. Both platforms predicted equally well the observed simvastatin (lactone and acid) pharmacokinetics (20–80 mg), BCRP and OATP1B1 drug–gene interactions (DGIs), and drug–drug interactions (DDIs) when co‐administered with CYP3A4 and OATP1B1 inhibitors/inducers. This study illustrates that in‐depth knowledge of established PBPK platforms is needed to enable an assessment of the consequences of PBPK platform selection. Specifically, this work provides insights on software differences and potential implications when bridging PBPK knowledge between Simcyp and PK‐Sim users. Finally, it provides a simvastatin model implemented in both platforms for risk assessment of metabolism‐ and transporter‐mediated DGIs and DDIs.

Development and Evaluation of a Physiologically Based Pharmacokinetic Model for Predicting Haloperidol Exposure in Healthy and Disease Populations

The physiologically based pharmacokinetic (PBPK) approach can be used to develop mathematical models for predicting the absorption, distribution, metabolism, and elimination (ADME) of administered drugs in virtual human populations. Haloperidol is a typical antipsychotic drug with a narrow therapeutic index and is commonly used in the management of several medical conditions, including psychotic disorders. Due to the large interindividual variability among patients taking haloperidol, it is very likely for them to experience either toxic or subtherapeutic effects. We intend to develop a haloperidol PBPK model for identifying the potential sources of pharmacokinetic (PK) variability after intravenous and oral administration by using the population-based simulator, PK-Sim. The model was initially developed and evaluated to predict the PK of haloperidol and its reduced metabolite in adult healthy population after intravenous and oral administration. After evaluating the developed PBPK model in healthy adults, it was used to predict haloperidol–rifampicin drug–drug interaction and was extended to tuberculosis patients. The model evaluation was performed using visual assessments, prediction error, and mean fold error of the ratio of the observed-to-predicted values of the PK parameters. The predicted PK values were in good agreement with the corresponding reported values. The effects of the pathophysiological changes and enzyme induction associated with tuberculosis and its treatment, respectively, on haloperidol PK, have been predicted precisely. For all clinical scenarios that were evaluated, the predicted values were within the acceptable two-fold error range.

Physiologically based pharmacokinetic (PBPK) modeling of flurbiprofen in different CYP2C9 genotypes

The aim of this study was to establish the physiologically based pharmacokinetic (PBPK) model of flurbiprofen related to CYP2C9 genetic polymorphism and describe the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes. PK-Sim® software was used for the model development and validation. A total of 16 clinical pharmacokinetic data for flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups were used for the PBPK modeling. Turnover number (k_cat) of CYP2C9 values were optimized to capture the observed profiles in different CYP2C9 genotypes. In the simulation, predicted fraction metabolized by CYP2C9, fraction excreted to urine, bioavailability, and volume of distribution were similar to previously reported values. Predicted plasma concentration-time profiles in different CYP2C9 genotypes were visually similar to the observed profiles. Predicted AUC_inf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.44-, 2.05-, and 3.67-fold higher than the CYP2C9*1/*1 genotype. The ranges of fold errors for AUC_inf, C_max, and t_1/2 were 0.84-1.00, 0.61-1.22, and 0.74-0.94 in development and 0.59-0.98, 0.52-0.97, and 0.61-1.52 in validation, respectively, which were within the acceptance criterion. Thus, the PBPK model was successfully established and described the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups. The present model could guide the decision-making of tailored drug administration strategy by predicting the pharmacokinetics of flurbiprofen in various clinical scenarios.

Non-clinical studies indicating lack of interactions between iron/calcium ions and linzagolix, an orally available GnRH antagonist

Linzagolix is an orally available gonadotropin-releasing hormone antagonist used to treat sex-hormone-dependent diseases in women. This study aimed to investigate drug-drug interactions between linzagolix and iron/calcium ions in the intended clinical setting by conducting pharmacokinetic studies in vitro and in rats.Insoluble precipitate formation with metal ions was evaluated by measuring linzagolix concentrations in four types of bio-relevant dissolution media (fasted/fed state simulated gastric fluid and fasted/fed state simulated gastric fluid version 2), and chelate complex formation with metal ions was evaluated by release of linzagolix from a cellulose membrane sac. In these in vitro studies, linzagolix showed no potential for insoluble precipitate formation under fasted/fed conditions and no chelate complex formation in the presence of metal ions.In rats, the plasma concentration-time profiles of linzagolix and iron ion were similar regardless of whether they were administered with or without ferrous sulfate and linzagolix choline at clinically relevant doses. Thus, linzagolix and iron ion had no effect on each other's absorption in vivo.In conclusion, linzagolix is unlikely to cause clinically relevant drug-drug interactions by chelating metal ions according to the results of in vitro and in vivo studies.

Physiologically based pharmacokinetic modelling to predict the pharmacokinetics of metoprolol in different CYP2D6 genotypes

Metoprolol, a selective β₁-adrenoreceptor blocking agent used in the treatment of hypertension, angina, and heart failure, is primarily metabolized by the CYP2D6 enzyme, which catalyzes α-hydroxylation and O-desmethylation. As CYP2D6 is genetically highly polymorphic and the enzymatic activity differs greatly depending on the presence of the mutant allele(s), the pharmacokinetic profile of metoprolol is highly variable depending on the genotype of CYP2D6. The aim of study was to develop the physiologically based pharmacokinetic (PBPK) model of metoprolol related to CYP2D6 genetic polymorphism for personalized therapy with metoprolol. For PBPK modelling, our previous pharmacogenomic data were used. To obtain kinetic parameters (K_m, V_max, and CL_int) of each genotype, the recombinant CYP enzyme of each genotype was incubated with metoprolol and metabolic rates were assayed. Based on these data, the PBPK model of metoprolol was developed and validated in different CYP2D6 genotypes using PK-Sim^® software. As a result, the input values for various parameters for the PBPK model were presented and the PBPK model successfully described the pharmacokinetics of metoprolol in each genotype group. The simulated values were within the acceptance criterion (99.998% confidence intervals) compared with observed values. The PBPK model developed in this study can be used for personalized pharmacotherapy with metoprolol in individuals of various races, ages, and CYP2D6 genotypes.

Physiologically based pharmacokinetic (PBPK) modeling of piroxicam with regard to CYP2C9 genetic polymorphism

Piroxicam is a non-steroidal anti-inflammatory drug used to alleviate symptoms of osteoarthritis and rheumatoid arthritis. CYP2C9 genetic polymorphism significantly influences the pharmacokinetics of piroxicam. The objective of this study was to develop and validate the piroxicam physiologically based pharmacokinetic (PBPK) model related to CYP2C9 genetic polymorphism. PK-Sim^® version 10.0 was used for the PBPK modeling. The PBPK model was evaluated by predicted and observed plasma concentration-time profiles, fold errors of predicted to observed pharmacokinetic parameters, and a goodness-of-fit plot. The turnover number (k_cat) of CYP2C9 was adjusted to capture the pharmacokinetics of piroxicam in different CYP2C9 genotypes. The population PBPK model overall accurately described and predicted the plasma concentration-time profiles in different CYP2C9 genotypes. In our simulations, predicted AUC_inf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.83-, 2.07-, and 6.43-fold higher than CYP2C9*1/*1 genotype, respectively. All fold error values for AUC, C_max, and t_1/2 were included in the acceptance criterion with the ranges of 0.57-1.59, 0.63-1.39, and 0.65-1.51, respectively. The range of fold error values for predicted versus observed plasma concentrations was 0.11-3.13. 93.9% of fold error values were within the two-fold range. Average fold error, absolute average fold error, and root mean square error were 0.93, 1.27, and 0.72, respectively. Our model accurately captured the pharmacokinetic alterations of piroxicam according to CYP2C9 genetic polymorphism.

Use of Physiologically-Based Kinetics Modelling to Reliably Predict Internal Concentrations of the UV Filter, Homosalate, After Repeated Oral and Topical Application

Ethical and legal considerations have led to increased use of non-animal methods to evaluate the safety of chemicals for human use. We describe the development and qualification of a physiologically-based kinetics (PBK) model for the cosmetic UV filter ingredient, homosalate, to support its safety without the need of generating further animal data. The intravenous (IV) rat PBK model, using PK-Sim®, was developed and validated using legacy in vivo data generated prior to the 2013 EU animal-testing ban. Input data included literature or predicted physicochemical and pharmacokinetic properties. The refined IV rat PBK model was subject to sensitivity analysis to identify homosalate-specific sensitive parameters impacting the prediction of Cmax (more sensitive than AUC(0-∞)). These were then considered, together with population modeling, to calculate the confidence interval (CI) 95% Cmax and AUC(0-∞). Final model parameters were established by visual inspection of the simulations and biological plausibility. The IV rat model was extrapolated to oral administration, and used to estimate internal exposures to doses tested in an oral repeated dose toxicity study. Next, a human PBK dermal model was developed using measured human in vitro ADME data and a module to represent the dermal route. Model performance was confirmed by comparing predicted and measured values from a US-FDA clinical trial (Identifier: NCT03582215, https://clinicaltrials.gov/). Final exposure estimations were obtained in a virtual population and considering the in vitro and input parameter uncertainty. This model was then used to estimate the Cmax and AUC(0-24 h) of homosalate according to consumer use in a sunscreen. The developed rat and human PBK models had a good biological basis and reproduced in vivo legacy rat and human clinical kinetics data. They also complied with the most recent WHO and OECD recommendations for assessing the confidence level. In conclusion, we have developed a PBK model which predicted reasonably well the internal exposure of homosalate according to different exposure scenarios with a medium to high level of confidence. In the absence of in vivo data, such human PBK models will be the heart of future completely non-animal risk assessments; therefore, valid approaches will be key in gaining their regulatory acceptance. Clinical Trial Registration: https://clinicaltrials.gov/, identifier, NCT03582215.

In-Silico Prediction of Oral Drug Bioavailability: A multi-boluses approach

This work focuses on a new mathematical model able to describe in a simple manner the intestinal physiology, in order to better study drug absorption and bioavailability. The aim of our model is to overcome the limitations of physiological pharmacokinetics models of the literature, introducing a different modelling approach. The core of the new proposed model is a Discrete-Continuous Approach (DCA): a sequence of boluses travels in the investigated portion of the intestine, in counter-current with blood that flows in continuous mode. No empirical equations are implemented in this model. Simulation results show an excellent correlation between the predicted and experimental concentration profile used to validate our model. Our new approach provides a simple tool, with a good reliability, to analyze a very complex phenomenon, using only few parameters.

Physiologically based pharmacokinetic modeling of candesartan related to CYP2C9 genetic polymorphism in adult and pediatric patients

Candesartan cilexetil is an angiotensin II receptor blocker and it is widely used to treat hypertension and heart failure. This drug is a prodrug that rapidly converts to candesartan after oral administration. Candesartan is metabolized by cytochrome P450 2C9 (CYP2C9) enzyme or uridine diphosphate glucurinosyltransferase 1A3, or excreted in an unchanged form through urine, biliary tract and feces. We investigated the effect of genetic polymorphism of CYP2C9 enzyme on drug pharmacokinetics using physiologically based pharmacokinetic (PBPK) modeling. In addition, by introducing the age and ethnicity into the model, we developed a model that can propose an appropriate dosage regimen taking into account the individual characteristics of each patient. To evaluate the suitability of the model, the results of a clinical trial on twenty-two healthy Korean subjects and their CYP2C9 genetic polymorphism data was applied. In this study, PK-Sim® was used to develop the PBPK model of candesartan.

Physiologically based pharmacokinetic (PBPK) modeling of meloxicam in different CYP2C9 genotypes

Meloxicam, a non-steroidal anti-inflammatory drug, is used for the treatment of rheumatoid arthritis and osteoarthritis. Cytochrome P450 (CYP) 2C9 and CYP3A4 are major and minor enzymes involved in the metabolism of meloxicam. Impaired enzyme activity of CYP2C9 variants increases the plasma exposures of meloxicam and the risk of adverse events. The objective of our study is to develop and validate the physiologically based pharmacokinetic (PBPK) model of meloxicam related to CYP2C9 genetic polymorphism using the PK-Sim^® software. In vitro k_cat of CYP2C9 was optimized in different CYP2C9 genotypes. The demographic and pharmacokinetic dataset for the development of the PBPK model was extracted from two previous clinical pharmacokinetic studies. Thirty-one clinical datasets, representing different dose regimens and demographic characteristics, were utilized to validate the PBPK model. The shapes of simulated plasma concentration-time profiles in each CYP2C9 genotype were visually similar to observed profiles. The predicted exposures (AUC_inf) of meloxicam in CYP2C9*1/*3, CYP2C9*1/*13, and CYP2C9*3/*3 genotypes were increased by 1.77-, 2.91-, and 8.35-fold compared to CYP2C9*1/*1 genotype, respectively. In all datasets for the development and validations, fold errors between predicted and observed pharmacokinetic parameters were within the two-fold error criteria. As a result, the PBPK model was appropriately established and properly described the pharmacokinetics of meloxicam in different CYP2C9 genotypes. This study is expected to contribute to reducing the risk of adverse events of meloxicam through optimization of meloxicam dosing in different CYP2C9 genotypes.

A Physiological Approach to Pharmacokinetics in Chronic Kidney Disease

The conventional approach to approximating the pharmacokinetics of drugs in patients with chronic kidney disease (CKD) only accounts for changes in the estimated glomerular filtration rate. However, CKD is a systemic and multifaceted disease that alters many body systems. Therefore, the objective of this exercise was to develop and evaluate a whole-body mechanistic approach to predicting pharmacokinetics in patients with CKD. Physiologically based pharmacokinetic models were developed in PK-Sim v8.0 (www.open-systems-pharmacology.org) to mechanistically represent the disposition of 7 compounds in healthy human adults. The 7 compounds selected were eliminated by glomerular filtration and active tubular secretion by the organic cation transport system to varying degrees. After a literature search, the healthy adult models were adapted to patients with CKD by numerically accounting for changes in glomerular filtration rate, kidney volume, renal perfusion, hematocrit, plasma protein concentrations, and gastrointestinal transit. Literature-informed interindividual variability was applied to the physiological parameters to facilitate a population approach. Model performance in CKD was evaluated against pharmacokinetic data from 8 clinical trials in the literature. Overall, integration of the CKD parameterization enabled exposure predictions that were within 1.5-fold error across all compounds and patients with varying stages of renal impairment. Notable improvement was observed over the conventional approach to scaling exposure, which failed in all but 1 scenario in patients with advanced CKD. Further research is required to qualify its use for first-in-CKD dose selection and clinical trial planning for a wider selection of renally eliminated compounds, including those subject to anion transport.

PBPK Modeling as a Tool for Predicting and Understanding Intestinal Metabolism of Uridine 5'-Diphospho-glucuronosyltransferase Substrates

Uridine 5′-diphospho-glucuronosyltransferases (UGTs) are expressed in the small intestines, but prediction of first-pass extraction from the related metabolism is not well studied. This work assesses physiologically based pharmacokinetic (PBPK) modeling as a tool for predicting intestinal metabolism due to UGTs in the human gastrointestinal tract. Available data for intestinal UGT expression levels and in vitro approaches that can be used to predict intestinal metabolism of UGT substrates are reviewed. Human PBPK models for UGT substrates with varying extents of UGT-mediated intestinal metabolism (lorazepam, oxazepam, naloxone, zidovudine, cabotegravir, raltegravir, and dolutegravir) have demonstrated utility for predicting the extent of intestinal metabolism. Drug–drug interactions (DDIs) of UGT1A1 substrates dolutegravir and raltegravir with UGT1A1 inhibitor atazanavir have been simulated, and the role of intestinal metabolism in these clinical DDIs examined. Utility of an in silico tool for predicting substrate specificity for UGTs is discussed. Improved in vitro tools to study metabolism for UGT compounds, such as coculture models for low clearance compounds and better understanding of optimal conditions for in vitro studies, may provide an opportunity for improved in vitro–in vivo extrapolation (IVIVE) and prospective predictions. PBPK modeling shows promise as a useful tool for predicting intestinal metabolism for UGT substrates.

Physiologically based pharmacokinetic (PBPK) modeling for prediction of celecoxib pharmacokinetics according to CYP2C9 genetic polymorphism

Celecoxib is a non-steroidal anti-inflammatory drug (NSAID) and a representative selective cyclooxygenase (COX)-2 inhibitor, which is commonly prescribed for osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, acute pain, and primary dysmenorrhea. It is mainly metabolized by CYP2C9 and partly by CYP3A4 after oral administration. Many studies reported that CYP2C9 genetic polymorphism has significant effects on the pharmacokinetics of celecoxib and the occurrence of adverse drug reactions. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model of celecoxib according to CYP2C9 genetic polymorphism for personalized pharmacotherapy. Initially, a clinical pharmacokinetic study was conducted where a single dose (200 mg) of celecoxib was administered to 39 healthy Korean subjects with CYP2C9*1/*1 or CYP2C9*1/*3 genotypes to obtain data for PBPK development. Based on the conducted pharmacokinetic study and a previous pharmacokinetic study involving subjects with CYP2C9*1/*13 and CYP2C9*3/*3 genotype, PBPK model for celecoxib was developed. A PBPK model for CYP2C9*1/*1 genotype group was developed and then scaled to other genotype groups (CYP2C9*1/*3, CYP2C9*1/*13 and CYP2C9*3/*3). After model development, model validation was performed with comparison of five pharmacokinetic studies. As a result, the developed PBPK model of celecoxib successfully described the pharmacokinetics of each CYP2C9 genotype group and its predicted values were within the acceptance criterion. Additionally, all the predicted values were within two-fold error range in comparison to the previous pharmacokinetic studies. This study demonstrates the possibility of determining the appropriate dosage of celecoxib for each individual through the PBPK modeling with CYP2C9 genomic information. This approach could contribute to the reduction of adverse drug reactions of celecoxib and enable precision medicine.

Evaluation of Physiologically Based Pharmacokinetic Models to Predict the Absorption of BCS Class I Drugs in Different Pediatric Age Groups

Age-related changes in many parameters affecting drug absorption remain poorly characterized. The objective of this study was to apply physiologically based pharmacokinetic (PBPK) models in pediatric patients to investigate the absorption and pharmacokinetics of 4 drugs belonging to the Biopharmaceutics Classification System (BCS) class I administered as oral liquid formulations. Pediatric PBPK models built with PK-Sim/MoBi were used to predict the pharmacokinetics of acetaminophen, emtricitabine, theophylline, and zolpidem in different pediatric populations. The model performance for predicting drug absorption and pharmacokinetics was assessed by comparing the predicted absorption profile with the deconvoluted dose fraction absorbed over time and predicted with observed plasma concentration-time profiles. Sensitivity analyses were performed to analyze the effects of changes in relevant input parameters on the model output. Overall, most pharmacokinetic parameters were predicted within a 2-fold error range. The absorption profiles were generally reasonably predicted, but relatively large differences were observed for acetaminophen. Sensitivity analyses showed that the predicted absorption profile was most sensitive to changes in the gastric emptying time (GET) and the specific intestinal permeability. The drug's solubility played only a minor role. These findings confirm that gastric emptying time, more than intestinal permeability or solubility, is a key factor affecting BCS class I drug absorption in children. As gastric emptying time is prolonged in the fed state, a better understanding of the interplay between food intake and gastric emptying time in children is needed, especially in the very young in whom the (semi)fed condition is the prevailing prandial state, and hence prolonged gastric emptying time seems more plausible than the fasting state.

Predictive Performance of Physiology-Based Pharmacokinetic Dose Estimates for Pediatric Trials: Evaluation With 10 Bayer Small-Molecule Compounds in Children

Development and guidance of dosing schemes in children have been supported by physiology-based pharmacokinetic (PBPK) modeling for many years. PBPK models are built on a generic basis, where compound- and system-specific parameters are separated and can be exchanged, allowing the translation of these models from adults to children by accounting for physiological differences. Owing to these features, PBPK modeling is a valuable approach to support clinical decision making for dosing in children. In this analysis, we evaluate pediatric PBPK models for 10 small-molecule compounds that were applied to support clinical decision processes at Bayer for their predictive power in different age groups. Ratios of PBPK-predicted to observed PK parameters for the evaluated drugs in different pediatric age groups were estimated. Predictive performance was analyzed on the basis of a 2-fold error range and the bioequivalence range (ie, 0.8 ≤ predicted/observed ≤ 1.25). For all 10 compounds, all predicted-to-observed PK ratios were within a 2-fold error range (n = 27), with two-thirds of the ratios within the bioequivalence range (n = 18). The findings demonstrate that the pharmacokinetics of these compounds was successfully and adequately predicted in different pediatric age groups. This illustrates the applicability of PBPK for guiding dosing schemes in the pediatric population.

The impact of the intestinal microbiota and the mucosal permeability on three different antibiotic drugs

BackgroundThe totality of bacteria, protozoa, viruses and fungi that lives in the human body is called microbiota. Human microbiota specifically colonizes the skin, the respiratory and urinary tract, the urogenital tract and the gastrointestinal system. This study focuses on the intestinal microbiota to explore the drug-microbiota relationship and, therefore, how the drug bioavailability changes in relation to the microbiota biodiversity to identify more personalized therapies, with the minimum risk of side effects. MethodsTo achieve this goal, we developed a new mathematical model with two compartments, the intestine and the blood, which takes into account the colonic mucosal permeability variation - measured by Ussing chamber system on human colonic mucosal biopsies - and the fecal microbiota composition, determined through microbiota 16S rRNA sequencing analysis. Both of the clinical parameters were evaluated in a group of Irritable Bowel Syndrome patients compared to a group of healthy controls. Key ResultsThe results show that plasma drug concentration increases as bacterial concentration decreases, while it decreases as intestinal length decreases too. ConclusionsThe study provides interesting data since in literature there are not yet mathematical models with these features, in which the importance of intestinal microbiota, the "forgotten organ", is considered both for the subject health state and in the nutrients and drugs metabolism.

A Physiologically-Based Pharmacokinetic Framework for Prediction of Drug Exposure in Malnourished Children

Malnutrition in children is a global health problem, particularly in developing countries. The effects of an insufficient supply of nutrients on body composition and physiological functions may have implications for drug disposition and ultimately affect the clinical outcome in this vulnerable population. Physiologically-based pharmacokinetic (PBPK) modeling can be used to predict the effect of malnutrition as it links physiological changes to pharmacokinetic (PK) consequences. However, the absence of detailed information on body composition and the limited availability of controlled clinical trials in malnourished children complicates the establishment and evaluation of a generic PBPK model in this population. In this manuscript we describe the creation of physiologically-based bridge to a malnourished pediatric population, by combining information on (a) the differences in body composition between healthy and malnourished adults and (b) the differences in physiology between healthy adults and children. Model performance was confirmed using clinical reference data. This study presents a physiologically-based translational framework for prediction of drug disposition in malnourished children. The model is readily applicable for dose recommendation strategies to address the urgent medicinal needs of this vulnerable population.

Physiologically Based Pharmacokinetic Modeling of Metoprolol Enantiomers and α-Hydroxymetoprolol to Describe CYP2D6 Drug-Gene Interactions

The beta-blocker metoprolol (the sixth most commonly prescribed drug in the USA in 2017) is subject to considerable drug–gene interaction (DGI) effects caused by genetic variations of the CYP2D6 gene. CYP2D6 poor metabolizers (5.7% of US population) show approximately five-fold higher metoprolol exposure compared to CYP2D6 normal metabolizers. This study aimed to develop a whole-body physiologically based pharmacokinetic (PBPK) model to predict CYP2D6 DGIs with metoprolol. The metoprolol (R)- and (S)-enantiomers as well as the active metabolite α-hydroxymetoprolol were implemented as model compounds, employing data of 48 different clinical studies (dosing range 5–200 mg). To mechanistically describe the effect of CYP2D6 polymorphisms, two separate metabolic CYP2D6 pathways (α-hydroxylation and O-demethylation) were incorporated for both metoprolol enantiomers. The good model performance is demonstrated in predicted plasma concentration–time profiles compared to observed data, goodness-of-fit plots, and low geometric mean fold errors of the predicted AUC_last (1.27) and C_max values (1.23) over all studies. For DGI predictions, 18 out of 18 DGI AUC_last ratios and 18 out of 18 DGI C_max ratios were within two-fold of the observed ratios. The newly developed and carefully validated model was applied to calculate dose recommendations for CYP2D6 polymorphic patients and will be freely available in the Open Systems Pharmacology repository.

Evaluation of Voriconazole CYP2C19 Phenotype-Guided Dose Adjustments by Physiologically Based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVES

Controversy exists regarding dose adjustment in patients treated with voriconazole due to the severity of the infections for which it is prescribed. The Dutch Pharmacogenetics Working Group (DPWG) recommends a 50% dose increase or decrease for cytochrome P450 (CYP) 2C19 ultrarapid (UM) or poor (PM) metabolizers, respectively. In contrast, for the previous phenotypes, the Clinical Pharmacogenetics Implementation Consortium (CPIC) voriconazole guideline only recommends a change of treatment. Based on observed data from single-dose bioequivalence studies and steady-state observed concentrations, we aimed to investigate voriconazole dose adjustments by means of physiologically based pharmacokinetic (PBPK) modeling.

METHODS

PBPK modeling was used to optimize voriconazole single-dose models for each CYP2C19 phenotype, which were extrapolated to steady state and evaluated for concordance with the therapeutic range of voriconazole. Based on optimized models, dose adjustments were evaluated for better adjustment to the therapeutic range.

RESULTS

Our models suggest that the standard dose may only be appropriate for normal metabolizers (NM), although they would benefit from a 50-100% loading dose increase. Intermediate metabolizers (IMs) and PMs required a daily dose reduction of 50 and 75%, respectively. Rapid metabolizers (RMs) and UMs required a daily dose increase of 100% and 300%, respectively.

CONCLUSION

The prescription of voriconazole in clinical practice should be personalized according to the CYP2C19 phenotype, followed by therapeutic drug monitoring of plasma concentrations to guide dose adjustment.

Ethambutol disposition in humans: Challenges and limitations of whole-body physiologically-based pharmacokinetic modelling in early drug development

Whole-body physiologically based pharmacokinetic (WB-PBPK) models have become an important tool in drug development, as they enable characterization of pharmacokinetic profiles across different organs based on physiological (systems-specific) and physicochemical (drug-specific) properties. However, it remains unclear which data are needed for accurate predictions when applying the approach to novel candidate molecules progressing into the clinic. In this work, as case study, we investigated the predictive performance of WB-PBPK models both for prospective and retrospective evaluation of the pharmacokinetics of ethambutol, considering scenarios that reflect different stages of development, including settings in which the data are limited to in vitro experiments, in vivo preclinical data, and when some clinical data are available. Overall, the accuracy of PBPK model-predicted systemic and tissue exposure was heavily dependant on prior knowledge about the eliminating organs. Whilst these findings may be specific to ethambutol, the challenges and potential limitations identified here may be relevant to a variety of drugs, raising questions about (1) the minimum requirements for prospective use of WB-PBPK models during the characterization of drug disposition and (2) implication of uncertainty for dose selection in humans.

Insights from pharmacokinetic models of host-microbiome drug metabolism

Increasing evidence suggests a role of the gut microbiota in patients' response to medicinal drugs. In our recent study, we combined genomics of human gut commensals and gnotobiotic animal experiments to quantify microbiota and host contributions to drug metabolism. Informed by experimental data, we built a physiology-based pharmacokinetic model of drug metabolism that includes intestinal compartments with microbiome drug-metabolizing activity. This model successfully predicted serum levels of metabolites of three different drugs, quantified microbial contribution to systemic drug metabolite exposure, and simulated the effect of different parameters on host and microbiota drug metabolism. In this addendum, we expand these simulations to assess the effect of microbiota on the systemic drug and metabolite levels under conditions of altered host physiology, microbiota drug-metabolizing activity or physico-chemical properties of drugs. This work illustrates how and under which circumstances the gut microbiome may influence drug pharmacokinetics, and discusses broader implications of expanded pharmacokinetic models.

A multiscale absorption and transit model for oral drug delivery: Formulation and applications during fasting conditions

Most Food and Drug Administration (FDA)-approved drugs are administered orally, despite the complex process of oral drug absorption that is difficult to analyze experimentally. Oral bioavailability is dependent on the drug compound as well as the physiological and anatomical states of the user. Thus, computational models have emerged to mechanistically capture and predict the oral absorption process. The current models are generally 0D compartmental models and are limited by (a) simplified physiological characteristics of the gastrointestinal tract (GIT), (b) semiempirical/analytical dissolution profiles of the tested drugs, (c) incorrect absorption for some drug BCS classes (class IIa, for example), (d) GITs size variability among population, (e) incorrectly predicting the absorption of drugs that are GIT target specific, and (f) erroneous mixing in the domain. In this study, we have developed a multiscale absorption and transit (MAT) toolkit to simulate the dissolution, transport, absorption, distribution, metabolism, and elimination of orally administered drugs in the human GIT at multiple levels. MAT was constructed by integrating the spatially accurate first-principles driven high-fidelity drug transport, dissolution, and absorption model in the human stomach and GIT using our recently published quasi-3D (Q3D) framework. The process integrated the multilayer intestine physiologically based pharmacokinetics models with the whole-body compartmental models to predict the systemic pharmacokinetics of oral drugs. The computational results showed that this multiscale tool was able to match the experimental concentration results (individual and population) better than the traditional compartmental models. Ultimately, MAT will be developed into a commercial product to meet urgent demands from pharmaceutical and biomedical industries.

Physiologically Based Pharmacokinetic Models to Predict Maternal Pharmacokinetics and Fetal Exposure to Emtricitabine and Acyclovir

Pregnancy is associated with physiological changes that may impact drug pharmacokinetics (PK). The goals of this study were to build maternal-fetal physiologically based pharmacokinetic (PBPK) models for acyclovir and emtricitabine, 2 anti(retro)viral drugs with active renal net secretion, and to (1) evaluate the predicted maternal PK at different stages of pregnancy; (2) predict the changes in PK target parameters following the current dosing regimen of these drugs throughout pregnancy; (3) evaluate the predicted concentrations of these drugs in the umbilical vein at delivery; (4) compare the model performance for predicting maternal PK of emtricitabine in the third trimester with that of previously published PBPK models; and (5) compare different previously published approaches for estimating the placental permeability of these 2 drugs. Results showed that the pregnancy PBPK model for acyclovir predicted all maternal concentrations within a 2-fold error range, whereas the model for emtricitabine predicted 79% of the maternal concentrations values within that range. Extrapolation of these models to earlier stages of pregnancy indicated that the change in the median PK target parameters remained well above the target threshold. Concentrations of acyclovir and emtricitabine in the umbilical vein were overall adequately predicted. The comparison of different emtricitabine PBPK models suggested an overall similar predictive performance in the third trimester, but the comparison of different approaches for estimating placental drug permeability revealed large differences. These models can enhance the understanding of the PK behavior of renally excreted drugs, which may ultimately inform pharmacotherapeutic decision making in pregnant women and their fetuses.

Combination of Gemcitabine with Cell-Penetrating Peptides: A Pharmacokinetic Approach Using In Silico Tools

Gemcitabine is an anticancer drug used to treat a wide range of solid tumors and is a first line treatment for pancreatic cancer. Our group has previously developed novel conjugates of gemcitabine with cell-penetrating peptides (CPP), and here we report some preliminary data regarding the pharmacokinetics of gemcitabine, two gemcitabine-CPP conjugates and respective CPP gathered from GastroPlus™, and analyze these results considering our previous evaluation of gemcitabine release and conjugates’ bioactivity. Additionally, seeking to shed some light on the relation between the penetration ability of CPP and their physicochemical properties, chemical descriptors for the 20 natural amino acids were calculated, a new principal property scale (z-scale) was created and CPP prediction models were developed, establishing quantitative structure-activity relationships (QSAR). The z-scores of the peptides conjugated with gemcitabine are presented and analyzed with the aforementioned data.

A Physiologically-Based Pharmacokinetic Model to Describe Ciprofloxacin Pharmacokinetics Over the Entire Span of Life

Background

Physiologically-based pharmacokinetic (PBPK) modeling has received growing interest as a useful tool for the assessment of drug pharmacokinetics by continuous knowledge integration.

Objective

The objective of this study was to build a ciprofloxacin PBPK model for intravenous and oral dosing based on a comprehensive literature review, and evaluate the predictive performance towards pediatric and geriatric patients.

Methods

The aim of this report was to establish confidence in simulations of the ciprofloxacin PBPK model along the development process to facilitate reliable predictions outside of the tested adult age range towards the extremes of ages. Therefore, mean data of 69 published clinical trials were identified and integrated into the model building, simulation and verification process. The predictive performance on both ends of the age scale was assessed using individual data of 258 subjects observed in own clinical trials.

Results

Ciprofloxacin model verification demonstrated no concentration-related bias and accurate simulations for the adult age range, with only 4.8% of the mean observed data points for intravenous administration and 12.1% for oral administration being outside the simulated twofold range. Predictions towards the extremes of ages for the area under the plasma concentration–time curve (AUC) and the maximum plasma concentration (C_max) over the entire span of life revealed a reliable estimation, with only two pediatric AUC observations outside the 90% prediction interval.

Conclusion

Overall, this ciprofloxacin PBPK modeling approach demonstrated the predictive power of a thoroughly informed middle-out approach towards age groups of interest to potentially support the decision-making process.

Electronic supplementary material

The online version of this article (10.1007/s40262-018-0661-6) contains supplementary material, which is available to authorized users.

A Physiology-Based Model of Human Bile Acid Metabolism for Predicting Bile Acid Tissue Levels After Drug Administration in Healthy Subjects and BRIC Type 2 Patients

Drug-induced liver injury (DILI) is a matter of concern in the course of drug development and patient safety, often leading to discontinuation of drug-development programs or early withdrawal of drugs from market. Hepatocellular toxicity or impairment of bile acid (BA) metabolism, known as cholestasis, are the two clinical forms of DILI. Whole-body physiology-based modelling allows a mechanistic investigation of the physiological processes leading to cholestasis in man. Objectives of the present study were: (1) the development of a physiology-based model of the human BA metabolism, (2) population-based model validation and characterisation, and (3) the prediction and quantification of altered BA levels in special genotype subgroups and after drug administration. The developed physiology-based bile acid (PBBA) model describes the systemic BA circulation in humans and includes mechanistically relevant active and passive processes such as the hepatic synthesis, gallbladder emptying, transition through the gastrointestinal tract, reabsorption into the liver, distribution within the whole body, and excretion via urine and faeces. The kinetics of active processes were determined for the exemplary BA glycochenodeoxycholic acid (GCDCA) based on blood plasma concentration-time profiles. The robustness of our PBBA model was verified with population simulations of healthy individuals. In addition to plasma levels, the possibility to estimate BA concentrations in relevant tissues like the intracellular space of the liver enhance the mechanistic understanding of cholestasis. We analysed BA levels in various tissues of Benign Recurrent Intrahepatic Cholestasis type 2 (BRIC2) patients and our simulations suggest a higher susceptibility of BRIC2 patients toward cholestatic DILI due to BA accumulation in the liver. The effect of drugs on systemic BA levels were simulated for cyclosporine A (CsA). Our results confirmed the higher risk of DILI after CsA administration in healthy and BRIC2 patients. The presented PBBA model enhances our mechanistic understanding underlying cholestasis and drug-induced alterations of BA levels in blood and organs. The developed PBBA model might be applied in the future to anticipate potential risk of cholestasis in patients.