An integrated pharmacokinetic/pharmacogenomic analysis of ABCB1 and SLCO1B1 polymorphisms on edoxaban exposure.

Edoxaban and its low-abundance, active metabolite M4 are substrates of P-glycoprotein (P-gp; MDR1) and organic anion transporter protein 1B1 (OATP1B1), respectively, and pharmacological inhibitors of P-gp and OATP1B1 can affect edoxaban and M4 pharmacokinetics (PK). In this integrated pharmacogenomic analysis, genotype and concentration-time data from 458 healthy volunteers in 14 completed phase 1 studies were pooled to examine the impact on edoxaban PK parameters of allelic variants of ABCB1 (rs1045642: C3435T) and SLCO1B1 (rs4149056: T521C), which encode for P-gp and OATP1B1. Although some pharmacologic inhibitors of P-gp and OATP1B1 increase edoxaban exposure, neither the ABCB1 C3435T nor the SLCO1B1 T521C polymorphism affected edoxaban PK. A slight elevation in M4 exposure was observed among SLCO1B1 C-allele carriers; however, this elevation is unlikely to be clinically significant as plasma M4 concentrations comprise <10% of total edoxaban levels.

Introduction

Edoxaban is a once-daily, orally administered, selective and direct inhibitor of factor Xa, a key mediator of both the intrinsic and extrinsic pathways of the clotting cascade.[1, 2] By binding to free Xa as well as Xa bound to the prothrombinase complex, edoxaban attenuates thrombin generation, and consequently, thrombus formation.[1, 2] Edoxaban is approved in Japan for the prevention of venous thromboembolism (VTE) following major orthopedic surgery.[3] Based on the results of 2 large clinical trials,[4, 5] edoxaban recently received approval for the prevention of stroke and systemic embolism in patients with nonvalvular atrial fibrillation and the treatment of venous thromboembolism in Japan,[6] the United States,[7] Switzerland[8] and the European Union.[9]

Edoxaban is a substrate of the efflux pump P-glycoprotein (P-gp; also known as multidrug resistance protein 1 or MDR1) and its metabolite M4 is a substrate of the hepatic uptake transporter organic anion transporter protein 1B1 (OATP1B1), as demonstrated in in vitro studies.[10] Drug–drug interaction studies have demonstrated that changes in transporter function may influence edoxaban and M4 exposure in a clinically significant manner. P-gp inhibitors, including quinidine and verapamil, increased edoxaban exposure in healthy subjects and in nonvalvular atrial fibrillation patients.[11, 12] Cyclosporine, which inhibits several transporters including OATP1B1, and rifampin increase M4 levels in healthy volunteers.[13, 14]

Genetic variations in ABCB1 (ATP-binding cassette, sub-family B, member 1) and SLCO1B1 (solute carrier organic anion transporter family, member 1B1), the genes that encode for P-gp and OATP1B1, respectively, are known to affect transporter function and thus substrate exposure. In particular, studies have indicated that a single-nucleotide polymorphism (SNP) in ABCB1, rs1045642—a C to T substitution at position 3435 of the coding sequence—alters P-gp expression and may affect the pharmacokinetics (PK) of some P-gp substrates, although reports regarding the nature of the effect are conflicting.[15, 16] Similarly, the SLCO1B1 polymorphism, rs4149056—a T to C substitution at position 521 of the coding sequence—is associated with reduced transporter activity; augments exposure of OATP1B1 substrates such as atorvastatin, rosuvastatin and pravastatin;[17, 18] and is associated with increased risk of statin-induced myalgia and rhabdomyolysis.[19]

Edoxaban undergoes minimal metabolism as nearly 73% of the absorbed drug is eliminated unmodified.[20] M4 (also known as D21-2393), which has anticoagulant properties similar to edoxaban (data on file), is a low-abundance, active metabolite of edoxaban, accounting for <10% of the total anticoagulant exposure.[20] The metabolic enzyme, carboxylesterase 1 (CES1), catalyzes the formation of M4 from edoxaban.[20] The terminal half-life of edoxaban ranges from 10 to 14 h,[21] and ~50% of the absorbed dose is renally eliminated.[20, 22] The in vitro plasma protein binding of edoxaban and M4 are ~55% and 80%, respectively.[21, 23]

Although increased edoxaban and M4 exposures are observed upon administration of the P-gp inhibitors quinidine and verapamil or the OATP1B1 inhibitor cyclosporine, the effects of common genetic variants that alter the function of these transporters on edoxaban and M4 PK are unknown. To answer these questions, a prospective-retrospective integrated analysis using pooled genotype and PK data from 14 completed edoxaban phase 1 studies was conducted.

Materials and methods

Analysis design and population selection

Twenty-four phase 1 edoxaban studies were considered for inclusion in the present integrated analysis. Completed studies were included in the current analysis provided they enrolled healthy subjects, included a treatment arm in which a single 60-mg edoxaban dose was administered under fasting conditions, and collected DNA samples that were available for genotyping. Only data from the edoxaban 60 mg treatment groups, in which edoxaban was administered alone and not with any potentially interacting drugs or food, were included in the analysis. Based on these criteria, DNA samples and PK data from 458 healthy subjects (365 males and 93 females) from 14 completed edoxaban phase 1 clinical studies were pooled for this analysis (Supplementary Table 1). All subjects provided written informed consent to participate in the clinical studies and supply genetic material. Each study was approved by an institutional review board and was carried out in compliance with the Declaration of Helsinki and the International Conference on Harmonisation.

For each of the included phase 1 studies, serial blood samples were collected at predose and up to 96 h postdose to determine edoxaban and M4 plasma concentrations in accordance with the specific study protocols.

Bioanalytic methods

Edoxaban and M4 plasma concentrations were measured by a validated liquid chromatography tandem mass spectrometry (LC–MS/MS) method for all studies by Advion BioServices (Ithaca, New York, USA). The lower limit of quantitation for edoxaban and M4 was 0.764 ng ml−1 and 0.0792 ng ml−1, respectively.

Genotype analysis

Validated commercial Taqman assays (Life Technologies, Waltham, MA, USA) were performed to determine ABCB1 rs1045642 and SLCO1B1 rs4149056 genotypes by Gentris Corporation (Morrisville, North Carolina, USA). Assays were carried out in duplicate. If results were not in agreement, the sample was genotyped using Sanger sequencing. Results were further confirmed by reanalyzing ~25 (~5%) randomly selected samples for both genotypes using Sanger sequencing in both the forward and reverse direction. Using this approach, the genotype call rate was 100%.

Pharmacokinetic and statistical analysis

An integrated database consisting of PK parameters and genotype data from each of the 458 subjects included in this analysis was created. The PK parameters AUCinf, AUC0–24, AUClast, Cmax, and C24 were computed using a noncompartmental approach with WinNonlin Professional software Version 4.0 or 5.2 (Pharsight Corporation, Montreal, Canada). PK parameters were statistically compared between genotypes using analysis of variance (ANOVA) in SAS Proc MIXED (SAS Institute, Cary, North Carolina, USA), with genotype and study ID as fixed effects. Point estimates and 90% confidence intervals (CI) for the ratios of the PK parameters between different allele variations of ABCB1 rs1045642 (CC vs CT, CC vs TT, and CT vs TT) and SLCO1B1 rs4149056 (C-carriers vs TT) were calculated using multiple pairwise comparisons by applying exponential transformation to the difference in least squares (LS) means calculated using ln-transformed values and CIs obtained from the ANOVA model.

Assuming the expected mean ratio was within the range of 95% to 105% with Bonferroni adjustment for multiple comparisons among the three rs1045642 variant subgroups, a genotype subgroup of at least 37, 105 and 98 subjects for edoxaban AUCinf, Cmax and C24, respectively, was required in order to have at least 80% power so that the 90% CI for the mean ratio would fall within the range of 80–125%. Using similar criteria, rs4149056 genotype subgroups required a minimum of 26, 75 and 70 subjects for edoxaban AUCinf, Cmax and C24, respectively, while the corresponding PK measures for M4 required 57, 103 and 96 subjects, respectively. PK parameters were not considered to be different between genotype groups if the 90% CIs for the geometric LS means ratios of the analysis groups were within the range of 80–125%.

The metabolite-to-parent ratios (MPRs) were calculated for AUClast, AUCinf and Cmax with adjustments for molecular weights and summarized by SLCO1B1 rs4149056 genotype (TT and C-carriers). The molecular weights used in the MPR calculations were 548.07 g mol−1 for edoxaban (parent) and 520.99 g mol−1 for the metabolite M4.

Results

Demographics and baseline characteristics

The demographics and baseline characteristics of subjects included in this analysis are listed by ABCB1 and SLCO1B1 genotype in Table 1. The mean age, body mass index, and renal function were similar across genotype groups. The ABCB1 and SLCO1B1 SNPs were in Hardy–Weinberg equilibrium for all race groups examined. The observed frequency of ABCB1 genotypes varied across ethnic groups, consistent with previously published findings.[24]

Table 1

Demographic summary of ABCB1 (rs1045642) and SLCO1B1 (rs4149056) genotyped population

	ABCB1 (rs1045642)			N=458	SLCO1B1 (rs4149056)
	CC (n=206)	CT (n=193)	TT (n=59)		TT (n=384)	CT (n=71)	CC (n=3)
Sex, n (%)
Male	167 (81.1)	150 (77.7)	48 (81.4)	365 (79.7)	311 (81.0)	53 (74.6)	1 (33.3)
Female	39 (18.9)	43 (22.3)	11 (18.6)	93 (20.3)	73 (19.0)	18 (25.4)	2 (66.7)
Race, n (%)
American Indian/Alaskan	0 (0.0)	2 (1.0)	0 (0.0)	2 (0.4)	1 (0.3)	1 (1.4)	0 (0.0)
Asian	1 (0.5)	6 (3.1)	3 (5.1)	10 (2.2)	9 (2.3)	1 (1.4)	0 (0.0)
Black or African American	150 (72.8)	91 (47.2)	9 (15.3)	250 (54.6)	226 (58.9)	24 (33.8)	0 (0.0)
White	43 (20.9)	89 (46.1)	45 (76.3)	177 (38.6)	134 (34.9)	40 (56.3)	3 (100.0)
Other	12 (5.8)	5 (2.6)	2 (3.4)	19 (4.1)	14 (3.6)	5 (7.0)	0 (0.0)
Ethnicity, n (%)
Hispanic/Latino	38 (18.4)	51 (26.4)	16 (27.1)	105 (22.9)	84 (21.9)	20 (28.2)	1 (33.3)
Not Hispanic/Latino	168 (81.6)	142 (73.6)	43 (72.9)	353 (77.1)	300 (78.1)	51 (71.8)	2 (66.7)
Age (years)
Mean±s.d.	32.1±7.4	30.2±6.9	30.9±8.5	31.1±7.4	31.4±7.4	29.6±6.9	29.3±12.1
Weight (kg)
Mean±s.d.	81.5±12.1	79.0±12.1	76.6±12.0	79.8±12.2	80.1±12.1	78.8±12.3	64.0±4.0
BMI (kg m−2)
Mean±s.d.	26.6±3.1	25.9±3.0	25.2±3.0	26.1±3.1	26.1±3.0	26.2±3.2	24.9±3.9
CLcr (ml min−1)
Mean±s.d.	134.0±26.7	137.0±25.2	138.0±25.3	136.0±25.9	134.5±26.3	141.0±22.7	146.7±25.2

Abbreviations: BMI, body mass index; CLcr, creatinine clearance; SD, standard deviation.

Due to the infrequent occurrence of SLCO1B1 CC homozygotes (0.66%) in the overall population, subjects with the CT or CC genotypes were combined together and analyzed as C-carriers. Thus, overall, 16.2% of subjects were C-carriers, and 83.8% of subjects were TT homozygotes. As previously reported in the literature, the frequency of SLCO1B1 T521C variants differed between races.[25]

Pharmacokinetic parameters

As shown in Figure 1a, ABCB1 genotype did not significantly alter mean plasma edoxaban concentration–time curves. Edoxaban PK parameters including total plasma exposure (AUCinf), peak (Cmax) and trough plasma concentrations (C24) were not affected by ABCB1 C3435T genotype (Figures 1b–d and Table 2). In addition, plasma exposure from the time of dosing to the last measurable concentration (AUClast) and plasma exposure up to 24 h after dosing (AUC0–24) exhibited a similar range of values across genotypes (Table 2). The values of ratios of geometric LS means for AUCinf, Cmax and C24 were 103.5 (90% CI, 98.2–109.1), 109.8 (90% CI, 100.3–120.2), and 93.4 (90% CI, 85.7–101.9), respectively, for a CT vs CC comparison, and 104.9 (90% CI, 96.8–113.6), 106.2 (90% CI, 92.6–121.8), and 92.3 (90% CI, 80.8–105.4), respectively, for a TT vs CC comparison. These data reveal that the ABCB1 C3435T SNP does not significantly affect edoxaban exposure.

Figure 1

The effect of the ABCB1 rs1045642 variant on edoxaban pharmacokinetics. (a) Mean concentration–time profiles by genotype. Box and whisker plots demonstrating the effect of ABCB1 rs1045642 genotypes on (b) AUCinf, (c) Cmax and (d) C24. Inset figure in (a) shows profile of interval from initial dosing to 12 h postdose. Error bars in (a) represent the standard deviation. Whiskers in (b, c and d) indicate the 10th and 90th percentiles; individual dots represent points outside of the 10th and 90th percentiles. Median values are represented numerically and as lines in box plots.

Table 2

Summary of pharmacokinetic parameters by ABCB1 rs1045642 genotype

Parameter	ABCB1 (rs1045642) genotypes
	CC (n=206)	n	CT (n=193)	n	TT (n=59)	n
AUClast (ng•h ml−1)	1760±447.2	206	1814±417.0	193	1820±434.8	59
AUC0–24 (ng•h ml−1)	1608±441.2	206	1671±420.3	193	1680±441.5	59
AUCinf (ng•h ml−1)	1789±450.2	205	1845±417.8	190	1863±436.5	58
Cmax (ng ml−1)	246±96.6	206	268±95.2	193	261±101.3	59
C24 (ng ml−1)	12.5±5.11	205	11.8±5.07	193	12.3±6.46	58

Abbreviations: AUC, area under the curve; AUC0–24, AUC from the time of dosing over a 0–24 h postdose time interval; C24, observed plasma concentration at 24 h postdose; AUCinf, AUC from the time of dosing to the last measurable concentration and extrapolated to infinity; AUClast, area under the plasma drug concentration–time curve from the time of dosing to the last measurable concentration; Cmax, maximum observed plasma concentration.

Data presented as mean±s.d.

The mean plasma concentration vs time profiles of edoxaban were also comparable among subjects independent of SLCO1B1 genotype (Figure 2a, left panel). Consistent with this observation, AUCinf, Cmax, and C24 of edoxaban were similar between C-carriers and TT homozygotes (Figure 2b and Table 3). Furthermore, the range of values for AUClast and AUC0–24 did not differ between C-carriers and TT homozygotes (Table 3). In comparing C-carriers vs TT homozygotes, the geometric LS means ratios for edoxaban AUCinf, Cmax and C24 were 101.0 (90% CI, 95.9–106.3), 105.4 (90% CI, 96.5–115.1) and 95.3 (90% CI, 87.5–103.7), respectively.

Figure 2

The effect of the SLCO1B1 rs4149056 variant on edoxaban and M4 pharmacokinetics. Mean concentration–time profiles by genotype of (a) edoxaban (left panel) and M4 (right panel). Box and whisker plots demonstrating the effect of SLCO1B1 rs4149056 genotypes on AUCinf, Cmax and C24 for (b) edoxaban and (c) M4. Inset figure in (a) shows profile of interval from initial dosing to 12 h postdose. Error bars in (a) represent the standard deviation. Whiskers in (b and c) indicate the 10th and 90th percentiles; individual dots represent points outside of the 10th and 90th percentiles. Median values are represented numerically and as lines in box plots.

Table 3

Summary of pharmacokinetic parameters of edoxaban and M4 by SLCO1B1 rs4149056 genotype

SLCO1B1 (rs4149056) genotypes
PK parameter	TT	n	C-carriers	n
Edoxaban
AUClast (ng • h ml−1)	1792.9±442.8	384	1777.2±380.8	74
AUC0–24 (ng • h ml−1)	1644.1±440.6	384	1643.8±392.9	74
AUCinf (ng • h ml−1)	1823.5±444.9	381	1815.5±382.4	72
Cmax (ng ml−1)	255.1±97.4	384	268.0±94.7	74
C24 (ng ml−1)	12.3±5.4	382	11.6±4.8	74
M4
AUClast (ng • h ml−1)	147.1±60.7	384	184.6±75.6	74
AUC0–24 (ng • h ml−1)	134.7±57.2	384	170.0±73.0	74
AUCinf (ng • h ml−1)	150.2±61.4	379	191.9±76.4	72
Cmax (ng ml−1)	21.2±10.7	384	26.7±13.6	74
C24 (ng ml−1)	1.1±0.57	382	1.3±0.74	74
MPR AUClast	8.7±2.9	384	10.9±3.6	74
MPR AUCinf	8.7±2.9	378	11.1±3.5	72
MPR Cmax	8.8±3.1	384	10.5±3.3	74

Abbreviations: AUC0–24, AUC from the time of dosing over a 0–24 h postdose time interval; C24, observed plasma concentration at 24 h postdose; AUCinf, AUC from the time of dosing to the last measurable concentration and extrapolated to infinity; AUClast, area under the plasma drug concentration–time curve from the time of dosing to the last measurable concentration; Cmax, maximum observed plasma concentration; MPR, metabolite to parent ratio of relevant PK parameter.

Data shown as mean±s.d.

Plasma M4 exposure was slightly higher among SLCO1B1 C-carriers compared to TT homozygotes, as assessed by the mean plasma concentration vs time curves of M4 (Figure 2a, right panel). Accordingly, C-carriers exhibited modestly higher levels of median total exposure and peak plasma concentration (Figure 2c). Median M4 trough levels were comparable between the 2 genotypes (Figure 2c), while mean AUClast and AUC0–24 values were greater among C-carriers than TT homozygotes (Table 3). The geometric LS means ratios comparing C-carriers with TT homozygotes for AUCinf, Cmax and C24 were 126.9 (90% CI, 117.4–137.1), 124.1 (90% CI, 111.9–137.5) and 114.9 (90% CI, 104.0–127.0), respectively; the 90% CI of the ratios of geometric LS means of C-carriers vs TT homozygotes exceeded the predefined bioequivalence range (90% CI, 80–125%) for each of these PK measures. Taken together, these data confirm that while the SLCO1B1 T521C SNP does not influence edoxaban exposure, M4 exposure is elevated among C-carriers.

Edoxaban and M4 have similar anticoagulant activities on a molar basis (data on file). MPRs for AUClast, AUCinf and Cmax of M4 to edoxaban ranged from 9% to 11%, indicating that M4 contributes a small amount to the overall anticoagulant effect (Table 3 and Figure 3). Together, these findings demonstrate that the SLCO1B1 T521C polymorphism increases M4 exposure, but not in a clinically meaningful way, as M4 plasma concentration is approximately one-tenth that of edoxaban.

Figure 3

Summary pharmacokinetics of edoxaban and M4 in plasma by SLCO1B1 rs4149056 genotype. Bar graphs represent mean values, and the error bars indicate standard deviation.

Safety

The single 60-mg dose of edoxaban was generally well tolerated by the healthy subjects in these 14 phase 1 studies.[13, 14, 26, 27, 28, 29]

Discussion

Common genetic variants and pharmacologic inhibition can result in altered transporter function, and thereby affect the PK of certain substrates. Allelic variants of ABCB1, including C3435T, which codes for a synonymous substitution in exon 26, are reported to alter the disposition of some P-gp substrates.[15, 24, 30, 31] For example, some reports indicate that exposure to digoxin, a P-gp substrate, is altered upon the administration of the P-gp inhibitor quinidine[32] and is affected by the ABCB1 C3435T SNP.[15, 30, 31] The SLCO1B1 T521C SNP results in a valine-to-alanine substitution in the uptake transporter OATP1B1, leading to decreased membrane expression[33] and transport activity of OATP1B1.[19, 34, 35] Concentrations of the OATP1B1 substrate rosuvastatin increase upon treatment with the OATP1B1 inhibitor cyclosporine[36] and are also higher among SLCO1B1 T521C CC homozygotes than TT homozygotes.[18]

Edoxaban is a substrate of the P-gp, and exposure to edoxaban increases upon coadministration of the P-gp inhibitors ketoconazole, quinidine, erythromycin, verapamil or dronedarone in both healthy subjects and nonvalvular atrial fibrillation patients.[11, 12] M4, an active metabolite of edoxaban, is a substrate of OATP1B1,[10] and coadministration of edoxaban with the OATP1B1 inhibitor cyclosporine increases exposure to M4.[13]

The present integrated analysis pooled PK and genotype data from 14 completed single-dose, phase 1 studies involving 458 healthy subjects to examine the impact of polymorphisms in two important transporter genes on edoxaban exposure. Results from this analysis demonstrate for the first time that common genetic variants of ABCB1 (C3435T) and SLCO1B1 (T521C) had minimal impact on edoxaban exposure. Edoxaban concentration–time profiles were comparable regardless of ABCB1 C3435T genotype. In addition, mean and median AUCinf, Cmax and C24 of edoxaban did not vary between the three ABCB1 genotypes. With regard to the SLCO1B1 T521C genotype, edoxaban concentration–time profiles and median PK parameters including AUCinf, Cmax and C24 were similar between TT homozygotes and C-carriers. This was expected as edoxaban is not known to be an OATP1B1 substrate. While a modest increase in M4 plasma exposure is observed among C-allele carriers, this association is unlikely to be clinically relevant given that M4, at concentrations reaching <10% of edoxaban levels, does not significantly contribute to the total anticoagulant activity. Therefore, while some pharmacologic inhibitors of P-gp and OATP1B1 may increase edoxaban and M4 exposure in a clinically relevant way,[11, 12, 13] these two common polymorphisms do not.

As concurrent administration of P-gp or OATP1B1 inhibitors with edoxaban increases edoxaban or M4 exposure, respectively, the effects of functional polymorphisms of ABCB1 and SLCO1B1 on edoxaban exposure were tested in this study. Neither ABCB1 C3435T nor SLCO1B1 T521C genotype status affected the exposure of edoxaban among healthy subjects, suggesting that these are not clinically significant polymorphisms. As such, the data reported here suggest that prescribers can follow approved dosing recommendations[21] for edoxaban without concern for the presence of these two genetic variants.

This study was limited to candidate SNPs and thus did not capture the effects of other known ABCB1 and SLCO1B1 polymorphisms or of haplotypes[24, 25] on edoxaban and M4 PK. The polymorphisms chosen for this analysis are common, well characterized, and known to be functionally active[24, 25] and thus were selected for initial pharmacogenomic signal detection. Because of the lack of an observed effect with these polymorphisms, additional SNPs and haplotypes within ABCB1 and SLCO1B1 were not characterized, as they are unlikely to have a substantial impact on edoxaban or M4 PK.

In pharmacogenomics, integrated multistudy analyses can be especially useful to answer important questions that otherwise cannot be resolved within small individual studies. Some phase 1 studies—which characterize the PK parameters of new drugs and investigate the effects of food and potential interacting drugs on those PK parameters—include pharmacogenomics sample collection, although their typically small sample size limits the investigation of pharmacogenomic effects. Although phases 2 and 3 studies are often sufficiently powered for pharmacogenomics analyses, these studies usually lack the comprehensive concentration–time data needed to explore potential pharmacogenomic effects on PK parameters. The current integrated analysis circumvents these limitations by pooling genetic and concentration–time data from multiple phase 1 studies. The resulting data set, containing genotypes, covariates, and full PK profiles, is large enough to conduct an adequately powered pharmacogenomic analysis of PK parameters. This approach has been successfully employed in previous studies to investigate the effects of CYP2C19 polymorphisms on the outcomes of clopidogrel treatment[37] and UGT1A1 polymorphisms on dolutegravir PK.[38]

Variability can be introduced into pooled analyses when different laboratories are used to collect genotype or concentration–time data. For the present study, all genetic data were generated by a single genotyping laboratory and all concentration–time data were generated by a single bioanalysis laboratory. However, some limitations remain. First, SNPs in CES1, which is involved in the metabolism of edoxaban to M4, have previously been reported to affect the disposition of other drugs[39, 40] but were not investigated in this study. In addition, while this analysis controlled for race, the effect of ABCB1 C3435T and SLCO1B1 T521C on edoxaban and M4 PK were not analyzed by race. The present study was adequately powered; however, its sample size limits further evaluation of the effect of race, additional polymorphisms, and haplotypes. Larger follow-up replication studies will be required to address the limitations of this analysis.