- Research
- Open access
- Published:
Low heritability in pharmacokinetics of talinolol: a pharmacogenetic twin study on the heritability of the pharmacokinetics of talinolol, a putative probe drug of MDR1 and other membrane transporters
Genome Medicine volume 8, Article number: 119 (2016)
Abstract
Background
Efflux transporters like MDR1 and MRP2 may modulate the pharmacokinetics of about 50 % of all drugs. It is currently unknown how much of the variation in the activities of important drug membrane transporters like MDR1 or MRP2 is determined by genetic or by environmental factors. In this study we assessed the heritability of the pharmacokinetics of talinolol as a putative probe drug for MDR1 and possibly other membrane transporters.
Methods
Talinolol pharmacokinetics were investigated in a repeated dose study in 42 monozygotic and 13 same-sex dizygotic twin pairs. The oral clearance of talinolol was predefined as the primary parameter. Heritability was analyzed by structural equation modeling and by within- and between-subject variance and talinolol clearance was correlated with polymorphisms in MDR1, MRP2, BCRP, MDR5, OATP1B1, and OCT1.
Results
Talinolol clearance varied approximately ninefold in the studied sample of healthy volunteers. The correlation of clearances between siblings was not significantly different for the monozygotic and dizygotic pairs. All data analyses consistently showed that variation of talinolol pharmacokinetics was mainly determined by environmental effects. Structural equation modeling attributed 53.5 % of the variation of oral clearance to common environmental effects influencing both siblings to the same extent and 46.5 % to unique environmental effects randomly affecting individual subjects. Talinolol pharmacokinetics were significantly dependent on sex, body mass index, total protein consumption, and vegetable consumption.
Conclusions
The twin study revealed that environmental factors explained much more of the variation in pharmacokinetics of talinolol than genetic factors.
Trial registration
European clinical trials database number: EUDRA-CT 2008-006223-31. Registered 26 September 2008. ClinicalTrials.gov number: NCT01845194.
Background
Efflux transmembrane transport of drugs and other molecules is highly variable within and between humans due to inherited and environmental factors. This can have profound consequences for the intended and adverse effects of drugs and for the biological effects of many other endogenous and exogenous molecules. Efflux transporters first came into focus as mediators of multi-drug resistance in cancer chemotherapy [1, 2]. Among the efflux transporters, multi-drug resistance protein 1 (MDR1, P-glycoprotein), encoded by the ABCB1 gene, has been the most extensively studied [3–9]. MDR1 is expressed on the apical membrane of enterocytes, hepatocytes, and renal tubular cells and at the blood–brain barrier and the blood–testes barrier and protects cells and organs against overload with numerous substances [10, 11]. Some animals completely lack MDR1 but we do not know of any humans completely lacking it [10]. Nevertheless, MDR1 expression and activity vary substantially between individuals. MDR1 expression is regulated by several mechanisms and substantial up-regulation can be induced by several substances and is mediated by PXR, RXR, and other nuclear transcription factors [12]. On the other hand, numerous environmental substances and drugs, such as cyclosporine or verapamil, can substantially inhibit P-glycoprotein [13, 14].
Inherited polymorphisms may contribute to variation in MDR1 activity. In particular, the polymorphisms C1236T (silent), G2677T/A (Ala893Ser/Thr), and C3435T (silent) have been extensively studied [15–24]. Understanding of the functional role of these variants in humans is complicated by the fact that these variants are in strong linkage disequilibrium [25, 26]. Genetic polymorphisms may play some role in the variation of MDR1 expression and activity. However, the genotype–phenotype correlation data obtained by different investigators is not entirely consistent. In conclusion, it is currently unknown how much of the variation in MDR1 activity as tested in vivo with putative probe drugs like talinolol or digoxin is explained by environmental or heritable factors.
Another efflux transporter of particular interest for the disposition of drugs, other foreign compounds, and endogenous substrates is the multidrug resistance-associated protein 2 (MRP2; also known as ABCC2 (ATP binding cassette subfamily C member 2)). MRP2 is expressed in a broad range of tissues and also contributes to multi-drug resistance in cancer therapy [27]. It is expressed on the apical membrane of cells in the liver, intestine, and kidneys [28, 29]. MRP2 is an important efflux transporter of conjugated bilirubin and inherited lack of activity is known as Dubin–Johnson syndrome [30]. Furthermore, MRP2-mediated drug interactions may cause liver toxicity [31]. Activity of MRP2 can also be greatly modified by transcriptional regulation and by inhibition of its transport function, but also by genetic polymorphisms. Among others, the impact of two genetic polymorphisms (-C24T and G1249A in MRP2) on, for example, the bioavailability of mycophenolate mofetil [32, 33] and tacrolimus [34, 35], have been extensively studied, with controversial results.
A number of probe drugs for the in vivo study of MDR1 activity have been suggested but none of them is perfect [36, 37]. Because of the low substrate specificity of many drug membrane transporters, many probe drugs may reflect the activity of several transporters. For instance, a typical probe drug for MDR1, fexofenadine, is also a substrate of several OATP influx transporters [37]. In this study, talinolol was chosen as a probe drug for drug membrane transporters. Talinolol appears suitable as a membrane transporter probe drug since less than 1 % of an orally administered dose was recovered as metabolites in urine [38, 39]; thus, metabolism plays only a minor role in its elimination.
Evidence for which transporters are possibly relevant in the pharmacokinetics of talinolol is based on multiple data sources. Talinolol is widely considered as an MDR1 probe drug both in vitro and in vivo [36]. Several investigators used talinolol in caco-2 cells as a model substrate of MDR1. Plasma concentrations after oral administration of talinolol were substantially (4.5-fold) higher in a MDR1 knockout rat compared with the wild-type rat [40]. Clinically, net intestinal secretion of talinolol was inhibited by the MDR1 inhibitor verapamil [41] and carbamazepin [42]. Additionally, rifampin-induced MDR1 and MRP2 expression resulted in decreased oral bioavailability of talinolol [43] and corresponding effects were seen with St. John’s wort [44]. Unexpectedly, two inhibitors of MDR1, grapefruit juice and verapamil, decreased oral bioavailability of talinolol, which was explained by preferential inhibition of intestinal influx transporters by polytropic inhibitors like verapamil [45, 46]. According to these data, it is very likely that MDR1 modulates the pharmacokinetics of talinolol in humans to a relevant extent.
MRP2 may also modulate the pharmacokinetics of talinolol. Expression of both MDR1 and MPR2 is greatly increased by typical inducers like carbamazepin and rifampin and the effects of these inducers on talinolol pharmacokinetics may thus also be explained by an involvement of MRP2. This is also supported by the low talinolol oral bioavailability found in carriers of the MPR2 exon 10 Val417Ile genotype (G1249A) [47].
A comprehensive in vitro screen of all influx and efflux transporters possibly relevant in talinolol pharmacokinetics has not been performed to our knowledge. Preliminary data indicated that talinolol may be a substrate of OATP1B1 [48] and with its cationic structure, talinolol may be a substrate of OCT1. In addition, other efflux transporters, including the polymorphic MDR5 (ABCB5) and BCRP (ABCG2) [49], may modulate the pharmacokinetics of talinolol; therefore, we also explore here the impact of polymorphisms known to be functional in these transporters.
The aim of the present study was to reveal how much of the variation in talinolol pharmacokinetics is due to heritable factors. Secondary questions were whether we can reproduce the effects of genetic polymorphisms which were associated with talinolol pharmacokinetics in earlier publications and how suitable a much lower test dose of talinolol (2.5 mg) might be for in vivo phenotyping. A low dose is apparently more safe and the lower dose was also scientifically interesting because the impact of MDR1 on talinolol pharmacokinetics might depend on the dose [50].
Methods
Subjects and study design
Monozygotic (MZ) and same-sex dizygotic (DZ) twin pairs were included if they were healthy according to medical history, clinical examination, electrocardiogram, and standard clinical biochemistry and hematology tests. Except for oral contraceptives no other drugs were allowed for 1 week prior to and 48 h after administration of talinolol.
A single dose of 50 mg talinolol (Cordanum®, AWD.pharma, Dresden, Germany) was administered to 114 subjects. Among these subjects were 43 MZ and 13 DZ twin pairs; two single twins were excluded from the analyses since the second twin could not participate in the study. Additionally, one MZ pair was excluded from the analyses because of missing concentration data for the last blood sample 22 h after talinolol administration. Therefore, all analyses on 50 mg talinolol are based on data from 110 subjects (42 MZ and 13 DZ twin pairs). Per protocol, talinolol 50 mg was administered on three different occasions to each subject with an at least 7-day wash-out period in between. Blood sampling was done prior to and 0.5, 1, 2, 3, 4, 5, 6, and 22 h after talinolol administration. Talinolol was part of a phenotyping drug cocktail in which midazolam, torsemide, metoprolol (all intravenously) [51], and caffeine (orally) [52] were also administered in a sequential order in the first 2 h of each study day. All the substances in this cocktail are not known to have any pharmacokinetic interactions and talinolol was the last administered drug and was given 1 h after caffeine.
In an additional add-on study aiming to compare low-dose talinolol, which is presumed to lack relevant cardiovascular effects, 2.5 mg talinolol was administered to 73 subjects: 29 MZ and 7 DZ twin pairs. One single twin was excluded from the analyses. The capsules containing 2.5 mg of talinolol were prepared from Cordanum® 50 tablets, which were ground and used to refill gelatin capsules in 2.5 mg portions. In this add-on study, blood sampling was done prior to and 0.5, 1, 1.5, 2, 3, 4, 6, and 8 h after talinolol administration.
Blood samples were centrifuged for 10 minutes (4000 rpm, 4 °C) and plasma samples were stored at −20 °C until laboratory analysis.
Talinolol plasma concentrations
Talinolol plasma concentrations were determined by liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis using deuterium labeled talinolol-d5 (Toronto Research Chemicals, Nr. T005002) as internal standard. For the 50 mg talinolol study sample preparation with protein precipitation and HPLC separation were performed as described [51]. Samples from the low-dose study were worked up with solid phase extraction. LC-MS/MS analysis was carried out on an Agilent 6460 triple quadrupol mass spectrometer (Agilent, Waldbronn, Germany) coupled to an Agilent 1200 HPLC system. Ionization mode was electrospray (ESI), polarity positive. The mass spectrometer was operated in the multiple reaction monitoring (MRM) mode using the transitions m/z 364.3 to 308 and m/z 369.3 to 313 for talinolol and the internal standard, respectively. For the 50-mg dose study, the limit of quantification (LOQ) was 1.5 pmol/ml. Quality control samples over the whole concentration range could be determined with a coefficient of variation (CV) better than 2.1 %. For the low-dose study, the LOQ was 0.2 pmol/ml. Over the whole concentration range (from 0.2 to 30 pmol/ml), a CV better than 5 % was achieved.
Molecular genetic analyses
Genotyping was performed by single-base primer extension using fluorescence-labeled dideoxynucleotides followed by capillary electrophoresis. Genotype calling was performed using Gene Mapper v3.7 Software ® (Applied Biosystems ®, Foster City, CA, USA). The following variants of efflux transporters were genotyped: MDR1 C1236T (rs1128503), G2677T/A (rs2032582), and C3435T (rs1045642); MDR5 C392T (rs17143212); MRP2 G1249A (rs2273697) and -C24T (rs717620), BCRP (ABCG2) G34A (rs2231137) and C421A (rs2231142). Additionally, the following haplotypes of the influx transporters were determined: OATP1B1 *1a, *1b, *5, *14, *15, and OCT1 alleles *1 to *10. In total, variants in 23 polymorphic genes were genotyped for assessment of zygosity. Monozygosity was concluded when both siblings of each pair had completely identical genotypes for all 23 polymorphic loci.
Statistics
Pharmacokinetic parameters were estimated with WinNonlin (Pharsight Corporation, Mountain View, USA). Total clearance, the primary parameter, normalized to bioavailability (Cl/F) was calculated as dose/AUCinfinity and in this study is referred to as “clearance”. The area under the concentration time curve (AUC) was calculated by the linear/log trapezoidal rule. AUCinfinity values were extrapolated based on the last predicted concentration using the terminal elimination rate constant (lambda z). AUC7h reflects the AUC in the first 7 h after dosage. Maximum concentration (Cmax) and time of maximum concentration (tmax) are given as measured. The terminal half-life (t1/2) was calculated as ln(2)/lambda z. The volume of distribution based on the terminal phase (Vz) was calculated as Dose/(lambda z × AUCinfinity).
Heritability was estimated by structural equation modeling using the “mets” package for the statistical programming environment R [53, 54]. Variation in clearance was decomposed into additive genetic effects (A), dominant genetic effects (D), common environmental effects (C), and unique environmental effects (E). Fit of the models ACE, ADE, AE, and CE was tested by the Akaike-information criterion (AIC). In addition, the repeated administration of 50 mg talinolol to each subject was used to calculate the genetic component (rGC) as described by Kalow et al. [55, 56]. Variance within (Vw) the subjects and variance between (Vb) the subjects were calculated for all subjects who received 50 mg talinolol on three different occasions. The mean of 50 calculations of rGC is given. For each calculation of rGC only one randomly selected sibling of each pair was used, since the two siblings of each twin pair are not independent units.
The non-parametric Mann–Whitney U test was used for comparison of two groups and the Jonckheere–Terpstra test for comparison of more than two groups. In the analyses of the effects of genetic polymorphisms on talinolol clearance, an adjustment for multiple testing was done by using the Bonferroni correction method. Due to the high number of tests on genetic polymorphisms (ten in total), significance was determined if p < 0.005. Multiple linear regression analysis was performed for evaluating the effects of the genetic polymorphisms as well as the confounders sex, age, body mass index, consumption of protein, and consumption of vegetables. Protein and vegetable consumption was recorded using validated nutrition questionnaires [57] and calculated as the mean of items “meat”, “sausage”, “fish”, “dairy products”, and “eggs” and as the mean of items “salad and raw vegetable” and “cooked vegetable”, respectively.
Results
Heritability of the variation in pharmacokinetics of talinolol was studied in 110 volunteers (42 monozygotic (MZ) and 13 dizygotic (DZ) twin pairs). No serious adverse events related to the study medication were reported. No statistically significant differences in the demographic data between MZ and DZ twin pairs were observed (Additional file 1).
Talinolol total clearance (Cl/F) measured after 50 mg talinolol varied 3.8-fold (in DZ) and 8.8-fold (in MZ); after adjustment for body weight, variation in clearance remained high with 2.9-fold and 6.3-fold in DZ and MZ twins, respectively. Similarly, strong variation was observed in AUCinfinity and AUC7h (Table 1).
No sequence effects were observed, i.e., the mean pharmacokinetic parameters were not significantly different between the study days. Clearance varied within subjects up to threefold between the three different study days (Fig. 1a) and the between-subject coefficient of variation (CV) was about 42 % for the whole study population. Within-twin pair correlations of the clearance were not significantly different between DZ and MZ twin pairs. With Pearson’s correlation coefficients of 0.65 (95 % confidence interval (CI) 0.43–0.80) for MZ and 0.82 (95 % CI 0.49–0.94) for DZ correlation, the correlation was unexpectedly nominally even higher in DZ than in MZ twin pairs (Fig. 1b, c). This indicates that common environmental effects substantially influenced the variation of the clearance.
According to structural equation modeling, common environmental (C) and unique environmental effects (E) best explained the variation in the clearance of talinolol (Table 2). Common environmental effects seem to account for 53.5 % and unique environmental effects for a further 46.5 % of the variation in the clearance (Table 2).
An alternative model (AE) explained 50 % of the variation by additive genetic effects and 50 % by environmental effects. However, statistically this model was clearly inferior according to the AIC.
The pharmacokinetics of drugs usually depends on body weight and body weight is partially also genetically determined. Results were similar after adjustment of talinolol clearance for body weight with slightly more unique environmental effects (C = 45 % (95 % CI 24–66 %), E = 55 % (95 % CI 34–76 %)).
An alternative indicator for heritability can be estimated from the comparison of within- and between-subject variation. Based on the within- and between-subject variance with the three times repetition of 50 mg talinolol, the genetic component (rGC) was 0.38 (95 % CI 0.35–0.41) for talinolol total oral clearance and 0.53 (95 % CI 0.50–0.57) for talinolol AUC (Table 3). Within- and between-subject variance was similar for tmax and Cmax, which indicates there were no relevant genetic influences on these two parameters (Table 3).
Although our study indicated only low to moderate heritability in talinolol pharmacokinetic variation, there may be up to 38 % heritability in the clearance according to the repeated measurement approach. Therefore, we evaluated the influence of genetic polymorphisms on the efflux transporters MDR1, MDR5, MRP2, and BCRP as well as the most relevant polymorphisms in the influx transporters OCT1 and OATP1B1 on talinolol clearance. Considering the moderate sample size of our study, only those polymorphisms which have been associated with transport activity in earlier studies were analyzed. Heterozygous or homozygous carriers of MDR1 1236 T alleles and MDR1 2677 T alleles showed a significant increase (p = 0.0025 and 0.0015, respectively) in the total oral clearance of talinolol (Table 4). Clearance was increased by 44 and 48 % in homozygous carriers of MDR1 1236 T and 2677 T alleles compared to homozygous carriers of 1236C or 2677G alleles, respectively. All other polymorphisms in the genes coding for the efflux transporters MDR1, MRP2, MDR5, and BCRP were not significant predictors of talinolol total oral clearance (Table 4). Concerning the influx transporters OATP1B1 and OCT1, after adjustment for multiple testing, none of the variants known to be functional showed any significant association with talinolol clearance (Additional file 2).
Multiple linear regression analysis was performed to evaluate the possible influences of demographic and environmental factors potentially modulating the in vivo activity of drug membrane transporters (Table 5). In this multifactorial analysis none of the genetic polymorphisms remained significant, but sex, body mass index, protein, and vegetable consumption significantly explained variation in the talinolol clearance; these factors each explained only about 4 % of the total oral clearance of talinolol.
In addition to the 50 mg study, a low-dose study using 2.5 mg talinolol was performed. This study included 29 MZ and 7 DZ twin pairs. To enhance the feasibility of phenotyping, we restricted blood sampling to a shorter interval and therefore AUC7h was the primary parameter. The median of the AUC7h was 0.58 mg*min/l and varied between 0.36 and 1.42 mg*min/l. Between-twin pairs Pearson’s correlation coefficients of the AUC7h were similar for both the 50-mg and the 2.5-mg groups (for 50 mg MZ = 0.73, DZ = 0.71; for 2.5 mg MZ = 0.71, DZ = 0.74; Fig. 1d, e). Structural equation modeling revealed that variation of the AUC7h was best explained by common environmental effects (68.5 %, 95 % CI 51.1–85.8 %) and unique environmental effects (31.5 %, 95 % CI 14.2–48.9 %).
While significant inter-individual variation in the time of maximum blood concentration was observed after administration of 50 mg talinolol (Fig. 2a), variation in the time–concentration curves was apparently much lower in the 2.5 mg group (Fig. 2b). After 50 mg talinolol, the time of maximum plasma concentration (tmax) varied between 1 and 4.1 h whereas after 2.5 mg, tmax was 4 h after administration in all 72 subjects (Fig. 2a, b). The maximum plasma concentration (Cmax) after 2.5 mg talinolol administration was 3.94 μg/l (median) and varied between 1.67 and 13.2 μg/l. Accordingly, the dose adjusted Cmax was 1.57 μg/l/mg in the 2.5-mg study phase and 1.85 μg/l/mg in the 50-mg study phase.
The AUC7h for the 50-mg and the 2.5-mg study phases showed a significant (p = 0.035) but only weak correlation, with a Person’s correlation coefficient of 0.25 (95 % CI 0.18–0.45) (Fig. 2c). Thus, according to the coefficient of determination, only about 6 % of the variation observed after administration of either dose was predictive for the variation in the other dose.
Discussion
The impact of genetic variation on the activities of MDR1 and MRP2 has been extensively studied in the past 20 years but is still controversial [10, 11]. Conflicting results have been reported, particularly for the abundant and widely studied variants in MDR1 (C1236T, G2677T/A, and C3435T) as well as for the -C24T and G1249A variants in MRP2 [10]. Also, the genotype–phenotype correlations for many other frequent variants have mostly not been consistently reproduced, whereas those variants with strong penetrance and clear genotype–phenotype correlations are relatively rare. In this situation, we went back in the logical sequence of research to the important first question, namely, how much of the variation in drug membrane transport as tested with in vivo probe drugs is caused by heritable factors? One approach to answer this question and to differentiate between genetic and environmental factors is the use of twin studies. The results of twin studies may provide a unique basis to direct further research [58]. If a major part of the variation in the pharmacokinetics of probe drugs can be explained by genetic factors, then further approaches like whole-genome sequencing to search for the underlying factors may be promising. If results from twin studies reveal no or only moderate genetic effects, then such genome-wide screening might not be very promising or could require major effort to achieve adequate statistical power.
We chose talinolol as a probe drug because it is not metabolized to a relevant extent and in vitro and in vivo studies revealed that it may be a good probe drug for MDR1 as well as probably for MRP2 [37]. Generally, a probe drug has great heuristic potential allowing extrapolation from effects in one drug to numerous other drugs. However, this probe drug approach has limitations because in reality pharmacokinetics of every probe drug is not dependent on just one transporter or enzyme but probably on several, and many transporters or enzymes have multiple substrate binding sites with differential substrate affinities. With talinolol, we do not know quantitatively how much of the AUC or total clearance depends on MDR1 since data obtained in the MDR1 knockout rat [40] may be compromised by unspecific or secondary effects of the knockout; also, other factors limit extrapolation from rat to human.
More than usual with pharmacokinetic data, the mean pharmacokinetic parameters of talinolol differed between the studies. Its mean oral clearance was 0.93 l/min in our study and varied between 0.58 l/min [59] and 1.95 l/min (calculated from dose divided by given AUCinfinity) [45] in other studies. Although this variation may result from minor differences in the pharmaceutical preparations, data analysis, and study populations, it could also indicate that some not yet clearly defined environmental factors have a major influence on the pharmacokinetics of talinolol.
To the best of our knowledge this is the first twin study on the heritability of talinolol pharmacokinetics as a marker of variation in drug membrane transport activities. We observed unexpectedly small heritability. The correlations between the siblings were nominally even higher for DZ than for MZ twin pairs (Fig. 1b, c) and, consistent with this, structural equation modeling revealed that common and unique environmental effects best explained the variation in talinolol oral clearance. Mechanistically, unique environmental effects include everything that acted randomly and not in the same manner in both siblings, whereas common environmental effects may include everything that is common to the siblings arising during pregnancy, childhood, and later life. Such common environmental effects may include all types of epigenetic effects but also any other factors which may be common due to growing up in the same family and environment.
The twin-study independent approach of Kalow et al. [55, 56], which uses repeated drug administration, revealed a low genetic component (rGC) of 0.38 and is generally in good agreement with the result from structural equation modeling. However, this component of 38 % includes both genetic effects and common environmental effects. This means that less (probably even much less) than 38 % of the variation is due to heritable factors. In contrast to the low genetic component in this study, other studies reveal high genetic components for other pharmacological variability, e.g., between 0.95 and 0.98 for the activity of N-acetyltransferase 2 (NAT2) [52, 55], 0.96 for plasma clearance of midazolam [60], 0.94 for the ethinylestradiol serum AUC [60], and 0.91 for the plasma AUC of metoprolol [51].
Several genetic polymorphisms in drug membrane transporters have already been analyzed in relation to talinolol pharmacokinetics. In our study, the genetic variants of MDR1 C1236T and G2677T/A were significantly associated with talinolol oral clearance, but this association was weak and even not significant after including other factors in a multifactorial analysis. For both polymorphisms clearance was enhanced in subjects heterozygously or homozygously carrying the 1236 T and 2677 T alleles, indicating higher activity of MDR1 compared to homozygous carriers of the 1236C or 2677G alleles. In line with this, both polymorphisms showed a lower AUC in subjects carrying hetero- or homozygous variants (1236 T and 2677 T or A) in another independent study on the oral bioavailability of talinolol [44]. However, two further studies found no significant effects of C3435T and G2677T/A [45, 61] as well as C1236T [45] on the talinolol AUC, indicating only weak functional effects and/or big effects of other not yet fully understood factors.
In this study we did not observe effects of some frequently studied and apparently functional MRP2, MDR5, and ABCG2 polymorphisms on talinolol clearance. The MRP2 finding is in contrast to Haenisch et al. [47], who reported an association of a lower AUC (and a lower oral bioavailability) and the 1249A variant. However, confirmation of earlier genotype–phenotype associations was not the primary aim of the present study and would have required a substantially bigger sample size, particularly for the quantitatively moderate effects.
The effects of genetic variation in efflux transporters on talinolol pharmacokinetics may be complicated by additional influences of genetically polymorphic influx transporters. A study by Bernsdorf et al. reported a significantly decreased talinolol half-life in carriers of the *1b variant in organic anion transporter OATP1B1 [48]. The *1b variant seems to code for a protein with high activity [62] but, in the study by Bernsdorf et al., subjects carrying the low-activity allele *15 of OATP1B1 also showed a decreased half-life. In line with these results, our study revealed that heterozygous carriers of the *1b allele (n = 21) tended to have a higher clearance compared to carriers of the wild-type allele *1a (Additional file 2) but homozygous carriers of the *1b allele (n = 4) had a lower clearance. The lack of consistency and the lack of a consistent effect in our study of the low activity allele *5 do not support a relevant role for OATP1B1 in the pharmacokinetics of talinolol.
Talinolol is a relatively hydrophilic drug and with a pKa of 9.3, about 98 % of the molecules are positively charged at physiological pH [50]. It is therefore a possible substrate of organic cation transporters [44]. For this reason we studied the influence of genetic variants in the genetically highly polymorphic transporter OCT1 on talinolol clearance. Surprisingly, although not significant, talinolol oral clearance tended to be increased in carriers of one or two reduced activity OCT1 alleles. Further in vitro research is required to elucidate whether OATP1B1 or OCT1 may enhance cellular uptake of talinolol.
Interestingly, protein and vegetable consumption both statistically significantly increased talinolol clearance. Our study was not designed to elucidate the underlying substances and mechanisms but it is reasonable that the pharmacokinetics of the mostly not metabolized drug talinolol is strongly dependent on transport activities in the intestinal mucosa and expression in these cells should depend on diet.
In an additional study, a subgroup of the volunteers received a low oral dose of 2.5 mg talinolol. The aim of this low-dose study was to enable a comparison between it and the 50 mg study. The so-called microdosing design has repeatedly been shown to be a valuable tool in the drug development process [63]. Although microdosing is defined by administration of doses below 100 μg, and in this sense we did not perform microdosing, the observations with our low range millidosing are notable. Interestingly, we observed quite different results compared to the therapeutic dose. The dose-normalized AUC7h was only about 50 % in the low dose talinolol group compared with the AUC from the comparable time interval in the 50 mg talinolol group. Also, tmax was achieved significantly earlier in the 50 mg talinolol group. These differences were already noted by Wetterich et al. [50]. One explanation may be that, in the high dose group, efflux transporters like MDR1 may be saturated very early so drug elimination, especially due to intestinal efflux, is delayed [50]. The very low correlation of the AUC7h for 2.5 mg and 50 mg talinolol (Fig. 2c) may similarly indicate first that this may be one example of a substance for which microdosing is not highly predictive for pharmacokinetics in therapeutic doses and second that the relevant transport processes and proteins involved may be different depending on the doses. This is supported by the fact that this correlation (Fig. 2c) was much lower than the correlation between the different study days with the same dose and between the siblings (Fig. 1d, e). In contrast to our study, another study with another MDR1 probe drug, fexofenadine, showed more rapid absorption of the drug with the smaller dose [64]. In this other study, however, 100 μg of fexofenadine was administered in 10 ml physiological saline and quickly swallowed and the high dose of 120 mg was administered as a tablet. Therefore, the different formulations may explain the contradictory results to our study.
Conclusions
The absorption, distribution, and elimination of talinolol seem to be mainly determined by environmental effects. This may indicate that the in vivo activity of MDR1 may also be dominantly determined by environmental factors. Some data presented in this study, however, may raise doubts about whether or not talinolol is really a good in vivo probe drug for MDR1. Further experimental and clinical research on the role of OATP1B1 and OCT1 may be scientifically interesting, but from our data one might simply conclude that frequent genetic polymorphisms with dominant effects on talinolol pharmacokinetics do not exist.
Abbreviations
- ABCC2:
-
ATP binding cassette subfamily C member 2
- AIC:
-
Akaike information criterion
- AUC:
-
Area under the curve
- Cl/F:
-
Total plasma clearance/bioavailability
- CV:
-
Coefficient of variation
- DZ:
-
Dizygotic
- LC-MS/MS:
-
Liquid chromatography–tandem mass spectrometry
- LOQ:
-
Limit of quantification
- MDR1:
-
Multi-drug resistance protein 1 (P-glycoprotein)
- MRP2:
-
Multidrug resistance-associated protein 2
- MZ:
-
Monozygotic
- OATP1B1:
-
Organic anion-transporting polypeptide 1B1
- OCT1:
-
Organic cation transporter 1
- rGC:
-
Genetic component
References
Pastan I, Gottesman MM. Multidrug resistance. Annu Rev Med. 1991;42:277–86.
Norris MD, Bordow SB, Marshall GM, Haber PS, Cohn SL, Haber M. Expression of the gene for multidrug-resistance-associated protein and outcome in patients with neuroblastoma. N Engl J Med. 1996;334:231–8.
Juliano RL, Ling V. A surface glycoprotein modulating drug permeability in Chinese hamster ovary cell mutants. Biochim Biophys Acta. 1976;455:152–62.
Chen CJ, Chin JE, Ueda K, Clark DP, Pastan I, Gottesman MM, Roninson IB. Internal duplication and homology with bacterial transport proteins in the mdr1 (P-glycoprotein) gene from multidrug-resistant human cells. Cell. 1986;47:381–9.
Schinkel AH. The physiological function of drug-transporting P-glycoproteins. Semin Cancer Biol. 1997;8:161–70.
Greiner B, Eichelbaum M, Fritz P, Kreichgauer HP, von Richter O, Zundler J, Kroemer HK. The role of intestinal P-glycoprotein in the interaction of digoxin and rifampin. J Clin Invest. 1999;104:147–53.
Fromm MF. P-glycoprotein: a defense mechanism limiting oral bioavailability and CNS accumulation of drugs. Int J Clin Pharmacol Ther. 2000;38:69–74.
Fellay J, Marzolini C, Meaden ER, Back DJ, Buclin T, Chave JP, Decosterd LA, Furrer H, Opravil M, Pantaleo G, et al. Response to antiretroviral treatment in HIV-1-infected individuals with allelic variants of the multidrug resistance transporter 1: a pharmacogenetics study. Lancet. 2002;359:30–6.
Sadhasivam S, Chidambaran V, Zhang X, Meller J, Esslinger H, Zhang K, Martin LJ, McAuliffe J. Opioid-induced respiratory depression: ABCB1 transporter pharmacogenetics. Pharmacogenomics J. 2015;15:119–26.
Bruhn O, Cascorbi I. Polymorphisms of the drug transporters ABCB1, ABCG2, ABCC2 and ABCC3 and their impact on drug bioavailability and clinical relevance. Expert Opin Drug Metab Toxicol. 2014;10:1337–54.
Wolking S, Schaeffeler E, Lerche H, Schwab M, Nies AT. Impact of genetic polymorphisms of ABCB1 (MDR1, P-glycoprotein) on drug disposition and potential clinical implications: update of the literature. Clin Pharmacokinet. 2015;54:709–35.
Geick A, Eichelbaum M, Burk O. Nuclear receptor response elements mediate induction of intestinal MDR1 by rifampin. J Biol Chem. 2001;276:14581–7.
Slater LM, Sweet P, Stupecky M, Gupta S. Cyclosporine-A reverses vincristine and daunorubicin resistance in acute lymphatic-leukemia in vitro. J Clin Investig. 1986;77:1405–8.
Tsuruo T, Iida H, Tsukagoshi S, Sakurai Y. Overcoming of vincristine resistance in P388 leukemia in vivo and in vitro through enhanced cytotoxicity of vincristine and vinblastine by verapamil. Cancer Res. 1981;41:1967–72.
Hoffmeyer S, Burk O, von Richter O, Arnold HP, Brockmoller J, Johne A, Cascorbi I, Gerloff T, Roots I, Eichelbaum M, Brinkmann U. Functional polymorphisms of the human multidrug-resistance gene: multiple sequence variations and correlation of one allele with P-glycoprotein expression and activity in vivo. Proc Natl Acad Sci U S A. 2000;97:3473–8.
Sakaeda T, Nakamura T, Horinouchi M, Kakumoto M, Ohmoto N, Sakai T, Morita Y, Tamura T, Aoyama N, Hirai M, et al. MDR1 genotype-related pharmacokinetics of digoxin after single oral administration in healthy Japanese subjects. Pharm Res. 2001;18:1400–4.
Becquemont L, Verstuyft C, Kerb R, Brinkmann U, Lebot M, Jaillon P, Funck-Brentano C. Effect of grapefruit juice on digoxin pharmacokinetics in humans. Clin Pharmacol Ther. 2001;70:311–6.
von Ahsen N, Richter M, Grupp C, Ringe B, Oellerich M, Armstrong VW. No influence of the MDR-1 C3435T polymorphism or a CYP3A4 promoter polymorphism (CYP3A4-V allele) on dose-adjusted cyclosporin A trough concentrations or rejection incidence in stable renal transplant recipients. Clin Chem. 2001;47:1048–52.
Min DI, Ellingrod VL. C3435T mutation in exon 26 of the human MDR1 gene and cyclosporine pharmacokinetics in healthy subjects. Ther Drug Monit. 2002;24:400–4.
Drescher S, Schaeffeler E, Hitzl M, Hofmann U, Schwab M, Brinkmann U, Eichelbaum M, Fromm MF. MDR1 gene polymorphisms and disposition of the P-glycoprotein substrate fexofenadine. Br J Clin Pharmacol. 2002;53:526–34.
Verstuyft C, Schwab M, Schaeffeler E, Kerb R, Brinkmann U, Jaillon P, Funck-Brentano C, Becquemont L. Digoxin pharmacokinetics and MDR1 genetic polymorphisms. Eur J Clin Pharmacol. 2003;58:809–12.
Wang D, Sadee W. Searching for polymorphisms that affect gene expression and mRNA processing: example ABCB1 (MDR1). AAPS J. 2006;8:E515–20.
Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM. A “silent” polymorphism in the MDR1 gene changes substrate specificity. Science. 2007;315:525–8.
Fung KL, Gottesman MM. A synonymous polymorphism in a common MDR1 (ABCB1) haplotype shapes protein function. Biochim Biophys Acta. 2009;1794:860–71.
Kim RB, Leake BF, Choo EF, Dresser GK, Kubba SV, Schwarz UI, Taylor A, Xie HG, McKinsey J, Zhou S, et al. Identification of functionally variant MDR1 alleles among European Americans and African Americans. Clin Pharmacol Ther. 2001;70:189–99.
Kroetz DL, Pauli-Magnus C, Hodges LM, Huang CC, Kawamoto M, Johns SJ, Stryke D, Ferrin TE, DeYoung J, Taylor T, et al. Sequence diversity and haplotype structure in the human ABCB1 (MDR1, multidrug resistance transporter) gene. Pharmacogenetics. 2003;13:481–94.
Cui Y, Konig J, Buchholz JK, Spring H, Leier I, Keppler D. Drug resistance and ATP-dependent conjugate transport mediated by the apical multidrug resistance protein, MRP2, permanently expressed in human and canine cells. Mol Pharmacol. 1999;55:929–37.
Schaub TP, Kartenbeck J, Konig J, Spring H, Dorsam J, Staehler G, Storkel S, Thon WF, Keppler D. Expression of the MRP2 gene-encoded conjugate export pump in human kidney proximal tubules and in renal cell carcinoma. J Am Soc Nephrol. 1999;10:1159–69.
Kamisako T, Gabazza EC, Ishihara T, Adachi Y. Molecular aspects of organic compound transport across the plasma membrane of hepatocytes. J Gastroenterol Hepatol. 1999;14:405–12.
Keppler D. The roles of MRP2, MRP3, OATP1B1, and OATP1B3 in conjugated hyperbilirubinemia. Drug Metab Dispos. 2014;42:561–5.
Pedersen JM, Matsson P, Bergstrom CA, Norinder U, Hoogstraate J, Artursson P. Prediction and identification of drug interactions with the human ATP-binding cassette transporter multidrug-resistance associated protein 2 (MRP2; ABCC2). J Med Chem. 2008;51:3275–87.
Naesens M, Kuypers DR, Verbeke K, Vanrenterghem Y. Multidrug resistance protein 2 genetic polymorphisms influence mycophenolic acid exposure in renal allograft recipients. Transplantation. 2006;82:1074–84.
Levesque E, Benoit-Biancamano MO, Delage R, Couture F, Guillemette C. Pharmacokinetics of mycophenolate mofetil and its glucuronide metabolites in healthy volunteers. Pharmacogenomics. 2008;9:869–79.
Ogasawara K, Chitnis SD, Gohh RY, Christians U, Akhlaghi F. Multidrug resistance-associated protein 2 (MRP2/ABCC2) haplotypes significantly affect the pharmacokinetics of tacrolimus in kidney transplant recipients. Clin Pharmacokinet. 2013;52:751–62.
Renders L, Frisman M, Ufer M, Mosyagin I, Haenisch S, Ott U, Caliebe A, Dechant M, Braun F, Kunzendorf U, Cascorbi I. CYP3A5 genotype markedly influences the pharmacokinetics of tacrolimus and sirolimus in kidney transplant recipients. Clin Pharmacol Ther. 2007;81:228–34.
FDA. Drug interaction studies--study design, data analysis, implications for dosing, and labeling recommendations. 2012. http://www.fda.gov/downloads/Drugs/GuidanceComplianceRegulatoryInformation/Guidances/UCM292362.pdf.
Ma JD, Tsunoda SM, Bertino Jr JS, Trivedi M, Beale KK, Nafziger AN. Evaluation of in vivo P-glycoprotein phenotyping probes: a need for validation. Clin Pharmacokinet. 2010;49:223–37.
Oertel R, Richter K, Trausch B, Berndt A, Gramatte T, Kirch W. Elucidation of the structure of talinolol metabolites in man. Determination of talinolol and hydroxylated talinolol metabolites in urine and analysis of talinolol in serum. J Chromatogr B Biomed Appl. 1994;660:353–63.
Trausch B, Oertel R, Richter K, Gramatte T. Disposition and bioavailability of the beta 1-adrenoceptor antagonist talinolol in man. Biopharm Drug Dispos. 1995;16:403–14.
Chu X, Zhang Z, Yabut J, Horwitz S, Levorse J, Li XQ, Zhu L, Lederman H, Ortiga R, Strauss J, et al. Characterization of multidrug resistance 1a/P-glycoprotein knockout rats generated by zinc finger nucleases. Mol Pharmacol. 2012;81:220–7.
Gramatte T, Oertel R. Intestinal secretion of intravenous talinolol is inhibited by luminal R-verapamil. Clin Pharmacol Ther. 1999;66:239–45.
Giessmann T, May K, Modess C, Wegner D, Hecker U, Zschiesche M, Dazert P, Grube M, Schroeder E, Warzok R, et al. Carbamazepine regulates intestinal P-glycoprotein and multidrug resistance protein MRP2 and influences disposition of talinolol in humans. Clin Pharmacol Ther. 2004;76:192–200.
Westphal K, Weinbrenner A, Zschiesche M, Franke G, Knoke M, Oertel R, Fritz P, von Richter O, Warzok R, Hachenberg T, et al. Induction of P-glycoprotein by rifampin increases intestinal secretion of talinolol in human beings: a new type of drug/drug interaction. Clin Pharmacol Ther. 2000;68:345–55.
Schwarz UI, Hanso H, Oertel R, Miehlke S, Kuhlisch E, Glaeser H, Hitzl M, Dresser GK, Kim RB, Kirch W. Induction of intestinal P-glycoprotein by St John’s wort reduces the oral bioavailability of talinolol. Clin Pharmacol Ther. 2007;81:669–78.
Schwarz UI, Seemann D, Oertel R, Miehlke S, Kuhlisch E, Fromm MF, Kim RB, Bailey DG, Kirch W. Grapefruit juice ingestion significantly reduces talinolol bioavailability. Clin Pharmacol Ther. 2005;77:291–301.
Schwarz UI, Gramatte T, Krappweis J, Berndt A, Oertel R, von Richter O, Kirch W. Unexpected effect of verapamil on oral bioavailability of the beta-blocker talinolol in humans. Clin Pharmacol Ther. 1999;65:283–90.
Haenisch S, May K, Wegner D, Caliebe A, Cascorbi I, Siegmund W. Influence of genetic polymorphisms on intestinal expression and rifampicin-type induction of ABCC2 and on bioavailability of talinolol. Pharmacogenet Genomics. 2008;18:357–65.
Bernsdorf A, Giessmann T, Modess C, Wegner D, Igelbrink S, Hecker U, Haenisch S, Cascorbi I, Terhaag B, Siegmund W. Simvastatin does not influence the intestinal P-glycoprotein and MPR2, and the disposition of talinolol after chronic medication in healthy subjects genotyped for the ABCB1, ABCC2 and SLCO1B1 polymorphisms. Br J Clin Pharmacol. 2006;61:440–50.
Zheng M, Zhang H, Dill DL, Clark JD, Tu S, Yablonovitch AL, Tan MH, Zhang R, Rujescu D, Wu M, et al. The role of Abcb5 alleles in susceptibility to haloperidol-induced toxicity in mice and humans. PLoS Med. 2015;12:e1001782.
Wetterich U, Spahn-Langguth H, Mutschler E, Terhaag B, Rosch W, Langguth P. Evidence for intestinal secretion as an additional clearance pathway of talinolol enantiomers: concentration- and dose-dependent absorption in vitro and in vivo. Pharm Res. 1996;13:514–22.
Matthaei J, Brockmoller J, Tzvetkov MV, Sehrt D, Sachse-Seeboth C, Hjelmborg JB, Moller S, Halekoh U, Hofmann U, Schwab M, Kerb R. Heritability of metoprolol and torsemide pharmacokinetics. Clin Pharmacol Ther. 2015;98:611–21.
Matthaei J, Tzvetkov MV, Strube J, Sehrt D, Sachse-Seeboth C, Hjelmborg JB, Moller S, Halekoh U, Hofmann U, Schwab M, et al. Heritability of caffeine metabolism: Environmental effects masking genetic effects on CYP1A2 activity but not on NAT2. Clin Pharmacol Ther. 2016. doi:10.1002/cpt.444 [Epub ahead of print].
Holst K, Scheike T. mets: analysis of multivariate event times, R package version 1.1.1. 2015. http://CRAN.R-project.org/package = mets.
R Core Team. R: A language and environment for statistical computing. Vienna: R foundation for statistical computing; 2015. http://www.R-project.org.
Kalow W, Tang BK, Endrenyi L. Hypothesis: comparisons of inter- and intra-individual variations can substitute for twin studies in drug research. Pharmacogenetics. 1998;8:283–9.
Kalow W, Endrenyi L, Tang B. Repeat administration of drugs as a means to assess the genetic component in pharmacological variability. Pharmacology. 1999;58:281–4.
Winkler G, Doring A. Validation of a short qualitative food frequency list used in several German large scale surveys. Z Ernahrungswiss. 1998;37:234–41.
Vesell ES. Twin studies in pharmacogenetics. Hum Genet Suppl. 1978;1:19–30.
Westphal K, Weinbrenner A, Giessmann T, Stuhr M, Franke G, Zschiesche M, Oertel R, Terhaag B, Kroemer HK, Siegmund W. Oral bioavailability of digoxin is enhanced by talinolol: evidence for involvement of intestinal P-glycoprotein. Clin Pharmacol Ther. 2000;68:6–12.
Ozdemir V, Kalow W, Tang BK, Paterson AD, Walker SE, Endrenyi L, Kashuba ADM. Evaluation of the genetic component of variability in CYP3A4 activity: a repeated drug administration method. Pharmacogenetics. 2000;10:373–88.
Siegmund W, Ludwig K, Giessmann T, Dazert P, Schroeder E, Sperker B, Warzok R, Kroemer HK, Cascorbi I. The effects of the human MDR1 genotype on the expression of duodenal P-glycoprotein and disposition of the probe drug talinolol. Clin Pharmacol Ther. 2002;72:572–83.
Mwinyi J, Johne A, Bauer S, Roots I, Gerloff T. Evidence for inverse effects of OATP-C (SLC21A6) 5 and 1b haplotypes on pravastatin kinetics. Clin Pharmacol Ther. 2004;75:415–21.
Burt T, Yoshida K, Lappin G, Vuong L, John C, de Wildt SN, Sugiyama Y, Rowland M. Microdosing and other phase-0 clinical trials: facilitating translation in drug development. Clin Transl Sci. 2016;9:74-88.
Lappin G, Shishikura Y, Jochemsen R, Weaver RJ, Gesson C, Houston B, Oosterhuis B, Bjerrum OJ, Rowland M, Garner C. Pharmacokinetics of fexofenadine: evaluation of a microdose and assessment of absolute oral bioavailability. Eur J Pharm Sci. 2010;40:125–31.
Acknowledgements
The contributions of Jakob Strube, Josephine Salaw, Stefanie Brunke, and Ulrike Schöning to the clinical study and Karoline Jobst, Robert Karandi, and Markus König to the laboratory analyses are gratefully acknowledged.
Funding
Financial support for the study was obtained from the German Ministry for Education and Research, BMBF, project “Virtual Liver” grant 0315755, from the Robert-Bosch Foundation (Stuttgart, Germany), and additional support from AstraZeneca.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are not publicly available due to data and privacy protection considerations but may be available on justified request.
Authors’ contributions
JB, DS, MS, and RK designed research. JM, JB, MVT, VG, CSS, DS, and UH performed research. JM, JB, and JBH analyzed data. JM, JB, MVT, and UH wrote the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the University of Göttingen, Germany (reference number 13/10/08) and by the German Federal Drug Administration (BfArM), (reference number 4034688). All subjects gave their written informed consent before participation.
Author information
Authors and Affiliations
Corresponding author
Additional files
Additional file 1:
Demographic data of the studied individuals. (DOCX 15 kb)
Additional file 2:
Influence of OCT1 and OATP1B1 on talinolol clearance. (DOCX 15 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Matthaei, J., Tzvetkov, M.V., Gal, V. et al. Low heritability in pharmacokinetics of talinolol: a pharmacogenetic twin study on the heritability of the pharmacokinetics of talinolol, a putative probe drug of MDR1 and other membrane transporters. Genome Med 8, 119 (2016). https://doi.org/10.1186/s13073-016-0372-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13073-016-0372-2