Metabolomic and Genome-wide Association Studies Reveal Potential Endogenous Biomarkers for OATP1B1

Utilization of OATP1B Biomarker Coproporphyrin-I to Guide Drug-Drug Interaction Risk Assessment: Evaluation by the Pharmaceutical Industry

Drug-drug interactions (DDIs) involving hepatic organic anion transporting polypeptides 1B1/1B3 (OATP1B) can be substantial, however, challenges remain for predicting interaction risk. Emerging evidence suggests that endogenous biomarkers, particularly coproporphyrin-I (CP-I), can be used to assess in vivo OATP1B activity. The present work under the International Consortium for Innovation and Quality in Pharmaceutical Development was aimed primarily at assessing CP-I as a biomarker for informing OATP1B DDI risk. Literature and unpublished CP-I data along with pertinent in vitro and clinical DDI information were collected to identify DDIs primarily involving OATP1B inhibition and assess the relationship between OATP1B substrate drug and CP-I exposure changes. Static models to predict changes in exposure of CP-I, as a selective OATP1B substrate, were also evaluated. Significant correlations were observed between CP-I area under the curve ratio (AUCR) or maximum concentration ratio (Cmax R) and AUCR of substrate drugs. In general, the CP-I Cmax R was equal to or greater than the CP-I AUCR. CP-I Cmax R < 1.25 was associated with absence of OATP1B-mediated DDIs (AUCR < 1.25) with no false negative predictions. CP-I Cmax R < 2 was associated with weak OATP1B-mediated DDIs (AUCR < 2). A correlation was identified between CP-I exposure changes and OATP1B1 static DDI predictions. Recommendations for collecting and interpreting CP-I data are discussed, including a decision tree for guiding DDI risk assessment. In conclusion, measurement of CP-I is recommended to inform OATP1B inhibition potential. The current analysis identified changes in CP-I exposure that may be used to prioritize, delay, or replace clinical DDI studies.

Progress in the Quantitative Assessment of Transporter-Mediated Drug-Drug Interactions Using Endogenous Substrates in Clinical Studies

Variations in drug transporter activities, caused by genetic polymorphism and drug-drug interactions (DDIs), alter the systemic exposure of substrate drugs, leading to differences in drug responses. Recently, some endogenous substrates of drug transporters, particularly the solute carrier family transporters such as OATP1B1, OATP1B3, OAT1, OAT3, OCT1, OCT2, MATE1, and MATE2-K, have been identified to investigate variations in drug transporters in humans. Clinical data obtained support their performance as surrogate probes in terms of specificity and reproducibility. Pharmacokinetic parameters of the endogenous biomarkers depend on the genotypes of drug transporters and the systemic exposure to perpetrator drugs. Furthermore, the development of physiologically based pharmacokinetic models for the endogenous biomarkers has enabled a top-down approach to obtain insights into the effect of perpetrators on drug transporters and to more precisely simulate the DDI with victim drugs, including probe drugs. The endogenous biomarkers can address the uncertainty in the DDI prediction in the preclinical and early phases of clinical development and have the potential to fulfill regulatory requirements. Therefore, the endogenous biomarkers should be able to predict disease effects on the variations in drug transporter activities observed in patients. This mini-review focuses on recent progress in the identification and use of the endogenous drug transporter substrate biomarkers and their application in drug development. SIGNIFICANCE STATEMENT: Advances in analytical methods have enabled the identification of endogenous substrates of drug transporters. Changes in the pharmacokinetic parameters (Cmax, AUC, or CLR) of these endogenous biomarkers relative to baseline values can serve as a quantitative index to assess variations in drug transporter activities during clinical studies and thereby provide more precise DDI predictions.

Characterization of Bile Acid Sulfate Conjugates as Substrates of Human Organic Anion Transporting Polypeptides

Drug interactions involving the inhibition of hepatic organic anion transporting polypeptides (OATPs) 1B1 and OATP1B3 are considered important. Therefore, we sought to study various sulfated bile acids (BA-S) as potential clinical OATP1B1/3 biomarkers. It was determined that BA-S [e.g., glycochenodeoxycholic acid 3-O-sulfate (GCDCA-S) and glycodeoxycholic acid 3-O-sulfate (GDCA-S)] are substrates of OATP1B1, OATP1B3, and sodium-dependent taurocholic acid cotransporting polypeptide (NTCP) transfected into human embryonic kidney 293 cells, with minimal uptake evident for other solute carriers (SLCs) like OATP2B1, organic anion transporter 2, and organic cation transporter 1. It was also shown that BA-S uptake by plated human hepatocytes (PHH) was inhibited (≥96%) by a pan-SLC inhibitor (rifamycin SV), and there was greater inhibition (≥77% versus ≤12%) with rifampicin (OATP1B1/3-selective inhibitor) than a hepatitis B virus myristoylated-preS1 peptide (NTCP-selective inhibitor). Estrone 3-sulfate was also used as an OATP1B1-selective inhibitor. In this instance, greater inhibition was observed with GDCA-S (76%) than GCDCA-S (52%). The study was expanded to encompass the measurement of GCDCA-S and GDCA-S in plasma of SLCO1B1 genotyped subjects. The geometric mean GDCA-S concentration was 2.6-fold (90% confidence interval 1.6, 4.3; P = 2.1 × 10^-4) and 1.3-fold (1.1, 1.7; P = 0.001) higher in individuals homozygous and heterozygous for the SLCO1B1 c.521T > C loss-of-function allele, respectively. For GCDCA-S, no significant difference was noted [1.2-fold (0.8, 1.7; P = 0.384) and 0.9-fold (0.8, 1.1; P = 0.190), respectively]. This supported the in vitro data indicating that GDCA-S is a more OATP1B1-selective substrate (versus GCDCA-S). It is concluded that GCDCA-S and GDCA-S are viable plasma-based OATP1B1/3 biomarkers, but they are both less OATP1B1-selective when compared to their corresponding 3-O-glucuronides (GCDCA-3G and GDCA-3G). Additional studies are needed to determine their utility versus more established biomarkers, such as coproporphyrin I, for assessing inhibitors with different OATP1B1 (versus OATP1B3) inhibition signatures.

Regulation of Human Endogenous Metabolites by Drug Transporters and Drug Metabolizing Enzymes: An Analysis of Targeted SNP-Metabolite Associations

Drug transporters and drug-metabolizing enzymes are primarily known for their role in the absorption, distribution, metabolism, and excretion (ADME) of small molecule drugs, but they also play a key role in handling endogenous metabolites. Recent cross-tissue co-expression network analyses have revealed a "Remote Sensing and Signaling Network" of multispecific, oligo-specific, and monospecific transporters and enzymes involved in endogenous metabolism. This includes many proteins from families involved in ADME (e.g., SLC22, SLCO, ABCC, CYP, UGT). Focusing on the gut-liver-kidney axis, we identified the endogenous metabolites potentially regulated by this network of ~1000 proteins by associating SNPs in these genes with the circulating levels of thousands of small, polar, bioactive metabolites, including free fatty acids, eicosanoids, bile acids, and other signaling metabolites that act in part via G-protein coupled receptors (GPCRs), nuclear receptors, and kinases. We identified 77 genomic loci associated with 7236 unique metabolites. This included metabolites that were associated with multiple, distinct loci, indicating coordinated regulation between multiple genes (including drug transporters and drug-metabolizing enzymes) of specific metabolites. We analyzed existing pharmacogenomic data and noted SNPs implicated in endogenous metabolite handling (e.g., rs4149056 in SLCO1B1) also affecting drug ADME. The overall results support the existence of close relationships, via interactions with signaling metabolites, between drug transporters and drug-metabolizing enzymes that are part of the Remote Sensing and Signaling Network, and with GPCRs and nuclear receptors. These analyses highlight the potential for drug-metabolite interactions at the interfaces of the Remote Sensing and Signaling Network and the ADME protein network.

OAT, OATP, and MRP Drug Transporters and the Remote Sensing and Signaling Theory

The coordinated movement of organic anions (e.g., drugs, metabolites, signaling molecules, nutrients, antioxidants, gut microbiome products) between tissues and body fluids depends, in large part, on organic anion transporters (OATs) [solute carrier 22 (SLC22)], organic anion transporting polypeptides (OATPs) [solute carrier organic (SLCO)], and multidrug resistance proteins (MRPs) [ATP-binding cassette, subfamily C (ABCC)]. Depending on the range of substrates, transporters in these families can be considered multispecific, oligospecific, or (relatively) monospecific. Systems biology analyses of these transporters in the context of expression patterns reveal they are hubs in networks involved in interorgan and interorganismal communication. The remote sensing and signaling theory explains how the coordinated functions of drug transporters, drug-metabolizing enzymes, and regulatory proteins play a role in optimizing systemic and local levels of important endogenous small molecules. We focus on the role of OATs, OATPs, and MRPs in endogenous metabolism and how their substrates (e.g., bile acids, short chain fatty acids, urate, uremic toxins) mediate interorgan and interorganismal communication and help maintain and restore homeostasis in healthy and disease states.

Role of Hepatocyte Transporters in Drug-Induced Liver Injury (DILI)-In Vitro Testing

Bile acids and bile salts (BA/BS) are substrates of both influx and efflux transporters on hepatocytes. Canalicular efflux transporters, such as BSEP and MRP2, are crucial for the removal of BA/BS to the bile. Basolateral influx transporters, such as NTCP, OATP1B1/1B3, and OSTα/β, cooperate with canalicular transporters in the transcellular vectorial flux of BA/BS from the sinusoids to the bile. The blockage of canalicular transporters not only impairs the bile flow but also causes the intracellular accumulation of BA/BS in hepatocytes that contributes to, or even triggers, liver injury. In the case of BA/BS overload, the efflux of these toxic substances back to the blood via MRP3, MRP4, and OST α/β is considered a relief function. FXR, a key regulator of defense against BA/BS toxicity suppresses de novo bile acid synthesis and bile acid uptake, and promotes bile acid removal via increased efflux. In drug development, the early testing of the inhibition of these transporters, BSEP in particular, is important to flag compounds that could potentially inflict drug-induced liver injury (DILI). In vitro test systems for efflux transporters employ membrane vesicles, whereas those for influx transporters employ whole cells. Additional in vitro pharmaceutical testing panels usually include cellular toxicity tests using hepatocytes, as well as assessments of the mitochondrial toxicity and accumulation of reactive oxygen species (ROS). Primary hepatocytes are the cells of choice for toxicity testing, with HepaRG cells emerging as an alternative. Inhibition of the FXR function is also included in some testing panels. The molecular weight and hydrophobicity of the drug, as well as the steady-state total plasma levels, may positively correlate with the DILI potential. Depending on the phase of drug development, the physicochemical properties, dosing, and cut-off values of BSEP IC₅₀ ≤ 25-50 µM or total C_ss,plasma/BSEP IC₅₀ ≥ 0.1 may be an indication for further testing to minimize the risk of DILI liability.

Effects of SLCO1B1 Genetic Variant on Metabolite Profile in Participants on Simvastatin Treatment

Organic-anion-transporting polypeptide 1B1 (OATP1B1), encoded by the solute carrier organic anion transporter family member 1B1 gene (SLCO1B1), is highly expressed in the liver and transports several endogenous metabolites into the liver, including statins. Previous studies have not investigated the association of SLCO1B1 rs4149056 variant with the risk of type 2 diabetes (T2D) or determined the metabolite signature of the C allele of SLCO1B1 rs4149056 (SLCO1B1 rs4149056-C allele) in a large randomly selected population. SLCO1B1 rs4149056-C inhibits OATP1B1 transporter and is associated with increased levels of blood simvastatin concentrations. Our study is to first to show that SLCO1B1 rs4149056 variant is not significantly associated with the risk of T2D, suggesting that simvastatin has a direct effect on the risk of T2D. Additionally, we investigated the effects of SLCO1B1 rs4149056-C on plasma metabolite concentrations in 1373 participants on simvastatin treatment and in 1368 age- and body-mass index (BMI)-matched participants without any statin treatment. We found 31 novel metabolites significantly associated with SLCO1B1 rs4149056-C in the participants on simvastatin treatment and in the participants without statin treatment. Simvastatin decreased concentrations of dicarboxylic acids, such as docosadioate and dodecanedioate, that may increase beta- and peroxisomal oxidation and increased the turnover of cholesterol into bile acids, resulting in a decrease in steroidogenesis due to limited availability of cholesterol for steroid synthesis. Our findings suggest that simvastatin exerts its effects on the lowering of low-density lipoprotein (LDL) cholesterol concentrations through several distinct pathways in the carriers of SLCO1B1 rs4149056-C, including dicarboxylic acids, bile acids, steroids, and glycerophospholipids.

A Metabolomics Approach for Predicting OATP1B-Type Transporter-Mediated Drug–Drug Interaction Liabilities

In recent years, various endogenous compounds have been proposed as putative biomarkers for the hepatic uptake transporters OATP1B1 and OATP1B3 that have the potential to predict transporter-mediated drug–drug interactions (DDIs). However, these compounds have often been identified from top–down strategies and have not been fully utilized as a substitute for traditional DDI studies. In an attempt to eliminate observer bias in biomarker selection, we applied a bottom–up, untargeted metabolomics screening approach in mice and found that plasma levels of the conjugated bile acid chenodeoxycholate-24-glucuronide (CDCA-24G) are particularly sensitive to deletion of the orthologous murine transporter Oatp1b2 (31-fold increase vs. wild type) or the entire Oatp1a/1b(−/−)cluster (83-fold increased), whereas the humanized transgenic overexpression of hepatic OATP1B1 or OATP1B3 resulted in the partial restoration of transport function. Validation studies with the OATP1B1/OATP1B3 inhibitors rifampin and paclitaxel in vitro as well as in mice and human subjects confirmed that CDCA-24G is a sensitive and rapid response biomarker to dose-dependent transporter inhibition. Collectively, our study confirmed the ability of CDCA-24G to serve as a sensitive and selective endogenous biomarker of OATP1B-type transport function and suggests a template for the future development of biomarkers for other clinically important xenobiotic transporters.

Effect of Cyclosporin A and Impact of Dose Staggering on OATP1B1/1B3 Endogenous Substrates and Drug Probes for Assessing Clinical Drug Interactions

This study was designed to assess the quantitative performance of endogenous biomarkers for organic anion transporting polypeptide (OATP) 1B1/1B3‐mediated drug‐drug interactions (DDIs). Ten healthy volunteers orally received OATP1B1/1B3 probe cocktail (0.2 mg pitavastatin, 1 mg rosuvastatin, and 2 mg valsartan) and an oral dose of cyclosporin A (CysA, 20 mg and 75 mg) separated by a 1‐hour interval (20 mg (−1 hour), and 75 mg (−1 hour)). CysA 75 mg was also given with a 3‐hour interval (75 mg (−3 hours)) to examine the persistence of OATP1B1/1B3 inhibition. The area under the plasma concentration‐time curve ratios (AUCRs) were 1.63, 3.46, and 2.38 (pitavastatin), 1.39, 2.16, and 1.81 (rosuvastatin), and 1.42, 1.77, and 1.85 (valsartan), at 20 mg, 75 mg (−1 hour) and 75 mg (−3 hours) of CysA, respectively. CysA effect on OATP1B1/1B3 was unlikely to persist at the dose examined. Among 26 putative OATP1B1/1B3 biomarkers evaluated, AUCR and maximum concentration ratio (C_maxR) of CP‐I showed the highest Pearson’s correlation coefficient with CysA AUC (0.94 and 0.93, respectively). Correlation between AUCR of pitavastatin, and C_maxR or AUCR of CP‐I were consistent between this study and our previous study using rifampicin as an OATP1B1/1B3 inhibitor. Nonlinear regression analysis of AUCR⁻¹ of pitavastatin and CP‐I against CysA C_max yielded K_{i,OATP1B1/1B3,app} (109 ± 35 and 176 ± 42 nM, respectively), similar to the K_{i,OATP1B1/1B3} estimated by our physiologically‐based pharmacokinetic model analysis described previously (107 nM). The endogenous OATP1B1/1B3 biomarkers, particularly C_maxR and AUCR of CP‐I, corroborates OATP1B1/1B3 inhibition and yields valuable information that improve accurate DDI predictions in drug development, and enhance our understanding of interindividual variability in the magnitude of DDIs.

Genome-Wide Meta-analysis Identifies Genetic Variants Associated With Glycemic Response to Sulfonylureas

OBJECTIVE

Sulfonylureas, the first available drugs for the management of type 2 diabetes, remain widely prescribed today. However, there exists significant variability in glycemic response to treatment. We aimed to establish heritability of sulfonylurea response and identify genetic variants and interacting treatments associated with HbA1c reduction.

RESEARCH DESIGN AND METHODS

As an initiative of the Metformin Genetics Plus Consortium (MetGen Plus) and the DIabetes REsearCh on patient straTification (DIRECT) consortium, 5,485 White Europeans with type 2 diabetes treated with sulfonylureas were recruited from six referral centers in Europe and North America. We first estimated heritability using the generalized restricted maximum likelihood approach and then undertook genome-wide association studies of glycemic response to sulfonylureas measured as HbA1c reduction after 12 months of therapy followed by meta-analysis. These results were supported by acute glipizide challenge in humans who were naïve to type 2 diabetes medications, cis expression quantitative trait loci (eQTL), and functional validation in cellular models. Finally, we examined for possible drug-drug-gene interactions.

RESULTS

After establishing that sulfonylurea response is heritable (mean ± SEM 37 ± 11%), we identified two independent loci near the GXYLT1 and SLCO1B1 genes associated with HbA1c reduction at a genome-wide scale (P < 5 × 10-8). The C allele at rs1234032, near GXYLT1, was associated with 0.14% (1.5 mmol/mol), P = 2.39 × 10-8), lower reduction in HbA1c. Similarly, the C allele was associated with higher glucose trough levels (β = 1.61, P = 0.005) in healthy volunteers in the SUGAR-MGH given glipizide (N = 857). In 3,029 human whole blood samples, the C allele is a cis eQTL for increased expression of GXYLT1 (β = 0.21, P = 2.04 × 10-58). The C allele of rs10770791, in an intronic region of SLCO1B1, was associated with 0.11% (1.2 mmol/mol) greater reduction in HbA1c (P = 4.80 × 10-8). In 1,183 human liver samples, the C allele at rs10770791 is a cis eQTL for reduced SLCO1B1 expression (P = 1.61 × 10-7), which, together with functional studies in cells expressing SLCO1B1, supports a key role for hepatic SLCO1B1 (encoding OATP1B1) in regulation of sulfonylurea transport. Further, a significant interaction between statin use and SLCO1B1 genotype was observed (P = 0.001). In statin nonusers, C allele homozygotes at rs10770791 had a large absolute reduction in HbA1c (0.48 ± 0.12% [5.2 ± 1.26 mmol/mol]), equivalent to that associated with initiation of a dipeptidyl peptidase 4 inhibitor.

CONCLUSIONS

We have identified clinically important genetic effects at genome-wide levels of significance, and important drug-drug-gene interactions, which include commonly prescribed statins. With increasing availability of genetic data embedded in clinical records these findings will be important in prescribing glucose-lowering drugs.

Performance of Plasma Coproporphyrin I and III as OATP1B1 Biomarkers in Humans

A previous study in 356 healthy Finnish volunteers showed that glycochenodeoxycholate 3‐O‐glucuronide (GCDCA‐3G) and glycodeoxycholate 3‐O‐glucuronide (GDCA‐3G) are promising biomarkers of organic anion transporting polypeptide 1B1 (OATP1B1). In the same cohort, we now evaluated the performances of two other OATP1B1 biomarkers, coproporphyrin I (CPI) and III (CPIII), and compared them with GCDCA‐3G and GDCA‐3G. Based on decreased (*5 and *15) and increased (*14 and *20) function SLCO1B1 haplotypes, we stratified the participants to poor, decreased, normal, increased, and highly increased OATP1B1 function groups. Fasting plasma CPI concentration was 68% higher in the poor (95% confidence interval, 44%, 97%; P = 1.74 × 10⁻¹⁰), 7% higher in the decreased (0%, 15%; P = 0.0385), 10% lower in the increased (3%, 18%; P = 0.0087), and 23% lower in the highly increased (1%, 40%; P = 0.0387) function group than in the normal function group. CPIII concentration was 27% higher (7%, 51%; P = 0.0071) in the poor function group than in the normal function group. CPI and CPIII detected poor OATP1B1 function with areas under the precision‐recall curve (AUPRC) of 0.388 (95% confidence interval, 0.197, 0.689) and 0.0798 (0.0485, 0.203), and receiver operating characteristic curve (AUROC) of 0.888 (0.851, 0.919) and 0.731 (0.682, 0.776). The AUPRC and AUROC of GCDCA‐3G were, however, 0.389 (0.258, 0.563) and 0.100 (−0.0046, 0.204; P = 0.0610) larger than those of CPI, and 0.697 (0.555, 0.831) and 0.257 (0.141, 0.373; P < 0.0001) larger than those of CPIII. In conclusion, these data indicate that plasma CPI outperforms CPIII in detecting altered OATP1B1 function, but GCDCA‐3G is an even more sensitive OATP1B1 biomarker than CPI.

Endogenous Biomarkers for SLC Transporter-Mediated Drug-Drug Interaction Evaluation

Membrane transporters play an important role in the absorption, distribution, metabolism, and excretion of xenobiotic substrates, as well as endogenous compounds. The evaluation of transporter-mediated drug-drug interactions (DDIs) is an important consideration during the drug development process and can guide the safe use of polypharmacy regimens in clinical practice. In recent years, several endogenous substrates of drug transporters have been identified as potential biomarkers for predicting changes in drug transport function and the potential for DDIs associated with drug candidates in early phases of drug development. These biomarker-driven investigations have been applied in both preclinical and clinical studies and proposed as a predictive strategy that can be supplanted in order to conduct prospective DDIs trials. Here we provide an overview of this rapidly emerging field, with particular emphasis on endogenous biomarkers recently proposed for clinically relevant uptake transporters.

Transporter-mediated drug-drug interactions: advancement in models, analytical tools, and regulatory perspective

Drug-drug interactions mediated by transporters are a serious clinical concern hence a tremendous amount of work has been done on the characterization of the transporter-mediated proteins in humans and animals. The underlying mechanism for the transporter-mediated drug-drug interaction is the induction or inhibition of the transporter which is involved in the cellular uptake and efflux of drugs. Transporter of the brain, liver, kidney, and intestine are major determinants that alter the absorption, distribution, metabolism, excretion profile of drugs, and considerably influence the pharmacokinetic profile of drugs. As a consequence, transporter proteins may affect the therapeutic activity and safety of drugs. However, mounting evidence suggests that many drugs change the activity and/or expression of the transporter protein. Accordingly, evaluation of drug interaction during the drug development process is an integral part of risk assessment and regulatory requirements. Therefore, this review will highlight the clinical significance of the transporter, their role in disease, possible cause underlying the drug-drug interactions using analytical tools, and update on the regulatory requirement. The recent in-silico approaches which emphasize the advancement in the discovery of drug-drug interactions are also highlighted in this review. Besides, we discuss several endogenous biomarkers that have shown to act as substrates for many transporters, which could be potent determinants to find the drug-drug interactions mediated by transporters. Transporter-mediated drug-drug interactions are taken into consideration in the drug approval process therefore we also provided the extrapolated decision trees from in-vitro to in-vivo, which may trigger the follow-up to clinical studies.

Metabolomic and transcriptomic analysis reveals endogenous substrates and metabolic adaptation in rats lacking Abcg2 and Abcb1a transporters

Abcg2/Bcrp and Abcb1a/Pgp are xenobiotic efflux transporters limiting substrate permeability in the gastrointestinal system and brain, and increasing renal and hepatic drug clearance. The systemic impact of Bcrp and Pgp ablation on metabolic homeostasis of endogenous substrates is incompletely understood. We performed untargeted metabolomics of cerebrospinal fluid (CSF) and plasma, transcriptomics of brain, liver and kidney from male Sprague Dawley rats (WT) and Bcrp/Pgp double knock-out (dKO) rats, and integrated metabolomic/transcriptomic analysis to identify putative substrates and perturbations in canonical metabolic pathways. A predictive Bayesian machine learning model was used to predict in silico those metabolites with greater substrate-like features for either transporters. The CSF and plasma levels of 169 metabolites, nutrients, signaling molecules, antioxidants and lipids were significantly altered in dKO rats, compared to WT rats. These metabolite changes suggested alterations in histidine, branched chain amino acid, purine and pyrimidine metabolism in the dKO rats. Levels of methylated and sulfated metabolites and some primary bile acids were increased in dKO CSF or plasma. Elevated uric acid levels appeared to be a primary driver of changes in purine and pyrimidine biosynthesis. Alterations in Bcrp/Pgp dKO CSF levels of antioxidants, precursors of neurotransmitters, and uric acid suggests the transporters may contribute to the regulation of a healthy central nervous system in rats. Microbiome-generated metabolites were found to be elevated in dKO rat plasma and CSF. The altered dKO metabolome appeared to cause compensatory transcriptional change in urate biosynthesis and response to lipopolysaccharide in brain, oxidation-reduction processes and response to oxidative stress and porphyrin biosynthesis in kidney, and circadian rhythm genes in liver. These findings present insight into endogenous functions of Bcrp and Pgp, the impact that transporter substrates, inhibitors or polymorphisms may have on metabolism, how transporter inhibition could rewire drug sensitivity indirectly through metabolic changes, and identify functional Bcrp biomarkers.

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

Metabolomics and metabolites in ischemic stroke

Stroke is a major reason for disability and the second highest cause of death in the world. When a patient is admitted to a hospital, it is necessary to identify the type of stroke, and the likelihood for development of a recurrent stroke, vascular dementia, and depression. These factors could be determined using different biomarkers. Metabolomics is a very promising strategy for identification of biomarkers. The advantage of metabolomics, in contrast to other analytical techniques, resides in providing low molecular weight metabolite profiles, rather than individual molecule profiles. Technically, this approach is based on mass spectrometry and nuclear magnetic resonance. Furthermore, variations in metabolite concentrations during brain ischemia could alter the principal neuronal functions. Different markers associated with ischemic stroke in the brain have been identified including those contributing to risk, acute onset, and severity of this pathology. In the brain, experimental studies using the ischemia/reperfusion model (IRI) have shown an impaired energy and amino acid metabolism and confirmed their principal roles. Literature data provide a good basis for identifying markers of ischemic stroke and hemorrhagic stroke and understanding metabolic mechanisms of these diseases. This opens an avenue for the successful use of identified markers along with metabolomics technologies to develop fast and reliable diagnostic tools for ischemic and hemorrhagic stroke.

Isobutyrylcarnitine as a Biomarker of OCT1 Activity and Interspecies Differences in its Membrane Transport

Genome-wide association studies have identified an association between isobutyrylcarnitine (IBC) and organic cation transporter 1 (OCT1) genotypes. Higher IBC blood concentrations in humans with active OCT1 genotypes and experimental studies with mouse OCT1 suggested an OCT1-mediated efflux of IBC. In this study, we wanted to confirm the suggested use of IBC as an endogenous biomarker of OCT1 activity and contribute to a better understanding of the mechanisms behind the association between blood concentrations of carnitine derivatives and OCT1 genotype. Blood and urine IBC concentrations were quantified in healthy volunteers regarding intra- and interindividual variation and correlation with OCT1 genotype and with pharmacokinetics of known OCT1 substrates. Furthermore, IBC formation and transport were studied in cell lines overexpressing OCT1 and its naturally occurring variants. Carriers of high-activity OCT1 genotypes had about 3-fold higher IBC blood concentrations and 2-fold higher amounts of IBC excreted in urine compared to deficient OCT1. This was likely due to OCT1 function, as indicated by the fact that IBC correlated with the pharmacokinetics of known OCT1 substrates, like fenoterol, and blood IBC concentrations declined with a 1 h time delay following peak concentrations of the OCT1 substrate sumatriptan. Thus, IBC is a suitable endogenous biomarker reflecting both, human OCT1 (hOCT1) genotype and activity. While murine OCT1 (mOCT1) was an efflux transporter of IBC, hOCT1 exhibited no IBC efflux activity. Inhibition experiments confirmed this data showing that IBC and other acylcarnitines, like butyrylcarnitine, 2-methylbutyrylcarnitine, and hexanoylcarnitine, showed reduced efflux upon inhibition of mOCT1 but not of hOCT1. IBC and other carnitine derivatives are endogenous biomarkers of hOCT1 genotype and phenotype. However, in contrast to mice, the mechanisms underlying the IBC-OCT1 correlation in humans is apparently not directly the OCT1-mediated efflux of IBC. A plausible explanation could be that hOCT1 mediates cellular concentrations of specific regulators or co-substrates in lipid and energy metabolism, which is supported by our in vitro finding that at baseline intracellular IBC concentration is about 6-fold lower alone by OCT1 overexpression.

A key role for the transporter OAT1 in systemic lipid metabolism

Organic anion transporter 1 (OAT1/SLC22A6) is a drug transporter with numerous xenobiotic and endogenous substrates. The Remote Sensing and Signaling Theory suggests that drug transporters with compatible ligand preferences can play a role in “organ crosstalk,” mediating overall organismal communication. Other drug transporters are well known to transport lipids, but surprisingly little is known about the role of OAT1 in lipid metabolism. To explore this subject, we constructed a genome-scale metabolic model using omics data from the Oat1 knockout mouse. The model implicated OAT1 in the regulation of many classes of lipids, including fatty acids, bile acids, and prostaglandins. Accordingly, serum metabolomics of Oat1 knockout mice revealed increased polyunsaturated fatty acids, diacylglycerols, and long-chain fatty acids and decreased ceramides and bile acids when compared with wildtype controls. Some aged knockout mice also displayed increased lipid droplets in the liver when compared with wildtype mice. Chemoinformatics and machine learning analyses of these altered lipids defined molecular properties that form the structural basis for lipid-transporter interactions, including the number of rings, positive charge/volume, and complexity of the lipids. Finally, we obtained targeted serum metabolomics data after short-term treatment of rodents with the OAT-inhibiting drug probenecid to identify potential drug–metabolite interactions. The treatment resulted in alterations in eicosanoids and fatty acids, further supporting our metabolic reconstruction predictions. Consistent with the Remote Sensing and Signaling Theory, the data support a role of OAT1 in systemic lipid metabolism.

Effects of Probenecid on Hepatic and Renal Disposition of Hexadecanedioate, an Endogenous Substrate of Organic Anion Transporting Polypeptide 1B in Rats

The aim of the present study was to investigate changes in plasma concentrations and tissue distribution of endogenous substrates of organic anion transporting polypeptide (OATP) 1B, hexadecanedioate (HDA), octadecanedioate (ODA), tetradecanedioate (TDA), and coproporphyrin-III, induced by its weak inhibitor, probenecid (PBD), in rats. PBD increased the plasma concentrations of these four compounds regardless of bile duct cannulation, whereas liver-to-plasma (K_p,liver) and kidney-to-plasma concentration ratios of HDA and TDA were reduced. Similar effects of PBD on plasma concentrations and K_p,liver of HDA, ODA, and TDA were observed in kidney-ligated rats, suggesting a minor contribution of renal disposition to the overall distribution of these three compounds. Tissue uptake clearance of deuterium-labeled HDA (d-HDA) in liver was 16-fold higher than that in kidney, and was reduced by 80% by PBD. This was compatible with inhibition by PBD of d-HDA uptake in isolated rat hepatocytes. Such inhibitory effects of PBD were also observed in the human OATP1B1-mediated uptake of d-HDA. Overall, the disposition of HDA is mainly determined by hepatic OATP-mediated uptake, which is inhibited by PBD. HDA might, thus, be a biomarker for OATPs minimally affected by urinary and biliary elimination in rats.

Multi-factorial pharmacokinetic interactions: unraveling complexities in precision drug therapy

Introduction: Precision drug therapy requires accounting for pertinent factors in pharmacokinetic (PK) inter-individual variability (i.e., pharmacogenetics, diseases, polypharmacy, and natural product use) that can cause sub-therapeutic or adverse effects. Although each of these individual factors can alter victim drug PK, multi-factorial interactions can cause additive, synergistic, or opposing effects. Determining the magnitude and direction of these complex multi-factorial effects requires understanding the rate-limiting redundant and/or sequential PK processes for each drug.Areas covered: Perturbations in drug-metabolizing enzymes and/or transporters are integral to single- and multi-factorial PK interactions. Examples of single factor PK interactions presented include gene-drug (pharmacogenetic), disease-drug, drug-drug, and natural product-drug interactions. Examples of multi-factorial PK interactions presented include drug-gene-drug, natural product-gene-drug, gene-gene-drug, disease-natural product-drug, and disease-gene-drug interactions. Clear interpretation of multi-factorial interactions can be complicated by study design, complexity in victim drug PK, and incomplete mechanistic understanding of victim drug PK.Expert opinion: Incorporation of complex multi-factorial PK interactions into precision drug therapy requires advances in clinical decision tools, intentional PK study designs, drug-metabolizing enzyme and transporter fractional contribution determinations, systems and computational approaches (e.g., physiologically-based pharmacokinetic modeling), and PK phenotyping of progressive diseases.

Identification of Glycochenodeoxycholate 3-O-Glucuronide and Glycodeoxycholate 3-O-Glucuronide as Highly Sensitive and Specific OATP1B1 Biomarkers

The aim of this study was to investigate the sensitivity and specificity of endogenous glycochenodeoxycholate and glycodeoxycholate 3-O-glucuronides (GCDCA-3G and GDCA-3G) as substrates for organic anion transporting polypeptide 1B1 (OATP1B1) in humans. We measured fasting levels of plasma GCDCA-3G and GDCA-3G using liquid chromatography-tandem mass spectrometry in 356 healthy volunteers. The mean plasma levels of both compounds were ~ 50% lower in women than in men (P = 2.25 × 10^-18 and P = 4.73 × 10^-9 ). In a microarray-based genome-wide association study, the SLCO1B1 rs4149056 (c.521T>C, p.Val174Ala) variation showed the strongest association with the plasma GCDCA-3G (P = 3.09 × 10^-30 ) and GDCA-3G (P = 1.60 × 10^-17 ) concentrations. The mean plasma concentration of GCDCA-3G was 9.2-fold (P = 8.77 × 10^-31 ) and that of GDCA-3G was 6.4-fold (P = 2.45x10^-13 ) higher in individuals with the SLCO1B1 c.521C/C genotype than in those with the c.521T/T genotype. No other variants showed independent genome-wide significant associations with GCDCA-3G or GDCA-3G. GCDCA-3G was highly efficacious in detecting the SLCO1B1 c.521C/C genotype with an area under the receiver operating characteristic curve of 0.996 (P < 0.0001). The sensitivity (98-99%) and specificity (100%) peaked at a cutoff value of 180 ng/mL for men and 90 ng/mL for women. In a haplotype-based analysis, SLCO1B1*5 and *15 were associated with reduced, and SLCO1B1*1B, *14, and *35 with increased OATP1B1 function. In vitro, both GCDCA-3G and GDCA-3G showed at least 6 times higher uptake by OATP1B1 than OATP1B3 or OATP2B1. These data indicate that the hepatic uptake of GCDCA-3G and GDCA-3G is predominantly mediated by OATP1B1. GCDCA-3G, in particular, is a highly sensitive and specific OATP1B1 biomarker in humans.

Current progress in identifying endogenous biomarker candidates for drug transporter phenotyping and their potential application to drug development

Drug transporters play important roles in the elimination of various compounds from the blood. Genetic variation and drug-drug interactions underlie the pharmacokinetic differences for the substrates of drug transporters. Some endogenous substrates of drug transporters have emerged as biomarkers to assess differences in drug transporter activity-not only in animals, but also in humans. Metabolomic analysis is a promising approach for identifying such endogenous substrates through their metabolites. The appropriateness of metabolites is supported by studies in vitro and in vivo, both in animals and through pharmacogenomic or drug-drug interaction studies in humans. This review summarizes current progress in identifying such endogenous biomarkers and applying them to drug transporter phenotyping.

Alteration in the Plasma Concentrations of Endogenous Organic Anion-Transporting Polypeptide 1B Biomarkers in Patients with Non-Small Cell Lung Cancer Treated with Paclitaxel

Paclitaxel has been considered to cause OATP1B-mediated drug-drug interactions at therapeutic doses; however, its clinical relevance has not been demonstrated. This study aimed to elucidate in vivo inhibition potency of paclitaxel against OATP1B1 and OATP1B3 using endogenous OATP1B biomarkers. Paclitaxel is an inhibitor of OATP1B1 and OATP1B3, with K_i of 0.579 ± 0.107 and 5.29 ± 3.87 μM, respectively. Preincubation potentiated its inhibitory effect on both OATP1B1 and OATP1B3, with K_i of 0.154 ± 0.031 and 0.624 ± 0.183 μM, respectively. Ten patients with non-small cell lung cancer who received 200 mg/m² of paclitaxel by a 3-hour infusion were recruited. Plasma concentrations of 10 endogenous OATP1B biomarkers-namely, coproporphyrin I, coproporphyrin III, glycochenodeoxycholate-3-sulfate, glycochenodeoxycholate-3-glucuronide, glycodeoxycholate-3-sulfate, glycodeoxycholate-3-glucuronide, lithocholate-3-sulfate, glycolithocholate-3-sulfate, taurolithocholate-3-sulfate, and chenodeoxycholate-24-glucuronide-were determined in the patients with non-small cell lung cancer on the day before paclitaxel administration and after the end of paclitaxel infusion for 7 hours. Paclitaxel increased the area under the plasma concentration-time curve (AUC) of the endogenous biomarkers 2- to 4-fold, although a few patients did not show any increment in the AUC ratios of lithocholate-3-sulfate, glycolithocholate-3-sulfate, and taurolithocholate-3-sulfate. Therapeutic doses of paclitaxel for the treatment of non-small cell lung cancer (200 mg/m²) will cause significant OATP1B1 inhibition during and at the end of the infusion. This is the first demonstration that endogenous OATP1B biomarkers could serve as surrogate biomarkers in patients. SIGNIFICANCE STATEMENT: Endogenous biomarkers can address practical and ethical issues in elucidating transporter-mediated drug-drug interaction (DDI) risks of anticancer drugs clinically. We could elucidate a significant increment of the plasma concentrations of endogenous OATP1B biomarkers after a 3-hour infusion (200 mg/m²) of paclitaxel, a time-dependent inhibitor of OATP1B, in patients with non-small cell lung cancer. The endogenous OATP1B biomarkers are useful to assess the possibility of OATP1B-mediated DDIs in patients and help in appropriately designing a dosing schedule to avoid the DDIs.

Recent progress in in vivo phenotyping technologies for better prediction of transporter-mediated drug-drug interactions

Clinical reports on transporter-mediated drug-drug interactions (TP-DDIs) have rapidly accumulated and regulatory guidance/guidelines recommend that sponsors consider performing quantitative prediction of TP-DDI risks in the process of drug development. In vitro experiments for characterizing the function of drug transporters have been established and various parameters such as the inhibition constant (K_i) of drugs and the intrinsic uptake/efflux clearance for a certain transporter can be obtained. However, many reports have indicated large discrepancies between the parameters estimated from in vitro experiments and those rationally explaining drug pharmacokinetics. Thus, it is essential to evaluate directly the function of each transporter isoform in vivo in humans. At present, several transporter substrate drugs and endogenous compounds have been recognized as probe substrates for a specific transporter and transporter function was evaluated by monitoring the plasma and urine concentration of those probes; however, few compounds specifically transported via a single transporter isoform have been found. For monitoring the intraorgan concentration of drugs, positron emission tomography can be a powerful tool and clinical examples for quantification of in vivo transporter function have been published. In this review, novel methodologies for in vivo phenotyping of transporter function are summarized.

Dose-Dependent Inhibition of OATP1B by Rifampicin in Healthy Volunteers: Comprehensive Evaluation of Candidate Biomarkers and OATP1B Probe Drugs

To address the most appropriate endogenous biomarker for drug–drug interaction risk assessment, eight healthy subjects received an organic anion transporting polypeptide 1B (OATP1B) inhibitor (rifampicin, 150, 300, and 600 mg), and a probe drug cocktail (atorvastatin, pitavastatin, rosuvastatin, and valsartan). In addition to coproporphyrin I, a widely studied OATP1B biomarker, we identified at least 4 out of 28 compounds (direct bilirubin, glycochenodeoxycholate‐3‐glucuronide, glycochenodeoxycholate‐3‐sulfate, and hexadecanedioate) that presented good sensitivity and dynamic range in terms of the rifampicin dose‐dependent change in area under the plasma concentration‐time curve ratio (AUCR). Their suitability as OATP1B biomarkers was also supported by the good correlation of AUC_0‐24h between the endogenous compounds and the probe drugs, and by nonlinear regression analysis (AUCR⁻¹ vs. rifampicin plasma C_max (maximum total concentration in plasma)) to yield an estimate of the inhibition constant of rifampicin. These endogenous substrates can complement existing OATP1B‐mediated drug–drug interaction risk assessment approaches based on agency guidelines in early clinical trials.

Associations between G6PD, OATP1B1 and BLVRA variants and susceptibility to neonatal hyperbilirubinaemia in a Chinese Han population

AIM

Hyperbilirubinaemia is a common disorder in newborns. The aim of this study was to investigate the associations between G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 variants and the risk of neonatal hyperbilirubinaemia in a Chinese neonate population.

METHODS

A total of 447 Chinese neonates with hyperbilirubinaemia were selected as the study group and 544 healthy subjects were recruited as the control group matched by baseline sex, age, feeding pattern and delivery mode. About 2 mL of peripheral venous blood was taken from all subjects. The single nucleotide polymorphisms (SNPs) of G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 loci were examined by the polymerase chain reaction and Sanger sequencing technique in the peripheral blood of all subjects.

RESULTS

For the G6PD 1388 G>A SNP, individuals carrying the A-allele were associated with a significantly increased risk of neonatal hyperbilirubinaemia (adjusted odds ratio (OR) = 1.49, P < 0.001, 95% confidence interval (CI): 1.31-1.67). This risk increased significantly in the CC genotype carriers at the rs4149056 locus of the SLCO1B1 gene (OR = 2.17, 95% CI: 1.87-2.33), whereas it decreased significantly in individuals carrying the G-allele at the rs699512 locus of the BLVRA gene (adjusted OR = 0.84, P = 0.01, 95% CI: 0.75-0.95). The G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 SNPs had a significant impact on serum total bilirubin levels. Moreover, individuals carrying the A-allele of G6PD 1388 G>A and BLVRA rs699512 had a significantly increased risk of developing neonatal hyperbilirubinaemia (OR = 5.01, P < 0.001, 95% CI: 3.42-7.85).

CONCLUSION

Genetic variants of bilirubin metabolism genes, including G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512, are associated with the risk of neonatal hyperbilirubinaemia, and are potential markers for predicting the disorder.

Unraveling the functional role of the orphan solute carrier, SLC22A24 in the transport of steroid conjugates through metabolomic and genome-wide association studies

Variation in steroid hormone levels has wide implications for health and disease. The genes encoding the proteins involved in steroid disposition represent key determinants of interindividual variation in steroid levels and ultimately, their effects. Beginning with metabolomic data from genome-wide association studies (GWAS), we observed that genetic variants in the orphan transporter, SLC22A24 were significantly associated with levels of androsterone glucuronide and etiocholanolone glucuronide (sentinel SNPs p-value <1x10-30). In cells over-expressing human or various mammalian orthologs of SLC22A24, we showed that steroid conjugates and bile acids were substrates of the transporter. Phylogenetic, genomic, and transcriptomic analyses suggested that SLC22A24 has a specialized role in the kidney and appears to function in the reabsorption of organic anions, and in particular, anionic steroids. Phenome-wide analysis showed that functional variants of SLC22A24 are associated with human disease such as cardiovascular diseases and acne, which have been linked to dysregulated steroid metabolism. Collectively, these functional genomic studies reveal a previously uncharacterized protein involved in steroid homeostasis, opening up new possibilities for SLC22A24 as a pharmacological target for regulating steroid levels.

Membrane transporter data to support kinetically-informed chemical risk assessment using non-animal methods: Scientific and regulatory perspectives

Humans are continuously exposed to low levels of thousands of industrial chemicals, most of which are poorly characterised in terms of their potential toxicity. The new paradigm in chemical risk assessment (CRA) aims to rely on animal-free testing, with kinetics being a key determinant of toxicity when moving from traditional animal studies to integrated in vitro-in silico approaches. In a kinetically informed CRA, membrane transporters, which have been intensively studied during drug development, are an essential piece of information. However, how existing knowledge on transporters gained in the drug field can be applied to CRA is not yet fully understood. This review outlines the opportunities, challenges and existing tools for investigating chemical-transporter interactions in kinetically informed CRA without animal studies. Various environmental chemicals acting as substrates, inhibitors or modulators of transporter activity or expression have been shown to impact TK, just as drugs do. However, because pollutant concentrations are often lower in humans than drugs and because exposure levels and internal chemical doses are not usually known in contrast to drugs, new approaches are required to translate transporter data and reasoning from the drug sector to CRA. Here, the generation of in vitro chemical-transporter interaction data and the development of transporter databases and classification systems trained on chemical datasets (and not only drugs) are proposed. Furtheremore, improving the use of human biomonitoring data to evaluate the in vitro-in silico transporter-related predicted values and developing means to assess uncertainties could also lead to increase confidence of scientists and regulators in animal-free CRA. Finally, a systematic characterisation of the transportome (quantitative monitoring of transporter abundance, activity and maintenance over time) would reinforce confidence in the use of experimental transporter/barrier systems as well as in established cell-based toxicological assays currently used for CRA.

Organic Anion Transporter Polypeptide 1B1 Polymorphism Modulates the Extent of Drug-Drug Interaction and Associated Biomarker Levels in Healthy Volunteers

Understanding transporter‐mediated drug–drug interactions is an integral part of risk assessment in drug development. Recent studies support the use of hexadecanedioate (HDA), tetradecanedioate (TDA), coproporphyrin (CP)‐I, and CP‐III as clinical biomarkers for evaluating organic anion‐transporting polypeptide (OATP)1B1 (SLCO1B1) inhibition. The current study investigated the effect of OATP1B1 genotype c.521T>C (OATP1B1‐Val174Ala) on the extent of interaction between cyclosporin A (CsA) and pravastatin, and associated endogenous biomarkers of the transporter (HDA, TDA, CP‐I, and CP‐III), in 20 healthy volunteers. The results show that the levels of each clinical biomarker and pravastatin were significantly increased in plasma samples of the volunteers following administration of pravastatin plus CsA compared with pravastatin plus placebo. The overall fold change in the area under the concentration–time curve (AUC) and maximum plasma concentration (C_max) was similar among the four biomarkers (1.8–2.5‐fold, paired t‐test P value < 0.05) in individuals who were homozygotes or heterozygotes of the major allele, c.521T. However, the fold change in AUC and C_max for HDA and TDA was significantly abolished in the subjects who were c.521‐CC, whereas the respective fold change in AUC and C_max for pravastatin and CP‐I and CP‐III were slightly weaker in individuals who were c.521‐CC compared with c.521‐TT/TC genotypes. In addition, this study provides the first evidence that SLCO1B1 c.521T>C genotype is significantly associated with CP‐I but not CP‐III levels. Overall, these results suggest that OATP1B1 genotype can modulate the effects of CsA on biomarker levels; the extent of modulation differs among the biomarkers.

Precision Medicine and Personalized Management of Lipoprotein and Lipid Disorders in Chronic and End-Stage Kidney Disease

Precision medicine is an emerging field that calls for individualization of treatment strategies based on characteristics unique to each patient. In lipid management, current guidelines are driven mainly by clinical trial results that presently indicate that patients with non-dialysis-dependent chronic kidney disease (CKD) should be treated with a β-hydroxy β-methylglutaryl-CoA reductase inhibitor, also known as statin therapy. For patients with end-stage kidney disease (ESKD) being treated with hemodialysis, statin therapy has not been shown to successfully reduce poor outcomes in trials and therefore is not recommended. The two major guidelines dissent on whether statin therapy should be of moderate or high intensity in non-dialysis-dependent CKD patients, but often leave the prescribing clinician to make that decision. These decisions often are complicated by the increased concerns for adverse events such as myopathies in patients with advanced kidney disease and ESKD. In the future, there may be an opportunity to further identify CKD and ESKD patients who are more likely to benefit from lipid-modifying therapy as opposed to those who likely will suffer from its side effects using precision medicine tools. For now, data from genetics studies and subgroup analyses may provide insight for future research directions in this field and we review some of the work that has been published in this regard.

Effect of Rifampicin on the Plasma Concentrations of Bile Acid-O-Sulfates in Monkeys and Human Liver-Transplanted Chimeric Mice With or Without Bile Flow Diversion

The present study examined the significance of enterohepatic circulation and the effect of rifampicin [an inhibitor of organic anion-transporting polypeptide 1B (OATP1B)] on the plasma concentrations of bile acid-O-sulfates (glycochenodeoxycholate-O-sulfate, lithocholate-O-sulfate, glycolithocholate-O-sulfate, and taurolithocholate-O-sulfate) in monkeys and human liver-transplanted chimeric mice (PXB mouse). Rifampicin significantly increased the area under the curve of bile acid-O-sulfates in monkeys (13-69 times) and PXB mice (13-25 times) without bile flow diversion. Bile flow diversion reduced the concentration of plasma bile acid-O-sulfates under control conditions in monkeys and the concentration of plasma glycochenodeoxycholate-O-sulfate in PXB mice. It also diminished diurnal variation of plasma lithocholate-O-sulfate, glycolithocholate-O-sulfate, and taurolithocholate-O-sulfate in PXB mice under control conditions. Bile flow diversion did not affect the plasma concentration of bile acid-O-sulfates in monkeys and PXB mice treated with rifampicin. Plasma coproporphyrin I and III levels were constant in monkeys throughout the study, even with bile flow diversion. This study demonstrated that bile acid-O-sulfates are endogenous OATP1B biomarkers in monkeys and PXB mice. Enterohepatic circulation can affect the baseline levels of plasma bile acid-O-sulfates and modify the effect of OATP1B inhibition.

A prospective study of serum metabolites and risk of ischemic stroke

OBJECTIVE

To identify promising blood-based biomarkers and novel etiologic pathways of disease risk, we applied an untargeted serum metabolomics profiling in a community-based prospective study of ischemic stroke (IS).

METHODS

In 3,904 men and women from the Atherosclerosis Risk In Communities study, Cox proportional hazard models were used to estimate the association of incident IS with the standardized level of 245 fasting serum metabolites individually, adjusting for age, sex, race, field center, batch, diabetes, hypertension, current smoking status, body mass index, and estimated glomerular filtration rate. Validation of results was carried out in an independent sample of 114 IS cases and 112 healthy controls.

RESULTS

Serum levels of 2 long-chain dicarboxylic acids, tetradecanedioate and hexadecanedioate, were strongly correlated (r = 0.88) and were associated with incident IS after adjusting for covariates (hazard ratio [95% confidence interval (CI)] 1.11 [1.06-1.16] and 1.12 [1.07-1.17], respectively; p < 0.0001). Analyses by IS subtypes suggested that these associations were specific to cardioembolic stroke (CES). Associations of tetradecanedioate and hexadecanedioate with IS were independently confirmed (odds ratio [95% CI] 1.76 [1.21; 2.56] and 1.60 [1.11; 2.32], respectively).

CONCLUSION

Two serum long-chain dicarboxylic acids, metabolic products of ω-oxidation of fatty acids, were associated with IS and CES independently of known risk factors. Pathways related to intracellular hexadecanedioate synthesis or those involved in its clearance from the circulation may mediate IS risk. These results highlight the potential of metabolomics to discover novel circulating biomarkers for stroke and to unravel novel pathways for IS and its subtypes.

[Quantitative assessment of the risk of OATP1B1/1B3-mediated drug-drug interactions]

Drug transporters play important roles in determining drug pharmacokinetics. Organic anion transporting polypeptides 1B1/1B3 (OATP1B1/1B3) are transporters mediating hepatic uptake of various anionic drugs. OATP1B1/1B3 activities are changed by genetic mutation and drug-drug interaction (DDI) that could lead to severe adverse reactions. Methods to address the precise DDI risk assessment have been developed in addition to the translational assessment from the results of in vitro studies. Using endogenous substrates as probes is an emerging approach that allows clinical assessment of the DDI risk in the early phase of drug development. Then, the clinical data will be subjected to the pharmacokinetic analysis using physiologically-based pharmacokinetic models to perform the more realistic DDI risk assessment with OATP1B1/1B3 substrate drugs. When drug targets are located inside the hepatocytes, DDI impact on the intrahepatic concentration is critical for their pharmacological actions. Positron emission tomography (PET) allows researchers to determine tissue concentration time profiles of the PET probe upon the inhibition of OATP1B1/1B3, and to estimate the change in kinetic parameter for each intrinsic process of hepatic elimination of PET probe. Integration of the clinical data into the PBPK model realizes more precise prediction of DDI impact on the pharmacokinetics of drugs, and their therapeutic effects.

Biomarkers for In Vivo Assessment of Transporter Function

Drug-drug interactions are a major concern not only during clinical practice, but also in drug development. Due to limitations of in vitro-in vivo predictions of transporter-mediated drug-drug interactions, multiple clinical Phase I drug-drug interaction studies may become necessary for a new molecular entity to assess potential drug interaction liabilities. This is a resource-intensive process and exposes study participants, who frequently are healthy volunteers without benefit from study treatment, to the potential risks of a new drug in development. Therefore, there is currently a major interest in new approaches for better prediction of transporter-mediated drug-drug interactions. In particular, researchers in the field attempt to identify endogenous compounds as biomarkers for transporter function, such as hexadecanedioate, tetradecanedioate, coproporphyrins I and III, or glycochenodeoxycholate sulfate for hepatic uptake via organic anion transporting polypeptide 1B or N¹-methylnicotinamide for multidrug and toxin extrusion protein-mediated renal secretion. We summarize in this review the currently proposed biomarkers and potential limitations of the substances identified to date. Moreover, we suggest criteria based on current experiences, which may be used to assess the suitability of a biomarker for transporter function. Finally, further alternatives and supplemental approaches to classic drug-drug interaction studies are discussed.

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

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

PBPK Modeling of Coproporphyrin I as an Endogenous Biomarker for Drug Interactions Involving Inhibition of Hepatic OATP1B1 and OATP1B3

The aim of the present study was to establish a physiologically based pharmacokinetic (PBPK) model for coproporphyrin I (CP-I), a biomarker supporting the prediction of drug-drug interactions (DDIs) involving hepatic organic anion transporting polypeptide 1B (OATP1B), using clinical DDI data with an OATP1B inhibitor rifampicin (300 and 600 mg, orally). The in vivo inhibition constants of rifampicin used as initial input parameters for OATP1Bs (K_i,u,OATP1Bs ) and multidrug resistance-associated protein two-mediated biliary excretion were estimated as 0.23 and 0.87 μM, respectively, from previous reports. Sensitivity analysis demonstrated that the K_i,u,OATP1Bs and biosynthesis rate of CP-I affected the magnitude of the interaction. K_i,u,OATP1Bs values optimized by nonlinear least-squares fitting were ~0.5-fold of the initial value. It was determined that the blood concentration-time profiles of four statins were well-predicted using corrected individual K_i,u,OATP1B values (ratio of in vitro K_i,u(statin) /in vitro K_i,u(CP-I) ). In conclusion, PBPK modeling of CP-I supports dynamic prediction of OATP1B-mediated DDIs.

LC-MS/MS bioanalysis of plasma 1, 14-tetradecanedioic acid and 1, 16-hexadecanedioic acid as candidate biomarkers for organic anion-transporting polypeptide mediated drug-drug interactions

AIM

A robust LC-MS/MS assay was developed to quantify endogenous 1, 14-tetradecanedioic acid (TDA) and 1, 16-hexadecanedioic acid (HDA) in human plasma as potential biomarkers for evaluating drug-drug interactions mediated by the hepatic drug transporters, organic anion-transporting polypeptides.

RESULTS

This assay was validated using fit-for-purpose approach over standard curve range of 2.5-1000 nM for TDA and HDA using analyte-free charcoal-stripped human plasma as the surrogate matrix. Chromatographic separation condition was successfully optimized to separate TDA from an interference peak while maintaining both analytes in neutral forms to minimize carryover issue.

CONCLUSION

The described assay is currently applied to a clinical study for evaluating TDA/HDA as potential substitute biomarkers for drug-drug interaction studies.

Quantitative Prediction of OATP-Mediated Drug-Drug Interactions With Model-Based Analysis of Endogenous Biomarker Kinetics

Quantitative prediction of the magnitude of transporter‐mediated clinical drug‐drug interactions (DDIs) solely from in vitro inhibition data remains challenging. The objective of the present work was to analyze the kinetic profile of an endogenous biomarker for organic anion‐transporting polypeptides 1B (OATP1B), coproporphyrin I (CPI), and to predict clinical DDIs with a probe OATP1B substrate (pravastatin) based on “in vivo” inhibition constants (K_i). The CPI kinetics in the presence and absence of strong and weak OATP1B inhibitors (rifampin and GDC‐0810) were described well with a one‐compartment model, and in vivo K_i were estimated. Clinical DDIs between pravastatin and these inhibitors were predicted using physiologically based pharmacokinetic (PBPK) models coupled with the estimated in vivo K_i and predicted magnitude matched well with the observed DDIs. In conclusion, model‐based analysis of the CPI profile has the potential to quantitatively predict liability of a new molecular entity (NME) as an OATP1B inhibitor early in drug development.

Genetic variation of SORBS1 gene is associated with glucose homeostasis and age at onset of diabetes: A SAPPHIRe Cohort Study

The SORBS1 gene plays an important role in insulin signaling. We aimed to examine whether common single-nucleotide polymorphisms (SNPs) of SORBS1 are associated with prevalence and incidence of diabetes, age at onset of diabetes, and the related traits of glucose homeostasis. A total of 1135 siblings from 492 ethnic Chinese families were recruited at baseline, and 630 were followed up for 5.19 ± 0.96 years. Nine SNPs including rs7081076, rs2281939, rs3818540, rs2274490, rs61739184, rs726176, rs2296966, rs17849148, and rs3193970 were genotyped and examined. To deal with correlated data of subjects within the same families, the generalized estimating equations approach was applied throughout all association analyses. The GG genotype of rs2281939 was associated with a higher risk of diabetes at baseline, an earlier onset of diabetes, and higher steady-state plasma glucose levels in the modified insulin suppression test. The minor allele T of rs2296966 was associated with higher prevalence and incidence of diabetes, an earlier onset of diabetes, and higher 2-h glucose during oral glucose tolerance test. These two SNPs revealed independent associations with age of diabetes onset as well as risk of diabetes at baseline. These findings supported that SORBS1 gene participates in the pathogenesis of diabetes.

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

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

Comparative Study of the Dose-Dependence of OATP1B Inhibition by Rifampicin Using Probe Drugs and Endogenous Substrates in Healthy Volunteers

PURPOSE

To evaluate association of the dose-dependent effect of rifampicin, an OATP1B inhibitor, on the plasma concentration-time profiles among OATP1B substrates drugs and endogenous substrates.

METHODS

Eight healthy volunteers received atorvastatin (1 mg), pitavastatin (0.2 mg), rosuvastatin (0.5 mg), and fluvastatin (2 mg) alone or with rifampicin (300 or 600 mg) in a crossover fashion. The plasma concentrations of these OATP1B probe drugs, total and direct bilirubin, glycochenodeoxycholate-3-sulfate (GCDCA-S), and coproporphyrin I, were determined.

RESULTS

The most striking effect of 600 mg rifampicin was on atorvastatin (6.0-times increase) and GCDCA-S (10-times increase). The AUC_0-24h of atorvastatin was reasonably correlated with that of pitavastatin (r² = 0.73) and with the AUC_0-4h of fluvastatin (r² = 0.62) and sufficiently with the AUC_0-24h of rosuvastatin (r² = 0.32). The AUC_0-24h of GCDCA-S was reasonably correlated with those of direct bilirubin (r² = 0.74) and coproporphyrin I (r² = 0.78), and sufficiently with that of total bilirubin (r² = 0.30). The AUC_0-24h of GCDCA-S, direct bilirubin, and coproporphyrin I were reasonably correlated with that of atorvastatin (r² = 0.48-0.70) [corrected].

CONCLUSION

These results suggest that direct bilirubin, GCDCA-S, and coproporphyrin I are promising surrogate probes for the quantitative assessment of potential OATP1B-mediated DDI.