Impact of Direction of Transport on the Evaluation of Inhibition Potencies of Multidrug and Toxin Extrusion Protein 1 Inhibitors

Prediction of Inhibitory Activity against the MATE1 Transporter via Combined Fingerprint- and Physics-Based Machine Learning Models

Renal secretion plays an important role in excretion of drug from the kidney. Two major transporters known to be highly involved in renal secretion are MATE1/2 K and OCT2, the former of which is highly related to drug-drug interactions. Among published in silico models for MATE inhibitors, a previous model obtained a ROC-AUC value of 0.78 using high throughput percentage inhibition data [J. Med. Chem. 2013, 56(3), 781-795] which we aimed to improve upon here using a combined fingerprint and physics-based approach. To this end, we collected 225 publicly available compounds with pIC50 values against MATE1. Subsequently, on the one hand, we performed a physics-based approach using an Alpha-Fold protein structure, from which we obtained MM-GB/SA scores for those compounds. On the other hand, we built Random Forest (RF) and message passing neural network models using extended-connectivity fingerprints with radius 4 (ECFP4) and chemical structures as graphs, respectively, which also included MM-GB/SA scores as input variables. In a five-fold cross-validation with a separate test set, we found that the best predictivity for the hold-out test was observed in the RF model (including ECFP4 and MM-GB/SA data) with an ROC-AUC of 0.833 ± 0.036; while that of the MM-GB/SA regression model was 0.742. However, the MM-GB/SA model did not show a dependency of the performance on the particular chemical space being predicted. Additionally, via structural interaction fingerprint analysis, we identified interacting residues with inhibitor as identical for those with noninhibitors, including substrates, such as Gln49, Trp274, Tyr277, Tyr299, Ile303, and Tyr306. The similar binding modes are consistent with the observed similar IC50 value inhibitor when using different substrates experimentally, and practically, this can release the experimental scientists from bothering of selecting substrates for MATE1. Hence, we were able to build highly predictive classification models for MATE1 inhibitory activity with both ECFP4 and MM-GB/SA score as input features, which is fit-for-purpose for use in the drug discovery process.

Extracellular and intracellular productions of lysophosphatidic acids and cyclic phosphatidic acids by lysophospholipase D from exogenously added lysophosphatidylcholines to cultured NRK52E cells

Lysophosphatidic acid (LPA) is a bioactive lysophospholipid that is a notable biomarker of kidney injury. However, it is not clear how LPA is produced in renal cells. In this study, we explored LPA generation and its enzymatic pathway in a rat kidney-derived cell, NRK52E cells. Culturing of NRK52E cells with acyl lysophosphatidylcholine (acyl LPC), or lyso-platelet activating factor (lysoPAF, alkyl LPC) was resulted in increased extracellular level of choline, co-product with LPA by lysophospholipase D (lysoPLD). Their activities were enhanced by addition of calcium ions to the cell culture medium, but failed to be inhibited by S32826, an autotaxin (ATX)-specific inhibitor. Liquid chromatography-tandem mass spectrometric analysis revealed the small, but significant extracellular production of acyl LPA/cyclic phosphatidic acid (cPA) and alkyl LPA/cPA. The mRNA expression of glycerophosphodiesterase (GDE) 7 with lysoPLD activity was elevated in confluent NRK52E cells cultured over 3 days. GDE7 plasmid-transfection of NRK52E cells augmented both extracellular and intracellular productions of LPAs (acyl and alkyl) as well as extracellular productions of cPAs (acyl and alkyl) from exogenous LPCs (acyl and alkyl). These results suggest that intact NRK52E cells are able to produce choline and LPA/cPA from exogenous LPCs through the enzymatic action of GDE7 that is located on the plasma membranes and intracellular membranes.

Reimagining the Framework Supporting the Static Analysis of Transporter Drug Interaction Risk; Integrated Use of Biomarkers to Generate Pan‐Transporter Inhibition Signatures

Solute carrier (SLC) transporters present as the loci of important drug-drug interactions (DDIs). Therefore, sponsors generate in vitro half-maximal inhibitory concentration (IC₅₀ ) data and apply regulatory agency-guided "static" methods to assess DDI risk and the need for a formal clinical DDI study. Because such methods are conservative and high false-positive rates are likely (e.g., DDI study triggered when liver SLC R value ≥ 1.04 and renal SLC maximal unbound plasma (C_max,u )/IC₅₀ ratio ≥ 0.02), investigators have attempted to deploy plasma- and urine-based SLC biomarkers in phase I studies to de-risk DDI and obviate the need for drug probe-based studies. In this regard, it was possible to generate in-house in vitro SLC IC₅₀ data for various clinically (biomarker)-qualified perpetrator drugs, under standard assay conditions, and then estimate "% inhibition" for each SLC and relate it empirically to published clinical biomarker data (area under the plasma concentration vs. time curve (AUC) ratio (AUCR, AUC_inhibitor /AUC_reference ) and % decrease in renal clearance (ΔCL_renal )). After such a "calibration" exercise, it was determined that only compounds with high R values (> 1.5) and C_max,u /IC₅₀ ratios (> 0.5) are likely to significantly modulate liver (AUCR > 1.25) and renal (ΔCL_renal  > 25%) biomarkers and evoke DDI risk. The % inhibition approach supports integration of liver and renal SLC data and allows one to generate pan-SLC inhibition signatures for different test perpetrators (e.g., SLC % inhibition ranking). In turn, such signatures can guide the selection of the most appropriate individual (or combinations of) biomarkers for testing in phase I studies.

Challenges and Opportunities for Improved Drug-Drug Interaction Predictions for Renal OCT2 and MATE1/2-K Transporters

Transporters contribute to renal elimination of drugs; therefore drug disposition can be impacted if transporters are inhibited by comedicant drugs. Regulatory agencies have provided guidelines to assess potential drug-drug interaction (DDI) risk for renal organic cation transporter 2 (OCT2) and multidrug and toxin extrusion 1 and 2-K (MATE1/2-K) transporters. Despite this, there are challenges with translating in vitro data using currently available tools to obtain a quantitative assessment of DDI risk in the clinic. Given the high number of drugs and new molecular entities showing in vitro inhibition toward OCT2 and/or MATE1/2-K and the lack of translation to clinically significant effects, it is reasonable to question whether the current in vitro assay design and modeling practice has led to unnecessary clinical evaluation. The aim of this review is to assess and discuss available in vitro and clinical data along with prediction models intended to provide clinical context of risk, including static models proposed by regulatory agencies and physiologically-based pharmacokinetic models, in order to identify best practices and areas of future opportunity. This analysis highlights that different in vitro assay designs, including substrate and cell systems used, strongly influence the derived concentration of drug producing 50% inhibition values and contribute to high variability observed across laboratories. Furthermore, the lack of sensitive index substrates coupled with specific inhibitors for individual transporters necessitates the use of complex models to evaluate clinical DDI risk.