Semi-mechanistic kidney model incorporating physiologically-relevant fluid reabsorption and transporter-mediated renal reabsorption: pharmacokinetics of γ-hydroxybutyric acid and L-lactate in rats

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.

Effect of Sex and Cross-Sex Hormone Treatment on Renal Monocarboxylate-Transporter Expression in Rats

Proton- and sodium-dependent monocarboxylate transporters (MCTs/SMCTs) are determinants of renal clearance through the renal reabsorption of monocarboxylate substrates. Prior studies with intact females and males, ovariectomized females and castrated males have revealed the hormonal regulation of renal monocarboxylate-transporter expression, prompting investigation into the regulatory role of individual hormones. The aim of the present study is to evaluate the effect of exogenous sex and cross-sex hormones on renal MCT1, MCT4, CD147 and SMCT1 mRNA and membrane-bound protein expression. Ovariectomized (OVX) females and castrated (CST) male Sprague Dawley rats received estrogen and/or progesterone, testosterone, or a corresponding placebo treatment for 21 days prior to kidney collection. The quantitative measurement of mRNA and membrane-protein levels were conducted using qPCR and Western blot. Quantitative analysis revealed the combination estrogen/progesterone treatment reduced membrane MCT1 and 4 expression and increased SMCT1 expression, while testosterone administration increased MCT1 membrane-protein expression. Correlation analysis indicated that plasma 17β-estradiol was negatively correlated with MCT1 and MCT4 membrane expression, while testosterone was positively correlated. In contrast, SMCT1 membrane expression was positively correlated with 17β-estradiol and progesterone concentrations. MCT1, MCT4, CD147 and SMCT1 renal expression are significantly altered in response to female and male sex hormones following sex and cross-sex hormone treatment in OVX and CST rats. Further studies are needed to understand the complex role of sex hormones, sex hormone receptors and the impact of puberty on MCT/SMCT regulation.

Effects of 1α,25-Dihydroxyvitamin D3 on the Pharmacokinetics of Procainamide and Its Metabolite N-Acetylprocainamide, Organic Cation Transporter Substrates, in Rats with PBPK Modeling Approach

In this study, possible changes in the expression of rat organic cationic transporters (rOCTs) and rat multidrug and toxin extrusion proteins (rMATEs) following treatment with 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) were investigated. Rats received intraperitoneal administrations of 1,25(OH)₂D₃ for four consecutive days, and the tissues of interest were collected. The mRNA expression of rOCT1 in the kidneys was significantly increased in 1,25(OH)₂D₃-treated rats compared with the control rats, while the mRNA expressions of rOCT2 and rMATE1 in the kidneys, rOCT1 and N-acetyltransferase-II (NAT-II) in the liver, and rOCT3 in the heart were significantly decreased. Changes in the protein expression of hepatic rOCT1 and renal rOCT2 and rMATE1 were confirmed by western blot analysis. We further evaluated the pharmacokinetics of procainamide (PA) hydrochloride and its major metabolite N-acetyl procainamide (NAPA) in the presence of 1,25(OH)₂D₃. When PA hydrochloride was administered intravenously at a dose 10 mg/kg to 1,25(OH)₂D₃-treated rats, a significant decrease in renal and/or non-renal clearance of PA and NAPA was observed. A physiological model for the pharmacokinetics of PA and NAPA in rats was useful for linking changes in the transcriptional and translational expressions of rOCTs and rMATE1 transporters to the altered pharmacokinetics of the drugs.

Physiologically-Based Pharmacokinetic Modeling for Drug-Drug Interactions of Procainamide and N-Acetylprocainamide with Cimetidine, an Inhibitor of rOCT2 and rMATE1, in Rats

Previous observations demonstrated that cimetidine decreased the clearance of procainamide (PA) and/or N-acetylprocainamide (NAPA; the primary metabolite of PA) resulting in the increased systemic exposure and the decrease of urinary excretion. Despite an abundance of in vitro and in vivo data regarding pharmacokinetic interactions between PA/NAPA and cimetidine, however, a mechanistic approach to elucidate these interactions has not been reported yet. The primary objective of this study was to construct a physiological model that describes pharmacokinetic interactions between PA/NAPA and cimetidine, an inhibitor of rat organic cation transporter 2 (rOCT2) and rat multidrug and toxin extrusion proteins (rMATE1), by performing extensive in vivo and in vitro pharmacokinetic studies for PA and NAPA performed in the absence or presence of cimetidine in rats. When a single intravenous injection of PA HCl (10 mg/kg) was administered to rats, co-administration of cimetidine (100 mg/kg) significantly increased systemic exposure and decreased the systemic (CL) and renal (CL_R) clearance of PA, and reduced its tissue distribution. Similarly, cimetidine significantly decreased the CL_R of NAPA formed by the metabolism of PA and increased the AUC of NAPA. Considering that these drugs could share similar renal secretory pathways (e.g., via rOCT2 and rMATE1), a physiologically-based pharmacokinetic (PBPK) model incorporating semi-mechanistic kidney compartments was devised to predict drug-drug interactions (DDIs). Using our proposed PBPK model, DDIs between PA/NAPA and cimetidine were successfully predicted for the plasma concentrations and urinary excretion profiles of PA and NAPA observed in rats. Moreover, sensitivity analyses of the pharmacokinetics of PA and NAPA showed the inhibitory effects of cimetidine via rMATE1 were probably important for the renal elimination of PA and NAPA in rats. The proposed PBPK model may be useful for understanding the mechanisms of interactions between PA/NAPA and cimetidine in vivo.

Effects of renal impairment on transporter-mediated renal reabsorption of drugs and renal drug-drug interactions: A simulation-based study

Renal impairment (RI) significantly impacts the clearance of drugs through changes in the glomerular filtration rate, protein binding and alterations in the expression of renal drug transport proteins and hepatic metabolizing enzymes. The objectives of this study were to evaluate quantitatively the effects of renal impairment on the pharmacokinetics of drugs undergoing renal transporter-mediated reabsorption. A previously published semi-mechanistic kidney model incorporating physiologically relevant fluid reabsorption and transporter-mediated active renal reabsorption (PMID: 26341876) was utilized in this study. The probe drug/transporter pair utilized was γ-hydroxybutyric acid (GHB) and monocarboxylate transporter 1 (SCL16A1, MCT1). γ-Hydroxybutyric acid concentrations in the blood and amount excreted into urine were simulated using ADAPT 5 for the i.v. dose range of 200-1500 mg/kg in rats and the impact of renal impairment on CL_R and AUC was evaluated. A 90% decrease in GFR resulted in a > 100-fold decrease in GHB CL_R . When expression of reabsorptive transporters was decreased and f_u was increased, CL_R approached GFR. The effect of renal impairment on CL_R was reduced when the expression of drug metabolizing enzymes (DME) was increased as a result of increased metabolic clearance; the converse held true when the DME expression was decreased. In conclusion, this study quantitatively demonstrated that the effects of renal insufficiency on the clearance of drugs is modulated by transporter expression, contribution of renal clearance to overall clearance, expression of drug metabolizing enzymes, fraction unbound and drug-drug interactions with inhibitors of renal transporters that may be increased in the presence of renal impairment.

Toxicologic/transport properties of NCS-382, a γ-hydroxybutyrate (GHB) receptor ligand, in neuronal and epithelial cells: Therapeutic implications for SSADH deficiency, a GABA metabolic disorder

We report the in vitro assessment of pharmacotoxicity for the high-affinity GHB receptor ligand, NCS-382, using neuronal stem cells derived from mice with a targeted deletion of the aldehyde dehydrogenase 5a1 gene (succinic semialdehyde dehydrogenase(SSADH)-deficient mice). These animals represent a phenocopy of the human disorder of GABA metabolism, SSADH deficiency, that metabolically features accumulation of both GABA and the GABA-analog γ-hydroxybutyric acid in conjunction with a nonspecific neurological phenotype. We demonstrate for the first time using MDCK cells that NCS-382 is actively transported and capable of inhibiting GHB transport. Following these in vitro assays with in vivo studies in aldh5a1^-/- mice, we found the ratio of brain/liver GHB to be unaffected by chronic NCS-382 administration (300mg/kg; 7 consecutive days). Employing a variety of cellular parameters (reactive oxygen and superoxide species, ATP production and decay, mitochondrial and lysosomal number, cellular viability and necrosis), we demonstrate that up to 1mM NCS-382 shows minimal evidence of pharmacotoxicity. As well, studies at the molecular level indicate that the effects of NCS-382 at 0.5mM are minimally toxic as evaluated using gene expression assay. The cumulative data provides increasing confidence that NCS-382 could eventually be considered in the therapeutic armament for heritable SSADH deficiency.

γ-Hydroxybutyric Acid (GHB) Pharmacokinetics and Pharmacodynamics: Semi-Mechanistic and Physiologically Relevant PK/PD Model

An overdose of γ-hydroxybutyric acid (GHB), a drug of abuse, results in fatality caused by severe respiratory depression. In this study, a semi-mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model was developed to characterize monocarboxylate transporter 1 (MCT1)-mediated transport of GHB, as well as effects of GHB on respiration frequency, for IV doses of 200, 600, and 1500 mg/kg in rats. The proposed PK/PD model for GHB consists of nonlinear metabolism of GHB in the liver, MCT1-mediated renal reabsorption with physiologically relevant concurrent fluid reabsorption, MCT1-mediated uptake into the brain, and direct effects of binding of GHB to GABA_B receptors on the PD parameter, respiration frequency. Michaelis-Menten affinity constants for metabolism, renal reabsorption, and uptake into and efflux from the brain were fixed to the observed in vitro values. The IC ₅₀ value for the effect of GHB on respiration frequency was fixed to a reported value for binding of GHB to GABA_B receptors. All physiological parameters were fixed to the reported values for a 300-g rat. The model successfully captured the GHB PK/PD data and was further validated using the data for a 600-mg/kg dose of GHB after IV bolus administration. Unbound GHB brain ECF/blood partition coefficient (Kp _u,u ) values obtained from the model agreed well with values calculated using experimental ECF concentrations obtained with brain microdialysis, demonstrating the physiological relevance of this model. Sensitivity analysis indicated that the PK/PD model was stable. In conclusion, we developed a semi-mechanistic and physiologically relevant PK/PD model of GHB using in vitro drug-transporter kinetics and in vivo PK/PD data in rats.

Novel minimal physiologically-based model for the prediction of passive tubular reabsorption and renal excretion clearance

PURPOSE

Develop a minimal mechanistic model based on in vitro-in vivo extrapolation (IVIVE) principles to predict extent of passive tubular reabsorption. Assess the ability of the model developed to predict extent of passive tubular reabsorption (Freab) and renal excretion clearance (CLR) from in vitro permeability data and tubular physiological parameters.

METHODS

Model system parameters were informed by physiological data collated following extensive literature analysis. A database of clinical CLR was collated for 157 drugs. A subset of 45 drugs was selected for model validation; for those, Caco-2 permeability (Papp) data were measured under pH6.5-7.4 gradient conditions and used to predict Freab and subsequently CLR. An empirical calibration approach was proposed to account for the effect of inter-assay/laboratory variation in Papp on the IVIVE of Freab.

RESULTS

The 5-compartmental model accounted for regional differences in tubular surface area and flow rates and successfully predicted the extent of tubular reabsorption of 45 drugs for which filtration and reabsorption were contributing to renal excretion. Subsequently, predicted CLR was within 3-fold of the observed values for 87% of drugs in this dataset, with an overall gmfe of 1.96. Consideration of the empirical calibration method improved overall prediction of CLR (gmfe=1.73 for 34 drugs in the internal validation dataset), in particular for basic drugs and drugs with low extent of tubular reabsorption.

CONCLUSIONS

The novel 5-compartment model represents an important addition to the IVIVE toolbox for physiologically-based prediction of renal tubular reabsorption and CLR. Physiological basis of the model proposed allows its application in future mechanistic kidney models in preclinical species and human.

Key to Opening Kidney for In Vitro-In Vivo Extrapolation Entrance in Health and Disease: Part II: Mechanistic Models and In Vitro-In Vivo Extrapolation

It is envisaged that application of mechanistic models will improve prediction of changes in renal disposition due to drug-drug interactions, genetic polymorphism in enzymes and transporters and/or renal impairment. However, developing and validating mechanistic kidney models is challenging due to the number of processes that may occur (filtration, secretion, reabsorption and metabolism) in this complex organ. Prediction of human renal drug disposition from preclinical species may be hampered by species differences in the expression and activity of drug metabolising enzymes and transporters. A proposed solution is bottom-up prediction of pharmacokinetic parameters based on in vitro-in vivo extrapolation (IVIVE), mediated by recent advances in in vitro experimental techniques and application of relevant scaling factors. This review is a follow-up to the Part I of the report from the 2015 AAPS Annual Meeting and Exhibition (Orlando, FL; 25th-29th October 2015) which focuses on IVIVE and mechanistic prediction of renal drug disposition. It describes the various mechanistic kidney models that may be used to investigate renal drug disposition. Particular attention is given to efforts that have attempted to incorporate elements of IVIVE. In addition, the use of mechanistic models in prediction of renal drug-drug interactions and potential for application in determining suitable adjustment of dose in kidney disease are discussed. The need for suitable clinical pharmacokinetics data for the purposes of delineating mechanistic aspects of kidney models in various scenarios is highlighted.