Non-targeted metabolomics for the identification of plasma metabolites associated with organic anion transporting polypeptide 1B1 function.

Our aim was to evaluate biomarkers for organic anion transporting polypeptide 1B1 (OATP1B1) function using a hypothesis-free metabolomics approach. We analyzed fasting plasma samples from 356 healthy volunteers using non-targeted metabolite profiling by liquid chromatography high-resolution mass spectrometry. Based on SLCO1B1 genotypes, we stratified the volunteers to poor, decreased, normal, increased, and highly increased OATP1B1 function groups. Linear regression analysis, and random forest (RF) and gradient boosted decision tree (GBDT) regressors were used to investigate associations of plasma metabolite features with OATP1B1 function. Of the 9152 molecular features found, 39 associated with OATP1B1 function either in the linear regression analysis (p < 10-5) or the RF or GBDT regressors (Gini impurity decrease > 0.01). Linear regression analysis showed the strongest associations with two features identified as glycodeoxycholate 3-O-glucuronide (GDCA-3G; p = 1.2 × 10-20 for negative and p = 1.7 × 10-19 for positive electrospray ionization) and one identified as glycochenodeoxycholate 3-O-glucuronide (GCDCA-3G; p = 2.7 × 10-16). In both the RF and GBDT models, the GCDCA-3G feature showed the strongest association with OATP1B1 function, with Gini impurity decreases of 0.40 and 0.17. In RF, this was followed by one GDCA-3G feature, an unidentified feature with a molecular weight of 809.3521, and the second GDCA-3G feature. In GBDT, the second and third strongest associations were observed with the GDCA-3G features. Of the other associated features, we identified with confidence two representing lysophosphatidylethanolamine 22:5. In addition, one feature was putatively identified as pregnanolone sulfate and one as pregnenolone sulfate. These results confirm GCDCA-3G and GDCA-3G as robust OATP1B1 biomarkers in human plasma.

Genome-Wide Association Study of Atorvastatin Pharmacokinetics: Associations With SLCO1B1, UGT1A3, and LPP.

In a genome-wide association study of atorvastatin pharmacokinetics in 158 healthy volunteers, the SLCO1B1 c.521T>C (rs4149056) variant associated with increased area under the plasma concentration-time curve from time zero to infinity (AUC0-∞) of atorvastatin (P = 1.2 × 10-10), 2-hydroxy atorvastatin (P = 4.0 × 10-8), and 4-hydroxy atorvastatin (P = 2.9 × 10-8). An intronic LPP variant, rs1975991, associated with reduced atorvastatin lactone AUC0-∞ (P = 3.8 × 10-8). Three UGT1A variants linked with UGT1A3*2 associated with increased 2-hydroxy atorvastatin lactone AUC0-∞ (P = 3.9 × 10-8). Furthermore, a candidate gene analysis including 243 participants suggested that increased function SLCO1B1 variants and decreased activity CYP3A4 variants affect atorvastatin pharmacokinetics. Compared with individuals with normal function SLCO1B1 genotype, atorvastatin AUC0-∞ was 145% (90% confidence interval: 98-203%; P = 5.6 × 10-11) larger in individuals with poor function, 24% (9-41%; P = 0.0053) larger in those with decreased function, and 41% (16-59%; P = 0.016) smaller in those with highly increased function SLCO1B1 genotype. Individuals with intermediate metabolizer CYP3A4 genotype (CYP3A4*2 or CYP3A4*22 heterozygotes) had 33% (14-55%; P = 0.022) larger atorvastatin AUC0-∞ than those with normal metabolizer genotype. UGT1A3*2 heterozygotes had 16% (5-25%; P = 0.017) smaller and LPP rs1975991 homozygotes had 34% (22-44%; P = 4.8 × 10-5) smaller atorvastatin AUC0-∞ than noncarriers. These data demonstrate that genetic variation in SLCO1B1, UGT1A3, LPP, and CYP3A4 affects atorvastatin pharmacokinetics. This is the first study to suggest that LPP rs1975991 may reduce atorvastatin exposure. [Correction added on 6 April, after first online publication: An incomplete sentence ("= 0.017) smaller in heterozygotes for UGT1A3*2 and 34% (22%, 44%; P × 10-5) smaller in homozygotes for LPP noncarriers.") has been corrected in this version.].

Solanidine is a sensitive and specific dietary biomarker for CYP2D6 activity.

BACKGROUND: Individual assessment of CYP enzyme activities can be challenging. Recently, the potato alkaloid solanidine was suggested as a biomarker for CYP2D6 activity. Here, we aimed to characterize the sensitivity and specificity of solanidine as a CYP2D6 biomarker among Finnish volunteers with known CYP2D6 genotypes. RESULTS: Using non-targeted metabolomics analysis, we identified 9152 metabolite features in the fasting plasma samples of 356 healthy volunteers. Machine learning models suggested strong association between CYP2D6 genotype-based phenotype classes with a metabolite feature identified as solanidine. Plasma solanidine concentration was 1887% higher in genetically poor CYP2D6 metabolizers (gPM) (n = 9; 95% confidence interval 755%, 4515%; P = 1.88 × 10-11), 74% higher in intermediate CYP2D6 metabolizers (gIM) (n = 89; 27%, 138%; P = 6.40 × 10-4), and 35% lower in ultrarapid CYP2D6 metabolizers (gUM) (n = 20; 64%, - 17%; P = 0.151) than in genetically normal CYP2D6 metabolizers (gNM; n = 196). The solanidine metabolites m/z 444 and 430 to solanidine concentration ratios showed even stronger associations with CYP2D6 phenotypes. Furthermore, the areas under the receiver operating characteristic and precision-recall curves for these metabolic ratios showed equal or better performances for identifying the gPM, gIM, and gUM phenotype groups than the other metabolites, their ratios to solanidine, or solanidine alone. In vitro studies with human recombinant CYP enzymes showed that solanidine was metabolized mainly by CYP2D6, with a minor contribution from CYP3A4/5. In human liver microsomes, the CYP2D6 inhibitor paroxetine nearly completely (95%) inhibited the metabolism of solanidine. In a genome-wide association study, several variants near the CYP2D6 gene associated with plasma solanidine metabolite ratios. CONCLUSIONS: These results are in line with earlier studies and further indicate that solanidine and its metabolites are sensitive and specific biomarkers for measuring CYP2D6 activity. Since potato consumption is common worldwide, this biomarker could be useful for evaluating CYP2D6-mediated drug-drug interactions and to improve prediction of CYP2D6 activity in addition to genotyping.

Screening of 16 major drug glucuronides for time-dependent inhibition of nine drug-metabolizing CYP enzymes - detailed studies on CYP3A inhibitors.

Time-dependent inhibition of cytochrome P450 (CYP) enzymes has been observed for a few glucuronide metabolites of clinically used drugs. Here, we investigated the inhibitory potential of 16 glucuronide metabolites towards nine major CYP enzymes in vitro. Automated substrate cocktail methods were used to screen time-dependent inhibition of CYP1A2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D6, 2J2 and 3A in human liver microsomes. Seven glucuronides (carvedilol ß-D-glucuronide, diclofenac acyl-ß-D-glucuronide, 4-hydroxyduloxetine ß-D-glucuronide, ezetimibe phenoxy-ß-D-glucuronide, raloxifene 4'-glucuronide, repaglinide acyl-ß-D-glucuronide and valproic acid ß-D-glucuronide) caused NADPH- and time-dependent inhibition of at least one of the CYPs investigated, including CYP2A6, CYP2C19 and CYP3A. In more detailed experiments, we focused on the glucuronides of carvedilol and diclofenac, which inhibited CYP3A. Carvedilol ß-D-glucuronide showed weak time-dependent inhibition of CYP3A, but the parent drug carvedilol was found to be a more potent inhibitor of CYP3A, with the half-maximal inhibitor concentration (IC50) decreasing from 7.0 to 1.1 µM after a 30-min preincubation with NADPH. The maximal inactivation constant (kinact) and the inhibitor concentration causing half of kinact (KI) for CYP3A inactivation by carvedilol were 0.051 1/min and 1.8 µM, respectively. Diclofenac acyl-ß-D-glucuronide caused time-dependent inactivation of CYP3A at high concentrations, with a 4-fold IC50 shift (from 400 to 98 µM after a 30-min preincubation with NADPH) and KI and kinact values of >2,000 µM and >0.16 1/min. In static predictions, carvedilol caused significant (>1.25-fold) increase in the exposure of the CYP3A substrates midazolam and simvastatin. In conclusion, we identified several glucuronide metabolites with CYP inhibitory properties. Based on detailed experiments, the inactivation of CYP3A by carvedilol may cause clinically significant drug-drug interactions.

Ticagrelor Increases Exposure to the Breast Cancer Resistance Protein Substrate Rosuvastatin.

Ticagrelor and rosuvastatin are often used concomitantly after atherothrombotic events. Several cases of rhabdomyolysis during concomitant ticagrelor and rosuvastatin have been reported, suggesting a drug-drug interaction. We showed recently that ticagrelor inhibits breast cancer resistance protein (BCRP) and organic anion transporting polypeptide (OATP) 1B1, 1B3, and 2B1-mediated rosuvastatin transport in vitro. The aim of this study was to investigate the effects of ticagrelor on rosuvastatin pharmacokinetics in humans. In a randomized, crossover study, 9 healthy volunteers ingested a single dose of 90 mg ticagrelor or placebo, followed by a single 10 mg dose of rosuvastatin 1 hour later. Ticagrelor 90 mg or placebo were additionally administered 12, 24, and 36 hours after their first dose. Ticagrelor increased rosuvastatin area under the plasma concentration-time curve (AUC) and peak plasma concentration 2.6-fold (90% confidence intervals: 1.8-3.8 and 1.7-4.0, P = 0.001 and P = 0.003), and prolonged its half-life from 3.1 to 6.6 hours (P = 0.009). Ticagrelor also decreased the renal clearance of rosuvastatin by 11% (3%-19%, P = 0.032). The N-desmethylrosuvastatin:rosuvastatin AUC0-10h ratio remained unaffected by ticagrelor. Ticagrelor had no effect on the plasma concentrations of the endogenous OATP1B substrates glycodeoxycholate 3-O-glucuronide, glycochenodeoxycholate 3-O-glucuronide, glycodeoxycholate 3-O-sulfate, and glycochenodeoxycholate 3-O-sulfate, or the sodium-taurocholate cotransporting polypeptide substrate taurocholic acid. These data indicate that ticagrelor increases rosuvastatin concentrations more than twofold in humans, probably mainly by inhibiting intestinal BCRP. Because the risk for rosuvastatin-induced myotoxicity increases along with rosuvastatin plasma concentrations, using ticagrelor concomitantly with high doses of rosuvastatin should be avoided.

Posaconazole-ibrutinib interaction cannot be avoided by staggered dosing: How to optimize ibrutinib dose during posaconazole treatment.

AIMS: Ibrutinib is used in the treatment of certain B-cell malignancies. Due to its CYP3A4-mediated metabolism and highly variable pharmacokinetics, it is prone to potentially harmful drug-drug interactions. METHODS: In a randomized, placebo-controlled, three-phase crossover study, we examined the effect of the CYP3A4-inhibiting antifungal posaconazole on ibrutinib pharmacokinetics. Eleven healthy participants ingested repeated doses of 300 mg of posaconazole either in the morning or in the evening, or placebo. A single dose of ibrutinib (30, 70 or 140 mg, respectively) was administered at 9 AM, 1 or 12 h after the preceding posaconazole/placebo dose. RESULTS: On average, morning posaconazole increased the dose-adjusted geometric mean area under the plasma concentration-time curve from zero to infinity (AUC0-∞ ) and peak plasma concentration (Cmax ) of ibrutinib 9.5-fold (90% confidence interval [CI] 6.3-14.3, P < 0.001) and 8.5-fold (90% CI 5.7-12.8, P < 0.001), respectively, while evening posaconazole increased those 10.3-fold (90% CI 6.7-16.0, P < 0.001) and 8.2-fold (90% CI 5.2-13.2, P < 0.001), respectively. Posaconazole had no significant effect on the half-life of ibrutinib, but substantially reduced the metabolite PCI-45227 to ibrutinib AUC0-∞ ratio. There were no significant differences in ibrutinib pharmacokinetics between morning and evening posaconazole phases. CONCLUSIONS: Posaconazole increases ibrutinib exposure substantially, by about 10-fold. This interaction cannot be avoided by dosing the drugs 12 h apart. In general, a 70-mg daily dose of ibrutinib should not be exceeded during posaconazole treatment to avoid potentially toxic systemic ibrutinib concentrations.

Inhibition of efflux transporters by poly ADP-ribose polymerase inhibitors.

Poly ADP-ribose polymerase (PARP) inhibitors have been approved for the treatment of various cancers. They share a similar mechanism of action but have differences in pharmacokinetic characteristics and potential for drug-drug interactions (DDI). This study evaluated the potential ATP-binding cassette transporter-mediated interactions between PARP inhibitors (niraparib, olaparib and rucaparib) and statins (atorvastatin and rosuvastatin). We studied the inhibitory activity of PARP inhibitors on breast cancer resistance protein (BCRP), multidrug resistance-associated protein 3 (MRP3) and P-glycoprotein (P-gp) using vesicular transport assays and determined the concentrations required for 50% inhibition (IC50 ). Then, we predicted the increase of statin exposure followed by the administration of PARP inhibitors using a mechanistic static model. Rucaparib was the strongest inhibitor of BCRP-mediated rosuvastatin transport (IC50 13.7 µM), followed by niraparib (42.6 µM) and olaparib (216 µM). PARP inhibitors did not affect MRP3. While niraparib appeared to inhibit P-gp, the inhibition showed large variability. The inhibition of intestinal BCRP by rucaparib, niraparib and olaparib was predicted to elevate rosuvastatin exposure by 52%, 37% and 24%, respectively. The interactions between PARP inhibitors and rosuvastatin are probably of minor clinical significance alone, but combined with other predisposing factors, they may increase the risk of rosuvastatin-associated adverse effects.

The effect of hydroxychloroquine on cholesterol metabolism in statin treated patients after myocardial infarction.

Background and aims: To evaluate the effect of hydroxychloroquine (HCQ) on serum and lipoprotein lipids and serum biomarkers of cholesterol synthesis and absorption in myocardial infarction patients with a high-dose statin. Methods: Myocardial infarction patients (n = 59) with a constant statin dose were randomized to receive hydroxychloroquine 300 mg (n = 31) or placebo (n = 28) daily for six months and followed up for one year. Results: Statin reduced total-c (-26 ± 22% in hydroxychloroquine and -28 ± 19% in placebo group, P = 0.931), LDL-c (-38 ± 26% vs. -44 ± 23%, respectively, P = 0.299), and cholesterol synthesis biomarkers zymostenol, desmosterol, and lathosterol ratios from baseline to one year (e.g., serum lathosterol ratio -17 ± 45% vs. -15 ± 41%, respectively, P < 0.001 for both, P = 0.623 between groups). Compensatorily, cholesterol absorption increased during the intervention (e.g., serum campesterol ratio 125 ± 90% vs. 113 ± 72%, respectively, P < 0.001 for both, P = 0.488 between groups). Hydroxychloroquine did not affect cholesterol concentrations or cholesterol absorption. It prevented the statin-induced increase in cholesterol precursor, desmosterol ratio, from six months to one year in the hydroxychloroquine group (P = 0.007 at one year compared to placebo). Conclusions: Combined with a high-dose statin, hydroxychloroquine had no additional effect on serum cholesterol concentration or cholesterol absorption. However, the findings suggest that hydroxychloroquine interferes with lanosterol synthesis, and thereafter, it temporarily interferes with the cholesterol synthesis pathway, best seen in halting the increase of the desmosterol ratio.Trial Registration ClinicalTrials.gov Identifier: NCT02648464.

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.

Extraesophageal reflux and reflux aspiration in dogs with respiratory diseases and in healthy dogs.

BACKGROUND: Salivary bile acids are used to diagnose extraesophageal reflux (EER) and to evaluate the risk of reflux aspiration that is associated with respiratory diseases in dogs. OBJECTIVES: To study total bile acid (TBA) concentrations in saliva and in bronchoalveolar lavage fluid (BALF) to investigate EER and reflux aspiration in dogs with respiratory diseases and in healthy dogs. ANIMALS:  Thirty-one West Highland White Terriers (WHWTs) with idiopathic pulmonary fibrosis (IPF), 12 dogs with inflammatory airway disease (IAD), 6 dogs with recurrent pneumonia (RP), 26 brachycephalic dogs (BD), 27 healthy WHWTs (HW), 52 healthy dogs (HD). All privately-owned dogs. METHODS: Saliva and BALF were collected from dogs in each group. RESULTS: Salivary TBA concentrations were higher in IPF (median 0.1692 µM, interquartile range [IQR] 0.1115-0.2925 µM, Cohen's d 3.4, 95% confidence interval [CI] 2.2-4.0, P < .001) and BD (0.0256 µM, IQR 0.0086-0.0417 µM, d 0.5, CI -0.1 to 1.1, P = .003) compared to HD (0 µM, IQR not quantifiable [n.q.]-0.0131 µM). Bronchoalveolar lavage fluid TBA concentrations were higher in IPF (0.0117 µM, IQR 0.0048-0.0361 µM, d 0.5, CI 0-1.1, P < .001) compared to HD (0 µM, IQR n.q.-0.0074 µM). CONCLUSION AND CLINICAL IMPORTANCE: Extraesophageal reflux and reflux aspiration occur in healthy dogs and those with respiratory diseases.

Novel inhibitors of breast cancer resistance protein (BCRP, ABCG2) among marketed drugs.

Drug-drug interactions (DDIs) are a major concern for the safe use of medications. Breast cancer resistance protein (BCRP) is a clinically relevant ATP-binding cassette (ABC) transporter for drug disposition. Inhibition of BCRP increases the plasma concentrations of BCRP substrate drugs, which potentially could lead to adverse drug reactions. The aim of the present study was to identify BCRP inhibitors amongst a library of 232 commonly used drugs and anticancer drugs approved by the United States Food and Drug Administration (FDA). BCRP inhibition studies were carried out using the vesicular transport assay. We found 75 drugs that reduced the relative transport activity of BCRP to less than 25% of the vehicle control and were categorized as strong inhibitors. The concentration required for 50% inhibition (IC50) was determined for 13 strong inhibitors that were previously poorly characterized for BCRP inhibition. The IC50 ranged from 1.1 to 11 µM, with vemurafenib, dabigatran etexilate and everolimus being the strongest inhibitors. According to the drug interaction guidance documents from the FDA and the European Medicines Agency (EMA), in vivo DDI studies are warranted if the theoretical intestinal luminal concentration of a drug exceeds its IC50 by tenfold. Here, the IC50 values for eight of the drugs were 100-fold lower than their theoretical intestinal luminal concentration. Moreover, a mechanistic static model suggested that vemurafenib, bexarotene, dabigatran etexilate, rifapentine, aprepitant, and ivacaftor could almost fully inhibit intestinal BCRP, increasing the exposure of concomitantly administered rosuvastatin over 90%. Therefore, clinical studies are warranted to investigate whether these drugs cause BCRP-mediated DDIs in humans.

A comprehensive pharmacogenomic study indicates roles for SLCO1B1, ABCG2 and SLCO2B1 in rosuvastatin pharmacokinetics.

AIMS: The aim was to comprehensively investigate the effects of genetic variability on the pharmacokinetics of rosuvastatin. METHODS: We conducted a genome-wide association study and candidate gene analyses of single dose rosuvastatin pharmacokinetics in a prospective study (n = 159) and a cohort of previously published studies (n = 88). RESULTS: In a genome-wide association meta-analysis of the prospective study and the cohort of previously published studies, the SLCO1B1 c.521 T > C (rs4149056) single nucleotide variation (SNV) associated with increased area under the plasma concentration-time curve (AUC) and peak plasma concentration of rosuvastatin (P = 1.8 × 10-12 and P = 3.2 × 10-15 ). The candidate gene analysis suggested that the ABCG2 c.421C > A (rs2231142) SNV associates with increased rosuvastatin AUC (P = .0079), while the SLCO1B1 c.388A > G (rs2306283) and SLCO2B1 c.1457C > T (rs2306168) SNVs associate with decreased rosuvastatin AUC (P = .0041 and P = .0076). Based on SLCO1B1 genotypes, we stratified the participants into poor, decreased, normal, increased and highly increased organic anion transporting polypeptide (OATP) 1B1 function groups. The OATP1B1 poor function phenotype associated with 2.1-fold (90% confidence interval 1.6-2.8, P = 4.69 × 10-5 ) increased AUC of rosuvastatin, whereas the OATP1B1 highly increased function phenotype associated with a 44% (16-62%; P = .019) decreased rosuvastatin AUC. The ABCG2 c.421A/A genotype associated with 2.2-fold (1.5-3.0; P = 2.6 × 10-4 ) increased AUC of rosuvastatin. The SLCO2B1 c.1457C/T genotype associated with 28% decreased rosuvastatin AUC (11-42%; P = .01). CONCLUSION: These data suggest roles for SLCO1B1, ABCG2 and SLCO2B1 in rosuvastatin pharmacokinetics. Poor SLCO1B1 or ABCG2 function genotypes may increase the risk of rosuvastatin-induced myotoxicity. Reduced doses of rosuvastatin are advisable for patients with these genotypes.

Genomewide Association Study of Simvastatin Pharmacokinetics.

We investigated genetic determinants of single-dose simvastatin pharmacokinetics in a prospective study of 170 subjects and a retrospective cohort of 59 healthy volunteers. In a microarray-based genomewide association study with the prospective data, the SLCO1B1 c.521T>C (p.Val174Ala, rs4149056) single nucleotide variation showed the strongest, genomewide significant association with the area under the plasma simvastatin acid concentration-time curve (AUC; P = 6.0 × 10-10 ). Meta-analysis with the retrospective cohort strengthened the association (P = 1.6 × 10-17 ). In a stepwise linear regression candidate gene analysis among all 229 participants, SLCO1B1 c.521T>C (P = 1.9 × 10-13 ) and CYP3A4 c.664T>C (p.Ser222Pro, rs55785340, CYP3A4*2, P = 0.023) were associated with increased simvastatin acid AUC. Moreover, the SLCO1B1 c.463C>A (p.Pro155Thr, rs11045819, P = 7.2 × 10-6 ) and c.1929A>C (p.Leu643Phe, rs34671512, P = 5.3 × 10-4 ) variants associated with decreased simvastatin acid AUC. Based on these results and the literature, we classified the volunteers into genotype-predicted OATP1B1 and CYP3A4 phenotype groups. Compared with the normal OATP1B1 function group, simvastatin acid AUC was 273% larger in the poor (90% confidence interval (CI), 137%, 488%; P = 3.1 × 10-6 ), 40% larger in the decreased (90% CI, 8%, 83%; P = 0.036), and 67% smaller in the highly increased function group (90% CI, 46%, 80%; P = 2.4 × 10-4 ). Intermediate CYP3A4 metabolizers (i.e., heterozygous carriers of either CYP3A4*2 or CYP3A4*22 (rs35599367)), had 87% (90% CI, 39%, 152%, P = 6.4 × 10-4 ) larger simvastatin acid AUC than normal metabolizers. These data suggest that in addition to no function SLCO1B1 variants, increased function SLCO1B1 variants and reduced function CYP3A4 variants may affect the pharmacokinetics, efficacy, and safety of simvastatin. Care is warranted if simvastatin is prescribed to patients carrying decreased function SLCO1B1 or CYP3A4 alleles.

Functional in vitro characterization of SLCO1B1 variants and simulation of the clinical pharmacokinetic impact of impaired OATP1B1 function.

Organic Anion Transporting Polypeptide 1B1 is important to the hepatic elimination and distribution of many drugs. If OATP1B1 function is decreased, it can increase plasma exposure of e.g. several statins leading to increased risk of muscle toxicity. First, we examined the impact of three naturally occurring rare variants and the frequent SLCO1B1 c.388A>G variant on in vitro transport activity with cellular uptake assay using two substrates: 2', 7'-dichlorofluorescein (DCF) and rosuvastatin. Secondly, LC-MS/MS based quantitative targeted absolute proteomics measured the OATP1B1 protein abundance in crude membrane fractions of HEK293 cells over-expressing these single nucleotide variants. Additionally, we simulated the effect of impaired OATP1B1 function on rosuvastatin pharmacokinetics to estimate the need for genotype-guided dosing. R57Q impaired DCF and rosuvastatin transport significantly yet did not change protein expression considerably, while N130D and N151S did not alter activity but increased protein expression. R253Q did not change protein expression but reduced DCF uptake and increased rosuvastatin Km. Based on pharmacokinetic simulations, doses of 30 mg (with 50% OATP1B1 function) and 20 mg (with 0% OATP1B1 function) result in plasma exposure similar to 40 mg dose (with 100% OATP1B1 function). Therefore dose reductions might be considered to avoid increased plasma exposure caused by function-impairing OATP1B1 genetic variants, such as R57Q.

Circulating oxylipin and bile acid profiles of dexmedetomidine, propofol, sevoflurane, and S-ketamine: a randomised controlled trial using tandem mass spectrometry.

Background: This exploratory study aimed to investigate whether dexmedetomidine, propofol, sevoflurane, and S-ketamine affect oxylipins and bile acids, which are functionally diverse molecules with possible connections to cellular bioenergetics, immune modulation, and organ protection. Methods: In this randomised, open-label, controlled, parallel group, Phase IV clinical drug trial, healthy male subjects (n=160) received equipotent doses (EC50 for verbal command) of dexmedetomidine (1.5 ng ml-1; n=40), propofol (1.7 µg ml-1; n=40), sevoflurane (0.9% end-tidal; n=40), S-ketamine (0.75 µg ml-1; n=20), or placebo (n=20). Blood samples for tandem mass spectrometry were obtained at baseline, after study drug administration at 60 and 130 min from baseline; 40 metabolites were analysed. Results: Statistically significant changes vs placebo were observed in 62.5%, 12.5%, 5.0%, and 2.5% of analytes in dexmedetomidine, propofol, sevoflurane, and S-ketamine groups, respectively. Data are presented as standard deviation score, 95% confidence interval, and P-value. Dexmedetomidine induced wide-ranging decreases in oxylipins and bile acids. Amongst others, 9,10-dihydroxyoctadecenoic acid (DiHOME) -1.19 (-1.6; -0.78), P<0.001 and 12,13-DiHOME -1.22 (-1.66; -0.77), P<0.001 were affected. Propofol elevated 9,10-DiHOME 2.29 (1.62; 2.96), P<0.001 and 12,13-DiHOME 2.13 (1.42; 2.84), P<0.001. Analytes were mostly unaffected by S-ketamine. Sevoflurane decreased tauroursodeoxycholic acid (TUDCA) -2.7 (-3.84; -1.55), P=0.015. Conclusions: Dexmedetomidine-induced oxylipin alterations may be connected to pathways associated with organ protection. In contrast to dexmedetomidine, propofol emulsion elevated DiHOMEs, oxylipins associated with acute respiratory distress syndrome, and mitochondrial dysfunction in high concentrations. Further research is needed to establish the behaviour of DIHOMEs during prolonged propofol/dexmedetomidine infusions and to verify the sevoflurane-induced reduction in TUDCA, a suggested neuroprotective agent. Clinical trial registration: NCT02624401.

Pharmacogenomics of celiprolol - evidence for a role of P-glycoprotein and organic anion transporting polypeptide 1A2 in celiprolol pharmacokinetics.

The aim of this study was to search for associations of genetic variants with celiprolol pharmacokinetics in a large set of pharmacokinetic genes, and, more specifically, in a set of previously identified candidate genes ABCB1, SLCO1A2, and SLCO2B1. To this end, we determined celiprolol single-dose (200 mg) pharmacokinetics and sequenced 379 pharmacokinetic genes in 195 healthy volunteers. Analysis with 46,064 common sequence variants in the 379 genes did not identify any novel genes associated with celiprolol exposure. The candidate gene analysis showed that the ABCB1 c.3435T>C and c.2677T/G>A, and the SLCO1A2 c.516A>C variants were associated with reduced celiprolol area under the plasma concentration-time curve (AUC0-∞ ). An alternative analysis with ABCB1 haplotypes showed that, in addition to SLCO1A2 c.516A>C, three ABCB1 haplotypes were associated with reduced celiprolol AUC0-∞ . A genotype scoring system was developed based on these variants and applied to stratify the participants to low and high celiprolol exposure genotype groups. The mean AUC0-∞ of celiprolol in the low exposure genotype group was 55% of the mean AUC0-∞ in the high exposure group (p = 1.08 × 10-11 ). In addition, the results showed gene-gene interactions in the effects of SLCO1A2 and ABCB1 variants on celiprolol AUC0-∞ (p < 5 × 10-6 ) suggesting an interplay between organic anion transporting polypeptide 1A2 and P-glycoprotein in celiprolol absorption. Taken together, these data indicate that P-glycoprotein and organic anion transporting polypeptide 1A2 play a role in celiprolol pharmacokinetics. Furthermore, patients with ABCB1 and SLCO1A2 genotypes associated with low celiprolol exposure may have an increased risk of poor blood-pressure lowering response to celiprolol.

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-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.

Comparative Hepatic and Intestinal Efflux Transport of Statins.

Previous studies have shown that lipid-lowering statins are transported by various ATP-binding cassette (ABC) transporters. However, because of varying methods, it is difficult to compare the transport profiles of statins. Therefore, we investigated the transport of 10 statins or statin metabolites by six ABC transporters using human embryonic kidney cell-derived membrane vesicles. The transporter protein expression levels in the vesicles were quantified with liquid chromatography-tandem mass spectrometry and used to scale the measured clearances to tissue levels. In our study, apically expressed breast cancer resistance protein (BCRP) and P-glycoprotein (P-gp) transported atorvastatin, fluvastatin, pitavastatin, and rosuvastatin. Multidrug resistance-associated protein 3 (MRP3) transported atorvastatin, fluvastatin, pitavastatin, and, to a smaller extent, pravastatin. MRP4 transported fluvastatin and rosuvastatin. The scaled clearances suggest that BCRP contributes to 87%-91% and 84% of the total active efflux of rosuvastatin in the small intestine and the liver, respectively. For atorvastatin, the corresponding values for P-gp-mediated efflux were 43%-79% and 66%, respectively. MRP3, on the other hand, may contribute to 23%-26% and 25%-37% of total active efflux of atorvastatin, fluvastatin, and pitavastatin in jejunal enterocytes and liver hepatocytes, respectively. These data indicate that BCRP may play an important role in limiting the intestinal absorption and facilitating the biliary excretion of rosuvastatin and that P-gp may restrict the intestinal absorption and mediate the biliary excretion of atorvastatin. Moreover, the basolateral MRP3 may enhance the intestinal absorption and sinusoidal hepatic efflux of several statins. Taken together, the data show that statins differ considerably in their efflux transport profiles. SIGNIFICANCE STATEMENT: This study characterized and compared the transport of atorvastatin, fluvastatin, pitavastatin, pravastatin, rosuvastatin, and simvastatin acid and four atorvastatin metabolites by six ABC transporters (BCRP, MRP2, MRP3, MRP4, MRP8, P-gp). Based on in vitro findings and protein abundance data, the study concludes that BCRP, MRP3, and P-gp have a major impact in the efflux of various statins. Together with in vitro metabolism, uptake transport, and clinical data, our findings are applicable for use in comparative systems pharmacology modeling of statins.

Comparative Hepatic and Intestinal Metabolism and Pharmacodynamics of Statins.

This study aimed to comprehensively investigate the in vitro metabolism of statins. The metabolism of clinically relevant concentrations of atorvastatin, fluvastatin, pitavastatin, pravastatin, rosuvastatin, simvastatin, and their metabolites were investigated using human liver microsomes (HLMs), human intestine microsomes (HIMs), liver cytosol, and recombinant cytochrome P450 enzymes. We also determined the inhibitory effects of statin acids on their pharmacological target, 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase. In HLMs, statin lactones were metabolized to a much higher extent than their acid forms. Atorvastatin lactone and simvastatin (lactone) showed extensive metabolism [intrinsic clearance (CLint) values of 3700 and 7400 µl/min per milligram], whereas the metabolism of the lactones of 2-hydroxyatorvastatin, 4-hydroxyatorvastatin, and pitavastatin was slower (CLint 20-840 µl/min per milligram). The acids had CLint values in the range <0.1-80 µl/min per milligram. In HIMs, only atorvastatin lactone and simvastatin (lactone) exhibited notable metabolism, with CLint values corresponding to 20% of those observed in HLMs. CYP3A4/5 and CYP2C9 were the main statin-metabolizing enzymes. The majority of the acids inhibited HMG-CoA reductase, with 50% inhibitory concentrations of 4-20 nM. The present comparison of the metabolism and pharmacodynamics of the various statins using identical methods provides a strong basis for further application, e.g., comparative systems pharmacology modeling. SIGNIFICANCE STATEMENT: The present comparison of the in vitro metabolic and pharmacodynamic properties of atorvastatin, fluvastatin, pitavastatin, pravastatin, rosuvastatin, and simvastatin and their metabolites using unified methodology provides a strong basis for further application. Together with in vitro drug transporter and clinical data, the present findings are applicable for use in comparative systems pharmacology modeling to predict the pharmacokinetics and pharmacological effects of statins at different dosages.

An automated cocktail method for in vitro assessment of direct and time-dependent inhibition of nine major cytochrome P450 enzymes - application to establishing CYP2C8 inhibitor selectivity.

We developed an in vitro high-throughput cocktail assay with nine major drug-metabolizing CYP enzymes, optimized for screening of time-dependent inhibition. The method was applied to determine the selectivity of the time-dependent CYP2C8 inhibitors gemfibrozil 1-O-ß-glucuronide and clopidogrel acyl-ß-D-glucuronide. In vitro incubations with CYP selective probe substrates and pooled human liver microsomes were conducted in 96-well plates with automated liquid handler techniques and metabolite concentrations were measured with quantitative UHPLC-MS/MS analysis. After determination of inter-substrate interactions and Km values for each reaction, probe substrates were divided into cocktails I (tacrine/CYP1A2, bupropion/CYP2B6, amodiaquine/CYP2C8, tolbutamide/CYP2C9 and midazolam/CYP3A4/5) and II (coumarin/CYP2A6, S-mephenytoin/CYP2C19, dextromethorphan/CYP2D6 and astemizole/CYP2J2). Time-dependent inhibitors (furafylline/CYP1A2, selegiline/CYP2A6, clopidogrel/CYP2B6, gemfibrozil 1-O-ß-glucuronide/CYP2C8, tienilic acid/CYP2C9, ticlopidine/CYP2C19, paroxetine/CYP2D6 and ritonavir/CYP3A) and direct inhibitor (terfenadine/CYP2J2) showed similar inhibition with single substrate and cocktail methods. Established time-dependent inhibitors caused IC50 fold shifts ranging from 2.2 to 30 with the cocktail method. Under time-dependent inhibition conditions, gemfibrozil 1-O-ß-glucuronide was a strong (>90% inhibition) and selective (<< 20% inhibition of other CYPs) inhibitor of CYP2C8 at concentrations ranging from 60 to 300 µM, while the selectivity of clopidogrel acyl-ß-D-glucuronide was limited at concentrations above its IC80 for CYP2C8. The time-dependent IC50 values of these glucuronides for CYP2C8 were 8.1 and 38 µM, respectively. In conclusion, a reliable cocktail method including the nine most important drug-metabolizing CYP enzymes was developed, optimized and validated for detecting time-dependent inhibition. Moreover, gemfibrozil 1-O-ß-glucuronide was established as a selective inhibitor of CYP2C8 for use as a diagnostic inhibitor in in vitro studies.