Novel Mechanistic PBPK Model to Predict Renal Clearance in Varying Stages of CKD by Incorporating Tubular Adaptation and Dynamic Passive Reabsorption

3D human tissue models and microphysiological systems for HIV and related comorbidities

Three-dimensional (3D) human tissue models/microphysiological systems (e.g., organs-on-chips, organoids, and tissue explants) model HIV and related comorbidities and have potential to address critical questions, including characterization of viral reservoirs, insufficient innate and adaptive immune responses, biomarker discovery and evaluation, medical complexity with comorbidities (e.g., tuberculosis and SARS-CoV-2), and protection and transmission during pregnancy and birth. Composed of multiple primary or stem cell-derived cell types organized in a dedicated 3D space, these systems hold unique promise for better reproducing human physiology, advancing therapeutic development, and bridging the human-animal model translational gap. Here, we discuss the promises and achievements with 3D human tissue models in HIV and comorbidity research, along with remaining barriers with respect to cell biology, virology, immunology, and regulatory issues.

Maximum likelihood estimation of renal transporter ontogeny profiles for pediatric PBPK modeling

Optimal treatment of infants with many renally cleared drugs must account for maturational differences in renal transporter (RT) activity. Pediatric physiologically-based pharmacokinetic (PBPK) models may incorporate RT activity, but this requires ontogeny profiles for RT activity in children, especially neonates, to predict drug disposition. Therefore, RT expression measurements from human kidney postmortem cortical tissue samples were normalized to represent a fraction of mature RT activity. Using these data, maximum likelihood estimated the distributions of RT activity across the pediatric age spectrum, including preterm and term neonates. PBPK models of four RT substrates (acyclovir, ciprofloxacin, furosemide, and meropenem) were evaluated with and without ontogeny profiles using average fold error (AFE), absolute average fold error (AAFE), and proportion of observations within the 5-95% prediction interval. Novel maximum likelihood profiles estimated ontogeny distributions for the following RT: OAT1, OAT3, OCT2, P-gp, URAT1, BCRP, MATE1, MRP2, MRP4, and MATE-2 K. Profiles for OAT3, P-gp, and MATE1 improved infant furosemide and neonate meropenem PBPK model AFE from 0.08 to 0.70 and 0.53 to 1.34 and model AAFE from 12.08 to 1.44 and 2.09 to 1.36, respectively, and improved the percent of data within the 5-95% prediction interval from 48% to 98% for neonatal ciprofloxacin simulations, respectively. Even after accounting for other critical population-specific maturational differences, novel RT ontogeny profiles substantially improved neonatal PBPK model performance, providing validated estimates of maturational differences in RT activity for optimal dosing in children.

Predicting complex kidney drug handling using a physiologically-based pharmacokinetic model informed by biomarker-estimated secretory clearance and blood flow

Kidney function-adjusted drug dosing is currently based solely on the estimated glomerular filtration rate (GFR), however, kidney drug handling is accomplished by a combination of filtration, tubular secretion, and re-absorption. Mechanistic physiologically-based pharmacokinetic (PBPK) models recapitulate anatomic compartments to predict elimination from estimated perfusion, filtration, secretion, and re-absorption, but clinical applications are limited by a lack of empiric individual-level measurements of these functions. We adapted and validated a PBPK model to predict drug clearance from individual biomarker-based estimates of kidney perfusion and secretory clearance. We estimated organic anion transporter-mediated secretion via kynurenic acid clearance and kidney blood flow (KBF) via isovalerylglycine clearance in human participants, incorporating these measurements with GFR into the model to predict kidney drug clearance. We compared measured and model-predicted clearances of administered tenofovir and oseltamivir, which are cleared by both filtration and secretion. There were 27 outpatients (age 55 ± 15 years, mean iohexol-GFR [iGFR] 76 ± 31 mL/min/1.73 m2 ) in this drug clearance study. The mean observed and mechanistic model-predicted tenofovir clearances were 169 ± 102 mL/min and 163 ± 80 mL/min, respectively; estimated mean error of the mechanistic model was 37.1 mL/min (95% confidence interval [CI]: 24-52.9), compared to a mean error of 41.8 mL/min (95% CI: 25-61.6) from regression model. The mean observed and model-predicted oseltamivir carboxylate clearances were 183 ± 104 mL/min and 179 ± 89 mL/min, respectively; estimated mean error of the mechanistic model was 42.9 mL/min (95% CI: 29.7-56.4), versus error of 48.1 mL/min (95% CI: 31.2-67.3) from the regression model. Individualized estimates of tubular secretion and KBF improved the accuracy of PBPK model-predicted tenofovir and oseltamivir kidney clearances, suggesting the potential for biomarker-informed measures of kidney function to refine personalized drug dosing.

Chronic kidney disease and physiologically based pharmacokinetic modeling: a critical review of existing models

INTRODUCTION

Physiologically based pharmacokinetic (PBPK) modeling is a paradigm shift in this era for determining the exposure of drugs in pediatrics, geriatrics, and patients with chronic diseases where clinical trials are difficult to conduct.

AREAS COVERED

This review has collated data regarding published PBPK models on chronic kidney disease (CKD), including the drug and system-specific input model parameters and model evaluation criteria. Four databases were used from 13th June 2023 to 10th July 2023 for identifying the relevant studies that met the inclusion/exclusion criteria. Alterations in plasma protein (albumin/alpha-1 acid glycoprotein), gastric emptying time, hematocrit, small intestinal transit time, the abundance of cytochrome (CYP) 450 enzymes, glomerular filtration rate, and physicochemical parameters for different drugs were explicitly elaborated from earlier reported studies. Moreover, model evaluation depicted that models in CKD for most of the included drugs were within the allowed two-fold error range.

EXPERT OPINION

This review will provide insights for researchers on applying PBPK models in managing patients with different levels of CKD to prevent undesirable side effects and increase the effectiveness of drug therapy.

Challenges Related to Acquisition of Physiological Data for Physiologically Based Pharmacokinetic (PBPK) Models in Postpartum, Lactating Women and Breastfed Infants-A Contribution from the ConcePTION Project

Physiologically based pharmacokinetic (PBPK) modelling is a bottom-up approach to predict pharmacokinetics in specific populations based on population-specific and medicine-specific data. Using an illustrative approach, this review aims to highlight the challenges of incorporating physiological data to develop postpartum, lactating women and breastfed infant PBPK models. For instance, most women retain pregnancy weight during the postpartum period, especially after excessive gestational weight gain, while breastfeeding might be associated with lower postpartum weight retention and long-term weight control. Based on a structured search, an equation for human milk intake reported the maximum intake of 153 mL/kg/day in exclusively breastfed infants at 20 days, which correlates with a high risk for medicine reactions at 2-4 weeks in breastfed infants. Furthermore, the changing composition of human milk and its enzymatic activities could affect pharmacokinetics in breastfed infants. Growth in breastfed infants is slower and gastric emptying faster than in formula-fed infants, while a slower maturation of specific metabolizing enzymes in breastfed infants has been described. The currently available PBPK models for these populations lack structured systematic acquisition of population-specific data. Future directions include systematic searches to fully identify physiological data. Following data integration as mathematical equations, this holds the promise to improve postpartum, lactation and infant PBPK models.

The Trends and Future Prospective of In Silico Models from the Viewpoint of ADME Evaluation in Drug Discovery

Drug discovery and development are aimed at identifying new chemical molecular entities (NCEs) with desirable pharmacokinetic profiles for high therapeutic efficacy. The plasma concentrations of NCEs are a biomarker of their efficacy and are governed by pharmacokinetic processes such as absorption, distribution, metabolism, and excretion (ADME). Poor ADME properties of NCEs are a major cause of attrition in drug development. ADME screening is used to identify and optimize lead compounds in the drug discovery process. Computational models predicting ADME properties have been developed with evolving model-building technologies from a simplified relationship between ADME endpoints and physicochemical properties to machine learning, including support vector machines, random forests, and convolution neural networks. Recently, in the field of in silico ADME research, there has been a shift toward evaluating the in vivo parameters or plasma concentrations of NCEs instead of using predictive results to guide chemical structure design. Another research hotspot is the establishment of a computational prediction platform to strengthen academic drug discovery. Bioinformatics projects have produced a series of in silico ADME models using free software and open-access databases. In this review, we introduce prediction models for various ADME parameters and discuss the currently available academic drug discovery platforms.

Incorporating Uremic Solute-mediated Inhibition of OAT1/3 Improves PBPK Prediction of Tenofovir Renal and Systemic Disposition in Patients with Severe Kidney Disease

BACKGROUND

Dose modification of renally secreted drugs in patients with chronic kidney disease (CKD) has relied on serum creatinine concentration as a biomarker to estimate glomerular filtration (GFR) under the assumption that filtration and secretion decline in parallel. A discrepancy between actual renal clearance and predicted renal clearance based on GFR alone is observed in severe CKD patients with tenofovir, a compound secreted by renal OAT1/3. Uremic solutes that inhibit OAT1/3 may play a role in this divergence.

METHODS

To examine the impact of transporter inhibition by uremic solutes on tenofovir renal clearance, we determined the inhibitory potential of uremic solutes hippuric acid, indoxyl sulfate, and p-cresol sulfate. The inhibition parameters (IC50) were incorporated into a previously validated mechanistic kidney model; simulated renal clearance and plasma PK profile were compared to data from clinical studies.

RESULTS

Without the incorporation of uremic solute inhibition, the PBPK model failed to capture the observed data with an absolute average fold error (AAFE) > 2. However, when the inhibition of renal uptake transporters and uptake transporters in the slow distribution tissues were included, the AAFE value was within the pre-defined twofold model acceptance criterion, demonstrating successful model extrapolation to CKD patients.

CONCLUSION

A PBPK model that incorporates inhibition by uremic solutes has potential to better predict renal clearance and systemic disposition of secreted drugs in patients with CKD. Ongoing research is warranted to determine if the model can be expanded to include other OAT1/3 substrate drugs and to evaluate how these findings can be translated to clinical guidance for drug selection and dose optimization in patients with CKD.

Physiologically based pharmacokinetic-pharmacodynamic evaluation of meropenem in CKD and hemodialysis individuals

Objective: Chronic kidney disease (CKD) has significant effects on renal clearance of drugs. The application of antibiotics in CKD patients to achieve the desired therapeutic effect is challenging. This study aims to determine meropenem plasma exposure in the CKD population and further investigate optimal dosing regimens. Methods: A healthy adult PBPK model was established using the meropenem's physicochemical parameters, pharmacokinetic parameters, and available clinical data, and it was scaled to the populations with CKD and dialysis. The differences between the predicted concentration, C_max, and AUC_last predicted and observed model values were assessed by mean relative deviations (MRD) and geometric mean fold errors (GMFE) values and plotting the goodness of fit plot to evaluate the model's performance. Finally, dose recommendations for CKD and hemodialysis populations were performed by Monte Carlo simulations. Results: The PBPK models of meropenem in healthy, CKD, and hemodialysis populations were successfully established. The MRD values of the predicted concentration and the GMFE values of C_max and AUC_last were within 0.5-2.0-fold of the observed data. The simulation results of the PBPK model showed the increase in meropenem exposure with declining kidney function in CKD populations. The dosing regimen of meropenem needs to be further adjusted according to the renal function of CKD patients. In patients receiving hemodialysis, since meropenem declined more rapidly during the on-dialysis session than the off-dialysis session, pharmacodynamic evaluations were performed for two periods separately, and respective optimal dosing regimens were determined. Conclusion: The established PBPK model successfully predicted meropenem pharmacokinetics in patients with CKD and hemodialysis and could further be used to optimize dosing recommendations, providing a reference for personalized clinical medication.

Physiologically based pharmacokinetic modelling and simulation to predict the plasma concentration profile of schaftoside after oral administration of total flavonoids of Desmodium styracifolium

Introduction: The total flavonoids of Desmodium styracifolium (TFDS) are the flavonoid extracts purified from Desmodii Styracifolii Herba. The capsule of TFDS was approved for the treatment of urolithiasis by NMPA in 2022. Schaftoside is the representative compound of TFDS that possesses antilithic and antioxidant effects. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model of schaftoside to simulate its plasma concentration profile in rat and human after oral administration of the total flavonoids of Desmodium styracifolium. Methods: The physiologically based pharmacokinetic model of schaftoside was firstly developed and verified by the pharmacokinetic data in rats following intravenous injection and oral administration of the total flavonoids of Desmodium styracifolium. Then the PBPK model was extrapolated to human with PK-Sim^® software. In order to assess the accuracy of the extrapolation, a preliminary multiple-dose clinical study was performed in four healthy volunteers aged 18-45 years old. The predictive performance of PBPK model was mainly evaluated by visual predictive checks and fold error of C_max and AUC_0-t of schaftoside (the ratio of predicted to observed). Finally, the adult PBPK model was scaled to several subpopulations including elderly and renally impaired patients. Results: Schaftoside underwent poor metabolism in rat and human liver microsomes in vitro, and in vivo it was extensively excreted into urine and bile as an unchanged form. By utilizing literature and experimental data, the PBPK model of schaftoside was well established in rat and human. The predicted plasma concentration profiles of schaftoside were consistent with the corresponding observed data, and the fold error values were within the 2-fold acceptance criterion. No significant pharmacokinetic differences were observed after extrapolation from adult (18-40 years old) to elderly populations (71-80 years) in PK-Sim^®. However, the plasma concentration of schaftoside was predicted to be much higher in renally impaired patients. The maximum steady-state plasma concentrations in patients with chronic kidney disease stage 3, 4 and 5 were 3.41, 12.32 and 23.77 times higher, respectively, than those in healthy people. Conclusion: The established PBPK model of schaftoside provided useful insight for dose selection of the total flavonoids of Desmodium styracifolium in different populations. This study provided a feasible way for the assessment of efficacy and safety of herbal medicines.

Predicting transporter mediated drug-drug interactions via static and dynamic physiologically based pharmacokinetic modeling: A comprehensive insight on where we are now and the way forward

The greater utilization and acceptance of physiologically-based pharmacokinetic (PBPK) modeling to evaluate the potential metabolic drug-drug interactions is evident by the plethora of literature, guidance's, and regulatory dossiers available in the literature. In contrast, it is not widely used to predict transporter-mediated DDI (tDDI). This is attributed to the unavailability of accurate transporter tissue expression levels, the absence of accurate in vitro to in vivo extrapolations (IVIVE), enzyme-transporter interplay, and a lack of specific probe substrates. Additionally, poor understanding of the inhibition/induction mechanisms coupled with the inability to determine unbound concentrations at the interaction site made tDDI assessment challenging. Despite these challenges, continuous improvements in IVIVE approaches enabled accurate tDDI predictions. Furthermore, the necessity of extrapolating tDDI's to special (pediatrics, pregnant, geriatrics) and diseased (renal, hepatic impaired) populations is gaining impetus and is encouraged by regulatory authorities. This review aims to visit the current state-of-the-art and summarizes contemporary knowledge on tDDI predictions. The current understanding and ability of static and dynamic PBPK models to predict tDDI are portrayed in detail. Peer-reviewed transporter abundance data in special and diseased populations from recent publications were compiled, enabling direct input into modeling tools for accurate tDDI predictions. A compilation of regulatory guidance's for tDDI's assessment and success stories from regulatory submissions are presented. Future perspectives and challenges of predicting tDDI in terms of in vitro system considerations, endogenous biomarkers, the use of empirical scaling factors, enzyme-transporter interplay, and acceptance criteria for model validation to meet the regulatory expectations were discussed.

PBPK Simulation-Based Evaluation of Ganciclovir Crystalluria Risk Factors: Effect of Renal Impairment, Old Age, and Low Fluid Intake

Dosing guidance is often lacking for chronic kidney disease (CKD) due to exclusion of such patients from pivotal clinical trials. Physiologically based pharmacokinetic (PBPK) modelling supports model-informed dosing when clinical data are lacking, but application of these approaches to patients with impaired renal function is not yet at full maturity. In the current study, a ganciclovir PBPK model was developed for patients with normal renal function and extended to CKD population. CKD-related changes in tubular secretion were explored in the mechanistic kidney model and implemented either as proportional or non-proportional decline relative to GFR. Crystalluria risk was evaluated in different clinical settings (old age, severe CKD and low fluid intake) by simulating ganciclovir medullary collecting duct (MCD) concentrations. The ganciclovir PBPK model captured observed changes in systemic pharmacokinetic endpoints in mild-to-severe CKD; these trends were evident irrespective of assumed pathophysiological mechanism of altered active tubular secretion in the model. Minimal difference in simulated ganciclovir MCD concentrations was noted between young adult and geriatric populations with normal renal function and urine flow (1 mL/min), with lower concentrations predicted for severe CKD patients. High crystalluria risk was identified at reduced urine flow (0.1 mL/min) as simulated ganciclovir MCD concentrations exceeded its solubility (2.6–6 mg/mL), irrespective of underlying renal function. The analysis highlighted the importance of appropriate distribution of virtual subjects’ systems data in CKD populations. The ganciclovir PBPK model illustrates the ability of this translational tool to explore individual and combined effects of age, urine flow, and renal impairment on local drug renal exposure.

Bridging the gap between in silico and in vivo by modeling opioid disposition in a kidney proximal tubule microphysiological system

Opioid overdose, dependence, and addiction are a major public health crisis. Patients with chronic kidney disease (CKD) are at high risk of opioid overdose, therefore novel methods that provide accurate prediction of renal clearance (CL_r) and systemic disposition of opioids in CKD patients can facilitate the optimization of therapeutic regimens. The present study aimed to predict renal clearance and systemic disposition of morphine and its active metabolite morphine-6-glucuronide (M6G) in CKD patients using a vascularized human proximal tubule microphysiological system (VPT-MPS) coupled with a parent-metabolite full body physiologically-based pharmacokinetic (PBPK) model. The VPT-MPS, populated with a human umbilical vein endothelial cell (HUVEC) channel and an adjacent human primary proximal tubular epithelial cells (PTEC) channel, successfully demonstrated secretory transport of morphine and M6G from the HUVEC channel into the PTEC channel. The in vitro data generated by VPT-MPS were incorporated into a mechanistic kidney model and parent-metabolite full body PBPK model to predict CL_r and systemic disposition of morphine and M6G, resulting in successful prediction of CL_r and the plasma concentration–time profiles in both healthy subjects and CKD patients. A microphysiological system together with mathematical modeling successfully predicted renal clearance and systemic disposition of opioids in CKD patients and healthy subjects.

3D cell culture models: Drug pharmacokinetics, safety assessment, and regulatory consideration

Nonclinical testing has served as a foundation for evaluating potential risks and effectiveness of investigational new drugs in humans. However, the current two-dimensional (2D) in vitro cell culture systems cannot accurately depict and simulate the rich environment and complex processes observed in vivo, whereas animal studies present significant drawbacks with inherited species-specific differences and low throughput for increased demands. To improve the nonclinical prediction of drug safety and efficacy, researchers continue to develop novel models to evaluate and promote the use of improved cell- and organ-based assays for more accurate representation of human susceptibility to drug response. Among others, the three-dimensional (3D) cell culture models present physiologically relevant cellular microenvironment and offer great promise for assessing drug disposition and pharmacokinetics (PKs) that influence drug safety and efficacy from an early stage of drug development. Currently, there are numerous different types of 3D culture systems, from simple spheroids to more complicated organoids and organs-on-chips, and from single-cell type static 3D models to cell co-culture 3D models equipped with microfluidic flow control as well as hybrid 3D systems that combine 2D culture with biomedical microelectromechanical systems. This article reviews the current application and challenges of 3D culture systems in drug PKs, safety, and efficacy assessment, and provides a focused discussion and regulatory perspectives on the liver-, intestine-, kidney-, and neuron-based 3D cellular models.