Estimation of the minimum permeability coefficient in rats for perfusion-limited tissue distribution in whole-body physiologically-based pharmacokinetics

A Complete Extension of Classical Hepatic Clearance Models Using Fractional Distribution Parameter fd in Physiologically Based Pharmacokinetics

The classical organ clearance models have been proposed to relate the plasma clearance CLp to probable mechanism(s) of hepatic clearance. However, the classical models assume the intrinsic capability of drug elimination (CLu,int) that is physically segregated from the vascular blood but directly acts upon the unbound drug concentration in the blood (fubCavg), and do not handle the transit-time delay between the inlet/outlet concentrations in their closed-form clearance equations. Therefore, we propose unified model structures that can address the internal blood concentration patterns of clearance organs in a more mechanistic/physiological manner, based on the fractional distribution parameter fd operative in PBPK. The basic partial/ordinary differential equations for four classical models are revisited/modified to yield a more complete set of extended clearance models, i.e., the Rattle, Sieve, Tube, and Jar models, which are the counterparts of the dispersion, series-compartment, parallel-tube, and well-stirred models. We demonstrate the feasibility of applying the resulting extended models to isolated perfused rat liver data for 11 compounds and an example dataset for in vitro-in vivo extrapolation of the intrinsic to the systemic clearances. Based on their feasibilities to handle such real data, these models may serve as an improved basis for applying clearance models in the future.

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.

Physiologically Based Pharmacokinetic Modelling to Predict Pharmacokinetics of Enavogliflozin, a Sodium-Dependent Glucose Transporter 2 Inhibitor, in Humans

Enavogliflozin is a sodium-dependent glucose cotransporter 2 (SGLT2) inhibitor approved for clinical use in South Korea. As SGLT2 inhibitors are a treatment option for patients with diabetes, enavogliflozin is expected to be prescribed in various populations. Physiologically based pharmacokinetic (PBPK) modelling can rationally predict the concentration–time profiles under altered physiological conditions. In previous studies, one of the metabolites (M1) appeared to have a metabolic ratio between 0.20 and 0.25. In this study, PBPK models for enavogliflozin and M1 were developed using published clinical trial data. The PBPK model for enavogliflozin incorporated a non-linear urinary excretion in a mechanistically arranged kidney model and a non-linear formation of M1 in the liver. The PBPK model was evaluated, and the simulated pharmacokinetic characteristics were in a two-fold range from those of the observations. The pharmacokinetic parameters of enavogliflozin were predicted using the PBPK model under pathophysiological conditions. PBPK models for enavogliflozin and M1 were developed and validated, and they seemed useful for logical prediction.

Theoretical Examination Seeking Tangible Physical Meanings of Slopes and Intercepts of Plasma Concentration-Time Relationships in Minimal Physiologically Based Pharmacokinetic Models

In minimal physiologically based pharmacokinetic (mPBPK) models, physiological (e.g., cardiac output) and anatomical (e.g., blood/tissue volumes) variables are utilized in the domain of differential equations (DEs) for mechanistic understanding of the plasma concentration-time relationships [Formula: see text]. Although fundamental biopharmaceutical variables in terms of distribution (e.g., [Formula: see text] and [Formula: see text]) and elimination kinetics (e.g., [Formula: see text]) in mPBPK provide greater insights in comparison to classical compartment models, an absence of kinetic elucidation of slopes and intercepts in light of such DE model parameters hinders more intuitive appreciation of [Formula: see text]. Therefore, this study seeks the tangible physical meanings of slopes and intercepts of the plasma concentration-time relationships in one- and two-tissue mPBPK models (i.e., m2CM and m3CM), with respect to time parameters that are readily understandable in PK analyses, i.e., the mean residence ([Formula: see text]) and transit ([Formula: see text]) times. Utilizing the explicit equations (EEs) for the slopes, intercepts, and areas of each exponential phase in the m2CM and m3CM, we theoretically and numerically examined the limiting/boundary conditions of such kinetic properties, based on the ratio of the longest tissue [Formula: see text] to the [Formula: see text] in the body (i.e., [Formula: see text]) that is useful for dissecting complex PBPK systems. The kinetic contribution of the area of each exponential phase to the total drug exposure was assessed to identify the elimination phase between the terminal and non-terminal phases of the [Formula: see text] in the m2CM and m3CM. This assessment provides improved understanding of the complexities inherent in all PBPK profiles and models.

Determinants of Biological Half-Lives and Terminal Slopes in Physiologically Based Pharmacokinetic Systems: Assessment of Limiting Conditions

In pharmacokinetic (PK) analyses, the biological half-life T_1/2 is usually determined in the terminal phase after drug administration, which is readily calculated from the relationship T_1/2 = ln2/λ_z where λ_z is the terminal-phase slope obtainable from non-compartmental analysis (NCA). Since kinetic understanding of λ_z has been limited to the theory of a one-compartment model, this study seeks kinetic determinants of λ_z in more complex plasma concentration-time profiles. We utilized physiologically based pharmacokinetic (PBPK) systems that are consistent with the assumptions of NCA (e.g., linear PK and elimination occurring from plasma) to interrelate λ_z and disposition kinetic parameters of PBPK models. In a mammillary form of PBPK models, the two boundary conditions of λ_z are the inverses of the mean residence time in the body (1/MRT_B = CL/V_SS) and the mean transit time through the kinetically largest tissue (1/MTT_max = Q_Tf_dR_b/V_TK_p). Importantly, the limiting conditions of λ_z between 1/MRT_B and 1/MTT_max are dependent on a simple product MRT_Bλ_z (P_det) and a simple ratio MTT_max/MRT_B (K_det), leading to introduction of the unitless product-ratio plot for determination of the limiting condition of λ_z in linear PK. We found that the MRT_Bλ_z value of 0.5 serves as a practical threshold determining whether λ_z is more closely associated with 1/MRT_B or 1/MTT_max. The current theory was applied for assessment of the terminal slope λ_z for observed PK data of various compounds in man and rat.

Determination of the Number of Tissue Groups of Kinetically Distinct Transit Time in Whole-Body Physiologically Based Pharmacokinetic (PBPK) Models I: Theoretical Consideration of Bottom-Up Approach of Lumping Tissues in Whole-Body PBPK

Minimal physiologically based pharmacokinetic (mPBPK) models, consisting of system-specific (e.g., tissue volume and blood flow) and drug-related (e.g., tissue-to-plasma partition coefficient) parameters, are practically useful for pharmacokinetic analyses. However, biopharmaceutical principles were not clear on how peripheral tissues, adopted in whole-body physiologically based pharmacokinetic (WB-PBPK) models, could be kinetically consolidated into one or two tissue groups in the mPBPK models. In this theoretical examination, we studied the relationship between the progressive tissue lumping in the direction from the longest mean transit time (MTT_max) to the shorter one(s) and the slopes of the terminal (λ_ter)/distributional phases, assuming tissues with comparable MTTs are kinetically combined. The appropriateness of lumping was ascertained by evaluating the impact of difference in tissue MTTs during the lumping on the analytical solution of WB-PBPK models. We found that the ratio of MTT_max to the mean residence time in the body, viz., K_det, is related to the change in λ_ter by the progressive lumping and can serve as an index for the robustness of λ_ter. Calculations with two extreme cases revealed that, for caffeine at K_det < 0.03, the change in λ_ter was minimal even when all peripheral tissues were collectively lumped, whereas for artesunic acid at K_det > 50, the tissue of MTT_max could not be kinetically combined even with the tissue having the second-longest MTT without significantly affecting λ_ter. Therefore, we proposed K_det as an index for the robustness of λ_ter during tissue lumping and for the number of tissue groups with distinct transit times in WB-PBPK.

Determination of the Number of Tissue Groups of Kinetically Distinct Transit Time in Whole-Body Physiologically Based Pharmacokinetic (PBPK) Models II: Practical Application of Tissue Lumping Theories for Pharmacokinetics of Various Compounds

In our companion paper, we described the theoretical basis for tissue lumping in whole-body physiologically based pharmacokinetic (WB-PBPK) models and found that K_det, a coefficient for determining the number of tissue groups of distinct transit time in WB-PBPK models, was related to the fractional change in the terminal slope (FCT) when tissues were progressively lumped from the longest transit time to shorter ones. This study was conducted to identify the practical threshold of K_det by applying the lumping theory to plasma/blood concentration-time relationships of 113 model compounds collected from the literature. We found that drugs having K_det < 0.3 were associated with FCT < 0.1 even when all peripheral tissues were lumped, resulting in comparable plasma concentration-time profiles between one-tissue minimal PBPK (mPBPK) and WB-PBPK models. For drugs with K_det ≥ 1, WB-PBPK profiles appeared similar with two-tissue mPBPK models by applying the rule of FCT < 0.1 for lumping slowly equilibrating tissues. The two-tissue mPBPK model also appeared appropriate in terms of concentration-time profiles for drugs with 0.3 ≤ K_det < 1, although, some compounds (15.9% of the total cases), but not all, in this range showed a slight (maximum of 18.9% of the total AUC) deviation from WB-PBPK models, indicating that the two-tissue model, with caution, could still be used for those cases. Comparison of kinetic parameters between traditional (model-fitting) and current (theoretical calculation) mPBPK analyses revealed their significant correlations. Collectively, these observations suggest that the number of tissue groups could be determined based on the K_det/FCT criteria, and plasma concentration-time profiles from WB-PBPK could be calculated using equations significantly less complex.

A Permeability- and Perfusion-based PBPK model for Improved Prediction of Concentration-time Profiles

To improve predictions of concentration‐time (C‐t) profiles of drugs, a new physiologically based pharmacokinetic modeling framework (termed ‘PermQ’) has been developed. This model includes permeability into and out of capillaries, cell membranes, and intracellular lipids. New modeling components include (i) lumping of tissues into compartments based on both blood flow and capillary permeability, and (ii) parameterizing clearances in and out of membranes with apparent permeability and membrane partitioning values. Novel observations include the need for a shallow distribution compartment particularly for bases. C‐t profiles were modeled for 24 drugs (7 acidic, 5 neutral, and 12 basic) using the same experimental inputs for three different models: Rodgers and Rowland (RR), a perfusion‐limited membrane‐based model (K_p,mem), and PermQ. K_p,mem and PermQ can be directly compared since both models have identical tissue partition coefficient parameters. For the 24 molecules used for model development, errors in V_ss and t_1/2 were reduced by 37% and 43%, respectively, with the PermQ model. Errors in C‐t profiles were reduced (increased EOC) by 43%. The improvement was generally greater for bases than for acids and neutrals. Predictions were improved for all 3 models with the use of parameters optimized for the PermQ model. For five drugs in a test set, similar results were observed. These results suggest that prediction of C‐t profiles can be improved by including capillary and cellular permeability components for all tissues.

Consideration of Fractional Distribution Parameter fd in the Chen and Gross Method for Tissue-to-Plasma Partition Coefficients: Comparison of Several Methods

PURPOSE

The tissue-to-plasma partition coefficient (K_p) describes the extent of tissue distribution in physiologically-based pharmacokinetic (PBPK) models. Constant-rate infusion studies are common for experimental determination of the steady-state K_p,ss, while the tissue-plasma concentration ratio (C_T/C_p) in the terminal phase after intravenous doses is often utilized. The Chen and Gross (C&G) method converts a terminal slope C_T/C_p to K_p,ss based on assumptions of perfusion-limited distribution in tissue-plasma equilibration. However, considering blood flow (Q_T) and apparent tissue permeability (f_upPS_in) in the rate of tissue distribution, this report extends the C&G method by utilizing a fractional distribution parameter (f_d).

METHODS

Relevant PBPK equations for non-eliminating and eliminating organs along with lung and liver were derived for the conversion of C_T/C_p values to K_p,ss. The relationships were demonstrated in rats with measured C_T/C_p and K_p,ss values and the model-dependent f_d for 8 compounds with a range of permeability coefficients. Several methods of assessing K_p were compared.

RESULTS

Utilizing f_d in an extended C&G method, our estimations of K_p,ss from C_T/C_p were improved, particularly for lower permeability compounds. However, four in silico methods for estimating K_p performed poorly across tissues in comparison with measured K_p values. Mathematical relationships between K_p and K_p,ss that are generally applicable for eliminating organs with tissue permeability limitations necessitates inclusion of an extraction ratio (ER) and f_d.

CONCLUSION

Since many different types/sources of K_p are present in the literature and used in PBPK models, these perspectives and equations should provide better insights in measuring and interpreting K_p values in PBPK.

Assessment of the Kochak-Benet Equation for Hepatic Clearance for the Parallel-Tube Model: Relevance of Classic Clearance Concepts in PK and PBPK

This report reviews concepts related to operation of the classic parallel-tube model (PTM) for hepatic disposition and examines two recent proposals of a newly derived equation to describe hepatic clearance (CL_H). It is demonstrated that the proposed equation is identical to a re-arrangement of an earlier relationship from Pang and Rowland and provides a means of calculation of intrinsic clearance (CL_int,PTM) rather than CL_H as posed. We further demonstrate how classic hepatic clearance models with an assumed CL_int, while subject to numerous limitations, remain highly useful and necessary in both traditional pharmacokinetics (PK) and physiologically based pharmacokinetic (PBPK) modeling.

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.

Meta-Assessment of Metformin Absorption and Disposition Pharmacokinetics in Nine Species

The objective of this study was to systematically assess literature datasets and quantitatively analyze metformin PK in plasma and some tissues of nine species. The pharmacokinetic (PK) parameters and profiles of metformin in nine species were collected from the literature. Based on a simple allometric scaling, the systemic clearances (CL) of metformin in these species highly correlate with body weight (BW) (R2 = 0.85) and are comparable to renal plasma flow in most species except for rabbit and cat. Reported volumes of distribution (VSS) varied appreciably (0.32 to 10.1 L/kg) among species. Using the physiological and anatomical variables for each species, a minimal physiologically based pharmacokinetic (mPBPK) model consisting of blood and two tissue compartments (Tissues 1 and 2) was used for modeling metformin PK in the nine species. Permeability-limited distribution (low fd1 and fd2) and a single tissue-to-plasma partition coefficient (Kp) value for Tissues 1 and 2 were applied in the joint mPBPK fitting. Nonlinear regression analysis for common tissue distribution parameters along with species-specific CL values reasonably captured the plasma PK profiles of metformin across most species, except for rat and horse with later time deviations. In separate fittings of the mPBPK model to each species, Tissue 2 was considered as slowly-equilibrating compartment consisting of muscle and skin based on in silico calculations of the mean transit times through tissues. The well-fitted mPBPK model parameters for absorption and disposition PK of metformin for each species were compared with in vitro/in vivo results found in the literature with regard to the physiological details and physicochemical properties of metformin. Bioavailability and absorption rates decreased with the increased BW among the species. Tissues such as muscle dominate metformin distribution with low permeability and partitioning while actual tissue concentrations found in rats and mice show likely transporter-mediated uptake in liver, kidney, and gastrointestinal tissues. Metformin has diverse pharmacologic actions, and this assessment revealed allometric relationships in its absorption and renal clearance but considerable variability in actual and modeled tissue distribution probably caused by transporter differences.

Application of physiologically based pharmacokinetic models to promote the development of veterinary drugs with high efficacy and safety

Physiologically based pharmacokinetic (PBPK) models have become important tools for the development of novel human drugs. Food-producing animals and pets comprise an important part of human life, and the development of veterinary drugs (VDs) has greatly impacted human health. Owing to increased affordability of and demand for drug development, VD manufacturing companies should have more PBPK models required to reduce drug production costs. So far, little attention has been paid on applying PBPK models for the development of VDs. This review begins with the development processes of VDs; then summarizes case studies of PBPK models in human or VD development; and analyzes the application, potential, and advantages of PBPK in VD development, including candidate screening, formulation optimization, food effects, target-species safety, and dosing optimization. Then, the challenges of applying the PBPK model to VD development are discussed. Finally, future opportunities of PBPK models in designing dosing regimens for intracellular pathogenic infections and for efficient oral absorption of VDs are further forecasted. This review will be relevant to readers who are interested in using a PBPK model to develop new VDs.

Physiologically Based Pharmacokinetics of Lysosomotropic Chloroquine in Rat and Human

A semimechanistic physiologically based pharmacokinetic (PBPK) model for chloroquine (CQ), a highly lysosomotropic weak base, was applied to digitized rat and human concentration versus time data. The PBPK model in rat featured plasma and red blood cell (RBC) concentrations, extensive and apparent nonlinear tissue distribution, fitted hepatic and renal intrinsic clearances, and a plasma half-life of about 1 day. Tissue-to-plasma CQ ratios at 50 hours after dosing were highest in lung, kidney, liver, and spleen (182-318) and lower in heart, muscle, brain, eye, and skin (11-66). The RBC-to-plasma ratio of 11.6 was assumed to reflect cell lipid partitioning. A lysosome-based extended model was used to calculate subcellular CQ concentrations based on tissue mass balances, fitted plasma, interstitial and free cytosol concentrations, and literature-based pH and pKa values. The CQ tissue component concentrations ranked as follows: lysosome > > acidic phospholipid > plasma = interstitial = cytosol ≥ neutral lipids. The extensive lysosome sequestration appeared to change over time and was attributed to lowering pH values caused by proton pump influx of hydrogen ions. The human-to-rat volume of distribution (Vss) ratio of 7 used to scale rat tissue partitioning to human along with estimation of hepatic clearance allowed excellent fitting of oral-dose PK (150-600 mg) of CQ with a 50-day half-life in humans. The prolonged PK of chloroquine was well characterized for rat and human with this PBPK model. The calculated intratissue concentrations and lysosomal effects have therapeutic relevance for CQ and other cationic drugs. SIGNIFICANCE STATEMENT: Sequestration in lysosomes is a major factor controlling the pharmacokinetics and pharmacology of chloroquine and other cationic drugs. This report provides comprehensive physiologic modeling of chloroquine distribution in tissues and overall disposition in rat and human that reveals expected complexities and inferences related to its subcellular association with various tissue components.

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.

The conditional stimulation of rat organic cation transporter 2, but not its human ortholog, by mesoridazine: the possibility of the involvement of the high-affinity binding site of the transporter in the stimulation

Objectives

To study the functional consequences of the human and rat forms of OCT2 in the presence of phenothiazines.

Methods

MDCK cells expressing human or rat OCT2 were established, and MPP+ transport was determined by uptake assays. Concentration dependency was studied for the stimulatory/inhibitory effects of phenothiazines on MPP+ transport.

Key findings

Among the 11 phenothiazines examined, the majority were found to have comparable effects on transporter function between the orthologous forms, while three phenothiazines, particularly mesoridazine, had complex impacts on transporter function. For rOCT2, mesoridazine stimulated transport at 0.1 and 1 μmMPP+ with the mesoridazine concentration–uptake curve becoming bell-shaped. This conditional effect became less pronounced at 30 μmMPP+, resulting in an inhibition curve with a typical profile. For hOCT2, mesoridazine behaved as a typical inhibitor of transporter function at all MPP+ concentrations, although the kinetics of inhibition were still affected by the substrate concentration.

Conclusions

The conditional stimulation by mesoridazine in rOCT2, and the lack thereof in hOCT2, may be a manifestation of the interaction of phenothiazine with substrate binding at the high-affinity site of the OCT2. As OCT2 was previously indicated in some drug–drug interactions, the conditional stimulation of OCT2 and its potential species-differences may be of practical relevance.