Use of Physiologically Based Pharmacokinetic Modeling to Evaluate the Effect of Chronic Kidney Disease on the Disposition of Hepatic CYP2C8 and OATP1B Drug Substrates

Mechanistic Determinants of Daprodustat Drug-Drug Interactions and Pharmacokinetics in Hepatic Dysfunction and Chronic Kidney Disease: Significance of OATP1B-CYP2C8 Interplay

Daprodustat is the first oral hypoxia-inducible factor prolyl hydroxylase inhibitor approved recently for the treatment of anemia caused by chronic kidney disease (CKD) in adults receiving dialysis. We evaluated the role of organic anion transporting polypeptide (OATP)1B-mediated hepatic uptake transport in the pharmacokinetics (PKs) of daprodustat using in vitro and in vivo studies, and physiologically-based PK (PBPK) modeling of its drug-drug interactions (DDIs) with inhibitor drugs. In vitro, daprodustat showed specific transport by OATP1B1/1B3 in the transfected cell systems and primary human and monkey hepatocytes. A single-dose oral rifampin (OATP1B inhibitor) reduced daprodustat intravenous clearance by a notable 9.9 ± 1.2-fold (P < 0.05) in cynomolgus monkeys. Correspondingly, volume of distribution at steady-state was also reduced by 5.0 ± 1.1-fold, whereas the half-life change was minimal (1.5-fold), corroborating daprodustat hepatic uptake inhibition by rifampin. A PBPK model accounting for OATP1B-CYP2C8 interplay was developed, which well described daprodustat PK and DDIs with gemfibrozil (CYP2C8 and OATP1B inhibitor) and trimethoprim (weak CYP2C8 inhibitor) within 25% error of the observed data in healthy subjects. About 18-fold increase in daprodustat area under the curve (AUC) following gemfibrozil treatment was found to be associated with strong CYP2C8 inhibition and moderate OATP1B inhibition. Moreover, PK modulation in hepatic dysfunction and subjects with CKD, in comparison to healthy control, was well-captured by the model. CYP2C8 and/or OATP1B inhibitor drugs (e.g., gemfibrozil, clopidogrel, rifampin, and cyclosporine) were predicted to perpetrate moderate-to-strong DDIs in healthy subjects, as well as, in target CKD population. Daprodustat can be used as a sensitive CYP2C8 index substrate in the absence of OATP1B modulation.

Membrane transporters in drug development and as determinants of precision medicine

The effect of membrane transporters on drug disposition, efficacy and safety is now well recognized. Since the initial publication from the International Transporter Consortium, significant progress has been made in understanding the roles and functions of transporters, as well as in the development of tools and models to assess and predict transporter-mediated activity, toxicity and drug-drug interactions (DDIs). Notable advances include an increased understanding of the effects of intrinsic and extrinsic factors on transporter activity, the application of physiologically based pharmacokinetic modelling in predicting transporter-mediated drug disposition, the identification of endogenous biomarkers to assess transporter-mediated DDIs and the determination of the cryogenic electron microscopy structures of SLC and ABC transporters. This article provides an overview of these key developments, highlighting unanswered questions, regulatory considerations and future directions.

Use of physiologically-based pharmacokinetic modeling to understand the effect of omeprazole administration on the pharmacokinetics of oral extended-release nifedipine

Proton pump inhibitors (PPIs) can affect the release of drugs from their dosage forms in vivo by elevating the gastric pH. Our recent clinical study has demonstrated that drug-drug interactions (DDIs) exist between a PPI, omeprazole, and nifedipine extended-release formulations, where systemic exposure of nifedipine was increased in subjects after multiple-dose pretreatment of omeprazole. However, the mechanism of the observed DDIs between omeprazole and nifedipine has not been well-understood, as the DDI may also be mediated through CYP3A4 enzyme inhibition in addition to the elevated gastric pH caused by omeprazole. This study used physiologically-based pharmacokinetic (PBPK) modeling and simulations to investigate the underlying mechanism of these complex DDIs. A formulation exhibiting differences in in vitro dissolution across physiological pH range and another formulation where pH does not impact dissolution appreciably (e.g., an osmotic pump) were chosen to characterize the potential impact of pH. The PBPK models incorporated two-stage in vitro release profiles via US Pharmacopeia 2 apparatus. PBPK simulations suggest that the elevated gastric pH following multiple-dose administration of omeprazole has a minimal effect on nifedipine pharmacokinetics (PKs), whereas CYP3A4-mediated DDI is likely the main driver to the observed change of nifedipine PKs in the presence of omeprazole. Compared to the osmotic formulation, the slightly increased exposure of nifedipine can be accounted for by the enhanced drug release in the pH-dependent formulation. The reported model-based approach may be useful in DDI risk assessments, product formulation designs, and bioequivalence evaluations.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

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.

Development of 4-Pyridoxic Acid PBPK Model to Support Biomarker-Informed Evaluation of OAT1/3 Inhibition and Effect of Chronic Kidney Disease

Monitoring endogenous biomarkers is increasingly used to evaluate transporter-mediated drug-drug interactions (DDIs) in early drug development and may be applied to elucidate changes in transporter activity in disease. 4-pyridoxic acid (PDA) has been identified as the most sensitive plasma endogenous biomarker of renal organic anion transporters (OAT1/3). Increase in PDA baseline concentrations was observed after administration of probenecid, a strong clinical inhibitor of OAT1/3 and also in patients with chronic kidney disease (CKD). The aim of this study was to develop and verify a physiologically-based pharmacokinetic (PBPK) model of PDA, to predict the magnitude of probenecid DDI and predict the CKD-related changes in PDA baseline. The PBPK model for PDA was first developed in healthy population, building on from previous population pharmacokinetic modeling, and incorporating a mechanistic kidney model to consider OAT1/3-mediated renal secretion. Probenecid PBPK model was adapted from the Simcyp database and re-verified to capture its dose-dependent pharmacokinetics (n = 9 studies). The PBPK model successfully predicted the PDA plasma concentrations, area under the curve, and renal clearance in healthy subjects at baseline and after single/multiple probenecid doses. Prospective simulations in severe CKD predicted successfully the increase in PDA plasma concentration relative to healthy (within 2-fold of observed data) after accounting for 60% increase in fraction unbound in plasma and additional 50% decline in OAT1/3 activity beyond the decrease in glomerular filtration rate. The verified PDA PBPK model supports future robust evaluation of OAT1/3 DDI in drug development and increases our confidence in predicting exposure and renal secretion in patients with CKD.

Understanding Statin-Roxadustat Drug-Drug-Disease Interaction Using Physiologically-Based Pharmacokinetic Modeling

A different drug-drug interaction (DDI) scenario may exist in patients with chronic kidney disease (CKD) compared with healthy volunteers (HVs), depending on the interplay between drug-drug and disease (drug-drug-disease interaction (DDDI)). Physiologically-based pharmacokinetic (PBPK) modeling, in lieu of a clinical trial, is a promising tool for evaluating these complex DDDIs in patients. However, the prediction confidence of PBPK modeling in the severe CKD population is still low when nonrenal pathways are involved. More mechanistic virtual disease population and robust validation cases are needed. To this end, we aimed to: (i) understand the implications of severe CKD on statins (atorvastatin, simvastatin, and rosuvastatin) pharmacokinetics (PK) and DDI; and (ii) predict untested clinical scenarios of statin-roxadustat DDI risks in patients to guide suitable dose regimens. A novel virtual severe CKD population was developed incorporating the disease effect on both renal and nonrenal pathways. Drug and disease PBPK models underwent a four-way validation. The verified PBPK models successfully predicted the altered PKs in patients for substrates and inhibitors and recovered the observed statin-rifampicin DDIs in patients and the statin-roxadustat DDIs in HVs within 1.25- and 2-fold error. Further sensitivity analysis revealed that the severe CKD effect on statins PK is mainly mediated by hepatic BCRP for rosuvastatin and OATP1B1/3 for atorvastatin. The magnitude of statin-roxadustat DDI in patients with severe CKD was predicted to be similar to that in HVs. PBPK-guided suitable dose regimens were identified to minimize the risk of side effects or therapeutic failure of statins when co-administered with roxadustat.

Emerging Roles of Uremic Toxins and Inflammatory Cytokines in the Alteration of Hepatic Drug Disposition in Patients with Kidney Dysfunction

Patients with kidney dysfunction exhibit distinct pharmacokinetic profiles compared to those with normal kidney function. Hence, it is desirable to monitor the drug efficacy and toxicity caused by fluctuations in plasma drug concentrations associated with kidney dysfunction. Recently, pharmacokinetic information of drugs excreted mainly through the urine of patients with kidney dysfunction has been reported via drug-labeling information. Pharmacokinetic changes in drugs mainly eliminated by the liver cannot be overlooked as drug metabolism and/or transport activity in the liver may also be altered in patients with kidney dysfunction; however, the underlying mechanisms remain unclear. To plan an appropriate dosage regimen, it is necessary to clarify the underlying processes of functional changes in pharmacokinetic proteins. In recent years, uremic toxins have been shown to reduce the activity and/or expression of renal and hepatic transporters. This inhibitory effect has been reported to be time-dependent. In addition, inflammatory cytokines, such as interleukin-6, released from immune cells activated by uremic toxins and/or kidney injury can reduce the expression levels of drug-metabolizing enzymes and transporters in human hepatocytes. In this mini-review, we have summarized the renal and hepatic pharmacokinetic changes as well as the potential underlying mechanisms in kidney dysfunction, such as the chronic kidney disease and acute kidney injury. SIGNIFICANCE STATEMENT: Patients with kidney dysfunction exhibit distinct pharmacokinetic profiles compared to those with normal kidney function. Increased plasma concentrations of uremic toxins and inflammatory cytokines during kidney disease may potentially affect the activities and/or expression levels of drug-metabolizing enzymes and transporters in the liver and kidneys.

Revisiting the in-vitro and in-vivo considerations for in-silico modelling of complex injectable drug products

Complex injectable drug products (CIDPs) have often been developed to modulate the pharmacokinetics along with efficacy for therapeutic agents used for remediation of chronic disorders. The effective development of CIDPs has exhibited complex kinetics associated with multiphasic drug release from the prepared formulations. Consequently, predictability of pharmacokinetic modelling for such CIDPs has been difficult and there is need for advanced complex computational models for the establishment of accurate prediction models for in-vitro-in-vivo correlation (IVIVC). The computational modelling aims at supplementing the existing knowledge with mathematical equations to develop formulation strategies for generation of predictable and discriminatory IVIVC. Such an approach would help in reduction of the burden of effect of hidden factors on preclinical to clinical translations. Computational tools like physiologically based pharmacokinetics (PBPK) modelling have combined physicochemical and physiological properties along with IVIVC characteristics of clinically used formulations. Such techniques have helped in prediction and understanding of variability in pharmacodynamic parameters of potential generic products to clinically used formulations like Doxil®, Ambisome®, Abraxane® in healthy and diseased population using mathematical equations. The current review highlights the important formulation characteristics, in-vitro, preclinical in-vivo aspects which need to be considered while developing a stimulatory predictive PBPK model in establishment of an IVIVC and in-vitro-in-vivo relationship (IVIVR).

Physiologically Based Pharmacokinetic Modelling to Predict Imatinib Exposures in Cancer Patients with Renal Dysfunction: A Case Study

Imatinib is mainly metabolised by CYP3A4 and CYP2C8 and is extensively bound to α-acid glycoprotein (AAG). A physiologically based pharmacokinetic (PBPK) model for imatinib describing the CYP3A4-mediated autoinhibition during multiple dosing in gastrointestinal stromal tumor patients with normal renal function was previously reported. After performing additional verification, the PBPK model was applied to predict the exposure of imatinib after multiple dosing in cancer patients with varying degrees of renal impairment. In agreement with the clinical data, there was a positive correlation between AAG levels and imatinib exposure. A notable finding was that for recovery of the observed data in cancer patients with moderate RI (CrCL 20 to 39 mL/min), reductions of hepatic CYP3A4 and CYP2C8 abundances, which reflect the effects of RI, had to be included in the simulations. This was not the case for mild RI (CrCL 40 to 50 mL/min). The results support the finding of the clinical study, which demonstrated that both AAG levels and the degree of renal impairment are key components that contribute to the interpatient variability associated with imatinib exposure. As indicated in the 2020 FDA draft RI guidance, PBPK modelling could be used to support an expanded inclusion of patients with RI in clinical studies.

Single-dose Pharmacokinetics of Eluxadoline in Healthy Participants With Normal Renal Function and Participants With Renal Impairment

Eluxadoline is approved for the treatment of diarrhea-predominant irritable bowel syndrome in the United States. The impact of renal impairment on the pharmacokinetic (PK) parameters of eluxadoline is currently unknown. This phase 1, open-label, parallel-group study evaluated the PK and safety profile of eluxadoline in 8 participants with renal impairment and 8 matched healthy controls. Of the participants with renal impairment, 2 had severe renal impairment (estimated glomerular filtration rate [eGFR] <30 mL/min/1.73 m² ) and 6 had end-stage renal disease while not yet on dialysis (eGFR <15 mL/min/1.73 m² ). The primary objective was to assess plasma and urine PKs, and plasma protein binding of eluxadoline. In participants with renal impairment, mean plasma concentrations of eluxadoline were consistently higher compared with matched healthy controls: 1.4-fold higher for mean maximum plasma concentration (C_max ) and 2.2-fold higher for mean area under the plasma concentration-time curve from time 0 to time t. The median time to C_max was 2.5 hours in both groups. Although eluxadoline is a locally acting drug with low oral bioavailability, because of the increased systemic exposure in participants with renal impairment as a cautionary measure the lower approved dose of 75 mg twice daily is recommended for individuals with severe renal impairment and end-stage renal disease while not yet on dialysis. Eluxadoline 100 mg single dose was well tolerated in participants with renal impairment and matched healthy controls.

Effect of Hepatic Impairment on OATP1B Activity: Quantitative Pharmacokinetic Analysis of Endogenous Biomarker and Substrate Drugs

Hepatic impairment (HI) is known to modulate drug disposition and may lead to elevated plasma exposure. The aim of this study was to quantitate the in vivo OATP1B-mediated hepatic uptake activity in populations with varying degrees of HI. First, we measured baseline levels of plasma coproporphyrin-I, an endogenous OATP1B biomarker, in an open-label, parallel cohort study in adult subjects with normal liver function and mild, moderate, and severe HI (n = 24, 6/cohort). The geometric mean plasma concentrations of coproporphyrin-I were 1.66-fold, 2.81-fold (P < 0.05), and 7.78-fold (P < 0.0001) higher in mild, moderate, and severe impairment than those healthy controls. Second, we developed a dataset of 21 OATP1B substrate drugs with HI data extracted from literature. Median disease-to-healthy plasma area under the curve (AUC) ratios for substrate drugs were ~ 1.4, 3.0, and 6.4 for mild, moderate, and severe HI, respectively. Additionally, significant linear relationship was noted between AUC ratios of substrate drugs without and with co-administration of rifampin, a prototypic OATP1B inhibitor, and AUC ratios in moderate (P < 0.01) and severe (P < 0.001) HI. Third, a physiologically-based pharmacokinetic model analysis was conducted with 10 substrate drugs following estimation of relative OATP1B functional activity in virtual disease population models using coproporphyrin-I data and verification of substrate models with rifampin drug-drug interaction data. This approach adequately predicted plasma AUC change particularly in moderate (9 of 10 within 2-fold) and severe (5 of 5 within 2-fold) HI. Collective findings indicate progressive reduction, by as much as 90-92%, in OATP1B activity in the HI population.

Torsemide Pharmacometrics in Healthy Adult Populations Including CYP2C9 Genetic Polymorphisms and Various Patient Groups through Physiologically Based Pharmacokinetic-Pharmacodynamic Modeling

Torsemide is a widely used diuretic in clinical practice. In this study, pharmacokinetic (PK) and pharmacodynamic (PD) simulations of torsemide for various population groups and exposure scenarios were performed through human-scale physiologically-based PK-PD (PBPK-PD) modeling of torsemide. For PBPK-PD modeling of torsemide, invitro and clinical data of torsemide reported previously were used. After exposure to clinical doses of torsemide, observed plasma (or serum) concentration and urine torsemide excretion profiles were used as PK-data, and observed urinary sodium excretion rate was used as PD-data. The model was then extended to take into account physiological and biochemical factors according to different CYP2C9 phenotypes or patient populations. The established model captured various torsemide clinical results well. Differences in torsemide PKs and PDs between patient groups or CYP2C9 genetic polymorphisms were modelologically identified. It was confirmed that degrees of differences in torsemide PKs and PDs by disease groups were greater than those according to different CYP2C9 phenotypes. According to torsemide administration frequency or dose change, it was confirmed that although the difference in plasma PKs between groups (healthy adult and patient groups) could increase to 14.80 times, the difference in PDs was reduced to 1.01 times. Results of this study suggested that it is very important to consider disease groups in the setting of torsemide clinical therapy and that it is difficult to predict PD proportionally with only differences in PKs of torsemide between population groups. The PBPK-PD model established in this study is expected to be utilized for various clinical cases involving torsemide application in the future, enabling optimal drug therapy.

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.

Quantitative Proteomics in Translational Absorption, Distribution, Metabolism, and Excretion and Precision Medicine

A reliable translation of in vitro and preclinical data on drug absorption, distribution, metabolism, and excretion (ADME) to humans is important for safe and effective drug development. Precision medicine that is expected to provide the right clinical dose for the right patient at the right time requires a comprehensive understanding of population factors affecting drug disposition and response. Characterization of drug-metabolizing enzymes and transporters for the protein abundance and their interindividual as well as differential tissue and cross-species variabilities is important for translational ADME and precision medicine. This review first provides a brief overview of quantitative proteomics principles including liquid chromatography-tandem mass spectrometry tools, data acquisition approaches, proteomics sample preparation techniques, and quality controls for ensuring rigor and reproducibility in protein quantification data. Then, potential applications of quantitative proteomics in the translation of in vitro and preclinical data as well as prediction of interindividual variability are discussed in detail with tabulated examples. The applications of quantitative proteomics data in physiologically based pharmacokinetic modeling for ADME prediction are discussed with representative case examples. Finally, various considerations for reliable quantitative proteomics analysis for translational ADME and precision medicine and the future directions are discussed. SIGNIFICANCE STATEMENT: Quantitative proteomics analysis of drug-metabolizing enzymes and transporters in humans and preclinical species provides key physiological information that assists in the translation of in vitro and preclinical data to humans. This review provides the principles and applications of quantitative proteomics in characterizing in vitro, ex vivo, and preclinical models for translational research and interindividual variability prediction. Integration of these data into physiologically based pharmacokinetic modeling is proving to be critical for safe, effective, timely, and cost-effective drug development.

Population Pharmacokinetic Modelling for Nifedipine to Evaluate the Effect of Parathyroid Hormone on CYP3A in Patients with Chronic Kidney Disease

Purpose

Parathyroid hormone (PTH) can induce the downregulation of CYP3A in chronic kidney disease (CKD). Nevertheless, the effect of PTH on CYP3A-mediated clearance pathways from a clinical perspective remains unclear.

Methods

This study employed population pharmacokinetic (PopPK) modeling to delineate potential changes in CYP3A activity in patients with CKD. Pharmacokinetic data for nifedipine, a typical CYP3A substrate, as well as covariate information, were prospectively collected from 157 patients with a total of 612 concentrations. PopPK data analysis was performed using a nonlinear mixed-effects model.

Results

The pharmacokinetics of nifedipine were optimally described according to a one-compartment model with zero-order absorption and first-order elimination. The estimated population parameters (and interindividual variability) were apparent clearance (CL/F) 49.61 L/h (58.33%) and apparent volume of distribution (V/F) 2300.26 L (45.62%), and the PTH level negatively correlated with CL/F. In comparison with the reference level, it was observed that the dosage of nifedipine should be reduced with the maximum boundary value of PTH, after a Monte Carlo simulation.

Conclusion

This study provides insight into the effects of PTH on CYP3A-mediated clearance pathways. Moreover, PTH could be used as a guide for the appropriate administration of CYP3A eliminated drugs in patients with CKD.

Utilization of Pharmacokinetic/Pharmacodynamic Modeling in Pharmacoepidemiological Studies: A Systematic Review on Antiarrhythmic and Glucose-Lowering Medicines

Introduction: In human pharmacology, there are two important scientific branches: clinical pharmacology and pharmacoepidemiology. Pharmacokinetic/pharmacodynamic (PK/PD) modeling is important in preclinical studies and randomized control trials. However, it is rarely used in pharmacoepidemiological studies on the effectiveness and medication safety where the target population is heterogeneous and followed for longer periods. The objective of this literature review was to investigate how far PK/PD modeling is utilized in observational studies on glucose-lowering and antiarrhythmic drugs.

Method: A systematic literature search of MEDLINE, Embase, and Web of Science was conducted from January 2010 to 21 February 2020. To calculate the utilization of PK/PD modeling in observational studies, we followed two search strategies. In the first strategy, we screened a 1% random set from 95,672 studies on glucose-lowering and antiarrhythmic drugs on inclusion criteria. In the second strategy, we evaluated the percentage of studies in which PK/PD modeling techniques were utilized. Subsequently, we divided the total number of included studies in the second search strategy by the total number of eligible studies in the first search strategy.

Results: The comprehensive search of databases and the manual search of included references yielded a total of 29 studies included in the qualitative synthesis of our systematic review. Nearly all 29 studies had utilized a PK model, whereas only two studies developed a PD model to evaluate the effectiveness of medications. In total, 16 out of 29 studies (55.1%) used a PK/PD model in the observational setting to study effect modification. The utilization of PK/PD modeling in observational studies was calculated as 0.42%.

Conclusion: PK/PD modeling techniques were substantially underutilized in observational studies of antiarrhythmic and glucose-lowering drugs during the past decade.

Coproporphyrin I as an Endogenous Biomarker to Detect Reduced OATP1B Activity and Shift in Elimination Route in Chronic Kidney Disease

Coproporphyrin I (CPI) is an endogenous biomarker of organic anion transporting polypeptide 1B transporter (OATP1B). CPI plasma baseline was reported to increase with severity of chronic kidney disease (CKD). Further, ratio of CPI area under the plasma concentration-time curve (AUCR) in the presence/absence of OATP1B inhibitor rifampin was higher in patients with CKD compared with healthy participants, in contrast to pitavastatin (a clinical OATP1B probe). This study investigated mechanism(s) contributing to altered CPI baseline in patients with CKD by extending a previously developed physiologically-based pharmacokinetic (PBPK) model to this patient population. CKD-related covariates were evaluated in a stepwise manner on CPI fraction unbound in plasma (f_u,p ), OATP1B-mediated hepatic uptake clearance (CL_active ), renal clearance (CL_R ), and endogenous synthesis (k_syn ). The CPI model successfully recovered increased baseline and rifampin-mediated AUCR in patients with CKD by accounting for the following disease-related changes: 13% increase in f_u,p , 29% and 39% decrease in CL_active in mild and moderate to severe CKD, respectively, decrease in CL_R proportional to decline in glomerular filtration rate, and 27% decrease in k_syn in severe CKD. Almost complete decline in CPI renal elimination in severe CKD increased its fraction transported by OATP1B, rationalizing differences in the CPI-rifampin interaction observed between healthy participants and patients with CKD. In conclusion, mechanistic modeling performed here supports CKD-related decrease in OATP1B function to inform prospective PBPK modeling of OATP1B-mediated drug-drug interaction in these patients. Monitoring of CPI allows detection of CKD-drug interaction risk for OATP1B drugs with combined hepatic and renal elimination which may be underestimated by extrapolating the interaction risk based on pitavastatin data in healthy participants.

Effect of Chronic Kidney Disease on the Renal Secretion via Organic Anion Transporters 1/3: Implications for Physiologically-Based Pharmacokinetic Modeling and Dose Adjustment

There is growing evidence that active tubular secretory clearance (CL_s) may not decline proportionally with the glomerular filtration rate (GFR) in chronic kidney disease (CKD), leading to the overestimation of renal clearance (CL_r) when using solely GFR to approximate disease effect on renal elimination. The clinical pharmacokinetic data of 33 renally secreted OAT1/3 substrates were collated to investigate the impact of mild, moderate, and severe CKD on CL_r, tubular secretion and protein binding (f_u,p). The f_u,p of the collated substrates ranged from 0.0026 to 1.0 in healthy populations; observed CKD‐related increase in the f_u,p (up to 2.7‐fold) of 8 highly bound substrates (f_u,p ≤ 0.2) was accounted for in the analysis. Use of prediction equation based on disease‐related changes in albumin resulted in underprediction of the CKD‐related increase in f_u,p of highly bound substrates, highlighting the necessity to measure protein binding in severe CKD. The critical analysis of clinical data for 33 OAT1/3 probes established that decrease in OAT1/3 activity proportional to the changes in GFR was insufficient to recapitulate effects of severe CKD on unbound tubular secretion clearance. OAT1/3‐mediated CL_s was estimated to decline by an additional 50% relative to the GFR decline in severe CKD, whereas change in active secretion in mild and moderate CKD was proportional to GFR. Consideration of this additional 50% decline in OAT1/3‐mediated CL_s is recommended for physiologically‐based pharmacokinetic models and dose adjustment of OAT1/3 substrates in severe CKD, especially for substrates with high contribution of the active secretion to CL_r.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

Clinical Implications of Altered Drug Transporter Abundance/Function and PBPK Modeling in Specific Populations: An ITC Perspective

The role of membrane transporters on pharmacokinetics (PKs), drug-drug interactions (DDIs), pharmacodynamics (PDs), and toxicity of drugs has been broadly recognized. However, our knowledge of modulation of transporter expression and/or function in the diseased patient population or specific populations, such as pediatrics or pregnancy, is still emerging. This white paper highlights recent advances in studying the changes in transporter expression and activity in various diseases (i.e., renal and hepatic impairment and cancer) and some specific populations (i.e., pediatrics and pregnancy) with the focus on clinical implications. Proposed alterations in transporter abundance and/or activity in diseased and specific populations are based on (i) quantitative transporter proteomic data and relative abundance in specific populations vs. healthy adults, (ii) clinical PKs, and emerging transporter biomarker and/or pharmacogenomic data, and (iii) physiologically-based pharmacokinetic modeling and simulation. The potential for altered PK, PD, and toxicity in these populations needs to be considered for drugs and their active metabolites in which transporter-mediated uptake/efflux is a major contributor to their absorption, distribution, and elimination pathways and/or associated DDI risk. In addition to best practices, this white paper discusses current challenges and knowledge gaps to study and quantitatively predict the effects of modulation in transporter activity in these populations, together with the perspectives from the International Transporter Consortium (ITC) on future directions.

Determining the Effects of Chronic Kidney Disease on Organic Anion Transporter1/3 Activity Through Physiologically Based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The renal excretion of drugs via organic anion transporters 1 and 3 (OAT1/3) is significantly decreased in patients with renal impairment. This study uses physiologically based pharmacokinetic models to quantify the reduction in OAT1/3-mediated secretion of drugs throughout varying stages of chronic kidney disease.

METHODS

Physiologically based pharmacokinetic models were constructed for four OAT1/3 substrates in healthy individuals: acyclovir, meropenem, furosemide, and ciprofloxacin. Observed data from drug-drug interaction studies with probenecid, a potent OAT1/3 inhibitor, were used to parameterize the contribution of OAT1/3 to the renal elimination of each drug. The models were then translated to patients with chronic kidney disease by accounting for changes in glomerular filtration rate, kidney volume, renal blood flow, plasma protein binding, and hematocrit. Additionally, a relationship was derived between the estimated glomerular filtration rate and the reduction in OAT1/3-mediated secretion of drugs based on the renal extraction ratios of ƿ-aminohippuric acid in patients with varying degrees of renal impairment. The relationship was evaluated in silico by evaluating the predictive performance of each final model in describing the pharmacokinetics (PK) of drugs across stages of chronic kidney disease.

RESULTS

OAT1/3-mediated renal excretion of drugs was found to be decreased by 27-49%, 50-68%, and 70-96% in stage 3, stage 4, and stage 5 of chronic kidney disease, respectively. In support of the parameterization, physiologically based pharmacokinetic models of four OAT1/3 substrates were able to adequately characterize the PK in patients with different degrees of renal impairment. Total exposure after intravenous administration was predicted within a 1.5-fold error and 85% of the observed data points fell within a 1.5-fold prediction error. The models modestly under-predicted plasma concentrations in patients with end-stage renal disease undergoing intermittent hemodialysis. However, results should be interpreted with caution because of the limited number of molecules analyzed and the sparse sampling in observed chronic kidney disease pharmacokinetic studies.

CONCLUSIONS

A quantitative understanding of the reduction in OAT1/3-mediated excretion of drugs in differing stages of renal impairment will contribute to better predictive accuracy for physiologically based pharmacokinetic models in drug development, assisting with clinical trial planning and potentially sparing this population from unnecessary toxic exposures.

Clinical Effects of Exercise Rehabilitation Combined with Repaglinide in the Treatment of Diabetes

Objective

Diabetes, a common endocrine and metabolic disease in clinical practice, generally manifests a certain defect in insulin secretion due to several factors, thereafter leading to a metabolic disorder such as hyperglycemia. This study was conducted to explore the clinical effects of repaglinide combined with exercise rehabilitation on improving the blood glucose of patients with diabetes.

Methods

In this retrospective study, 100 patients with diabetes treated in our hospital from January 2018 to January 2020 were assessed for eligibility and recruited. They were assigned at a ratio of 1 : 1 to receive either repaglinide (control group) or repaglinide plus exercise rehabilitation (experimental group). Outcome measures include fasting blood glucose, 2 h postprandial blood glucose, glycosylated hemoglobin, time to normal blood glucose, blood glucose fluctuation, insulin dosage, adverse reactions, and blood glucose adequate rate.

Results

All eligible patients showed similar pretreatment fasting blood glucose, glycosylated hemoglobin, and 2 h postprandial blood glucose (P > 0.05). After treatment, repaglinide plus exercise rehabilitation resulted in lower levels of fasting blood glucose, glycosylated hemoglobin, and 2 h postprandial blood glucose versus repaglinide alone (P < 0.05). Repaglinide plus exercise rehabilitation was associated with a significantly shorter time to normal blood glucose and a milder fluctuation versus repaglinide (P < 0.05). The incidence of adverse reactions and blood glucose adequate rate was 6% and 94% in the experimental group and 50% and 52% in the control group, respectively (P < 0.05).

Conclusion

Repaglinide plus exercise rehabilitation results in effective blood glucose control and reduced incidence of adverse reactions and yields a promising efficacy, so it is worthy of clinical promotion and application.

Use of Physiologically Based Pharmacokinetic Modeling to Evaluate the Impact of Chronic Kidney Disease on CYP3A4-Mediated Metabolism of Saxagliptin

We characterized the impact of chronic kidney disease (CKD) on the cytochrome P450 (CYP) 3A4–mediated metabolism of saxagliptin to its metabolite, 5‐hydroxysaxagliptin, using a physiologically based pharmacokinetic (PBPK) model. A PBPK model of saxagliptin and its CYP3A4 metabolite, 5‐hydroxysaxagliptin, was constructed and validated for oral doses ranging from 5 to 100 mg. The observed ratios of area under the plasma concentration–time curve (AUC) and maximum plasma concentration (C_max) between healthy subjects and subjects with CKD were compared with those predicted using PBPK model simulations. Simulations were performed with virtual CKD populations having decreased CYP3A4 activity (ie, 64%‐75% of the healthy subjects’ CYP3A4 abundance) and preserved CYP3A4 activity (ie, 100% of the healthy subjects’ CYP3A4 abundance). We found that simulations using decreased CYP3A4 activity generally overpredicted the ratios of saxagliptin AUC and C_max in CKD compared with those using preserved CYP3A4 activity. Similarly, simulations using decreased CYP3A4 activity underpredicted the ratio of 5‐hydroxysaxagliptin AUC in moderate and severe CKD compared with simulations using preserved CYP3A4 activity. These findings suggest that decreased CYP3A4 activity in CKD underpredicts saxagliptin clearance compared with that observed clinically. Preserving CYP3A4 activity in CKD more closely estimates saxagliptin clearance and 5‐hydroxysaxagliptin exposure changes observed in vivo. Our findings suggest that there is no clinically meaningful impact of CKD on the metabolism of saxagliptin by CYP3A4. Since saxagliptin is not a highly sensitive substrate and validated probe for CYP3A4, this work represents a case study of a CYP3A4 substrate‐metabolite pair and is not a generalization for all CYP3A4 substrates.

Application of a Physiologically Based Pharmacokinetic Approach to Predict Theophylline Pharmacokinetics Using Virtual Non-Pregnant, Pregnant, Fetal, Breast-Feeding, and Neonatal Populations

Perinatal pharmacology is influenced by a myriad of physiological variables that are changing dynamically. The influence of these covariates has not been assessed systemically. The objective of this work was to use theophylline as a model drug and to predict its pharmacokinetics before, during (including prediction of the umbilical cord level), and after pregnancy as well as in milk (after single and multiple doses) and in neonates using a physiological-based pharmacokinetic (PBPK) model. Neonatal theophylline exposure from milk consumption was projected in both normal term and preterm subjects. Predicted infant daily doses were calculated using theophylline average and maximum concentration in the milk as well as an estimate of milk consumption. Predicted concentrations and parameters from the PBPK model were compared to the observed data. PBPK predicted theophylline concentrations in non-pregnant and pregnant populations at different gestational weeks were within 2-fold of the observations and the observed concentrations fell within the 5th-95th prediction interval from the PBPK simulations. The PBPK model predicted an average cord-to-maternal plasma ratio of 1.0, which also agrees well with experimental observations. Predicted postpartum theophylline concentration profiles in milk were also in good agreement with observations with a predicted milk-to-plasma ratio of 0.68. For an infant of 2 kg consuming 150 ml of milk per day, the lactation model predicted a relative infant dose (RID) of 12 and 17% using predicted average (C_avg,ss) and maximum (C_max,ss) concentration in milk at steady state. The maximum RID of 17% corresponds to an absolute infant daily dose of 1.4 ± 0.5 mg/kg/day. This dose, when administered as 0.233 mg/kg every 4 h, to resemble breastfeeding frequency, resulted in plasma concentrations as high as 3.9 (1.9-6.8) mg/L and 2.8 (1.3-5.3) (5th-95th percentiles) on day 7 in preterm (32 GW) and full-term neonatal populations.

Quantification of CYP3A and Drug Transporters Activity in Healthy Young, Healthy Elderly and Chronic Kidney Disease Elderly Patients by a Microdose Cocktail Approach

Background: Ageing and chronic kidney disease (CKD) affect pharmacokinetic (PK) parameters. Since mechanisms are related and remain unclear, cytochrome P450 (CYP) 3A and drug transporter activities were investigated in the elderly with or without CKD and compared to healthy adults using a microdose cocktail. Methods: Healthy young participants (n = 20), healthy elderly participants (n = 16) and elderly patients with CKD (n = 17) received, in study period 1, a single dose of microdose cocktail probe containing 30 µg midazolam, 750 µg dabigatran etexilate, 100 µg atorvastatin, 10 µg pitavastatin, and 50 µg rosuvastatin. After a 14-day wash-out period, healthy young participants continued to study period 2 with the microdose cocktail plus rifampicin. PK parameters including area under the plasma concentration-time curve (AUC), maximum plasma drug concentration (C_max), and half-life were estimated before making pairwise comparisons of geometric mean ratios (GMR) between groups. Results: AUC and C_max GMR (95% confidence interval; CI) of midazolam, a CYP3A probe substrate, were increased 2.30 (1.70-3.09) and 2.90 (2.16-3.88) fold in healthy elderly and elderly patients with CKD, respectively, together with a prolonged half-life. AUC and C_max GMR (95%CI) of atorvastatin, another CYP3A substrate, was increased 2.14 (1.52-3.02) fold in healthy elderly and 4.15 (2.98-5.79) fold in elderly patients with CKD, indicating decreased CYP3A activity related to ageing. Associated AUC changes in the probe drug whose activity could be modified by intestinal P-glycoprotein (P-gp) activity, dabigatran etexilate, were observed in patients with CKD. However, whether the activity of pitavastatin and rosuvastatin is modified by organic anion transporting polypeptide 1B (OATP1B) and of breast cancer resistance protein (BCRP), respectively, in elderly participants with or without CKD was inconclusive. Conclusions: CYP3A activity is reduced in ageing. Intestinal P-gp function might be affected by CKD, but further confirmation appears warranted. Clinical Trial Registration:http://www.thaiclinicaltrials.org/ (TCTR 20180312002 registered on March 07, 2018).

Assessing CYP2C8-Mediated Pharmaceutical Excipient-Drug Interaction Potential: A Case Study of Tween 80 and Cremophor EL-35

Pharmaceutical excipients (PEs) are substances included in drug formulations. Recent studies have revealed that some PEs can affect the activity of metabolic enzymes and drug transporters; however, the effects of PEs on CYP2C8 and its interaction potential with drugs remain unclear. In this study, we evaluated the effects of Tween 80 and EL−35 on CYP2C8 in vitro and further investigated their impacts on the PK of paclitaxel (PTX) in rats after single or multiple doses. The in vitro study indicated that Tween 80 and EL−35 inhibited CYP2C8 activity in human and rat liver microsomes. EL−35 also decreased the expression of CYP2C8 in HepG2 cells. In the in vivo study, Tween 80 did not alter the PK of PTX after single or multiple doses, whereas EL−35 administered for 14 days significantly increased the AUC and MRT of PTX. Further analysis indicated that multiple-dose EL−35 reduced the expression of Cyp2c22 and production of 6-OH-PTX in the rat liver. Our study suggested that short-term exposure to both PEs did not affect the PK of PTX in rats, but multiple doses of EL−35 increased the AUC and MRT of PTX by downregulating the hepatic expression of Cyp2c22. Such effects should be taken into consideration during drug formulation and administration.

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.

Application of Individualized PBPK Modeling of Rate Data to Evaluate the Effect of Hemodialysis on Nonrenal Clearance Pathways

The aim of this study was to apply individualized, physiologically based pharmacokinetic modeling of ¹⁴ CO₂ production rates (iPBPK-R) measured by the erythromycin breath test to characterize the effect of hemodialysis on the function of nonrenal clearance pathways in patients with end-stage renal disease. Twelve patients previously received ¹⁴ C-erythromycin intravenously pre- and post-hemodialysis. Serial breath samples were collected after each dose over 2 hours. Eight PBPK parameters were co-estimated across periods, whereas activity of cytochrome P450 (CYP) 3A4 clearance was independently estimated for each period. Inhibition coefficients for organic anion transporting polypeptide (OATP), P-glycoprotein, and multidrug resistance-associated protein 2 activities were also estimated. Nonrenal clearance parameter estimates were explored regarding sex differences and correlations with uremic toxins and were used in hierarchical cluster analysis (HCA). Relationships between the function of nonrenal clearance pathways and uremic toxin concentrations were explored. Mean CYP 3A4 clearance increased by 2.2% post-hemodialysis. Uptake transporter activity was highly intervariable across hemodialysis. Females had 22% and 30% higher median CYP3A4 activity than males pre- and post-hemodialysis, respectively. Exploratory HCA indicated that using both CYP3A4 and OATP activity parameters at pre- and post-hemodialysis best identifies heterogeneous patients. This is the first study to use the iPBPK-R approach to simultaneously estimate multiple in vivo nonrenal elimination pathways in individual patients with kidney disease and to assess the effect of hemodialysis.

Renohepatic crosstalk: a review of the effects of acute kidney injury on the liver

Several theories regarding acute kidney injury (AKI)-related mortality have been entertained, although mounting evidence supports the paradigm that impaired kidney function directly and adversely affects the function of several remote organs. The kidneys and liver are fundamental to human metabolism and detoxification, and it is therefore hardly surprising that critical illness complicated by hepatorenal dysfunction portends a poor prognosis. Several diseases can simultaneously impact the proper functioning of the liver and kidneys, although this review will address the impact of AKI on liver function. While evidence for this relationship in humans remains sparse, we present supportive studies and then discuss the most likely mechanisms by which AKI can cause liver dysfunction. These include 'traditional' complications of AKI (uremia, volume overload and acute metabolic acidosis, among others) as well as systemic inflammation, hepatic leukocyte infiltration, cytokine-mediated liver injury and hepatic oxidative stress. We conclude by addressing the therapeutic implications of these findings to clinical medicine.

Physiologically-Based Pharmacokinetic Modeling in Renal and Hepatic Impairment Populations: A Pharmaceutical Industry Perspective

The predictive performance of physiologically‐based pharmacokinetics (PBPK) models for pharmacokinetics (PK) in renal impairment (RI) and hepatic impairment (HI) populations was evaluated using clinical data from 29 compounds with 106 organ impairment study arms were collected from 19 member companies of the International Consortium for Innovation and Quality in Pharmaceutical Development. Fifty RI and 56 HI study arms with varying degrees of organ insufficiency along with control populations were evaluated. For RI, the area under the curve (AUC) ratios of RI to healthy control were predicted within twofold of the observed ratios for > 90% (N = 47/50 arms). For HI, > 70% (N = 43/56 arms) of the hepatically impaired to healthy control AUC ratios were predicted within twofold. Inaccuracies, typically overestimation of AUC ratios, occurred more in moderate and severe HI. PBPK predictions can help determine the need and timing of organ impairment study. It may be suitable for predicting the impact of RI on PK of drugs predominantly cleared by metabolism with varying contribution of renal clearance. PBPK modeling may be used to support mild impairment study waivers or clinical study design.

Physiologically-Based Pharmacokinetic Modelling of Creatinine-Drug Interactions in the Chronic Kidney Disease Population

Elevated serum creatinine (S_Cr ) caused by the inhibition of renal transporter(s) may be misinterpreted as kidney injury. The interpretation is more complicated in patients with chronic kidney disease (CKD) due to altered disposition of creatinine and renal transporter inhibitors. A clinical study was conducted in 17 patients with CKD (estimated glomerular filtration rate 15-59 mL/min/1.73 m² ); changes in S_Cr were monitored during trimethoprim treatment (100-200 mg/day), administered to prevent recurrent urinary infection, relative to the baseline level. Additional S_Cr -interaction data with trimethoprim, cimetidine, and famotidine in patients with CKD were collated from the literature. Our published physiologically-based creatinine model was extended to predict the effect of the CKD on S_Cr and creatinine-drug interaction. The creatinine-CKD model incorporated age/sex-related differences in creatinine synthesis, CKD-related glomerular filtration deterioration; change in transporter activity either proportional or disproportional to glomerular filtration rate (GFR) decline were explored. Optimized models successfully recovered baseline S_Cr from 64 patients with CKD (geometric mean fold-error of 1.1). Combined with pharmacokinetic models of inhibitors, the creatinine model was used to simulate transporter-mediated creatinine-drug interactions. Use of inhibitor unbound plasma concentrations resulted in 66% of simulated S_Cr interaction data within the prediction limits, with cimetidine interaction significantly underestimated. Assuming that transporter activity deteriorates disproportional to GFR decline resulted in higher predicted sensitivity to transporter inhibition in patients with CKD relative to healthy patients, consistent with sparse clinical data. For the first time, this novel modelling approach enables quantitative prediction of S_Cr in CKD and delineation of the effect of disease and renal transporter inhibition in this patient population.

Dose Optimization in Kidney Disease: Opportunities for PBPK Modeling and Simulation

Kidney disease affects pharmacokinetic (PK) profiles of not only renally cleared drugs but also nonrenally cleared drugs. The impact of kidney disease on drug disposition has not been fully elucidated, but describing the extent of such impact is essential for conducting dose optimization in kidney disease. Accurate evaluation of kidney function has been a clinical interest for dose optimization, and more scientists pay attention and conduct research for clarifying the role of drug transporters, metabolic enzymes, and their interplay in drug disposition as kidney disease progresses. Physiologically based pharmacokinetic (PBPK) modeling and simulation can provide valuable insights for dose optimization in kidney disease. It is a powerful tool to integrate discrete knowledge from preclinical and clinical research and mechanistically investigate system- and drug-dependent factors that may contribute to the changes in PK profiles. PBPK-based prediction of drug exposures may be used a priori to adjust dosing regimens and thereby minimize the likelihood of drug-related toxicity. With real-time clinical studies, parameter estimation may be performed with PBPK approaches that can facilitate identification of sources of interindividual variability. PBPK modeling may also facilitate biomarker research that aids dose optimization in kidney disease. U.S. Food and Drug Administration guidances related to conduction of PK studies in kidney impairment and PBPK documentation provide the foundation for facilitating model-based dose-finding research in kidney disease.

6-Hydroxyindole is an endogenous long-lasting OATP1B1 inhibitor elevated in renal failure patients

The hepatic uptake transporter organic anion transporting polypeptide (OATP) 1B1 is inhibited by some uremic toxins; however, direct inhibition can only partially explain the delayed systemic elimination of substrate drugs in renal failure patients. This study aimed to examine the long-lasting inhibition of OATP1B1 by uremic toxins and their metabolites. Preincubation of HEK293/OATP1B1 cells with 21 uremic toxins resulted in almost no change in the uptake of a typical substrate [³H]estrone-3-sulfate (E₁S), although some directly inhibited [³H]E₁S uptake. In contrast, preincubation with an indole metabolite, 6-hydroxyindole, reduced [³H]E₁S uptake, even after the inhibitor was washed out before [³H]E₁S incubation. Such long-lasting inhibition by 6-hydroxyindole was time-dependent and recovered after a 3-h incubation without 6-hydroxyindole. Preincubation with 6-hydroxyindole increased the Km for [³H]E₁S uptake with minimal change in Vmax. This was compatible with no change in the cell-surface expression of OATP1B1, as assessed by a biotinylation assay. Preincubation with 6-hydroxyindole reduced [³H]E₁S uptake in human hepatocytes without changes in OATP1B1 mRNA. Plasma concentration of 6-hydroxyindole in renal failure patients increased as renal function decreased, but might be insufficient to exhibit potent OATP1B1 inhibition. In conclusion, 6-hydroxyindole is an endogenous long-lasting OATP1B1 inhibitor with elevated plasma concentrations in renal failure patients.

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

Chronic kidney disease (CKD) has significant effects on renal clearance (CLr ) of drugs. Physiologically-based pharmacokinetic (PBPK) models have been used to predict CKD effects on transporter-mediated renal active secretion and CLr for hydrophilic nonpermeable compounds. However, no studies have shown systematic PBPK modeling of renal passive reabsorption or CLr for hydrophobic permeable drugs in CKD. The goal of this study was to expand our previously developed and verified mechanistic kidney model to develop a universal model to predict changes in CLr in CKD for permeable and nonpermeable drugs that accounts for the dramatic nonlinear effect of CKD on renal passive reabsorption of permeable drugs. The developed model incorporates physiologically-based tubular changes of reduced water reabsorption/increased tubular flow rate per remaining functional nephron in CKD. The final adaptive kidney model successfully (absolute fold error (AFE) all < 2) predicted renal passive reabsorption and CLr for 20 permeable and nonpermeable test compounds across the stages of CKD. In contrast, use of proportional glomerular filtration rate reduction approach without addressing tubular adaptation processes in CKD to predict CLr generated unacceptable CLr predictions (AFE = 2.61-7.35) for permeable compounds in severe CKD. Finally, the adaptive kidney model accurately predicted CLr of para-amino-hippuric acid and memantine, two secreted compounds, in CKD, suggesting successful integration of active secretion into the model, along with passive reabsorption. In conclusion, the developed adaptive kidney model enables mechanistic predictions of in vivo CLr through CKD progression without any empirical scaling factors and can be used for CLr predictions prior to assessment of drug disposition in renal impairment.

A Microdose Cocktail to Evaluate Drug Interactions in Patients with Renal Impairment

Renal impairment (RI) is known to influence the pharmacokinetics of nonrenally eliminated drugs, although the mechanism and clinical impact is poorly understood. We assessed the impact of RI and single dose oral rifampin (RIF) on the pharmacokinetics of CYP3A, OATP1B, P-gp, and BCRP substrates using a microdose cocktail and OATP1B endogenous biomarkers. RI alone had no impact on midazolam (MDZ), maximum plasma concentration (C_max ), and area under the curve (AUC), but a progressive increase in AUC with RI severity for dabigatran (DABI), and up to ~2-fold higher AUC for pitavastatin (PTV), rosuvastatin (RSV), and atorvastatin (ATV) for all degrees of RI was observed. RIF did not impact MDZ, had a progressively smaller DABI drug-drug interaction (DDI) with increasing RI severity, a similar 3.1-fold to 4.4-fold increase in PTV and RSV AUC in healthy volunteers and patients with RI, and a diminishing DDI with RI severity from 6.1-fold to 4.7-fold for ATV. Endogenous biomarkers of OATP1B (bilirubin, coproporphyrin I/III, and sulfated bile salts) were generally not impacted by RI, and RIF effects on these biomarkers in RI were comparable or larger than those in healthy volunteers. The lack of a trend with RI severity of PTV and several OATP1B biomarkers, suggests that mechanisms beyond RI directly impacting OATP1B activity could also be considered. The DABI, RSV, and ATV data suggest an impact of RI on intestinal P-gp, and potentially BCRP activity. Therefore, DDI data from healthy volunteers may represent a worst-case scenario for clinically derisking P-gp and BCRP substrates in the setting of RI.

Pharmacokinetics in children with chronic kidney disease

In children, the main causes of chronic kidney disease (CKD) are congenital diseases and glomerular disorders. CKD is associated with multiple physiological changes and may therefore influence various pharmacokinetic (PK) parameters. A well-known consequence of CKD on pharmacokinetics is a reduction in renal clearance due to a decrease in the glomerular filtration rate. The impact of renal impairment on pharmacokinetics is, however, not limited to a decreased elimination of drugs excreted by the kidney. In fact, renal dysfunction may lead to modifications in absorption, distribution, transport, and metabolism as well. Currently, insufficient evidence is available to guide dosing decisions on many commonly used drugs. Moreover, the impact of maturation on drug disposition and action should be taken into account when selecting and dosing drugs in the pediatric population. Clinicians should take PK changes into consideration when selecting and dosing drugs in pediatric CKD patients in order to avoid toxicity and increase efficiency of drugs in this population. The aim of this review is to summarize known PK changes in relation to CKD and to extrapolate available knowledge to the pediatric CKD population to provide guidance for clinical practice.

Utilization of Physiologically Based Pharmacokinetic Modeling in Clinical Pharmacology and Therapeutics: an Overview

The purpose of this review was to assess the advancement of applications for physiologically based pharmacokinetic (PBPK) modeling in various therapeutic areas. We conducted a PubMed search, and 166 articles published between 2012 and 2018 on FDA-approved drug products were selected for further review. Qualifying publications were summarized according to therapeutic area, medication(s) studied, pharmacokinetic model type utilized, simulator program used, and the applications of that modeling. The results showed a 13-fold increase in the number of papers published from 2012 to 2018, with the largest proportion of articles dedicated to the areas of infectious diseases, oncology, and neurology, and application extensions including prediction of drug-drug interactions due to metabolism and/or transporter-mediated effects and understanding drug kinetics in special populations. In addition, we profiled several high-impact studies whose results were used to guide package insert information and formulate dose recommendations. These results show that while utilization of PBPK modeling has drastically increased over the past several years, regulatory support, lack of easy-to-use systems for clinicians, and challenges with model validation remain major challenges for the widespread adoption of this practice in institutional and ambulatory settings. However, PBPK modeling will continue to be a useful tool in the future to assess therapeutic drug monitoring and the growing field of personalized medicine.

Physiologically-Based Pharmacokinetic Models for Evaluating Membrane Transporter Mediated Drug-Drug Interactions: Current Capabilities, Case Studies, Future Opportunities, and Recommendations

Physiologically-based pharmacokinetic (PBPK) modeling has been extensively used to quantitatively translate in vitro data and evaluate temporal effects from drug-drug interactions (DDIs), arising due to reversible enzyme and transporter inhibition, irreversible time-dependent inhibition, enzyme induction, and/or suppression. PBPK modeling has now gained reasonable acceptance with the regulatory authorities for the cytochrome-P450-mediated DDIs and is routinely used. However, the application of PBPK for transporter-mediated DDIs (tDDI) in drug development is relatively uncommon. Because the predictive performance of PBPK models for tDDI is not well established, here, we represent and discuss examples of PBPK analyses included in regulatory submission (the US Food and Drug Administration (FDA), the European Medicines Agency (EMA), and the Pharmaceuticals and Medical Devices Agency (PMDA)) across various tDDIs. The goal of this collaborative effort (involving scientists representing 17 pharmaceutical companies in the Consortium and from academia) is to reflect on the use of current databases and models to address tDDIs. This challenges the common perceptions on applications of PBPK for tDDIs and further delves into the requirements to improve such PBPK predictions. This review provides a reflection on the current trends in PBPK modeling for tDDIs and provides a framework to promote continuous use, verification, and improvement in industrialization of the transporter PBPK modeling.

Decreased Disposition of Anticancer Drugs Predominantly Eliminated via the Liver in Patients with Renal Failure

Background:

Evidence has revealed that renal impairment can affect the systemic exposure of drugs which are predominantly eliminated via the liver. The modulation of drug-metabolizing enzymes and transporters expressed in the liver and/or small intestine by diverse entities, including uremic toxins, in systemic circulation of patients with severe renal failure is considered as the cause of atypical pharmacokinetics, which sometimes induce undesirable adverse events that are especially critical for drugs with narrow therapeutic window such as anticancer drugs. A dosing strategy for anticancer drugs in these patients needs to be established.

Methods:

The effects of renal impairment on the systemic exposure and safety of anticancer drugs were summarized. The proposed mechanisms for the alterations in the pharmacokinetics of these anticancer drugs were also discussed.

Results:

Changes in pharmacokinetics and clinical response were reported in 9 out of 10 cytotoxic anticancer drugs investigated, although available information was limited and sometimes controversial. Systemic exposure of 3 out of 16 tyrosine kinase inhibitors was higher in patients with severe renal failure than that in patients with normal kidney function. An increase in systemic exposure of anticancer drugs in patients with renal impairment is likely to be observed for substrates of OATP1B1, despite the limited evidence.

Conclusion:

The molecular basis for the effect of uremia on non-renal drug elimination still needed to be clarified with further studies to generate generalizable concepts, which may provide insights into establishing better clinical usage of anticancer drugs, i.e. identifying patients at risk and dose adjustment.

Development and Evaluation of Physiologically Based Pharmacokinetic Drug-Disease Models for Predicting Rifampicin Exposure in Tuberculosis and Cirrhosis Populations

The physiologically based pharmacokinetic (PBPK) approach facilitates the construction of novel drug–disease models by allowing incorporation of relevant pathophysiological changes. The aim of the present work was to explore and identify the differences in rifampicin pharmacokinetics (PK) after the application of its single dose in healthy and diseased populations by using PBPK drug–disease models. The Simcyp^® simulator was used as a platform for modeling and simulation. The model development process was initiated by predicting rifampicin PK in healthy population after intravenous (i.v) and oral administration. Subsequent to successful evaluation in healthy population, the pathophysiological changes in tuberculosis and cirrhosis population were incorporated into the developed model for predicting rifampicin PK in these populations. The model evaluation was performed by using visual predictive checks and the comparison of mean observed/predicted ratios (ratio_(Obs/pred)) of the PK parameters. The predicted PK parameters in the healthy population were in adequate harmony with the reported clinical data. The incorporation of pathophysiological changes in albumin concentration in the tuberculosis population revealed improved prediction of clearance. The developed PBPK drug–disease models have efficiently described rifampicin PK in tuberculosis and cirrhosis populations after administering single drug dose, as the ratio_(Obs/pred) for all the PK parameters were within a two-fold error range. The mechanistic nature of the developed PBPK models may facilitate their extension to other diseases and drugs.

Quantitative Translation of Microfluidic Transporter in Vitro Data to in Vivo Reveals Impaired Albumin-Facilitated Indoxyl Sulfate Secretion in Chronic Kidney Disease

Indoxyl sulfate (IxS), a highly albumin-bound uremic solute, accumulates in chronic kidney disease (CKD) due to reduced renal clearance. This study was designed to specifically investigate the role of human serum albumin (HSA) in IxS renal secretion via organic anion transporter 1 (OAT1) in a microfluidic system and subsequently apply quantitative translation of in vitro data to predict extent of change in IxS renal clearance in CKD stage IV relative to healthy. Conditionally immortalized human proximal tubule epithelial cells overexpressing OAT1 were incubated with IxS (5-200 μM) in the HSA-free medium or in the presence of either HSA or CKD-modified HSA. IxS uptake in the presence of HSA resulted in more than 20-fold decrease in OAT1 affinity (K_m,u) and 37-fold greater in vitro unbound intrinsic clearance (CL_int,u) versus albumin-free condition. In the presence of CKD-modified albumin, K_m,u increased four-fold and IxS CL_int,u decreased almost seven-fold relative to HSA. Fold-change in parameters exceeded differences in IxS binding between albumin conditions, indicating additional mechanism and facilitating role of albumin in IxS OAT1-mediated uptake. Quantitative translation of IxS in vitro OAT1-mediated CL_int,u predicted a 60% decrease in IxS renal elimination as a result of CKD, in agreement with the observed data (80%). The findings of the current study emphasize the role of albumin in IxS transport via OAT1 and explored the impact of modifications in albumin on renal excretion via active secretion in CKD. For the first time, this study performed quantitative translation of transporter kinetic data generated in a novel microfluidic in vitro system to a clinically relevant setting. Knowledge gaps and future directions in quantitative translation of renal drug disposition from microphysiological systems are discussed.

Integration of Food Animal Residue Avoidance Databank (FARAD) empirical methods for drug withdrawal interval determination with a mechanistic population-based interactive physiologically based pharmacokinetic (iPBPK) modeling platform: example for flunixin meglumine administration

Violative chemical residues in animal-derived food products affect food safety globally and have impact on the trade of international agricultural products. The Food Animal Residue Avoidance Databank program has been developing scientific tools to provide appropriate withdrawal interval (WDI) estimations after extralabel drug use in food animals for the past three decades. One of the tools is physiologically based pharmacokinetic (PBPK) modeling, which is a mechanistic-based approach that can be used to predict tissue residues and WDIs. However, PBPK models are complicated and difficult to use by non-modelers. Therefore, a user-friendly PBPK modeling framework is needed to move this field forward. Flunixin was one of the top five violative drug residues identified in the United States from 2010 to 2016. The objective of this study was to establish a web-based user-friendly framework for the development of new PBPK models for drugs administered to food animals. Specifically, a new PBPK model for both cattle and swine after administration of flunixin meglumine was developed. Population analysis using Monte Carlo simulations was incorporated into the model to predict WDIs following extralabel administration of flunixin meglumine. The population PBPK model was converted to a web-based interactive PBPK (iPBPK) framework to facilitate its application. This iPBPK framework serves as a proof-of-concept for further improvements in the future and it can be applied to develop new models for other drugs in other food animal species, thereby facilitating the application of PBPK modeling in WDI estimation and food safety assessment.