Development of a Dynamic Physiologically Based Mechanistic Kidney Model to Predict Renal Clearance

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.

The impact of kidney function on plasma neurofilament light and phospho-tau 181 in a community-based cohort: the Shanghai Aging Study

BACKGROUND

The blood-based biomarkers are approaching the clinical practice of Alzheimer's disease (AD). Chronic kidney disease (CKD) has a potential confounding effect on peripheral protein levels. It is essential to characterize the impact of renal function on AD markers.

METHODS

Plasma phospho-tau181 (P-tau181), and neurofilament light (NfL) were assayed via the Simoa HD-X platform in 1189 dementia-free participants from the Shanghai Aging Study (SAS). The estimated glomerular filter rate (eGFR) was calculated. The association between renal function and blood NfL, P-tau181 was analyzed. An analysis of interactions between various demographic and comorbid factors and eGFR was conducted.

RESULTS

The eGFR levels were negatively associated with plasma concentrations of NfL and P-tau181 (B = - 0.19, 95% CI - 0.224 to - 0.156, P < 0.001; B = - 0.009, 95% CI - 0.013 to -0.005, P < 0.001, respectively). After adjusting for demographic characteristics and comorbid diseases, eGFR remained significantly correlated with plasma NfL (B = - 0.010, 95% CI - 0.133 to - 0.068, P < 0.001), but not with P-tau181 (B = - 0.003, 95% CI - 0.007 to 0.001, P = 0.194). A significant interaction between age and eGFR was found for plasma NfL (Pinteraction < 0.001). In participants ≥ 70 years and with eGFR < 60 ml/min/1.73 m2, the correlation between eGFR and plasma NfL was significantly remarkable (B = - 0.790, 95% CI - 1.026 to - 0,554, P < 0.001).

CONCLUSIONS

Considering renal function and age is crucial when interpreting AD biomarkers in the general aging population.

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.

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.

Utility of Minimal Physiologically Based Pharmacokinetic Models for Assessing Fractional Distribution, Oral Absorption, and Series-Compartment Models of Hepatic Clearance

Minimal physiologically based pharmacokinetic (mPBPK) models are physiologically relevant, require less information than full PBPK models, and offer flexibility in pharmacokinetics (PK). The well-stirred hepatic model (WSM) is commonly used in PBPK, whereas the more plausible dispersion model (DM) poses computational complexities. The series-compartment model (SCM) mimics the DM but is easier to operate. This work implements the SCM and mPBPK models for assessing fractional tissue distribution, oral absorption, and hepatic clearance using literature-reported blood and liver concentration-time data in rats for compounds mainly cleared by the liver. Further handled were various complexities, including nonlinear hepatic binding and metabolism, differing absorption kinetics, and sites of administration. The SCM containing one to five (n) liver subcompartments yields similar fittings and provides comparable estimates for hepatic extraction ratio (ER), prehepatic availability (Fg ), and first-order absorption rate constants (ka ). However, they produce decreased intrinsic clearances (CLint ) and liver-to-plasma partition coefficients (Kph ) with increasing n as expected. Model simulations demonstrated changes in intravenous and oral PK profiles with alterations in Kph and ka and with hepatic metabolic zonation. The permeability (PAMPA P) of the various compounds well explained the fitted fractional distribution (fd ) parameters. The SCM and mPBPK models offer advantages in distinguishing systemic, extrahepatic, and hepatic clearances. The SCM allows for incorporation of liver zonation and is useful in assessing changes in internal concentration gradients potentially masked by similar blood PK profiles. Improved assessment of intraorgan drug concentrations may offer insights into active moieties driving metabolism, biliary excretion, pharmacodynamics, and hepatic toxicity. SIGNIFICANCE STATEMENT: The minimal physiologically based pharmacokinetic model and the series-compartment model are useful in assessing oral absorption and hepatic clearance. They add flexibility in accounting for various drug- or system-specific complexities, including fractional distribution, nonlinear binding and saturable hepatic metabolism, and hepatic zonation. These models can offer improved insights into the intraorgan concentrations that reflect physiologically active moieties often driving disposition, pharmacodynamics, and toxicity.

The Role of "Physiologically Based Pharmacokinetic Model (PBPK)" New Approach Methodology (NAM) in Pharmaceuticals and Environmental Chemical Risk Assessment

Physiologically Based Pharmacokinetic (PBPK) models are mechanistic tools generally employed in the pharmaceutical industry and environmental health risk assessment. These models are recognized by regulatory authorities for predicting organ concentration-time profiles, pharmacokinetics and daily intake dose of xenobiotics. The extension of PBPK models to capture sensitive populations such as pediatric, geriatric, pregnant females, fetus, etc., and diseased populations such as those with renal impairment, liver cirrhosis, etc., is a must. However, the current modelling practices and existing models are not mature enough to confidently predict the risk in these populations. A multidisciplinary collaboration between clinicians, experimental and modeler scientist is vital to improve the physiology and calculation of biochemical parameters for integrating knowledge and refining existing PBPK models. Specific PBPK covering compartments such as cerebrospinal fluid and the hippocampus are required to gain mechanistic understanding about xenobiotic disposition in these sub-parts. The PBPK model assists in building quantitative adverse outcome pathways (qAOPs) for several endpoints such as developmental neurotoxicity (DNT), hepatotoxicity and cardiotoxicity. Machine learning algorithms can predict physicochemical parameters required to develop in silico models where experimental data are unavailable. Integrating machine learning with PBPK carries the potential to revolutionize the field of drug discovery and development and environmental risk. Overall, this review tried to summarize the recent developments in the in-silico models, building of qAOPs and use of machine learning for improving existing models, along with a regulatory perspective. This review can act as a guide for toxicologists who wish to build their careers in kinetic modeling.

Integration of Engineered Delivery with the Pharmacokinetics of Medical Candidates via Physiology-Based Pharmacokinetics

Physiologically based pharmacokinetic (PBPK) modeling is a mechanistic computational model that can be used to predict a drug product's ADME (absorption, distribution, metabolism, and excretion) and pharmacokinetics (PK). In recent years, PBPK modeling and simulation has been used increasingly to address many biopharmaceutics and clinical pharmacology questions, such as the effect of formulations, intrinsic factors (age, organ dysfunction, etc.), and extrinsic factors (comedications, food) on the PK of an investigational drug product. In this chapter, we will briefly introduce various PBPK models for ADME prediction and general procedures for PBPK modeling and simulations. The readers are encouraged to read updated literature on new applications of PBPK modeling and simulation which is still an emerging area in pharmaceutical development.

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.

Crosstalk of physiological pH and chemical pKa under the umbrella of physiologically based pharmacokinetic modeling of drug absorption, distribution, metabolism, excretion, and toxicity

Introduction: Physiological pH and chemical pKa are two sides of the same coin in defining the ionization of a drug in the human body. The Henderson-Hasselbalch equation and pH-partition hypothesis form the theoretical base to define the impact of pH-pKa crosstalk on drug ionization and thence its absorption, distribution, metabolism, excretion, and toxicity (ADMET).Areas covered: Human physiological pH is not constant, but a diverse, dynamic state regulated by various biological mechanisms, while the chemical pKa is generally a constant defining the acidic dissociation of the drug at various environmental pH. Works on pH-pKa crosstalk are scattered in the literature, despite its significant contributions to drug pharmacokinetics, pharmacodynamics, safety, and toxicity. In particular, its impacts on drug ADMET have not been effectively linked to the physiologically based pharmacokinetic (PBPK) modeling and simulation, a powerful tool increasingly used in model-informed drug development (MIDD).Expert opinion: Lacking a full consideration of the interactions of physiological pH and chemical pKa in a PBPK model limits scientists' capability in mechanistically describing the drug ADMET. This mini-review compiled literature knowledge on pH-pKa crosstalk and its impacts on drug ADMET, from the viewpoint of PBPK modeling, to pave the way to a systematic incorporation of pH-pKa crosstalk into PBPK modeling and simulation.

Incorporating renal excretion via the OCT2 transporter in physiologically based kinetic modelling to predict in vivo kinetics of mepiquat in rat

The present study aimed at incorporating active renal excretion via the organic cation transporter 2 (OCT2) into a generic rat physiologically based kinetic (PBK) model using an in vitro human renal proximal tubular epithelial cell line (SA7K) and mepiquat chloride (MQ) as the model compound. The Vmax (10.5 pmol/min/mg protein) and Km (20.6 μM) of OCT2 transport of MQ were determined by concentration-dependent uptake in SA7K cells using doxepin as inhibitor. PBK model predictions incorporating these values in the PBK model were 6.7-8.4-fold different from the reported in vivo data on the blood concentration of MQ in rat. Applying an overall scaling factor that also corrects for potential differences in OCT2 activity in the SA7K cells and in vivo kidney cortex and species differences resulted in adequate predictions for in vivo kinetics of MQ in rat (2.3-3.2-fold). The results indicate that using SA7K cells to define PBK parameters for active renal OCT2 mediated excretion with adequate scaling enables incorporation of renal excretion via the OCT2 transporter in PBK modelling to predict in vivo kinetics of mepiquat in rat. This study demonstrates a proof-of-principle on how to include active renal excretion into generic PBK models.

Diagnosing intake and rationalizing toxicities associated with 5F-MDMB-PINACA and 4F-MDMB-BINACA abuse

5F-MDMB-PINACA and 4F-MDMB-BINACA are synthetic cannabinoids (SCs) that elicit cannabinoid psychoactive effects. Defining pharmacokinetic-pharmacodynamic (PK-PD) relationships governing SCs and their metabolites are paramount to investigating their in vivo toxicological outcomes. However, the disposition kinetics and cannabinoid receptor (CB) activities of the primary metabolites of SCs are largely unknown. Additionally, reasons underlying the selection of ester hydrolysis metabolites (EHMs) as urinary biomarkers are often unclear. Here, metabolic reaction phenotyping was performed to identify key metabolizing enzymes of the parent SCs. Hepatic clearances of parent SCs and their EHMs were estimated from microsomal metabolic stability studies. Renal clearances were simulated using a mechanistic kidney model incorporating in vitro permeability and organic anionic transporter 3 (OAT3)-mediated uptake data. Overall clearances were considered in tandem with estimated volumes of distribution for in vivo biological half-lives (t_1/2) predictions. Interactions of the compounds with CB1 and CB2 were investigated using a G-protein coupled receptor activation assay. We demonstrated that similar enzymatic isoforms were implicated in the metabolism of 5F-MDMB-PINACA and 4F-MDMB-BINACA. Our in vivo t_1/2 determinations verified the rapid elimination of parent SCs and suggest prolonged circulation of their EHMs. The pronounced attenuation of the potencies and efficacies of the metabolites against CB1 and CB2 further suggests how toxic manifestations of SC abuse are likely precipitated by augmented exposure to parent SCs. Notably, basolateral OAT3-mediated uptake of the EHMs substantiates their higher urinary abundance. These novel insights underscore the importance of mechanistic, quantitative and systematic characterization of PK-PD relationships in rationalizing the toxicities of SCs.

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.

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.

Physiologically Based Pharmacokinetic Modelling for First-In-Human Predictions: An Updated Model Building Strategy Illustrated with Challenging Industry Case Studies

Physiologically based pharmacokinetic modelling is well established in the pharmaceutical industry and is accepted by regulatory agencies for the prediction of drug-drug interactions. However, physiologically based pharmacokinetic modelling is valuable to address a much wider range of pharmaceutical applications, and new regulatory impact is expected as its full power is leveraged. As one example, physiologically based pharmacokinetic modelling is already routinely used during drug discovery for in-vitro to in-vivo translation and pharmacokinetic modelling in preclinical species, and this leads to the application of verified models for first-in-human pharmacokinetic predictions. A consistent cross-industry strategy in this application area would increase confidence in the approach and facilitate further learning. With this in mind, this article aims to enhance a previously published first-in-human physiologically based pharmacokinetic model-building strategy. Based on the experience of scientists from multiple companies participating in the GastroPlus™ User Group Steering Committee, new Absorption, Distribution, Metabolism and Excretion knowledge is integrated and decision trees proposed for each essential component of a first-in-human prediction. We have reviewed many relevant scientific publications to identify new findings and highlight gaps that need to be addressed. Finally, four industry case studies for more challenging compounds illustrate and highlight key components of the strategy.

A Novel Physiologically Based Model of Creatinine Renal Disposition to Integrate Current Knowledge of Systems Parameters and Clinical Observations

Creatinine is the most common clinical biomarker of renal function. As a substrate for renal transporters, its secretion is susceptible to inhibition by drugs, resulting in transient increase in serum creatinine and false impression of damage to kidney. Novel physiologically based models for creatinine were developed here and (dis)qualified in a stepwise manner until consistency with clinical data. Data from a matrix of studies were integrated, including systems data (common to all models), proteomics-informed in vitro-in vivo extrapolation of all relevant transporter clearances, exogenous administration of creatinine (to estimate endogenous synthesis rate), and inhibition of different renal transporters (11 perpetrator drugs considered for qualification during creatinine model development and verification on independent data sets). The proteomics-informed bottom-up approach resulted in the underprediction of creatinine renal secretion. Subsequently, creatinine-trimethoprim clinical data were used to inform key model parameters in a reverse translation manner, highlighting best practices and challenges for middle-out optimization of mechanistic models.

Mechanistic PBPK Modeling of Urine pH Effect on Renal and Systemic Disposition of Methamphetamine and Amphetamine

The effect of urine pH on renal excretion and systemic disposition has been observed for many drugs and metabolites. When urine pH is altered, tubular ionization, passive reabsorption, renal clearance, and systemic exposure of drugs and metabolites may all change dramatically, raising clinically significant concerns. Surprisingly, the urine pH effect on drug disposition is not routinely explored in humans, and regulatory agencies have neither developed guidance on this issue nor required industry to conduct pertinent human trials. In this study, we hypothesized that physiologically based pharmacokinetic (PBPK) modeling could be used as a cost-effective method to examine potential urine pH effect on drug and metabolite disposition. Our previously developed and verified mechanistic kidney model was integrated with a full-body PBPK model to simulate renal clearance and area under the plasma concentration-time curve (AUC) with varying urine pH statuses using methamphetamine and amphetamine as model compounds. We first developed and verified drug models for methamphetamine and amphetamine under normal urine pH condition [absolute average fold error (AAFE) < 1.25 at study level]. Then, acidic and alkaline urine scenarios were simulated. Our simulation results show that the renal excretion and plasma concentration-time profiles for methamphetamine and amphetamine could be recapitulated under different urine pH (AAFE < 2 at individual level). The methamphetamine-amphetamine parent-metabolite full-body PBPK model also successfully simulated amphetamine plasma concentration-time profiles (AAFE < 1.25 at study level) and amphetamine/methamphetamine urinary concentration ratios (AAFE < 2 at individual level) after dosing methamphetamine. This demonstrates that our mechanistic PBPK model can predict urine pH effect on systemic and urinary disposition of drugs and metabolites. SIGNIFICANCE STATEMENT: Our study shows that integrating mechanistic kidney model with full-body physiologically based pharmacokinetic model can predict the magnitude of alteration in renal excretion and area under the plasma concentration-time curve (AUC) of drugs and metabolites when urine pH is changed. This provides a cost-effective method to evaluate the likelihood of renal and systemic disposition changes due to varying urine pH. This is important because multiple drugs and diseases can alter urine pH, leading to quantitatively and clinically significant changes in drug and metabolite disposition that may require adjustment of therapy.

A Physiologically Based in Silico Tool to Assess the Risk of Drug-Related Crystalluria

Drug precipitation in the nephrons of the kidney can cause drug-induced crystal nephropathy (DICN). To aid mitigation of this risk in early drug discovery, we developed a physiologically based in silico model to predict DICN in rats, dogs, and humans. At a minimum, the likelihood of DICN is determined by the level of systemic exposure to the molecule, the molecule's physicochemical properties and the unique physiology of the kidney. Accordingly, the proposed model accounts for these properties in order to predict drug exposure relative to solubility along the nephron. Key physiological parameters of the kidney were codified in a manner consistent with previous reports. Quantitative structure-activity relationship models and in vitro assays were used to estimate drug-specific physicochemical inputs to the model. The proposed model was calibrated against urinary excretion data for 42 drugs, and the utility for DICN prediction is demonstrated through application to 20 additional drugs.

Predicting tubular reabsorption with a human kidney proximal tubule tissue-on-a-chip and physiologically-based modeling

Kidney is a major route of xenobiotic excretion, but the accuracy of preclinical data for predicting in vivo clearance is limited by species differences and non-physiologic 2D culture conditions. Microphysiological systems can potentially increase predictive accuracy due to their more realistic 3D environment and incorporation of dynamic flow. We used a renal proximal tubule microphysiological device to predict renal reabsorption of five compounds: creatinine (negative control), perfluorooctanoic acid (positive control), cisplatin, gentamicin, and cadmium. We perfused compound-containing media to determine renal uptake/reabsorption, adjusted for non-specific binding. A physiologically-based parallel tube model was used to model reabsorption kinetics and make predictions of overall in vivo renal clearance. For all compounds tested, the kidney tubule chip combined with physiologically-based modeling reproduces qualitatively and quantitatively in vivo tubular reabsorption and clearance. However, because the in vitro device lacks filtration and tubular secretion components, additional information on protein binding and the importance of secretory transport is needed in order to make accurate predictions. These and other limitations, such as the presence of non-physiological compounds such as antibiotics and bovine serum albumin in media and the need to better characterize degree of expression of important transporters, highlight some of the challenges with using microphysiological devices to predict in vivo pharmacokinetics.

Development of an in silico prediction system of human renal excretion and clearance from chemical structure information incorporating fraction unbound in plasma as a descriptor

Prediction of pharmacokinetic profiles of new chemical entities is essential in drug development to minimize the risks of potential withdrawals. The excretion of unchanged compounds by the kidney constitutes a major route in drug elimination and plays an important role in pharmacokinetics. Herein, we created in silico prediction models of the fraction of drug excreted unchanged in the urine (f_e) and renal clearance (CL_r), with datasets of 411 and 401 compounds using freely available software; notably, all models require chemical structure information alone. The binary classification model for f_e demonstrated a balanced accuracy of 0.74. The two-step prediction system for CL_r was generated using a combination of the classification model to predict excretion-type compounds and regression models to predict the CL_r value for each excretion type. The accuracies of the regression models increased upon adding a descriptor, which was the observed and predicted fraction unbound in plasma (f_u,p); 78.6% of the samples in the higher range of renal clearance fell within 2-fold error with predicted f_u,p value. Our prediction system for renal excretion is freely available to the public and can be used as a practical tool for prioritization and optimization of compound synthesis in the early stage of drug discovery.

Sampling Site Has a Critical Impact on Physiologically Based Pharmacokinetic Modeling

It has been shown that arterial (central) and venous (peripheral) plasma drug concentrations can be very different. While pharmacokinetic studies typically measure drug concentrations from the peripheral vein such as the arm vein, physiologically based pharmacokinetic (PBPK) models generally output simulated concentrations from the central venous compartment that physiologically represents the right atrium, a merge of the superior and inferior vena cava. In this study, a physiologically based peripheral forearm sampling site model was developed and verified using nicotine, ketamine, lidocaine, and fentanyl as model drugs. This verified model allows output of simulated peripheral venous concentrations that can be meaningfully compared with observed pharmacokinetic data from the arm vein. The generalized effect of PBPK model sampling site on simulation output was investigated. Drugs and metabolites with large volumes of distribution showed considerable concentration discrepancy between the simulated central venous compartment and the peripheral arm vein after intravenous or oral administration, resulting in significant differences in values for C _max and time taken to reach C _max (t _max ) In addition, the simulated central venous metabolite profile showed an unexpected profile that was not observed in the peripheral arm vein. Using fentanyl as a model compound, we show that using the wrong sampling site in PBPK models can lead to biased model evaluation and subsequent erroneous model parameter optimization. Such an error in model parameters along with the discrepant sampling site could dramatically mislead the pharmacokinetic prediction in unstudied clinical scenarios, affecting the assessment of drug safety and efficacy. Overall, this study shows that PBPK model publications should specify the model sampling sites and match them with those employed in clinical studies. SIGNIFICANCE STATEMENT: Our study shows that sampling from the central venous compartment (right atrium) during physiologically based pharmacokinetic model development gives rise to biased model evaluation and erroneous model parameterization when observed data are collected from the peripheral arm vein. This can lead to a clinically significant error in predictions of plasma concentration-time profiles in unstudied scenarios. To address this error, we developed and verified a novel peripheral sampling site model to simulate arm vein drug concentrations that can be applied to different drug dosing scenarios.

Emerging Role of Organ-on-a-Chip Technologies in Quantitative Clinical Pharmacology Evaluation

The recently enacted Prescription Drug User Fee Act (PDUFA) VI includes in its performance goals "enhancing regulatory science and expediting drug development." The key elements in "enhancing regulatory decision tools to support drug development and review" include "advancing model-informed drug development (MIDD)." This paper describes (i) the US Food and Drug Administration (FDA) Office of Clinical Pharmacology's continuing efforts in developing quantitative clinical pharmacology models (disease, drug, and clinical trial models) to advance MIDD, (ii) how emerging novel tools, such as organ-on-a-chip technologies or microphysiological systems, can provide new insights into physiology and disease mechanisms, biomarker identification and evaluation, and elucidation of mechanisms of adverse drug reactions, and (iii) how the single organ or linked organ microphysiological systems can provide critical system parameters for improved physiologically-based pharmacokinetic and pharmacodynamic evaluations. Continuous public-private partnerships are critical to advance this field and in the application of these new technologies in drug development and regulatory review.