Abstract
Talniflumate is a phthalidyl ester of niflumic acid, which has potent analgesic and anti-inflammatory effects and is widely used to treat inflammatory disorders, such as rheumatoid arthritis. To screen the possible genetic factors affecting the pharmacokinetics (PK) of talniflumate, 23 male Korean volunteers were enrolled from two separate bioequivalence studies. All subjects received 740 mg (two tablets) talniflumate in a standard 2×2 cross-over model in a randomized order. For the genetic study, PK parameters of the reference drug were used. We used Illumina Human610Quad v1.0 DNA Analysis BeadChip for whole genome single nucleotide polymorphism (SNP) analysis and whole genome genotyping data were processed by linear regression analysis for PK parameters. Whole genome analysis revealed 1498 significant SNPs (P < 0.0001) for Cmax, 65 significant SNPs (P < 0.0001) for Tmax, and 1491 significant SNPs (P < 0.0001) for AUCinf. For clinical pharmacological purposes, we selected SNPs from drug metabolizing enzymes and transporters, and analyzed the PK parameters of various genotypes. Two SNPs (rs11165069 from ABCA4 (p=0.00002); rs17847036 from CYP2C9 (p=0.000001)) showed significant associations with talniflumate Cmax. In the Tmax group, two SNPs (rs3787555 from CYP24A1 (p=0.00035); rs2275034 from ABCA4 (p=0.000587)) showed significant associations with talniflumate Tmax. In the AUCinf group, two SNPs (rs11165069 from ABCA4 (p=0.00002); rs12461006 from SLC1A6 (p=0.00008)) exhibited significant associations with talniflumate absorption. These results show that genetic factors could affect the PK parameters, and provide information that may be used in the development of personalized talniflumate therapy.
Talniflumate, 3-pyridinecarboxylic acid 2-{[(3-trifluoromethyl) phenyl]amino}-1,3-dihydro-3-oxo-1-soenzofuranyl ester, is a phthalidyl ester of niflumic acid and is currently used to treat osteoarthritis and rheumatoid arthritis.[12] Niflumic acid shows rapid absorption followed by extensive metabolism, involving hydroxylation or glucuronidation.[3] However, the occurrence of side effects such as gastrointestinal irritation has been reported both in experimental animals and in clinical use. In order to lower ulcergenic activity, talniflumate was synthesized by esterification of the carboxyl group of niflumic acid.[4] Many genetic and non-genetic factors such as age, sex, life style, organ function, concomitant therapy, and the nature of the disease influence the effects of medications.[5] Following the Human Genome Project, it is understood that genetic factors can influence the response of an individual to a drug.[6] In clinical pharmacogenomics, drug-metabolizing enzymes.[78] drug transporters.[9] and other genes have been known to influence the individual differences in drug response.
Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation, and are proven to significantly affect drug absorption, distribution, biotransformation, and excretion.[10] Our previous studies showed that the pharmacokinetics (PK) of simvastatin and pioglitazone are affected by the genotype of each individual.[1112]
The purpose of this study was to screen possible candidate genes affecting the PK of talniflumate. A total of 23 male Korean volunteers were recruited from two bioequivalence studies after approval by the Institutional Review Board (IRB) of Kyung Hee University Hospital. In this study, DNA from each subject was analyzed using Illumina Human610Quad v1.0 DNA Analysis BeadChip. Linear regression analysis was performed for significant SNPs affecting PK parameters, such as maximum measured plasma concentration (Cmax), time of the maximum measured plasma concentration (Tmax), and area under the plasma concentration-time curve from zero to infinity (AUCinf). For clinical pharmacological purposes, further analysis was performed on significant SNPs encoding phase I and II drug metabolizing enzymes to examine their effect on the PK of talniflumate.
Volunteers were healthy Korean males who participated in two talniflumate bioequivalence tests at Kyung Hee Clinical Medical Research Institute of Kyung Hee University Hospital. The clinical protocol was approved by the IRB of Kyung Hee University Hospital and volunteers were recruited by direct call. Twenty-three out of 60 subjects from these previous tests participated in this pharmacogenetics study after providing written informed consent. The demographic characteristics of the volunteers are summarized in Table 1. Subjects ranged in age from 20 to 35 years (26 ± 3.6 years), in weight from 57.0 to 90.0 kg (68.7 ± 8.3 kg), and in height from 160.0 to 184.0 cm (173.3 ± 4.8 cm) (Table 1).
Previous studies were based on two talniflumate bioequivalence tests. Each bioequivalence study of two 370 mg talniflumate formulations (reference drug, Somalgen tablet 370 mg, Kunhwa Co., Ltd.) was conducted in healthy male Korean volunteers after a single dose administration in a randomized cross-over study with a washout period of at least one week. The subjects were hospitalized (Kyung Hee University Hospital, Seoul, Korea) and fasted overnight (10 h) and for 4 h after each drug administration. The doses were administered at 8.00 a.m. of each dosing day along with 240 ml of tap water. No food was allowed until 4 h after dose administration. Approximately 7 ml of blood for talniflumate and active metabolite assays were drawn into heparinized tubes through indwelling cannula before (0 h) and at 0.5, 1, 1.5, 2, 2.5, 3, 4, 6, 9, 12, 15 h after dosing. Blood samples were centrifuged at 3,000 rpm for 10 min; plasma was separated and kept frozen at −70℃ until assayed. Plasma was analyzed for talniflumate concentration using a validated high-performance liquid chromatography method.
PK parameters (Cmax, Tmax, AUCt and AUCinf) were calculated by noncompartmental models in WinNonlin v5.2 (Table 1).
DNA extraction was performed as previously described.[11] From May to August 2008, blood samples were obtained from 23 participants. After obtaining informed consent, approximately 10 ml of whole blood was collected from each participant. Blood samples were drawn into sterile tubes containing ethylenediaminetetraacetic acid (EDTA) and stored at −70℃ until the isolation of genomic DNA by a standard phenol chloroform extraction method.
SNPs of 23 healthy male volunteers were analyzed by Standard Illumina procedures using Illumina BeadStation 500G (Illumina Human610Quad v1.0 DNA Analysis BeadChip) as previously described.[1112] Intensity files (*.idat) were processed by BeadStudio GT module 3.3.4 with default analysis settings. Each SNP was analyzed independently to cluster and identify genotypes. Genotype calls were generated by comparing experimental data with those in the supplied cluster file (*.egt). The SNP set was filtered on the basis of genotype call rates (≥ 95%), and minor allele frequency (MAF ≥ 0.05).
Statistical analysis was performed as previously described.[1112] Hardy Weinberg equilibrium (HWE) was calculated for individual SNPs using an exact test. All of the SNPs reported in this manuscript have HWE P-values > 0.001. To calculate the degree (Beta), 95% confidence intervals (CI) and P-value while controlling for age, height and weight as covariants in a linear regression analysis with significant SNPs on individual PK parameters of talniflumate, unadjusted genotypic association with additive, dominant and recessive models were tested by calculating the Beta and P-value using PLINK version 1.06 (Shaun Purcell, USA). When the additive model was used, each genotype was independently coded as 0, 1, or 2 and analyzed. In the dominant model, a homozygote major allele and another two genotypes were coded as 0 and 1 and analyzed. In the recessive model, a homozygote minor allele and another two genotypes were coded a 0 and 1 and analyzed. After filtering, SNPs were analyzed on chromosome 1 through chromosome 22 for each group.
Significant SNPs were identified from each PK parameter (Cmax, Tmax and AUCinf). These included SNPs of phase I and II drug metabolizing enzymes, which were selected to further analyze their effect on the PK parameters (two SNPs in Cmax, two SNPs in Tmax and two SNPs in AUCinf).
The PK parameters of talniflumate are shown in Table 1. Cmax is 230.9 ± 113.5 µg/L, Tmax is 3.0 ± 1.2 h, AUCt is 1058.6 ± 557.1 µg·h/L, and AUCinf is 1100 ± 621.9 µg/L. These data were used for linear regression analysis.
Linear regression analysis was performed on the PK parameters of talniflumate (Cmax, Tmax, and AUCinf). In this analysis, 25,361 significant SNPs in the Cmax group, 24,067 significant SNPs in the Tmax group, and 23,604 significant SNPs in the AUCinf group showed significant associations (P<0.05). Among those significant SNPs, 1,498 SNPs from the Cmax group, 65 SNPs from the Tmax group, and 1,491 SNPs from the AUCinf group showed highly significant associations with each parameter (Table 2, 3 and 4, P<0.0001). All data were adjusted for age, height, and weight. For clinical pharmacological purposes, significant SNPs encoding phase I and II drug metabolizing enzymes were selected for further analysis.
We selected two SNPs for further analysis from the Cmax group (rs11165069 from ABCA4; rs17847036 from CYP2C9 Table 2), the Tmax group (rs3787555 from CYP24A1; rs2275034 from ABCA4, Table 3), and the AUCinf group (rs11165069 from ABCA4; rs12461006 from SLC1A6, Table 4).
ABCA4 (rs11165069, P = 0.000001 in Cmax group; rs2275034, P = 0.000587 in the Tmax group; rs11165069, P = 0.00002 in AUCinf group) showed statistically significant changes by genotype (Fig. 1, 2, and 3 and Table 5). ABCA4 gene polymorphisms affect talniflumate PK parameters. These data show that genetic polymorphism of ABCA4 has a profound effect on the PK of talniflumate.
In the Cmax group, CYP2C9 rs17847036 showed statistically significant changes according to genotype; the AG genotype exhibited a greater Cmax value (Fig. 1 and Table 5).
Many factors affect talniflumate blood concentration, such as diet, age, liver function and disease condition. The genes involved in the process of drug absorption, distribution, metabolism and excretion may affect the PK of talniflumate [101314] but the effect of pharmacogenetics factors has yet to be investigated. In this study, we performed whole genome association studies and regression analysis to investigate the genetic factors influencing the PK of talniflumate.
Regression analysis of significant SNPs and SNPs of candidate genes revealed several SNPs which significantly affect the PK of talniflumate. Interestingly, although SNPs of ABCA4 have low frequency, their effect was significant, resulting in increases in Cmax, Tmax and AUCinf. Although these results were only obtained via screening studies, they are the first data concerning the influence of genes, including those responsible for metabolizing enzymes and transporters, on the PK of talniflumate.
We identified the significant association of SNP rs11165069 (P=0.000001) (Fig. 1A) in intron 24 of the transporter gene ABCA4 with the Cmax and AUCinf of talniflumate. In our study, ABCA4 (rs1165069) with a TC increased the absorption of talniflumate by approximately four times compared to the CC genotype. Therefore, the TC genotype resulted in greater Cmax (0.000001) and AUCinf (P=0.00002) values. However, because of the low frequency of the TC genotype of ABCA4 rs11165069, further analysis is required. Furthermore, subjects in the Tmax group with the minor C allele of the intron 37 SNP rs2275034 (P=0.00058) (Fig. 2B) of ABCA4 showed a more than 2-fold delayed Tmax. From the result, it can be assumed that ABCA4 can affect the absorption of talniflumate and it is possible that the PK parameters can vary according to genotype.
Subjects with the minor A allele of the synonymous exon 2 SNP rs17847036 (P=0.000001) of the phase I enzyme gene CYP2C9 showed an increase in Cmax (Fig. 1B). It is well understood that CYP2C9 influences the metabolism of drug types [15] and these results indicate that the genotype of CYP2C9 could influence talniflumate plasma concentration.
Other genes, such as CYP24A1 (rs3787555) and SLC1A6 (rs12461006), showed significant associations with the PK of talniflumate, influencing Tmax and AUCinf, respectively.
To the best of our knowledge, this is the first whole genome association study regarding the PK of talniflumate. However, the limitations of this study must be acknowledged. As this was a screening study, we performed only univariate analysis between genotypes and PK parameters. Therefore, further confirmative study on the PK of talniflumate with subgroup analysis according to genotypes is required. In conclusion, this study is useful in explaining why the PK of talniflumate may vary depending on the individual. As several additional drugs and enzymes may also influence talniflumate metabolism, we propose that genetic differences are major contributing factors in the PK of this drug.
Notes
References
1. Lee HW, Won KJ, Cho SH, Ha YH, Park WS, Yim HT, et al. Quantitation of niflumic acid in human plasma by high-performance liquid chromatography with ultraviolet absorbance detection and its application to a bioequivalence study of talniflumate tablets. J Chromatogr B Analyt Technol Biomed Life Sci. 2005; 821:215–220.
2. Knight D. Talniflumate (Genaera). Curr Opin Investig Drugs. 2004; 5:557–562.
3. Lan SJ, Chando TJ, Weliky I, Schreiber EC. Metabolism of niflumic acid-14C: absorption, excretion and biotransformation by human and dog. J Pharmacol Exp Ther. 1973; 186:323–330. PMID: 4719785.
4. Los M, Boned JE, Piccininali C. New esters of substituted anilinonicotinic and phenylanthranilic acids. Farmaco Sci. 1981; 36:372–385. PMID: 7238853.
5. Evans WE, McLeod HL. Pharmacogenomics--drug disposition, drug targets, and side effects. N Engl J Med. 2003; 348:538–549. PMID: 12571262.
6. Voisey J, Morris CP. SNP technologies for drug discovery: a current review. Curr Drug Discov Technol. 2008; 5:230–235. PMID: 18690891.
7. Gardiner SJ, Begg EJ. Pharmacogenetics, drug-metabolizing enzymes, and clinical practice. Pharmacol Rev. 2006; 58:521–590. PMID: 16968950.
8. Iyer KR, Sinz MW. Characterization of Phase I and Phase II hepatic drug metabolism activities in a panel of human liver preparations. Chem Biol Interact. 1999; 118:151–169. PMID: 10359459.
9. McLean C, Wilson A, Kim RB. Impact of Transporter Polymorphisms on Drug Development: Is It Clinically Significant? J Clin Pharmacol. 2016; 56(Suppl 7):S40–S58. DOI: 10.1002/jcph.691. PMID: 27385178.
10. Daly TM, Dumaual CM, Miao X, Farmen MW, Njau RK, Fu DJ, et al. Multiplex assay for comprehensive genotyping of genes involved in drug metabolism, excretion, and transport. Clin Chem. 2007; 53:1222–1230. PMID: 17510302.
11. Im S, Kim BH, Lee K, Kwack K, Yim SV. Screening study for genetic polymorphisms affecting pharmacokinetics of simvastatin. Transl Clin Pharmacol. 2016; 24:43–54.
12. Yun JY, Kim BH, Lee JH, Lee K, Kwack K, Yim SV. Screening study for genetic polymorphisms affecting pharmacokinetics of pioglitazone. Transl Clin Pharmacol. 2016; 24:194–202.
13. Mizuno N, Niwa T, Yotsumoto Y, Sugiyama Y. Impact of drug transporter studies on drug discovery and development. Pharmacol Rev. 2003; 55:425–461. PMID: 12869659.
14. Müller M. Pharmacogenomics and drug response. Int J Clin Pharmacol Ther. 2003; 41:231–240. PMID: 12816175.
15. Topić E, Stefanović M, Samardzija M. Association between the CYP2C9 polymorphism and the drug metabolism phenotype. Clin Chem Lab Med. 2004; 42:72–78. PMID: 15061384.
Table 1
Abbreviations: Cmax, maximum measured plasma concentration; Tmax, time of the maximum measured plasma concentration; AUCt, area under the plasma concentration-time curve from time zero to time of last measurable concentration; AUCinf, area under the plasma concentration-time curve from zero to infinity; SD, standard deviation.