Advertisement
Research ArticleArticle

Discovery and Validation of Pyridoxic Acid and Homovanillic Acid as Novel Endogenous Plasma Biomarkers of Organic Anion Transporter (OAT) 1 and OAT3 in Cynomolgus Monkeys

Hong Shen, David M. Nelson, Regina V. Oliveira, Yueping Zhang, Colleen A. Mcnaney, Xiaomei Gu, Weiqi Chen, Ching Su, Michael D. Reily, Petia A. Shipkova, Jinping Gan, Yurong Lai, Punit Marathe and W. Griffith Humphreys
Drug Metabolism and Disposition February 2018, 46 (2) 178-188; DOI: https://doi.org/10.1124/dmd.117.077586

Abstract

Perturbation of organic anion transporter (OAT) 1- and OAT3-mediated transport can alter the exposure, efficacy, and safety of drugs. Although there have been reports of the endogenous biomarkers for OAT1/3, none of these have all of the characteristics required for a clinical useful biomarker. Cynomolgus monkeys were treated with intravenous probenecid (PROB) at a dose of 40 mg/kg in this study. As expected, PROB increased the area under the plasma concentration-time curve (AUC) of coadministered furosemide, a known substrate of OAT1 and OAT3, by 4.1-fold, consistent with the values reported in humans (3.1- to 3.7-fold). Of the 233 plasma metabolites analyzed using a liquid chromatography–tandem mass spectrometry (LC-MS/MS)–based metabolomics method, 29 metabolites, including pyridoxic acid (PDA) and homovanillic acid (HVA), were significantly increased after either 1 or 3 hours in plasma from the monkeys pretreated with PROB compared with the treated animals. The plasma of animals was then subjected to targeted LC-MS/MS analysis, which confirmed that the PDA and HVA AUCs increased by approximately 2- to 3-fold by PROB pretreatments. PROB also increased the plasma concentrations of hexadecanedioic acid (HDA) and tetradecanedioic acid (TDA), although the increases were not statistically significant. Moreover, transporter profiling assessed using stable cell lines constitutively expressing transporters demonstrated that PDA and HVA are substrates for human OAT1, OAT3, OAT2 (HVA), and OAT4 (PDA), but not OCT2, MATE1, MATE2K, OATP1B1, OATP1B3, and sodium taurocholate cotransporting polypeptide. Collectively, these findings suggest that PDA and HVA might serve as blood-based endogenous probes of cynomolgus monkey OAT1 and OAT3, and investigation of PDA and HVA as circulating endogenous biomarkers of human OAT1 and OAT3 function is warranted.

Introduction

The organic anion transporters (OATs) OAT1 (SLC22A6) and OAT3 (SLC22A8) function as influx transporters that are mainly expressed on the basolateral membrane of renal proximal tubular cells (RPTCs) and mediate cellular uptake of substrates from blood into the cells (Motohashi et al., 2002). OAT1 and OAT3 can serve as the loci of drug-drug interactions (DDIs), and such DDIs can lead to undesired changes in systemic and local exposures of victim drugs and toxins (Morrissey et al., 2013; Nigam et al., 2015). For example, furosemide (FSM) is taken up from the blood in RPTCs via OAT1 and OAT3 (Hasannejad et al., 2004), and probenecid (PROB), an inhibitor of OAT1 and OAT3, causes significant alterations in the pharmacokinetic parameters of FSM in humans through the increased area under the plasma drug concentration-time curve (AUC), decreased total plasma clearance (CLTOT) and renal clearance (CLR), and decreased fraction excreted in the urine (fR) (Chennavasin et al., 1979; Smith et al., 1980; Vree et al., 1995). Moreover, a few uremic toxins recently have been demonstrated to be inhibitors of OAT1 and OAT3, and have the potential to inhibit renal secretion clearance of drug substrates of OAT1 and OAT3 in patients with chronic kidney disease (Hsueh et al., 2016). On the other hand, concomitant use of PROB is recommended by the U.S. Food and Drug Administration for reducing the pronounced nephrotoxicity of cidofovir by inhibiting OAT1 and OAT3 to reduce the exposure of cidofovir in RPTCs to an extent that results in an acceptable benefit-risk balance (http://www.accessdata.fda.gov/drugsatfda_docs/label/1999/020638s003lbl.pdf).

Animal models have often been used to assess preclinical drug-OAT interactions. For example, the clinical famotidine-PROB interaction was reproduced in cynomolgus monkeys recently since PROB caused an approximate 90% reduction in the tubular secretion clearance of famotidine in monkeys (4.58 ± 1.25 vs. 0.38 ± 0.36 ml/min per kilogram) (Tahara et al., 2006). Similarly, PROB significantly decreased the renal tubular secretion clearance of famotidine from 2.80 ± 0.36 to 0.37 ± 0.07 ml/min per kilogram in humans. Furthermore, the protective effect of PROB treatment on the nephrotoxicity of cidofovir in humans was able to be recapitulated in monkeys (Lacy et al., 1998). These results suggest that monkey is the more suitable animal model to predict the clinical DDIs involving OAT1 and OAT3.

Recently, endogenous biomarkers have been envisioned as simple, fast, and cost-effective tools to monitor transporter activity in preclinical and clinical settings to facilitate development of a drug candidate (Ito et al., 2012; Müller et al., 2015; Lai et al., 2016; Yee et al., 2016; Bergagnini-Kolev et al., 2017). Therefore, the identification of a sensitive endogenous biomarker for OAT1 and OAT3 would be of great value. In this regard, Imamura et al. (2014) reported that 6β-hydroxycortisol could be an endogenous probe for OAT3 inhibition evaluation since PROB significantly changed the AUC and CLR of 6β-hydroxycortisol in healthy subjects. However, 6β-hydroxycortisol is formed from cortisol by hepatic CYP3A4, and many drugs known to affect this drug-metabolizing enzyme can also change exposure of 6β-hydroxycortisol (Peng et al., 2011). Very recently, Tsuruya et al. (2016) reported that taurine and glycochenodeoxycholate sulfate were endogenous biomarkers of OAT1 and OAT3, respectively. However, the plasma levels of taurine and glycochenodeoxycholate sulfate were not changed significantly by PROB treatment compared with control even though the renal secretion and CLR were significantly decreased (Tsuruya et al., 2016). Reduced activities of blood-facing OAT1 and OAT3 by PROB are supposed to increase plasma concentration of a sensitive and selective endogenous probe, mimicking systemic alterations of probe drugs such as FSM. To our knowledge, all published endogenous biomarkers of OAT1 and OAT3 could not recapitulate plasma drug concentration-time profiles in the presence of OAT inhibitors. Several features including specificity, sensitivity, predictability, reproducibility, acute response, and accessibility have been considered in the identification and validation of ideal endogenous biomarkers for drug transporters (Chu et al., 2017; Mariappan et al., 2017; Rodrigues et al., 2017). Given the aforementioned reasons, there is a need for novel plasma biomarkers of renal OATs.

In the present studies we provide direct experimental evidence that several organic anionic compounds, including pyridoxic acid (PDA), homovanilic acid (HVA), hexadecanedioic acid (HDA), and tetradecanedioic acid (TDA) are potential endogenous biomarkers of OAT1 and OAT3 in cynomolgus monkeys. Specifically, an untargeted metabolomics analysis was applied to plasma samples to screen endogenous compounds that were associated with OAT1 and OAT3 inhibition in monkeys. Follow-up quantitative liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis further characterized the time-plasma concentration profiles of selected endogenous compounds (i.e., PDA, HVA, had, and TDA) after PROB administration in cynomolgus monkeys. Moreover, transporter profiling assessed using human embryonic kidney (HEK) 293 cells stably transfected with major human renal and hepatic drug transporters demonstrated that PDA and HVA are OAT1 and OAT3 substrates and thus potential novel plasma endogenous biomarkers of OAT1 and OAT3 inhibition.

Materials and Methods

Materials.

[3H]Penciclovir (1.3 Ci/mmol) and [14C]metformin (98 mCi/mmol) were purchased from Moravek Biochemicals, Inc. (Brea, CA), and [3H]para-aminohippuric acid (4.5 Ci/mmol), [3H]estrone-3-sulfate (44.0 Ci/mmol), [3H]estradiol-17β-D-glucuronide (34.3 Ci/mmol), and [3H]cholecystokinin octapeptide (97.5 mCi/mmol) were purchased from PerkinElmer Life and Analytical Sciences (Waltham, MA). Nonradiolabeled FSM, PROB, and taurocholic acid, and the corresponding stable isotope-labeled internal standards were purchased from Toronto Research Chemicals Inc. (North York, Ontario, Canada). 4-PDA (≥98%) was purchased from Sigma-Aldrich (St. Louis, MO), and HVA (≥98%) was purchased from Acros Organics (Geel, Belgium). Hexadecanedioic acid-d28, tetradecanedioic acid-d24, and enalapril maleate-d5 (Enal-d5) were obtained from CDN Isotopes (Pointe-Claire, Quebec, Canada). Other nonradiolabeled compounds were purchased from either Sigma-Aldrich or Toronto Research Chemicals Inc. and were of analytical grade. Cell culture media and reagents were purchased from Invitrogen (Carlsbad, CA) or Mediatech, Inc., A Corning Subsidiary (Manassas, VA).

Cynomolgus Monkeys Pharmacokinetic FSM-PROB Interaction Study Protocol.

To identify and verify the endogenous biomarkers of OAT1 and OAT3, we conducted a series of experiments: cynomolgus monkey FSM-PROB interaction, metabolomics, targeted LC tandem MS (LC-MS/MS), and transporter profiling experiments. The experimental workflow is shown in Fig. 1.

Fig. 1.
Fig. 1.

Schematic overview of experimental workflow for identification and validation of endogenous biomarkers of OAT1 and OAT3. Plasma samples were collected from monkey FSM-PROB interaction study, and metabolites were quantified by LC-MS/MS. Raw data were extracted, and analyzed by various tools to identify associations between drug treatments and concentrations, determine significant correlations, and integrate results with transporter knowledge. Three selected potential endogenous biomarkers were further validated by LC-MS/MS with deuterated internal standards to confirm metabolomics observations. Transporter uptake studies were performed to determine whether the selected probes were selective substrates for human OAT1 and OAT3.

A single-dose, three-period crossover intravenous pharmacokinetic DDI study was carried out at the Bristol-Myers Squibb Company, and three male cynomolgus monkeys with body weights ranging from 5.3 to 6.0 kg during the study periods were included in this study. The experiment was performed in accordance with the National Institutes of Health guidelines (https://grants.nih.gov/grants/olaw/guide-for-the-care-and-use-of-laboratory-animals.pdf) and approved by the Bristol-Myers Squibb Animal Care and Use Committee. The animals were housed in a temperature- and humidity-controlled room with a 12-hour light/dark cycle.

In the first period, 40 mg/kg PROB dissolved in a sodium hydroxide solution (0.1 N) and titrated to neutral pH with hydrochloric acid was intravenously infused via a femoral vein catheter over 5 minutes (5 ml/kg). Venous blood samples (2 ml) were collected before and 0.08, 0.17, 0.25, 0.5, 0.75, 1, 2, 3, 5, 7, and 24 hours after administration in K2-EDTA-containing tubes. In the second period, after a washout period of 7 days, 2 mg/kg FSM dissolved in saline was intravenously infused via a femoral vein over 5 minutes (5 ml/kg), and blood samples were collected before and 0.08, 0.17, 0.25, 0.5, 0.75, 1, 2, 3, 5, 7, and 24 hours after administration. In the third period, 8 days after administrations of the previous dose of FSM, a FSM saline solution (2 mg/kg) was given by femoral vein infusion over 5 minutes to each monkey 30 minutes after PROB administration (40 mg/kg, i.v. infusion for 5 minutes). Blood was sampled before and 0.08, 0.17, 0.25, 0.5, 0.75, 1, 2, 3, 5, 7, and 24 hours after start of the FSM infusion. Blood samples were spun for 5 minutes at 13,000 rpm within 1 hour to obtain plasma. Urine samples were collected using metabolic cages for the following intervals in all three periods: 0–3, 3–7, and 7–24 hours after administration, and the volume of urine was recorded. The plasma and urine samples were stored at −80°C until LC-MS/MS and metabolomics analyses were conducted.

LC-MS/MS Analysis of Furosemide and Probenecid.

The LC-MS/MS analysis was performed on a SCIEX Triple Quad API-4000 system (AB SCIEX, Framingham, MA) coupled with a Shimadzu Nexera LD-30AD ultra-performance LC system (Shimadzu, Columbia, MD). The chromatographic separation was performed on an Agilent Zorbax RRHD SB-C8 column (2.0 × 100 mm, 1.7 µm) from Agilent Technologies (Santa Clara, CA) using mobile phases of 0.1% formic acid in water and 0.1% formic acid in acetonitrile (ACN). The flow rate was 0.7 ml/min and total run time was 3.7 minutes. The LC column was maintained at 60°C. The analytes were monitored using selected reaction monitoring in negative ionization mode with the optimized nebulizing and desolvation gases. The source temperature was set at 400°C, and the declustering potential and collision energy were optimized. FSM, furosemide-d5, and PROB were detected at the selected reaction monitoring transitions of m/z 329.1 → 285.0, 334.1 → 291.0, and 283.9 → 239.9, respectively.

Before the analysis, urine samples were diluted 100-fold into blank plasma and treated as previously described for plasma. The plasma and diluted urine samples were then extracted using protein precipitation. Specifically, the 100 µl plasma and diluted urine samples were mixed with 100 µl of ACN containing 400 nM of furosemide-d5 (internal standard), followed by vortex mixing on a mixer for 5 minutes at room temperature. The mixed solutions were then filtered using a MultiScreen hydrophilic flitration plate (Billerica, Millipore, MA) by centrifugation at 4000g for 10 minutes, and the filtered solution was then injected (5 µl) onto LC-MS/MS. The assay was qualified over the analytical range of 1–5000 nM using a linear 1/x2 weighed regression.

Plasma Metabolomic Profiling.

Frozen plasma samples were thawed and 50 µl aliquots were subjected to protein precipitation by the addition of 600 µl of methanol containing 0.1% formic acid and stable-labeled internal standards. The samples were mixed with a vortex mixer and subjected to centrifugation for 10 minutes at 2700g (4000 rpm). An aliquot (50 µl) of the resulting supernatants was transferred to a 96-well plate for further hydrophilic interaction LC (HILIC) LC-MS analysis and an additional 100 μl was transferred to a separate 96-well plate for reverse phase (RP) LC-MS analysis. Both plates were dried to completeness under a nitrogen stream at room temperature. For RP LC-MS analysis, the dried samples were reconstituted by the addition of 20 μl of methanol, followed by vigorous shaking and the addition of 180 μl of water. The samples were mixed with a vortex mixer and subjected to centrifugation for 10 minutes at 4000 rpm. The resulting supernatants were transferred to a 96-well plate, from which 10 μl was directly injected for analysis. For HILIC LC-MS analysis, the dried samples were reconstituted by the addition of 20 μl of water, followed by vigorous shaking and the addition of 180 μl of 50:50 MeOH:ACN. The samples were mixed with a vortex mixer and subjected to centrifugation for 10 minutes at 4000 rpm. The resulting supernatants were transferred to a 96-well plate, from which 10 μl was directly injected for analysis.

HILIC and RP LC-MS analyses were performed on a Nexera ×2 LC-30AD (Shimadzu, Somerset, NJ) ultra-high-performance LC system connected to an Exactive Plus (Thermo Fisher Scientific, Waltham, MA) mass spectrometer. The ultra-high-performance LC column for HILIC analyses was an Acquity BEH-NH2, 2.1 × 150 mm, 1.7 μ (Waters Corporation, Milford, MA) with mobile phases A (95:5 water:ACN, 10 mM NH4OAc, 0.05% NH4OH) and B (ACN, 0.05% NH4OH) at a flow rate of 300 µl/min with starting conditions of 95% B–37% B at 3.5 minutes, hold for 4 minutes, and down to starting conditions at 7 minutes, for a total run time of 11 minutes. The ultra-high-performance LC column for RP analyses was an Acquity BEH C18, 2.1 × 150 mm, 1.7 μ (Waters Corporation) with mobile phases A (water, 0.1% formic acid) and B (98:2ACN: water, 0.1% formic acid) at a flow rate of 600 µl/min with starting conditions of 100% A, to 80% A at 3 minutes, 40% A at 4 minutes, and 100% B by 7 minutes, and after a 2-minute hold, down to starting conditions at 9 minutes, for a total run time of 11 minutes. Both HILIC and RP LC-MS data were collected in positive and negative polarities (separate injections) at 35,000 resolution and expected mass accuracy of 5 ppm.

LC-MS data analysis was performed using software developed in-house, Expedient Data Mining, as described previously (Hnatyshyn et al., 2013). Direct data input from raw files was performed using the Thermo Fisher Scientific MSFileReader software. The resulting list of components was matched with accurate mass and retention time values using an in-house database of endogenous metabolites stored within Expedient Data Mining. The annotated table of component integrals was exported to Microsoft Excel for further statistical analysis using in-house Visual Basic scripts for Microsoft Excel. For these analyses, mean intensities for each treatment group were compared with the relevant concurrent control. For each component, fold change was calculated by dividing the treatment group value by the control group value and the P values were calculated as a pairwise comparison for two-tailed distribution, using the Student’s t test (Excel statistics package; Microsoft; Redmond, Washington).

Characterization of PDA and HVA Uptake in Stable Cell Lines Constitutively Expressing Major Renal Drug Transporters.

Uptake studies were performed as described previously (Shen et al., 2016b). Uptake of PDA and HVA was first measured at a single concentration for 5-minute incubation with OAT1-, OAT2-, OAT3-, OAT4-, OCT2-, MATE1-, MATE2K-, OATP1B1-, OATP1B3-, and sodium taurocholate cotransporting polypeptide (NTCP)-HEK cells, and then time- and concentration-dependent uptake was measured. The initial concentration used for PDA was 1 μM because the physiologic baseline level ranged from 0.67 to 2.5 µM in cynomolgus monkeys (Fig. 3A). The concentration studied for HVA was 5 μM because of limited bioanalytical sensitivity of this compound in Mock-HEK cells. In addition, the transport Michaelis-Menten constant value (Km) of 274 ± 100 μM (that corresponds to the substrate concentration at which the uptake rate is half that of the Vmax value) has been reported in rat Oat2-mediated transport for HVA (Mori et al., 2003), which is significantly greater than the selected testing concentrations (274 ± 100 vs. 5 μM). No corresponding data were available for PDA. Kinetic transport experiments were conducted under linear-uptake conditions or for shortest incubation duration with acceptable analytical sensitivity (Supplemental Fig. 1).

Cells were grown to confluence in 24-well poly-D-lysine-coated plates 2–3 days after seeding at a cell density of 500,000 cells per well (BD Biosciences, San Jose, CA). All experiments were conducted at 37°C using a working solution containing Hanks’ balanced salt solution (HBSS) supplemented with 10 mM HEPES (pH 7.4 for OAT1-, OAT2-, OAT3-, OAT4-, OCT2-, OATP1B1, OATP1B3-, and NTCP-HEK cells, and pH 8.4 for MATE1- and MATE2K-HEK cells, respectively), with a probe substrate, PDA, or HVA. The probe substrates used were [3H]para-aminohippuric acid (OAT1), [3H]penciclovir (OAT2), [3H]estrone-3-sulfate (OAT3 and OAT4), [14C]metformin (OCT2, MATE1, and MATE2K), [3H]estradiol-17β-D-glucuronide (OATP1B1), [3H] cholecystokinin octapeptide [OATP1B3, and taurocholic acid (NTCP)]. Compounds were dissolved in dimethyl sulfoxide and diluted in HBSS (maximum 0.2% dimethyl sulfoxide). Plating medium was removed, and cell monolayers were rinsed twice with prewarmed HBSS. Incubations were started by the addition of 200 μl of substrate prewarmed at 37°C. After incubation for desired time at 37°C, the cell monolayers were rinsed three times with 500 μl of ice-cold HBSS. The cells were then lysed with 300 μl buffer (0.1% Triton X-100 or methanol), and compound concentrations in the cell lysates were measured by either liquid scintillation counting (Tri-Carb 2910 TR Liquid Scintillation Analyzer; PerkinElmer Life and Analytical Sciences) or LC-MS/MS as described subsequently.

LC-MS/MS Analysis of PDA and HVA.

Stock PDA solution (390 µg/ml, 2.13 mM) was prepared by dissolving PDA in 0.5% NH4OH solution. Stock HVA (1.8 mg/ml, 9.5 mM) solution was prepared in water. Stock solutions were stored in darkness at −30°C and brought to room temperature before use. The highest calibration standard (5000 ng/ml) was prepared by diluting appropriate aliquots of PDA and HVA stock solutions to a final volume of 2 ml with 1% bovine serum albumin (BSA) in phosphate-buffered saline (PBS) (pH 7.4). Additional standard solutions were obtained by serial dilution from the 5000 ng/ml standard with 1% BSA in PBS (pH 7.4) to final concentrations of 2500, 1000, 500, 250, 100, 50, 25, 10, and 5 ng/ml. Quality control (QC) samples were also prepared with 1% BSA in PBS (pH 7.4) at three concentration levels: 7.5, 75, and 750 ng/ml. The internal standard stock solution of Enal-d5 was prepared at 1 mg/ml in methanol. The internal standard working solution containing 1000 ng/ml of Enal-d5 was prepared by dilution of internal standard stock solution with methanol and stored at 4°C.

Calibration curves for PDA and HVA were fitted by the 1/x weighted least-squares quadratic. The calibration curves for PDA and HVA ranged from 5 to 5000 ng/ml. All coefficients of determination (R2) of the calibration lines were ≥0.98. The mean accuracy (percentage of true value) of individual calibrators was ≥15%. The lower limit of quantitation, defined as the lowest concentration that could be determined with precision and accuracy of ± 20%, was 5 ng/ml for PDA and HVA.

Sample stability was also determined for QC samples (n = 5) that were extracted and stored in the instrument autosampler under refrigerated conditions (5°C). Samples were stable for 24 hours. Average recovery for QC samples (n = 5) prepared with 1% BSA in PBS (pH 7.4) was 101% and 102% for PDA and HVA, respectively.

Urine and plasma samples were diluted by factors of 1:1, 1:2, 1:10 or 1:100 with 1% BSA in PBS (pH 7.4) to ensure the levels were within the range of the calibration curve. Aliquots (50 µl) of diluted animal samples (urine and plasma), calibrator standards, or QCs were transferred to a 2 ml 96-well plate and mixed for 1 minute at 2000 rpm. The samples were extracted by adding 200 µl of Enal-d5 in methanol (1000 ng/ml) and vortex mixed for 5 minutes at 2000 rpm, and then centrifuged at 4°C for 10 minutes at 4000g. The supernatant (170 µl) was transferred to a 500 µl 96-well plate and dried under heated nitrogen (45°C). Samples were then reconstituted with 50 μl of water:methanol (98:2) containing 0.1% formic acid, vortex mixed for 1 minute at 2000 rpm, followed by 10 minutes of centrifugation at 4000g at 4°C before MS analysis. To avoid PDA and HVA degradation, all samples were protected from direct light exposure during sample preparation and analysis.

LC-MS/MS analyses were carried out on a Waters Acquity UPLC system, consisting of a Acquity binary solvent manager and Acquity sample manager with sample organizer (Waters) coupled to a SCIEX 6500 tandem quadrupole mass spectrometer (Applied Biosystems/MDS SCIEX, Toronto, Canada) equipped with an electrospray ionization source. The analytes (10 µl) were separated on an Acquity UPLC BEH130 C18 (2.1 × 100 mm; 1.7 µm particle size) column and eluted by a gradient program as follows: held 2% B for 0.5 minutes, 2% B–20% B for 2.5 minutes, 20% B–98% B for 1 minute, 98% B–2% B for 0.01 minute, and retained 1 minute for equilibration. The column was heated at 45°C, and the flow rate was 500 µl/min. The mobile phase consisted of an aqueous phase (A: 0.1% formic acid in water) and an organic phase (B: 0.1% formic acid in ACN).

The electrospray ionization source was operated in negative ion mode and its main working parameters were set as follows: ion spray voltage, −4.5 kV; ion source temperature, 550°C; declustering potential, −40 V for PDA and HVA and −150 V for Enal-d5; collision energy, −35 V; entrance potential, −10 V; and collision cell exit potential, −10 V. Multiple reactions monitoring measurements of PDA and HVA analytes were performed using individually optimized cone voltage and collision energy. The multiple reactions monitoring precursor/product ion transitions were as follows: m/z 182 > 138.0 for PDA, 181.0 > 122.0 for HVA, and 380.3 > 114.2 for the internal standard, Enal-d5. The dwell time established for each transition was 50 milliseconds. All peak integration and data processing methods were performed using SCIEX Analyst 1.6.2 (Applied Biosystems/MDS SCIEX).

Pharmacokinetic, Transport, and Statistical Analyses.

The AUC from zero to 24 hours (AUC0–24 h) was calculated using the mixed trapezoidal rule. The AUC from zero to infinity includes AUC0–24 h and one extrapolated to infinity from the last measured concentration. The volume of distribution at steady-state (VdSS) was determined by the noncompartmental method:Embedded Image(1)where AUMC is the area under the curve of the first moment of the concentration-time curve. The CLTOT was calculated from the following equation:

Embedded Image(2)

The pharmacokinetics parameters including AUC0–24 h, AUC, VdSS, and CLTOT for FSM, PROB, PDA, and HVA, following single intravenous administration of PROB and FSM with or without coadministration of PROB were analyzed with a mixed trapezoidal model using the Kinetica program (Thermo Electron, Philadelphia, PA). The CLR was estimated from the following equation:Embedded Image(3)where Xe0–24 h is the cumulative amount of unchanged FSM excreted in urine over 24 hours. The renal extraction ratio (ERR) of FSM was calculated by the following equation:Embedded Image(4)where fu is the fraction of unbound FSM in human plasma reported (i.e., 0.041) (Rane et al., 1978), and GFR is the glomerular filtration rate in cynomolgus monkeys (i.e., 10.4 ml/min) (Davies and Morris, 1993). The fraction of FSM excreted unchanged in the urine (fR) was calculated by dividing Xe0–24 h by the dose. The nonrenal plasma clearance (CLNR) was estimated as the difference between the total plasma and renal clearances. The fraction excreted by nonrenal routes (fNR) was calculated by dividing the nonrenal clearance by the total plasma clearance. Paired Student’s t test was performed to compare pharmacokinetic parameters between groups using GraphPad Prism version 7 (GraphPad Software Inc., San Diego, CA). A P value of less than 0.05 was considered to be statistically significant (*P < 0.05, **P < 0.01, and ***P < 0.001). To compare pharmacokinetic parameters between groups, the data are also reported as the geometric mean ratio with a two-sided 90% confidence interval.

Transport data represent the results from a single study run in triplicate and a minimum of two experiments on different days. The results are reported as mean ± S.D. (n = 3). To estimate the transport kinetics parameters of PDA and HVA into transporter-expressing HEK 293 cells, the transporter-mediated uptake was calculated by subtracting the uptake in MOCK-HEK cells from that in transporter-expressing HEK 293 cells. The following equation was used to estimate the parameters:Embedded Image(5)where V is the rate of uptake measured at the given concentration; Vmax is the maximal rate of uptake; Km represents the Michaelis-Menten constant at which the transport rate is half its maximal value; and [S] is the substrate concentration.

Statistical differences between cell lines or treatments were determined by an unpaired two-tailed Student t test (GraphPad Prism version 7; GraphPad Software Inc.). A P value of less than 0.05 was considered to be statistically significant (*P < 0.05, **P < 0.01, and ***P < 0.001).

Results

Effects of PROB on Pharmacokinetics of FSM in Cynomolgus Monkeys.

The mean FSM plasma concentration-time profiles after intravenous administration of 2 mg/kg of FSM in the absence and presence of PROB (40 mg/kg, i.v.) are illustrated in Fig. 2A. The plasma concentrations of FSM were higher in the presence than in the absence of PROB at all time points and in all animals. This statement was supported by the significant increase in AUC0–24 h (11.5 ± 0.7 and 47.7 ± 5.0 µM⋅h after intravenous administration of 2 mg/kg FSM alone and with 40 mg/kg PROB, respectively; P < 0.01) (Table 1). The geometric mean and 90% confidence interval of the FSM AUC0–24 h ratio was 4.1 (3.6–4.8). Consistently, the CLTOT of FSM was reduced significantly by PROB pretreatment (8.8 ± 0.6 vs. 2.1 ± 0.3 ml/min per kilogram; P < 0.001). However, the administration of PROB did not cause a significant difference in VdSS and terminal elimination half-life of FSM although the mean values were decreased by PROB pretreatment (0.33 ± 0.11 vs. 0.21 ± 0.03 l/kg and 5.9 ± 1.4 vs. 4.7 ± 0.8 hour, respectively; P > 0.05) (Table 1).

Fig. 2.
Fig. 2.

(A) Plasma concentration-time curves of FSM administered intravenously alone (2 mg/kg; open circles) or with PROB (40 mg/kg, i.v.; closed circles). (B) Plasma concentration time curves of PROB administered intravenously alone (40 mg/kg; open triangles) or with FSM (2 mg/kg, i.v.; closed triangles). (C) Effects of PROB on urinary exertion rate of FSM administered intravenously alone (2 mg/kg; open bar) or with PROB (40 mg/kg, i.v.; closed bar). Data are shown as mean ± S.D. (n = 3). **P < 0.01, significantly different from urinary excretion rate in the absence of PROB.

View this table:
TABLE 1

Pharmacokinetic parameters of FSM and PROB in cynomolgus monkeys (n = 3) following intravenous administration of PROB alone (40 mg/kg), FSM alone (2 mg/kg), and PROB with FSM

Pharmacokinetic parameters were determined as described in Materials and Methods. Data are shown as mean ± S.D. (n = 3). PROB was intravenously dosed 30 minutes prior to FSM administration in the coadministration treatment group.

The mean urinary excretion rates of FSM in cynomolgus monkeys are shown in Fig. 2C. A 2.1- to 2.4-fold reduction of FSM urinary excretion in the presence of PROB was observed (P < 0.01). The CLR and extraction ratio of FSM were reduced significantly with PROB pretreatment (3.7 ± 0.6 vs. 0.44 ± 0.12 ml/min per kilogram and 49.5 ± 8.8 vs. 5.9 ± 1.4, respectively; P < 0.01). Furthermore, the CLNR of FSM in the presence of PROB was also reduced markedly (5.0 ± 0.6 vs. 1.7 ± 0.2 ml/min per kilogram, respectively; P < 0.01). Therefore, the reduced CLTOT by PROB pretreatment was not solely due to either decreased CLR or CLNR.

Figure 2B illustrates the PROB concentration in plasma after intravenous PROB alone and concurrently with FSM. The PROB plasma concentrations 24 hours after PROB administration (C24 h) were not significantly different (7.3 ± 5.7 and 8.9 ± 6.5 μM, P > 0.05) (Table 1). Comparing other PROB pharmacokinetic parameters in the absence or presence of FSM indicates that the mean PROB AUC0–24 h, C24 h, CLTOT, and terminal elimination half-life values were almost identical (Table 1).

Identification of Potential Plasma Endogenous Biomarkers of OAT1 and OAT3 by Metabolomics.

To identify potential endogenous plasma probes of OAT1 and OAT3, LC-MS-based metabolomics was used to determine the alterations in plasma concentrations of endogenous compounds between treatments in cynomolgus monkeys. A total of 233 metabolites of known structural identity by matching accurate mass and retention time values with an in-house database of endogenous metabolites were measured, and the concentration values were used to determine statistical significance using the paired Student’s t test (Supplemental Table 1). Of those metabolites monitored, 29 endogenous molecules were identified to be present at concentrations at least 3-fold higher at 1 and/or 3 hours in PROB pretreatment groups (administered alone or with FSM) (Table 2), suggesting that the changes of endogenous metabolites in the plasma were associated with inhibition of monkey OAT1 and OAT3, with associated reduction of OAT-mediated renal clearance. However, only some metabolites were able to return to baseline at 24 hours, which were included as a criterion for transporter biomarker candidacy.

View this table:
TABLE 2

Comparison of plasma concentrations of the 29 selected endogenous metabolites in cynomolgus monkeys (n = 3) following intravenous administration of PROB alone (40 mg/kg), FSM alone (2 mg/kg), and PROB with FSM

Plasma concentrations of approximately 230 endogenous metabolites were determined by liquid chromatography-mass spectrometry (Supplemental Table 1). Each determination was made in samples collected at 1, 3, and 24 hours postdose from three cynomolgus monkeys, and values for compounds with ≥3.0-fold different plasma concentration at any of three time points in the monkeys treated with PROB and FSM are presented. The P values were determined by paired Student’s t test (Supplemental Table 1).

Based on the structure similarities, the identified metabolites can be categorized into three subset groups. The first subset of metabolites are long-chain dicarboxylic acids and derivatives (Table 2), which include TDA [HOOC(CH2)12COOH] and HDA [HOOC(CH2)14COOH]. The second subset consists of small acids that are produced by gut microflora. The metabolites include indole-3 acetic acid, cresol sulfate, phenyl sulfate, phenyllactic acid, indoxyl sulfate, and phenylacetylglycine. The third are amino acids and derivatives, which include hydroxy isovaleric acid, tyrosine, and aspartic acid. The metabolites whose plasma concentrations were significantly increased by PROB pretreatments include PDA, HVA, glucuronic acid, pantothenic acid, xanthurenic acid, and kynurenic acid (Table 2). The administration of PROB caused increases in plasma PDA and HVA concentrations at 1 and 3 hours (approximately 3- to 6-fold), and the concentrations returned to baseline at 24 hours (Table 2). PDA and HVA represent novel types of endogenous biomarkers not known previously to interact with monkey and human OAT1 and OAT3. Therefore, we conducted follow-up experiments to determine the time course of change in plasma PDA and HVA concentrations and transporter-expressing cell uptake to comprehend the sensitivity and specificity of the probes as described subsequently.

Time-Dependent Effects on Plasma PDA, HVA, HDA, and TDA Levels of PROB Pretreatment.

Administration of 40 mg/kg PROB intravenously caused increases in plasma PDA and HVA concentrations in monkeys. The plasma PDA concentration increased by approximately 3-fold over the first 24 hours (113.5 ± 23.3 and 120.0 ± 17.9 vs. 40.5 ± 2.2 µM⋅h, P < 0.05), and gradually declined to the basal level (1.5 ± 0.3 µM) at 24 hours after PROB pretreatment (Fig. 3A; Table 3). The PDA concentrations were greater at any time points after PROB pretreatment compared with those after FSM administration. Similarly, the increase in plasma HVA concentration was persistent. The plasma HVA concentration peaked at 3–4 hours, returning to the baseline (85 ± 32 nM) within 24 hours after PROB pretreatment with approximately 2-fold increase in AUC (6.5 ± 1.0 and 8.7 ± 1.9 vs. 4.1 ± 0.7 µM⋅h, P < 0.05) (Fig. 3B; Table 3). As shown in Table 3, PROB pretreatment significantly decreased CLR of PDA compared with FSM treatment (P < 0.05). The pretreatment also decreased the CLR value of HVA, although the reduction was not statistically significant.

Fig. 3.
Fig. 3.

Effect of PROB on plasma concentration-time curves of PDA (A), HVA (B), TDA (C), and HDA (D) in three cynomolgus monkeys after intravenous administration of PROB alone (40 mg/kg; open squares), FSM alone (2 mg/kg, close triangles), and PROB with FSM (open circles).

View this table:
TABLE 3

Pharmacokinetic parameters of PDA, HVA, TDA, and HDA in three cynomolgus monkeys following intravenous administration of PROB alone (40 mg/kg), FSM alone (2 mg/kg), and PROB with FSM

Pharmacokinetic parameters were determined as described in Materials and Methods. Data are shown as mean ± S.D. (n = 3). PROB was intravenously dosed 30 minutes prior to FSM administration in the coadministration treatment group.

We also detected the effects of PROB on plasma concentrations of TDA and HDA since they are reported as substrates of OAT1 and OAT3. PROB pretreatment increased plasma concentrations of TDA and other long-chain dicarboxylic acids in monkeys (Table 2). Figure 3, C and D show that PROB pretreatments increase plasma concentrations, although the increases were not statistically significant (Table 3).

Measurement of PDA and HVA Uptake in Stable Cell Lines Constitutively Expressing Renal Transporter.

To determine whether PDA and HVA are substrates for major renal drug transporters, the cellular uptake of the molecules was measured at a single concentration (1 and 5 µM for PDA and HVA, respectively) in human OAT1-, OAT2-, OAT3-, OAT4-, OCT2-, MATE1-, MATE2K-, OATP1B1, OATP1B3-, and NTCP-HEK cells after 5-minute incubation. The cellular uptake of PDA and HVA into the HEK cells stably transfected with the cynomolgus monkey renal OATs was not evaluated because in vitro models are not available in the United States (Tahara et al., 2005). The uptake of PDA in OAT1- and OAT3-HEK cells was approximately 55- and 52-fold higher than that in the control cells (1757 ± 703 vs. 30.6 ± 6.7 and 1548 ± 40.0 vs. 29.8 ± 5.9 pmol/mg after 5 minutes, respectively) (Fig. 4A). In addition, the PDA uptake in OAT1- and OAT3-HEK cells is significantly inhibited by PROB, a potent OAT1 and OAT3 inhibitor (Fig. 4A). Furthermore, PDA is a substrate for human OAT4 because the uptake of PDA in OAT4-HEK cells is significantly greater than that in Mock-HEK cells (1.8-fold; P < 0.01) and the uptake is reduced by 1 mM PROB (Fig. 5A). On the other hand, there was no significant OAT2-, OCT2-, MATE1-, and MATE2K-mediated uptake of PDA compared with the control (P > 0.05) (Figs. 4A and 5A). Similarly, the uptake of HVA in OAT1- and OAT3-expressing cells was approximately 45- and 6-fold higher than those in the control cells (8.33 ± 0.40 vs. 0.18 ± 0.02 and 1.09 ± 0.09 vs. 0.18 ± 0.02 pmol/mg after 5 minutes, respectively; P < 0.001) (Fig. 4B). Additionally, there was also significant OAT2-mediated uptake of HVA compared with the control (5-fold; P < 0.001). However, there was no significant OCT2-, MATE1-, and MATE2K-mediated uptake of HVA compared with the control (P > 0.05). Moreover, both PDA and HVA are not substrates for human hepatic transporters OATP1B1, OATP1B3, and NTCP because there are no significant differences in the uptake between the transporter-overexpressing cells and mock cells (P > 0.05) (Supplemental Fig. 2).

Fig. 4.
Fig. 4.

Profiling of the transport of PDA and HVA by major drug transporters expressed at the basolateral membrane of RPTC. Uptake in the HEK cells stably transfected with the control vector (Mock-HEK), OAT1, OAT2, OAT3, or OCT2 was measured after 5-minute incubation at 37°C with PDA (1 µM) (A), HVA (5 µM) (B), and radiolabeled probe substrates [1 µM [3H]para-aminohippuric acid (OAT1), 1 µM [3H]penciclovir (OAT2), 1 µM [3H]estrone-3-sulfate (OAT3), and 2 µM [14C]metformin (OAT2)] (C). Incubations were conducted in the absence and presence of 1 mM PROB (OAT1 and OAT3), 100 µM indomethacin (OAT2), or 100 μM pyrimethanmine (OCT2) to evaluate the effects of these inhibitors on PDA and HVA uptake. Each value represents the mean ± S.D. (n = 3). *P < 0.05; **P < 0.01; ***P < 0.001 statistically significantly different from uptake in Mock-HEK cells, and $P < 0.05; $$P < 0.01; $$$P < 0.001 statistically significantly different from uptake in the absence of an inhibitor.

Fig. 5.
Fig. 5.

Profiling of the transport of PDA and HVA by major drug transporters expressed at the apical membrane of RPTC. Uptake in the HEK cells stably transfected with the control vector (Mock-HEK), OAT4, MATE1, or MATE2K was measured after 5-minute incubation at 37°C with PDA (1 µM) (A), HVA (5 µM) (B), and [3H]estrone-3-sulfate (1 µM) and [14C]metformin (2 µM) (C). Incubations were conducted in the absence and presence of 1 mM PROB (OAT4) or 100 μM pyrimethanmine (MATE1 and MATE2K) to evaluate the effects of these inhibitors on PDA and HVA uptake. Each value represents the mean ± S.D. (n = 3). ***P < 0.001 statistically significantly different from uptake in Mock-HEK cells, and $$$P < 0.001 statistically significantly different from uptake in the absence of an inhibitor.

The affinities of OAT1- and OAT3-meditated transport for PDA and HVA were further determined over a range of PDA and HVA concentrations (0.2–500 µM) in OAT1- and OAT3-HEK cells after 2.5-minute incubation. The incubation time was set by the linear-uptake condition and the lower limit of bioanalytical methods (Supplemental Fig. 1). As shown in Fig. 6, the apparent Km values were 33.0 ± 5.3 and 52.1 ± 15.3 µM for OAT1- and OAT3-mediated uptake of PDA, respectively, while the apparent Km values of OAT1-, OAT-2, and OAT3-mediated uptake of HVA were 108 ± 6.2, 124 ± 10.8, and 438 ± 63 µM, respectively.

Fig. 6.
Fig. 6.

Concentration-dependent uptake of PDA (A) and HVA (B) by human OAT1 and OAT3. Cells were incubated with PDA or HVA (0.23–500 µM) for 2 minutes (linear range). The net OAT1- (open triangle) or OAT3-mediated uptake (open circle) was calculated by subtracting that in mock cells from that in OAT1- or OAT3-HEK cells. The curves represent the best fit of the Michaelis-Menten equation. Each value represents the mean ± S.D. (n = 3).

Discussion

OAT transporters are mainly expressed in RPTCs and mediate active renal secretion of their substrates. Inhibition of OAT transporter function could lead to decreased renal elimination and increased plasma exposure of xenobiotics and endogenous metabolites that are transported by these transporters. The general approach of using cynomolgus monkey as a model for study of human transporter function has become routine and takes advantage of the similarity and function of drug transporters between humans and monkeys (Tahara et al., 2006; Shen et al., 2013, 2016b). It is recognized that such approaches can identify endogenous biomarker candidates that are selective for a specific transporter (Chu et al., 2015; Shen et al., 2016a; Thakare et al., 2017). To demonstrate that cynomolgus monkey permits quantitative prediction of human OAT-mediated DDIs, we measured the effects of the clinically relevant OAT inhibitor PROB on the pharmacokinetics of the OAT1 and OAT3 probe substrate FSM in monkeys, and compared these results with the changes observed in human subjects. As shown in Fig. 2, the intravenous administration of PROB generated a plasma total C24 h of 7.3 ± 5.7 µM (Fig. 2; Table 1). Based on the IC50 values for both monkey OAT1 and OAT3 (<10 µM) (Tahara et al., 2005), sustained inhibition of OAT1 and OAT3 transport function is anticipated. In agreement with the in vitro inhibition, a pronounced increase in AUC of FSM, a known substrate of OAT1 and OAT3, resulting from PROB pretreatment was evident (Fig. 2; Table 1). At the relevant dose (i.e., 40 mg/kg, i.v.), PROB decreased the FSM CLR value by 8.3-fold in cynomolgus monkeys. The decreases in the CLR values of coadministered FSM observed in the clinic were 3.6- to 5.1-fold for PROB (Chennavasin et al., 1979; Smith et al., 1980). Furthermore, PROB decreased FSM CLNR in monkeys and humans by 2.9- and 1.5-fold, respectively, suggesting inhibition of FSM metabolism in both species. In the current study, we demonstrated that intravenous pretreatment of cynomolgus monkeys with PROB (40 mg/kg) increased in the AUC of FSM by 4.1-fold, which is consistent with the values (3.1- to 3.7-fold) reported in humans and indicates that DDIs associated with OAT transporter inhibitions could be reproduced in cynomoglous monkey. These results support previous findings in cynomolgus monkeys where it was found that OAT1 and OAT3 exhibit very similar transport kinetics compared with corresponding human orthologs, and cynomolgus monkeys are recommended as more appropriate alternative systems for predicting DDIs involving renal drug transporters in humans (Tahara et al., 2005, 2006).

In the present study, we applied a combination of fit-for-purpose untargeted metabolomics and quantitative multiple reactions monitoring LC-MS/MS methods (Fig. 1) to determine if OAT inhibition is associated with changes in plasma concentrations of endogenous metabolites that are reflective of OAT function. Of the 233 plasma metabolites examined using the metabolomics method, 29 metabolites including PDA and HVA were significantly increased at 1 or 3 hours in plasma from monkeys by PROB. Of those, a number of metabolites identified had not previously been reported as substrates for OAT1 and OAT3 (Table 2). LC-MS/MS methods were further optimized for determination of plasma PDA and HVA concentration-time profiles. We found PDA and HVA to be consistently elevated in plasma from cynomolgus monkeys pretreated with PROB alone and with FSM (AUC0–24 h of 2.8- to 2.9-fold and 1.6- to 2.1-fold, respectively), similar to the increase in plasma FSM concentration (4.1-fold) (Figs. 2 and 3; Tables 1 and 3). Furthermore, we demonstrated that PDA and HVA are substrates for human OAT1, OAT3, OAT2 (HVA), and OAT4 (PDA) but not OCT2, MATE1, MATE2K, OATP1B1, OATP1B3, and NTCP (Figs. 4 and 5; Supplemental Fig. 2). However, we did not study the transport of PDA and HVA by cynomolgus monkey OAT1 and OAT3 because in vitro cell models are not available in the United States. The species-dependent differences in OAT1- and OAT3-mediated transport of PDA and HVA cannot be excluded, although monkey OAT1 and OAT3 exhibit similar transport kinetics (i.e., Km and Vmax values) to human orthologs (Tahara et al., 2005). Collectively, the results from the present study demonstrate that circulating PDA and HVA can be potentially used as endogenous biomarkers of OAT1 and OAT3 inhibition.

A few years ago, an investigation was conducted examining the changes of plasma and urine metabolites with murine Oat1 (Slc22a6) deficiency. Using untargeted metabolomics analysis, the researchers identified several physiologically important metabolites, including PDA, that were increased in the plasma from Oat1 knockout mice compared with wild-type animals, suggesting they are substrates for mouse Oat1 (Wikoff et al., 2011). In agreement, PROB increased the plasma HVA levels in rhesus monkey (Bacopoulos et al., 1978). PDA has been shown to undergo active renal secretion in humans, and elevated plasma PDA was observed in patients with renal insufficiency (Coburn et al., 2002). It is worth noting that the present study shows a considerable variation in the urinary excretion and CLR of PDA and HVA (Fig. 3, B and D; Table 3). In addition, the volumes of urine samples collected in PROB groups during normal diuresis and forced diuresis (i.e., administration of PROB alone and with FSM) are significantly different. Despite large variation, PROB pretreatments, either alone or with FSM, significantly reduced the CLR values of PDA in monkeys (ratios of 0.10 and 0.33–0.46, respectively), similar to the decrease in FSM CLR [0.12 (0.11–0.13)] (Tables 1 and 3).

Although the plasma levels of HVA are significantly increased by PROB pretreatments compared with FSM treatment, the CLR value was only slightly reduced by PROB pretreatments (ratios of 0.8 and 0.49–0.50, respectively) and the changes were not statistically significant. It is unclear why PROB pretreatments decreased the CLR of HVA to a lesser extent compared with PDA. However, the following observations may shed some light. Based on the in vitro uptake experiments, the renal basolateral transporters OAT1, OAT2, and OAT3 are involved in the uptake of HVA while PDA is only transported by OAT1 and OAT3 (Figs. 4 and 6). OAT2 is less sensitive to the inhibition by PROB compared with OAT1 and OAT3 (IC50 of 393–766 µM vs. <10 µM) (Enomoto et al., 2002; Jia et al., 2015). In contrast to our findings, previous studies examining expression of human OATs in the choroid plexus and their interactions with neurotransmitter metabolites have reported that HVA is not a substrate for human OATs (OAT1 and OAT3) in mouse cells (second segment of the proximal tubule cells, S2 cells) stably expressing OAT1 or OAT3 (Alebouyeh et al., 2003); the difference in transport profiling is probably due to the use of different host cells (HEK 293 cells vs. S2 cells), and endogenous expression of mouse Oat1, Oat2, and Oat3 in S2 host cells would mask HVA uptake. Although we showed high-affinity uptake of PDA and HVA by basolateral transporters OAT1, OAT3, and/or OAT2 into stably transfected HEK cells (Figs. 4 and 6), the mechanism(s) of the export of PDA and HVA out of the proximal tubule cells into the urine is not fully understood. These molecules may be extruded out of the cell by the efflux transporters MRP2 and MRP4 since most MRP2 and MRP4 substrates are organic anions, which are localized in the apical membrane of RPTCs. In addition, the renal apical transporter OAT4 may play a role in the transport of PDA from RPTCs into the urine since OAT4 can operate in efflux mode in addition to influx fashion (Hagos et al., 2007).

PROB is a selective inhibitor of OAT1 and OAT3 since its IC50 values toward other transporters are greater than 25 µM (Supplemental Table 2), although for other monkey transporters there is still paucity of complete data for PROB inhibition. Therefore, even though the results presented herein showcase PDA and HVA as novel cynomolgus monkey plasma OAT biomarkers, it cannot be assumed that the results translate directly to human subjects. In addition, we cannot exclude that the formation of the biomarkers studied may be affected by PROB. Moreover, PDA is the dead-end catabolite of the B6 vitamins including pyridoxine, pyridoxamine, and pyridoxal (Merrill and Henderson, 1990). Although the plasma PDA baseline level exhibits acceptable interindividual variability in monkeys and humans (1.6 ± 0.4 µM and 43 ± 17 nM) (Coburn et al., 2002), the use of vitamin B6 supplements causes elevation of PDA levels far above the normal physiologic level (Zempleni, 1995; Coburn et al., 2002). Therefore, the use of vitamin B6 supplements should be excluded in a clinical study with the goal of assessing a change in PDA plasma level.

Emerging metabolomics and genome-wide association data revealed that HDA and TDA, two fatty acid dicarboxylates, were potential endogenous biomarkers of human OATP1B1 (Yee et al., 2016). Although HDA and TDA are substrates of OATP1B1, OAT1 and OAT3 were also involved in the disposition of HDA and TDA (Yee et al., 2016). Consistently, we have observed that administration of rifampin, a selective inhibitor of hepatic OATPs over renal OATs, increased the AUCs of HDA and TDA by approximately 2-fold in human subjects (Shen et al., 2017). Although these results indicated that HDA and TDA are endogenous biomarkers of OATP1B1, the in vitro transport profiling suggested the involvement of OAT1 and OAT3 in the disposition of HDA and TDA (Yee et al., 2016). Metabolomics analysis in this study demonstrated the significant elevations in plasma levels of at least 20 dicarboxylic acids and fatty acids in the monkeys pretreated with PROB (Supplemental Table 1; Table 2). A similar trend was more evident for TDA with significant increases of AUC in the LC-MS/MS data set (Fig. 3C; Table 3). A comparable degree of increase has been observed with HDA (Fig. 3D; Table 3). As a result, HDA and TDA may serve as dual hepatic OATP and renal OAT biomarkers although further studies are needed for the relevance of clinical interactions. The concentration-time profiles for TDA and HDA appear different in the PROB-alone and with FSM groups although the concentrations in both groups are greater than those in that of FSM alone (Fig. 3, C and D). This is likely due to the additional inhibitory effect of FSM toward OAT1- and OAT3-mediated transport of TDA and HDA (Hasannejad et al., 2004; Nieskens et al., 2016).

In conclusion, the use of monkey as a transporter-mediated DDI model along with metabolomics has been demonstrated to be a useful approach for transporter biomarker identification. Our investigations showed that PDA and HVA are novel blood biomarkers of monkey renal OAT since the changes in plasma exposures of PDA and HVA are similar to that of a probe substrate in monkeys. Although interspecies differences in the transport and disposition of PDA and HVA need be considered, the information obtained in this monkey study on the extent of these endogenous compounds after pretreatment with PROB along with the results from stably transporter-overexpressing cell lines suggest that PDA and HVA are candidate biomarkers of OAT1 and OAT3 for use in clinical settings.

Acknowledgments

We thank Joseph Cantone and Dr. Dieter Drexler for supporting the bioanalysis of PDA and HVA in the OAT4, OATP1B1, OATP1B3, and NTCP uptake experiments.

Authorship Contributions

Participated in research design: Shen, Nelson, Oliveira, Lai, Humphreys.

Conducted experiments: Shen, Oliveira, Zhang, Mcnaney, Gu, Chen, Su.

Contributed new reagents or analytic tools: Shen, Nelson, Oliveira.

Performed data analysis: Shen, Nelson, Oliveira, Zhang, Shipkova, Reily, Gan, Marathe, Humphreys.

Wrote or contributed to the writing of the manuscript: Shen, Nelson, Oliveira, Lai, Marathe, Humphreys.

Footnotes

Abbreviations

ACN
acetonitrile
AUC
area under the plasma concentration-time curve
BSA
bovine serum albumin
CLNR
nonrenal plasma clearance
CLR
renal clearance
CLTOT
total plasma clearance
DDI
drug-drug interaction
Enal-d5
enalapril maleate-d5
FSM
furosemide
HBSS
Hanks’ balanced salt solution
HDA
hexadecanedioic acid
HEK
human embryonic kidney
HILIC
hydrophilic interaction liquid chromatography
HVA
homovanillic acid
LC
liquid chromatography
LC-MS/MS
liquid chromatography–tandem mass spectrometry
MS
mass spectrometry
NTCP
sodium taurocholate cotransporting polypeptide
OAT
organic anion transporter
PBS
phosphate-buffered saline
PDA
pyridoxic acid
PROB
probenecid
QC
quality control
RP
reverse phase
RPTC
renal proximal tubular cell
TDA
tetradecanedioic acid
Vdss
volume of distribution at steady-state

References

  1. ↵
    1. Alebouyeh M,
    2. Takeda M,
    3. Onozato ML,
    4. Tojo A,
    5. Noshiro R,
    6. Hasannejad H,
    7. Inatomi J,
    8. Narikawa S,
    9. Huang XL,
    10. Khamdang S, et al.
    (2003) Expression of human organic anion transporters in the choroid plexus and their interactions with neurotransmitter metabolites. Journal of Pharmacological Sciences 93:430–436.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Bacopoulos NG,
    2. Heninger GR, and
    3. Roth RH
    (1978) Effects of haloperidol and probenecid on plasma and CSF dopamine metabolites in the rhesus monkey (Macacca mulatta). Life Sci 23:1805–1812.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Bergagnini-Kolev MC,
    2. Hebert MF,
    3. Easterling TR, and
    4. Lin YS
    (2017) Pregnancy increases the renal secretion of N(1)-methylnicotinamide, an endogenous probe for renal cation transporters, in patients prescribed metformin. Drug Metab Dispos 45:325–329.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Chennavasin P,
    2. Seiwell R,
    3. Brater DC, and
    4. Liang WM
    (1979) Pharmacodynamic analysis of the furosemide-probenecid interaction in man. Kidney Int 16:187–195.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Chu X,
    2. Chan GH, and
    3. Evers R
    (2017) Identification of endogenous biomarkers to predict the propensity of drug candidates to cause hepatic or renal transporter-mediated drug-drug interactions. J Pharm Sci 106:2357–2367.
    OpenUrl
  6. ↵
    1. Chu X,
    2. Shih SJ,
    3. Shaw R,
    4. Hentze H,
    5. Chan GH,
    6. Owens K,
    7. Wang S,
    8. Cai X,
    9. Newton D,
    10. Castro-Perez J, et al.
    (2015) Evaluation of cynomolgus monkeys for the identification of endogenous biomarkers for hepatic transporter inhibition and as a translatable model to predict pharmacokinetic interactions with statins in humans. Drug Metab Dispos 43:851–863.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Coburn SP,
    2. Reynolds RD,
    3. Mahuren JD,
    4. Schaltenbrand WE,
    5. Wang Y,
    6. Ericson KL,
    7. Whyte MP,
    8. Zubovic YM,
    9. Ziegler PJ,
    10. Costill DL, et al.
    (2002) Elevated plasma 4-pyridoxic acid in renal insufficiency. Am J Clin Nutr 75:57–64.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Davies B and
    2. Morris T
    (1993) Physiological parameters in laboratory animals and humans. Pharm Res 10:1093–1095.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Enomoto A,
    2. Takeda M,
    3. Shimoda M,
    4. Narikawa S,
    5. Kobayashi Y,
    6. Kobayashi Y,
    7. Yamamoto T,
    8. Sekine T,
    9. Cha SH,
    10. Niwa T, et al.
    (2002) Interaction of human organic anion transporters 2 and 4 with organic anion transport inhibitors. J Pharmacol Exp Ther 301:797–802.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Hagos Y,
    2. Stein D,
    3. Ugele B,
    4. Burckhardt G, and
    5. Bahn A
    (2007) Human renal organic anion transporter 4 operates as an asymmetric urate transporter. J Am Soc Nephrol 18:430–439.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Hasannejad H,
    2. Takeda M,
    3. Taki K,
    4. Shin HJ,
    5. Babu E,
    6. Jutabha P,
    7. Khamdang S,
    8. Aleboyeh M,
    9. Onozato ML,
    10. Tojo A, et al.
    (2004) Interactions of human organic anion transporters with diuretics. J Pharmacol Exp Ther 308:1021–1029.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Hnatyshyn S,
    2. Shipkova P, and
    3. Sanders M
    (2013) Expedient data mining for nontargeted high-resolution LC-MS profiles of biological samples. Bioanalysis 5:1195–1210.
    OpenUrl
  13. ↵
    1. Hsueh CH,
    2. Yoshida K,
    3. Zhao P,
    4. Meyer TW,
    5. Zhang L,
    6. Huang SM, and
    7. Giacomini KM
    (2016) Identification and quantitative assessment of uremic solutes as inhibitors of renal organic anion transporters, OAT1 and OAT3. Mol Pharm 13:3130–3140.
    OpenUrl
  14. ↵
    1. Imamura Y,
    2. Tsuruya Y,
    3. Damme K,
    4. Heer D,
    5. Kumagai Y,
    6. Maeda K,
    7. Murayama N,
    8. Okudaira N,
    9. Kurihara A,
    10. Izumi T, et al.
    (2014) 6β-Hydroxycortisol is an endogenous probe for evaluation of drug-drug interactions involving a multispecific renal organic anion transporter, OAT3/SLC22A8, in healthy subjects. Drug Metab Dispos 42:685–694.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Ito S,
    2. Kusuhara H,
    3. Kumagai Y,
    4. Moriyama Y,
    5. Inoue K,
    6. Kondo T,
    7. Nakayama H,
    8. Horita S,
    9. Tanabe K,
    10. Yuasa H, et al.
    (2012) N-methylnicotinamide is an endogenous probe for evaluation of drug-drug interactions involving multidrug and toxin extrusions (MATE1 and MATE2-K). Clin Pharmacol Ther 92:635–641.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Jia W,
    2. Du F,
    3. Liu X,
    4. Jiang R,
    5. Xu F,
    6. Yang J,
    7. Li L,
    8. Wang F,
    9. Olaleye OE,
    10. Dong J, et al.
    (2015) Renal tubular secretion of tanshinol: molecular mechanisms, impact on its systemic exposure, and propensity for dose-related nephrotoxicity and for renal herb-drug interactions. Drug Metab Dispos 43:669–678.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Lacy SA,
    2. Hitchcock MJ,
    3. Lee WA,
    4. Tellier P, and
    5. Cundy KC
    (1998) Effect of oral probenecid coadministration on the chronic toxicity and pharmacokinetics of intravenous cidofovir in cynomolgus monkeys. Toxicol Sci 44:97–106.
  18. ↵
    1. Lai Y,
    2. Mandlekar S,
    3. Shen H,
    4. Holenarsipur VK,
    5. Langish R,
    6. Rajanna P,
    7. Murugesan S,
    8. Gaud N,
    9. Selvam S,
    10. Date O, et al.
    (2016) Coproporphyrins in plasma and urine can be appropriate clinical biomarkers to recapitulate drug-drug interactions mediated by organic anion transporting polypeptide inhibition. J Pharmacol Exp Ther 358:397–404.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Mariappan TT,
    2. Shen H, and
    3. Marathe P
    (2017) Endogenous biomarkers to assess drug-drug interactions by drug transporters and enzymes. Curr Drug Metab 18:757–768.
    OpenUrl
  20. ↵
    1. Merrill AH Jr. and
    2. Henderson JM
    (1990) Vitamin B6 metabolism by human liver. Ann N Y Acad Sci 585:110–117.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Mori S,
    2. Takanaga H,
    3. Ohtsuki S,
    4. Deguchi T,
    5. Kang YS,
    6. Hosoya K, and
    7. Terasaki T
    (2003) Rat organic anion transporter 3 (rOAT3) is responsible for brain-to-blood efflux of homovanillic acid at the abluminal membrane of brain capillary endothelial cells. J Cereb Blood Flow Metab 23:432–440.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Morrissey KM,
    2. Stocker SL,
    3. Wittwer MB,
    4. Xu L, and
    5. Giacomini KM
    (2013) Renal transporters in drug development. Annu Rev Pharmacol Toxicol 53:503–529.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Motohashi H,
    2. Sakurai Y,
    3. Saito H,
    4. Masuda S,
    5. Urakami Y,
    6. Goto M,
    7. Fukatsu A,
    8. Ogawa O, and
    9. Inui K
    (2002) Gene expression levels and immunolocalization of organic ion transporters in the human kidney. J Am Soc Nephrol 13:866–874.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Müller F,
    2. Pontones CA,
    3. Renner B,
    4. Mieth M,
    5. Hoier E,
    6. Auge D,
    7. Maas R,
    8. Zolk O, and
    9. Fromm MF
    (2015) N(1)-methylnicotinamide as an endogenous probe for drug interactions by renal cation transporters: studies on the metformin-trimethoprim interaction. Eur J Clin Pharmacol 71:85–94.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Nieskens TT,
    2. Peters JG,
    3. Schreurs MJ,
    4. Smits N,
    5. Woestenenk R,
    6. Jansen K,
    7. van der Made TK,
    8. Röring M,
    9. Hilgendorf C,
    10. Wilmer MJ, et al.
    (2016) A human renal proximal tubule cell line with stable organic anion transporter 1 and 3 expression predictive for antiviral-induced toxicity. AAPS J 18:465–475.
    OpenUrl
  26. ↵
    1. Nigam SK,
    2. Wu W,
    3. Bush KT,
    4. Hoenig MP,
    5. Blantz RC, and
    6. Bhatnagar V
    (2015) Handling of drugs, metabolites, and uremic toxins by kidney proximal tubule drug transporters. Clin J Am Soc Nephrol 10:2039–2049.
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Peng CC,
    2. Templeton I,
    3. Thummel KE,
    4. Davis C,
    5. Kunze KL, and
    6. Isoherranen N
    (2011) Evaluation of 6β-hydroxycortisol, 6β-hydroxycortisone, and a combination of the two as endogenous probes for inhibition of CYP3A4 in vivo. Clin Pharmacol Ther 89:888–895.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Rane A,
    2. Villeneuve JP,
    3. Stone WJ,
    4. Nies AS,
    5. Wilkinson GR, and
    6. Branch RA
    (1978) Plasma binding and disposition of furosemide in the nephrotic syndrome and in uremia. Clin Pharmacol Ther 24:199–207.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Rodrigues AD,
    2. Taskar KS,
    3. Kusuhara H, and
    4. Sugiyama Y
    (2017) Endogenous probes for drug transporters: balancing vision with reality. Clin Pharmacol Ther [published ahead of print].
  30. ↵
    1. Shen H,
    2. Chen W,
    3. Drexler DM,
    4. Mandlekar S,
    5. Holenarsipur VK,
    6. Shields EE,
    7. Langish R,
    8. Sidik K,
    9. Gan J,
    10. Humphreys WG, et al.
    (2017) Comparative evaluation of plasma bile acids, dehydroepiandrosterone sulfate, hexadecanedioate, and tetradecanedioate with coproporphyrins I and III as markers of OATP inhibition in healthy subjects. Drug Metab Dispos 45:908–919.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Shen H,
    2. Dai J,
    3. Liu T,
    4. Cheng Y,
    5. Chen W,
    6. Freeden C,
    7. Zhang Y,
    8. Humphreys WG,
    9. Marathe P, and
    10. Lai Y
    (2016a) Coproporphyrins I and III as functional markers of OATP1B activity: in vitro and in vivo evaluation in preclinical species. J Pharmacol Exp Ther 357:382–393.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Shen H,
    2. Liu T,
    3. Jiang H,
    4. Titsch C,
    5. Taylor K,
    6. Kandoussi H,
    7. Qiu X,
    8. Chen C,
    9. Sukrutharaj S,
    10. Kuit K, et al.
    (2016b) Cynomolgus monkey as a clinically relevant model to study transport involving renal organic cation transporters: in vitro and in vivo evaluation. Drug Metab Dispos 44:238–249.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Shen H,
    2. Yang Z,
    3. Zhao W,
    4. Zhang Y, and
    5. Rodrigues AD
    (2013) Assessment of vandetanib as an inhibitor of various human renal transporters: inhibition of multidrug and toxin extrusion as a possible mechanism leading to decreased cisplatin and creatinine clearance. Drug Metab Dispos 41:2095–2103.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Smith DE,
    2. Gee WL,
    3. Brater DC,
    4. Lin ET, and
    5. Benet LZ
    (1980) Preliminary evaluation of furosemide-probenecid interaction in humans. J Pharm Sci 69:571–575.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Tahara H,
    2. Kusuhara H,
    3. Chida M,
    4. Fuse E, and
    5. Sugiyama Y
    (2006) Is the monkey an appropriate animal model to examine drug-drug interactions involving renal clearance? Effect of probenecid on the renal elimination of H2 receptor antagonists. J Pharmacol Exp Ther 316:1187–1194.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Tahara H,
    2. Shono M,
    3. Kusuhara H,
    4. Kinoshita H,
    5. Fuse E,
    6. Takadate A,
    7. Otagiri M, and
    8. Sugiyama Y
    (2005) Molecular cloning and functional analyses of OAT1 and OAT3 from cynomolgus monkey kidney. Pharm Res 22:647–660.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Thakare R,
    2. Gao H,
    3. Kosa RE,
    4. Bi YA,
    5. Varma MVS,
    6. Cerny MA,
    7. Sharma R,
    8. Kuhn M,
    9. Huang B,
    10. Liu Y, et al.
    (2017) Leveraging of rifampicin-dosed cynomolgus monkeys to identify bile acid 3-O-sulfate conjugates as potential novel biomarkers for organic anion-transporting polypeptides. Drug Metab Dispos 45:721–733.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Tsuruya Y,
    2. Kato K,
    3. Sano Y,
    4. Imamura Y,
    5. Maeda K,
    6. Kumagai Y,
    7. Sugiyama Y, and
    8. Kusuhara H
    (2016) Investigation of endogenous compounds applicable to drug-drug interaction studies involving the renal organic anion transporters, OAT1 and OAT3, in humans. Drug Metab Dispos 44:1925–1933.
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Vree TB,
    2. van den Biggelaar-Martea M, and
    3. Verwey-van Wissen CP
    (1995) Probenecid inhibits the renal clearance of frusemide and its acyl glucuronide. Br J Clin Pharmacol 39:692–695.
    OpenUrlPubMed
  40. ↵
    1. Wikoff WR,
    2. Nagle MA,
    3. Kouznetsova VL,
    4. Tsigelny IF, and
    5. Nigam SK
    (2011) Untargeted metabolomics identifies enterobiome metabolites and putative uremic toxins as substrates of organic anion transporter 1 (Oat1). J Proteome Res 10:2842–2851.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Yee SW,
    2. Giacomini MM,
    3. Hsueh CH,
    4. Weitz D,
    5. Liang X,
    6. Goswami S,
    7. Kinchen JM,
    8. Coelho A,
    9. Zur AA,
    10. Mertsch K, et al.
    (2016) Metabolomic and genome-wide association studies reveal potential endogenous biomarkers for OATP1B1. Clin Pharmacol Ther 100:524–536.
    OpenUrl
  42. ↵
    1. Zempleni J
    (1995) Pharmacokinetics of vitamin B6 supplements in humans. J Am Coll Nutr 14:579–586.
    OpenUrlPubMed
Back to top

In this issue

Download PDF
Article Alerts
Email Article
Citation Tools
Research ArticleArticle

Identification of Functional Plasma Markers of OAT1 and OAT3

Hong Shen, David M. Nelson, Regina V. Oliveira, Yueping Zhang, Colleen A. Mcnaney, Xiaomei Gu, Weiqi Chen, Ching Su, Michael D. Reily, Petia A. Shipkova, Jinping Gan, Yurong Lai, Punit Marathe and W. Griffith Humphreys
Drug Metabolism and Disposition February 1, 2018, 46 (2) 178-188; DOI: https://doi.org/10.1124/dmd.117.077586
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

Jump to section