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.

Intestinal P-gp and Putative Hepatic OATP1B Induction: International Transporter Consortium Perspective on Drug Development Implications

There is an increasing interest in transporter induction (i.e., decreased systemic drug exposure due to increased efflux-limited absorption or transporter-mediated clearance) as a mechanism of drug-drug interactions (DDIs), although evidence of clinical relevance is still evolving. Intestinal P-glycoprotein (P-gp) and hepatic organic anion transporting polypeptides 1B (OATP1B) can be important determinants of drug absorption and disposition, as well as targets for DDIs. Current data indicate that intestinal P-gp protein levels can be induced up to threefold to fourfold in humans primarily with pregnane X receptor (PXR) activators, and that this induction can decrease the systemic exposure of drugs with P-gp efflux-limited absorption (e.g., ≤ 67% decrease in the exposure of total dabigatran following rifampin multiple oral dosing). Evaluation of the clinical relevance of P-gp induction as a DDI mechanism must consider the induction potential of the perpetrator drug for P-gp and attenuation of exposure of the victim drug in the context of its therapeutic window. Practical drug development recommendations are provided herein. Reports are contradictory on OATP1B induction by PXR activators in human hepatocytes and liver biopsies. Some clinical investigations demonstrated that rifampin pretreatment decreased exposure of OATP1B substrates, while other studies found no differences, and the potential involvement of other mechanisms in these observed DDIs cannot be definitively ruled out. Thus, further studies are needed to understand hepatic OATP1B induction and potential involvement of other mechanisms contributing to reduced exposure of OATP1B substrates. This review critically summarizes the state-of-the-art on intestinal P-gp and hepatic OATP1B induction, and highlights implications for drug development.

Utilising Magnetically Isolated Lysosomes for Direct Quantification of Intralysosomal Drug Concentrations by LC-MS/MS Analysis: An Investigatory Study With Imipramine

Lysosomes are acidic intracellular organelles that can extensively sequester basic lipophilic drugs due to pH and membrane partitioning, and therefore may significantly influence subcellular drug concentrations. Current in vitro methods for lysosomal drug sequestration evaluation typically lack the ability to accurately and sensitively quantify drug concentrations directly within the lysosome. In the current study, magnetic lysosomal isolation was used in the lysosome rich rat NR8383 cell line and combined with LC-MS/MS analysis to quantify intralysosomal concentrations and lysosomal partitioning (Kp_Lysosome) values of imipramine. The purity of the isolated lysosomes was validated by enzymatic and electron microscopy analysis. Lysosomal imipramine accumulation was explored using 2 methods: addition of imipramine to cells followed by lysosomal isolation (Method 1), and direct addition of imipramine to isolated lysosomes (Method 2). This work highlighted that both experimental buffers and ATP influence intralysosomal drug concentrations, and non-specific drug binding and re-distribution limits the use of Method 1. Method 2 may benefit future lysosomal drug accumulation studies, as imipramine demonstrated high Kp_Lysosome values (3500), comparable to in silico predictions. This study reports a novel method for the direct quantification of intralysosomal drug concentrations that has the ability to be adapted to other cell types.

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 Models as Framework for Understanding Biomarker Disposition: Prediction of Creatinine-Drug Interactions

Creatinine is widely used as a biomarker of glomerular filtration, and, hence, renal function. However, transporter-mediated secretion also contributes to its renal clearance, albeit to a lesser degree. Inhibition of these transporters causes transient serum creatinine elevation, which can be mistaken as impaired renal function. The current study developed mechanistic models of creatinine kinetics within physiologically based framework accounting for multiple transporters involved in creatinine renal elimination, assuming either unidirectional or bidirectional-OCT2 transport (driven by electrochemical gradient). Robustness of creatinine models was assessed by predicting creatinine-drug interactions with 10 perpetrators; performance evaluation accounted for 5% intra-individual variability in serum creatinine. Models showed comparable predictive performances of the maximum steady-state effect regardless of OCT2 directionality assumptions. However, only the bidirectional-OCT2 model successfully predicted the minimal effect of ranitidine. The dynamic nature of models provides clear advantage to static approaches and most advanced framework for evaluating interplay between multiple processes in creatinine renal disposition.

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.

A population pharmacokinetic model for simvastatin and its metabolites in children and adolescents

Purpose

Poor adherence to dietary/behaviour modifications as interventions for hypercholesterolemia in paediatric patients often necessitates the initiation of statin therapy. The aim of this study was to develop a joint population pharmacokinetic model for simvastatin and four metabolites in children and adolescents to investigate sources of variability in simvastatin acid exposure in this patient population, in addition to SLCO1B1 genotype status.

Methods

Plasma concentrations of simvastatin and its four metabolites, demographic and polymorphism data for OATP1B1 and CYP3A5 were analysed utilising a population pharmacokinetic modelling approach from an existing single oral dose (10 mg < 17 years and 20 mg ≥ 18 years) pharmacokinetic dataset of 32 children and adolescents.

Results

The population PK model included a one compartment disposition model for simvastatin with irregular oral absorption described by two parallel absorption processes each consisting of sequential zero and first-order processes. The data for each metabolite were described by a one-compartment disposition model with the formation and elimination apparent parameters estimated. The model confirmed the statistically significant effect of c.521T>C (rs4149056) on the pharmacokinetics of the active metabolite simvastatin acid in children/adolescents, consistent with adult data. In addition, age was identified as a covariate affecting elimination clearances of 6-hydroxymethyl simvastatin acid and 3, 5 dihydrodiol simvastatin metabolites.

Conclusion

The model developed describes the pharmacokinetics of simvastatin and its metabolites in children/adolescents capturing the effects of both c.521T>C and age on variability in exposure in this patient population. This joint simvastatin metabolite model is envisaged to facilitate optimisation of simvastatin dosing in children/adolescents.

Electronic supplementary material

The online version of this article (10.1007/s00228-019-02697-y) contains supplementary material, which is available to authorized users.

Hepatic Organic Anion Transporting Polypeptide-Mediated Clearance in the Beagle Dog: Assessing In Vitro-In Vivo Relationships and Applying Cross-Species Empirical Scaling Factors to Improve Prediction of Human Clearance

In the present study, the beagle dog was evaluated as a preclinical model to investigate organic anion transporting polypeptide (OATP)-mediated hepatic clearance. In vitro studies were performed with nine OATP substrates in three lots of plated male dog hepatocytes ± OATP inhibitor cocktail to determine total uptake clearance (CL_uptake) and total and unbound cell-to-medium concentration ratio (Kp_uu). In vivo intrinsic hepatic clearances (CL_int,H) were determined following intravenous drug administration (0.1 mg/kg) in male beagle dogs. The in vitro parameters were compared with those previously reported in plated human, monkey, and rat hepatocytes; the ability of cross-species scaling factors to improve prediction of human in vivo clearance was assessed. CL_uptake in dog hepatocytes ranged from 9.4 to 135 µl/min/10⁶ cells for fexofenadine and telmisartan, respectively. Active process contributed >75% to CL_uptake for 5/9 drugs. Rosuvastatin and valsartan showed Kp_uu > 10, whereas cerivastatin, pitavastatin, repaglinide, and telmisartan had Kp_uu < 5. The extent of hepatocellular binding in dog was consistent with other preclinical species and humans. The bias (2.73-fold) obtained from comparison of predicted versus in vivo dog CL_int,H was applied as an average empirical scaling factor (ESF_av) for in vitro-in vivo extrapolation of human CL_int,H The ESF_av based on dog reduced underprediction of human CL_int,H for the same data set (geometric mean fold error = 2.1), highlighting its utility as a preclinical model to investigate OATP-mediated uptake. The ESF_av from all preclinical species resulted in comparable improvement of human clearance prediction, in contrast to drug-specific empirical scalars, rationalized by species differences in expression and/or relative contribution of particular transporters to drug hepatic uptake.

Advancing Predictions of Tissue and Intracellular Drug Concentrations Using In Vitro, Imaging and Physiologically Based Pharmacokinetic Modeling Approaches

This white paper examines recent progress, applications, and challenges in predicting unbound and total tissue and intra/subcellular drug concentrations using in vitro and preclinical models, imaging techniques, and physiologically based pharmacokinetic (PBPK) modeling. Published examples, regulatory submissions, and case studies illustrate the application of different types of data in drug development to support modeling and decision making for compounds with transporter-mediated disposition, and likely disconnects between tissue and systemic drug exposure. The goals of this article are to illustrate current best practices and outline practical strategies for selecting appropriate in vitro and in vivo experimental methods to estimate or predict tissue and plasma concentrations, and to use these data in the application of PBPK modeling for human pharmacokinetic (PK), efficacy, and safety assessment in drug development.

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

Chronic kidney disease (CKD) differentially affects the pharmacokinetics (PK) of nonrenally cleared drugs via certain pathways (e.g., cytochrome P450 (CYP)2D6); however, the effect on CYP2C8‐mediated clearance is not well understood because of overlapping substrate specificity with hepatic organic anion‐transporting polypeptides (OATPs). This study used physiologically based pharmacokinetic (PBPK) modeling to delineate potential changes in CYP2C8 or OATP1B activity in patients with CKD. Drugs analyzed are predominantly substrates of CYP2C8 (rosiglitazone and pioglitazone), OATP1B (pitavastatin), or both (repaglinide). Following initial model verification, pharmacokinetics (PK) of these drugs were simulated in patients with severe CKD considering changes in glomerular filtration rate (GFR), plasma protein binding, and activity of either CYP2C8 and/or OATP1B in a stepwise manner. The PBPK analysis suggests that OATP1B activity could be decreased up to 60% in severe CKD, whereas changes to CYP2C8 are negligible. This improved understanding of CKD effect on clearance pathways could be important to inform the optimal use of nonrenally eliminated drugs in patients with CKD.

Comprehensive Evaluation of the Utility of 20 Endogenous Molecules as Biomarkers of OATP1B Inhibition Compared with Rosuvastatin and Coproporphyrin I

Endogenous biomarkers can be clinically relevant tools for the assessment of transporter function in vivo and corresponding drug-drug interactions (DDIs). The aim of this study was to perform systematic evaluation of plasma data obtained for 20 endogenous molecules in the same healthy subjects (n = 8-12) in the absence and presence of organic anion transporting polypeptide (OATP) inhibitor rifampicin (600 mg, single dose). The extent of rifampicin DDI magnitude [the ratio of the plasma concentration-time area under the curve (AUCR)], estimated fraction transported (f_T), and baseline variability was compared across the biomarkers and relative to rosuvastatin and coproporphyrin I (CPI). Out of the 20 biomarkers investigated tetradecanedioate (TDA), hexadecanedioate (HDA), glycocholic acid, glycodeoxycholic acid (GDCA), taurodeoxycholic acid (TDCA), and coproporphyrin III (CPIII) showed the high AUCR (2.1-8.5) and f_T (0.5-0.76) values, indicative of substantial OATP1B-mediated transport. A significant positive correlation was observed between the individual GDCA and TDCA AUCRs and the magnitude of rosuvastatin-rifampicin interaction. The CPI and CPIII AUCRs were significantly correlated, but no clear trend was established with the rosuvastatin AUCR. Moderate interindividual variability (15%-62%) in baseline exposure and AUCR was observed for TDA, HDA, and CPIII. In contrast, bile acids demonstrated high interindividual variability (69%-113%) and significant decreases in baseline plasma concentrations during the first 4 hours. This comprehensive analysis in the same individuals confirms that none of the biomarkers supersede CPI in the evaluation of OATP1B-mediated DDI risk. Monitoring of CPI and GDCA/TDCA may be beneficial for dual OATP1B/sodium-taurocholate cotransporting polypeptide inhibitors with consideration of challenges associated with large inter- and intraindividual variability observed for bile acids. Benefit of monitoring combined biomarkers (CPI, one bile acid and one fatty acid) needs to be confirmed with larger data sets and against multiple OATP1B clinical probes and perpetrators.

Influence of Transporter Polymorphisms on Drug Disposition and Response: A Perspective From the International Transporter Consortium

Advances in genomic technologies have led to a wealth of information identifying genetic polymorphisms in membrane transporters, specifically how these polymorphisms affect drug disposition and response. This review describes the current perspective of the International Transporter Consortium (ITC) on clinically important polymorphisms in membrane transporters. ITC suggests that, in addition to previously recommended polymorphisms in ABCG2 (BCRP) and SLCO1B1 (OATP1B1), polymorphisms in the emerging transporter, SLC22A1 (OCT1), be considered during drug development. Collectively, polymorphisms in these transporters are important determinants of interindividual differences in the levels, toxicities, and response to many drugs.

ITC Commentary on Metformin Clinical Drug-Drug Interaction Study Design That Enables an Efficacy- and Safety-Based Dose Adjustment Decision

Metformin drug-drug interaction (DDI) studies are conducted during development of drugs that inhibit organic cation transporters and/or multidrug and toxin extrusion proteins (OCTs/MATEs). Monitoring solely changes in systemic exposure, the typical DDI study endpoint appears inadequate for metformin, which is metabolically stable, has poor passive membrane permeability, and undergoes transporter-mediated tissue distribution and clearance. Evaluation of renal clearance, antihyperglycemic effects, and potentially lactate as an exploratory safety marker, can support rational metformin dose adjustment. The proposed DDI study design aims to adequately inform metformin dosing during comedication.

In Vitro-In Vivo Extrapolation of OATP1B-Mediated Drug-Drug Interactions in Cynomolgus Monkey

Hepatic organic anion-transporting polypeptides (OATP) 1B1 and 1B3 are clinically relevant transporters associated with significant drug-drug interactions (DDIs) and safety concerns. Given that OATP1Bs in cynomolgus monkey share >90% degree of gene and amino acid sequence homology with human orthologs, we evaluated the in vitro-in vivo translation of OATP1B-mediated DDI risk using this preclinical model. In vitro studies using plated cynomolgus monkey hepatocytes showed active uptake K_m values of 2.0 and 3.9 µM for OATP1B probe substrates, pitavastatin and rosuvastatin, respectively. Rifampicin inhibited pitavastatin and rosuvastatin active uptake in monkey hepatocytes with IC₅₀ values of 3.0 and 0.54 µM, respectively, following preincubation with the inhibitor. Intravenous pharmacokinetics of ²H₄-pitavastatin and ²H₆-rosuvastatin (0.2 mg/kg) and the oral pharmacokinetics of cold probes (2 mg/kg) were studied in cynomolgus monkeys (n = 4) without or with coadministration of single oral ascending doses of rifampicin (1, 3, 10, and 30 mg/kg). A rifampicin dose-dependent reduction in i.v. clearance of statins was observed. Additionally, oral pitavastatin and rosuvastatin plasma exposure increased up to 19- and 15-fold at the highest dose of rifampicin, respectively. Use of in vitro IC₅₀ obtained following 1 hour preincubation with rifampicin (0.54 µM) predicted correctly the change in mean i.v. clearance and oral exposure of statins as a function of mean unbound maximum plasma concentration of rifampicin. This study demonstrates quantitative translation of in vitro OATP1B IC₅₀ to predict DDIs using cynomolgus monkey as a preclinical model and provides further confidence in application of in vitro hepatocyte data for the prediction of clinical OATP1B-mediated DDIs.