Next Article in Journal
Editorial for the Special Issue: Pharmacokinetics of Orally Administered Drugs, 2nd Edition
Next Article in Special Issue
Computational Analysis of S1PR1 SNPs Reveals Drug Binding Modes Relevant to Multiple Sclerosis Treatment
Previous Article in Journal
Recent Review on Biological Barriers and Host–Material Interfaces in Precision Drug Delivery: Advancement in Biomaterial Engineering for Better Treatment Therapies
Previous Article in Special Issue
Novel Genetic Variants Explaining Severe Adverse Drug Events after Clinical Implementation of DPYD Genotype-Guided Therapy with Fluoropyrimidines: An Observational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variation in CYP2D6, UGT1A4, SLC6A2 and SLCO1B1 Alters the Pharmacokinetics and Safety of Mirabegron

by
Paula Soria-Chacartegui
1,
Patricia Cendoya-Ramiro
1,
Eva González-Iglesias
1,
Samuel Martín-Vílchez
1,
Andrea Rodríguez-Lopez
1,
Gina Mejía-Abril
1,
Manuel Román
1,
Sergio Luquero-Bueno
1,
Dolores Ochoa
1 and
Francisco Abad-Santos
1,2,*
1
Clinical Pharmacology Department, Hospital Universitario de La Princesa, Faculty of Medicine, Instituto de Investigación Sanitaria La Princesa (IP), Universidad Autónoma de Madrid (UAM), 28006 Madrid, Spain
2
Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd), Instituto de Salud Carlos III, 28029 Madrid, Spain
*
Author to whom correspondence should be addressed.
Pharmaceutics 2024, 16(8), 1077; https://doi.org/10.3390/pharmaceutics16081077
Submission received: 23 July 2024 / Revised: 13 August 2024 / Accepted: 15 August 2024 / Published: 17 August 2024

Abstract

:
Mirabegron is a drug used in overactive bladder (OAB) treatment. Genetic variation in pharmacogenes might alter its pharmacokinetics, affecting its efficacy and safety. This research aimed to analyze the impact of genetic variation on mirabegron pharmacokinetics and safety. Volunteers from three bioequivalence trials (n = 79), treated with a single or a multiple dose of mirabegron 50 mg under fed or fasting conditions, were genotyped for 115 variants in pharmacogenes and their phenotypes were inferred. A statistical analysis was performed, searching for associations between genetics, pharmacokinetics and safety. CYP2D6 intermediate metabolizers showed a higher elimination half-life (t1/2) (univariate p-value (puv) = 0.018) and incidence of adverse reactions (ADRs) (puv = 0.008, multivariate p (pmv) = 0.010) than normal plus ultrarapid metabolizers. The UGT1A4 rs2011425 T/G genotype showed a higher t1/2 than the T/T genotype (puv = 0.002, pmv = 0.003). A lower dose/weight corrected area under the curve (AUC/DW) and higher clearance (CL/F) were observed in the SLC6A2 rs12708954 C/C genotype compared to the C/A genotype (puv = 0.015 and 0.016) and ADR incidence was higher when the SLCO1B1 function was decreased (puv = 0.007, pmv = 0.010). The lower elimination and higher ADR incidence when CYP2D6 activity is reduced suggest it might be a useful biomarker in mirabegron treatment. UGT1A4, SLC6A2 and SLCO1B1 might also be involved in mirabegron pharmacokinetics.

1. Introduction

Overactive bladder (OAB) syndrome is a condition characterized by unstable bladder contractions, which causes urinary urgency, with or without incontinence, frequent urination and nocturia [1]. Observed in approximately 12% of the population, its incidence and morbidity increase with age, particularly after the age of forty [1,2,3]. It has a highly negative impact on the quality of life, which emphasizes the need for an effective treatment. OAB management is typically based on a combination of behavioral therapies, such as weight loss and pelvic floor exercises, and pharmacological treatment, which includes antimuscarinics (also known as anticholinergics) and β3-agonists [4]. Mirabegron is a first-in-class drug that belongs to the family of β3-adrenoreceptor agonists [5]. It was approved by the European Medicines Agency (EMA) in 2012 for OAB treatment, after it was demonstrated to reduce urinary frequency, incontinence and urgency and to improve health-related quality of life (HR-QOL) [5].
The mechanism of action of mirabegron is the upregulation of the sympathetic nervous activity, which promotes the relaxation of the detrusor muscle of the bladder. Consequently, the storage period before feeling the need to urinate is prolonged [1]. Mirabegron represents a promising pharmacological alternative to anticholinergics, which are associated with several adverse drug reactions (ADRs) that reduce treatment adherence [1]. Furthermore, β3-agonists represent the only effective treatment for OAB when it is caused by the metabolic syndrome or by menopausal hormonal changes [6].
Mirabegron is administered orally following a posology of 25–50 mg every 24 h, with an absolute bioavailability of 29–35% [7]. The maximum plasma concentration (Cmax) is reached 3–4 h (tmax) after drug intake, showing a plasmatic protein binding of approximately 71%. It is metabolized to inactivate products, potentially by the action of enzymes from the cytochrome P450 (CYP) and the UDP-glucuronosyltransferase (UGT) families, such as CYP2D6, CYP3A4, UGT2B7, UGT1A3 or UGT1A8 [7,8]. Mainly excreted through urine (55%) and feces (34%), its elimination half-life (t1/2) is 50 h, while the total body clearance (CL) from plasma is 57 L/h [7].
The most frequently reported ADRs in mirabegron treatment, which are less common than with anticholinergics [1], are headache, tachycardia, urinary tract infection and dizziness [7]. Its occurrence is influenced by several factors, such as sex, age, concomitant medication and genetics [9]. Consequently, pharmacogenetics may provide physicians with effective clues to optimize drug efficacy and safety [9]. Pharmacogenetics is an approach to precision medicine that tailors drug selection and dosing to the patient’s genetic features [10] and it is already implemented in clinical practice for a variety of drugs [11,12,13,14,15,16,17]. However, no pharmacogenetic clinical guidelines are available to direct the treatment with mirabegron or with any antimuscarinic drug [18]. Thus, in order to enhance pharmacogenetic implementation in OAB treatment, it is essential to investigate the effect of genetic variation on these drugs’ behavior and response. Consequently, the main objective of this pharmacogenetic research was to study the impact of genetic variation on mirabegron pharmacokinetics and safety. This research forms part of “La Princesa Multidisciplinary Initiative for the Implementation of Pharmacogenetics” (PriME-PGx) [19].

2. Materials and Methods

2.1. Study Population

The population included in the research was composed of healthy volunteers from three bioequivalence clinical trials conducted at the Clinical Trials Unit of the Hospital Universitario de La Princesa (UECHUP), Madrid, Spain. Subjects self-reported their biogeographic origin, sex and age. A total of 104 volunteers provided written informed consent to participate in the bioequivalence trials, and 79 of them provided written informed consent to be included in the observational pharmacogenetic study. The protocol and the informed consent form of the three clinical trials and the pharmacogenetic study (code SFC-FG-2020-1, IRB/Code: 4176) were reviewed and approved by the Spanish Drug Agency (AEMPS) and by the Ethics Committee on Clinical Research with Drugs of the Hospital Universitario de La Princesa. Moreover, the studies were conducted in accordance with the guidelines of Good Clinical Practice [20], current Spanish legislation and the Revised Declaration of Helsinki [21].

2.2. Clinical Trial Design

The three trials included in this research were phase I, open-label, randomized, replicated and crossover bioequivalence clinical trials of mirabegron 50 mg prolonged-release tablets (Table 1). Bioequivalence trials constitute a specific type of phase I clinical trial in which two formulations of the same drug are compared to ascertain the similarity of their pharmacokinetic profile [22]. These studies are frequently conducted in order to approve the commercialization of a test formulation (T) by demonstrating its bioequivalence to a reference formulation (R). In clinical trials A and B, the administration of a single mirabegron dose under fasting and fed conditions was respectively examined. Fed conditions in clinical trial B were ensured by a high-fat 900 kcal breakfast, in accordance with EMA guidelines for bioequivalence trials [23]. The administration of 12 oral doses in fasting conditions was analyzed in trial C (Table 1).
The inclusion criteria for all the clinical trials included males and females aged from 18 to 55, with no clinically significant organic or psychic conditions and no abnormalities in medical records, in laboratory analysis and in vital signs and electrocardiogram. The exclusion criteria included the use of any prescribed medication in the past 15 days; body mass index (BMI) outside the range 18.5–30 kg/m2; a history of drug sensitivity; use of recreational drugs, alcohol or tobacco; blood donation in the past month; pregnant or breastfeeding women; participation in other clinical trials in the past three months; and a history of swallowing problems.
To ensure both formulations were administered twice to each volunteer, clinical trials A and B were structured in four single dose periods, with four different sequences: T-R-T-R, T-R-R-T, R-T-R-T and R-T-T-R. Hospitalization lasted from at least 10 h before to 12 h after the intake of the T or the R formulation (depending on the sequence) in each period. The washout period, i.e., the time between periods, necessary to ensure complete drug elimination, was 14 days. Clinical trial C was also structured in four periods. However, 11 doses (once daily) were administered in the first and the third period, whereas the second and fourth period consisted of the administration of a single dose 24 h after the last administration from the previous period. Thus, there was no washout time and two different sequences were followed: T-T-R-R and R-R-T-T. Subjects were hospitalized from at least 10 h prior to the start of the first and the third period to 12 h following the conclusion of the second and fourth period, respectively.

2.3. Pharmacokinetic and Safety Analysis

In clinical trials A and B, at least 16 blood samples were collected from pre-dose to 72 h after drug intake in the four periods. In clinical trial C, between 16 and 21 blood samples were extracted from pre-dose to 24 h after the last dose in each period. In the three clinical trials, blood samples were immediately centrifuged after extraction, and plasma was stored at −20 °C until its shipment to an external laboratory for mirabegron plasmatic quantification. The measurement of mirabegron plasma concentration was performed in accordance with the EMA guidelines by reversed phase ultra-high-performance liquid chromatography coupled to tandem mass spectrometry (UPLC-MS/MS), with a lower limit of quantification (LLOQ) of 100 pg/mL. The study was only blinded for the determination of mirabegron plasmatic concentration.
The estimation of the pharmacokinetic parameters was performed by non-compartmental methods, using the WinNonLin Professional version 8.3 software (Scientific Consulting, Inc., Cary, NC, USA). The maximum plasma concentration (Cmax; Cmax-ss in steady state for clinical trial C) and the time to reach it (tmax) were directly obtained from the plasma concentration–time curves. The total area under the curve (AUC) was calculated as the sum of the AUC from pre-dose to the last observed concentration point (AUC0–t), obtained by the trapezoidal method, and the AUC from the t point to infinity (AUCt–∞). The AUCt–∞ was calculated as the ratio of the last detectable concentration (Ct) divided by the elimination slope (Ke) of the log-linear portion from the concentration–time curve obtained by linear regression. In trial C, the AUC between dosing intervals (AUCtau) was calculated with the trapezoidal method. The t1/2 was calculated as ln2/Ke and the total CL adjusted for bioavailability (CL/F) was calculated as the dose divided by AUC and weight. Both formulations were deemed to be bioequivalent in the three trials; thus, the arithmetic mean of the parameters obtained in the four periods was calculated for each subject.
The occurrence of adverse effects (AEs) was documented through direct questioning of the subjects and spontaneous reporting. The algorithm of the Spanish Pharmacovigilance System was applied to assess AE causality [24]. The AEs with a definite, probable or possible causal relationship with mirabegron were considered ADRs.

2.4. Genotyping and Phenotyping

DNA was extracted from the blood samples collected during the clinical trials using a Maxwell® RSC Automated DNA extractor (Promega Biotech Iberica S.L, Madrid, Spain). The customized array PriME-PGx Very Important Pharmacogene Open Array (VIPOA), based on TaqMan® probes, was used for genotyping in a QuantStudio 12 k Flex real-time PCR System (Applied Biosystems, ThermoFisher, Waltham, MA, USA) with an Open Array thermal block. A total of 115 genetic variants from genes codifying for metabolizing enzymes and transporters were analyzed (Supplementary Table S1). Nucleotide variant annotation was performed according to the National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database [25]. The sample size available for each genetic variant differs due to the occurrence of low-frequent genotyping failures.
The presence of duplications, deletions or hybrid structures in the gene CYP2D6 was studied. A copy number variation assay for both ends of the gene (5′ untranslated region and exon 9) was performed in the QuantStudio 12 k Flex real-time PCR System (Applied Biosystems, ThermoFisher, Waltham, MA, USA) with a 96-well thermal block [26].
The Clinical Pharmacogenetics Implementation Consortium (CPIC) and the Pharmacogene Variation Consortium (PharmVar) database were used for allele definition [27,28]. The pharmacogenetic phenotype for CYP2B6 [13], CYP2C19 [11], CYP2C9 [12], CYP2D6 [14], CYP3A4 [29], CYP3A5 [15], NUDT15 [16], SLCO1B1 [12], TPMT [16] and UGT1A1 [17] was inferred from diplotypes based on their corresponding CPIC and Dutch Pharmacogenetic Working Group (DPWG) guidelines. For CYP2C8 and NAT2, phenotype was inferred as described in previous works [30,31]. The remaining genetic variants were analyzed individually.

2.5. Statistical Analysis

SPSS Statistics software version 29.0 (SPSS Inc., Chicago, IL, USA) was used for the statistical analysis. To correct the effect of dose and weight on exposure, the AUC (trials A and B), the AUCtau (trial C) and the Cmax were divided by the dose/weight ratio (DW). The normality of the dependent variables was analyzed by the Shapiro–Wilk test. Those variables not following a normal distribution were logarithmic-transformed, and normality was re-analyzed. The significance level was established at p-value (p) < 0.05. The demographic characteristics (age, height, weight and BMI) were analyzed according to sex, biogeographic origin and clinical trial design.
First, a univariate pharmacokinetic analysis was performed. Pharmacokinetic parameters were analyzed according to sex, biogeographic origin, clinical trial design, genotypes and phenotypes. For the analysis of normally distributed dependent variables according to independent variables with two or with three or more categories, t-tests and ANOVA tests followed by Bonferroni post hoc tests were respectively performed. For the analysis of not normally distributed dependent variables according to independent variables with two or with three or more categories, Mann–Whitney tests and Kruskal–Wallis test were conducted, respectively. Afterwards, a multivariate pharmacokinetic analysis was performed by means of linear regression. Those independent variables with significant differences in the univariate analysis (univariate p-value (puv) < 0.05) were included. Significant associations are shown as the multivariate p-value (pmv), the standardized β coefficient (β) and the R2.
For safety evaluation, differences in the AUC and Cmax were analyzed in relation to ADR incidence using the t-test and the Mann–Whitney test, respectively. In the univariate safety analysis, ADR incidence was analyzed according to sex, biogeographic origin, clinical trial, genotypes and phenotypes by a Pearson’s Chi2 test or a Fisher’s exact test when the expected frequency in at least 20% of cells was lower than 5. Subsequently, a multivariate safety analysis was conducted by means of logistic regression. The independent variables significantly associated with ADR incidence in the univariate analysis (puv < 0.05) were included. Significance is shown as pmv and odds ratio (OR).

3. Results

3.1. Demographic Characteristics

The height and weight of females and Latin Americans were lower than in males (p < 0.001 and p < 0.001, respectively) and Europeans (p = 0.005 and p = 0.01, respectively) (Table 2). No differences were identified in the demographic parameters according to the clinical trial (Table 2).

3.2. Pharmacokinetics

Females showed a higher AUC and Cmax (455.74 ± 152.50 ng·h/mL and 47.25 [20.95–57.58] ng/mL) compared to males (354.88 ± 120.42 ng·h/mL and 33.11 [18.63–44.88] ng/mL, p = 0.002 and p = 0.009, respectively); however, the differences disappeared after DW correction (Table 3).
The AUC/DW and the Cmax/DW were lower in clinical trial B compared to the other clinical trials (AUC/DW: puv < 0.001 compared to trials A and C, pmv < 0.001, β = 0.668, R2 = 0.438; and Cmax/DW: puv < 0.001 compared to trials A and C, pmv < 0.001, β = 0.739, R2 = 0.539). Conversely, the tmax and the CL/F were higher (tmax: puv = 0.014 compared to trial A, puv < 0.001 compared to trial C, pmv < 0.001, β = −0.426, R2 = 0.169; and CL/F: puv < 0.001 compared to trials A and C, pmv < 0.001, β = −0.691, R2 = 0.469). Furthermore, the t1/2 in clinical trial C was lower than in clinical trials A and B (puv < 0.001 compared to trials A and B, pmv < 0.001, β = 0.535, R2 = 0.381). No differences in the pharmacokinetic parameters according to biogeographic origin were observed (Table 3).

3.3. Pharmacogenetics

No CYP2D6 poor metabolizers (PMs) were available in this research. A higher t1/2 was observed in CYP2D6 intermediate metabolizers (IMs) compared to ultrarapid (UMs) and normal metabolizers (NMs) (puv = 0.028). Significance was transformed into a tendency after Bonferroni post hoc tests (NMs vs. IMs, puv = 0.074, UMs vs. IMs puv = 0.124 and NMs vs. UMs puv = 0.624). However, the t1/2 remained significantly higher in IMs compared to UMs and NMs merged in a unique group (UMs + NMs) (puv = 0.018). Differences in the t1/2 were accompanied by a not significantly lower AUC/DW and Cmax/DW in CYP2D6 UMs compared to NMs and IMs. The UGT1A3 rs2008584 A/A genotype was associated with lower tmax than the UGT1A3 rs2008584 A/G genotype (puv = 0.014). Also, the t1/2 was higher in the UGT1A4 rs2011425 T/G genotype compared to the T/T genotype (puv = 0.002, pmv = 0.003, β = 0.268, R2 = 0.381) (Table 4).
The Cmax/DW was significantly lower in NAT2 slow acetilators (SAs) compared to intermediate (IAs) and rapid acetilators (RAs) (puv = 0.046). Even though differences disappeared in pair comparison, the observed lower Cmax/DW in NAT2 SAs was maintained when IAs and RAs were merged in a unique group (puv = 0.016), accompanied by a tendency towards a lower AUC/DW (puv = 0.089) and higher tmax and CL/F in SAs (puv = 0.083 and puv = 0.076, respectively).
A lower AUC/DW and a higher CL/F were observed in the SLC6A2 rs12708954 C/C genotype compared to the C/A genotype (puv = 0.015 and puv = 0.016, respectively). Lastly, a lower Cmax/DW was observed in the SLC19A1 rs1051266 A/A genotype compared to the A/G genotype (puv = 0.016) (Table 4). For CYP3A4, the AUC/DW and Cmax/DW were not significantly higher in IMs compared to NMs (Supplementary Table S2). No further differences were observed in the pharmacokinetic parameters according to the remaining genetic variants (Supplementary Table S2).

3.4. Safety

ADRs were present in ten volunteers (12.66%) in total (two in clinical trials A and C, and six in clinical trial B, Table 5), five of whom suffered one ADR, three of whom suffered two ADRs and two of whom suffered three ADRs. Headache was the most common ADR (9 times), followed by nausea (3 times), vomiting (2 times), loose stool, palpitations and insomnia (1 time each). ADR incidence was similar between males and females (7.90% vs. 17.07%, p = 0.314), between Europeans and Latin Americans (5.60% vs. 14.80%, p = 0.414) and among clinical trials A, B and C (9.10% vs. 20.00% vs. 7.40%, respectively, p = 0.331). The AUC and Cmax in volunteers suffering ADRs and volunteers not experiencing them were not significantly different when the three trials were analyzed collectively and individually, with the exception of a tendency towards a higher AUC in the volunteers experiencing ADRs compared to those with no ADRs in clinical trial B (p = 0.070) (Table 5).
CYP2D6 IMs showed a higher ADR incidence (37.5%) than NMs (5.5%) and UMs (0%) (puv = 0.008, pmv = 0.010, OR = 9.83). Also, ADRs were more frequent in volunteers with decreased function (DF) for SLCO1B1 (30.4%) compared to those with normal function (NF) (5.7%) (puv = 0.007, pmv = 0.010, OR = 11.24). No further differences in ADR incidence according to genotypes or phenotypes were observed.

4. Discussion

The successful treatment of the OAB condition is of great relevance due to its negative impact on quality of life and healthcare systems and its indirect costs due to the loss of work capacity [32,33]. Using mirabegron for OAB treatment could potentially decrease the global treatment costs in comparison to antimuscarinic drugs [34]. However, the effect of genetic variation on mirabegron response is still uncertain. Thus, this research aimed to identify genetic factors explaining variability in mirabegron pharmacokinetics and safety.
The pharmacokinetic parameters observed in this research were consistent with those previously reported. The AUC, Cmax and tmax described in previous works with mirabegron ranged between 315–350 ng·h/mL, 16.1–39.2 ng/mL and 3–5 h, respectively, which are comparable to the results observed in this research [7,35,36,37]. Conversely, the t1/2 was lower than the 50 h indicated in the drug label [7]. However, a mirabegron t1/2 of 23–25 h and 40 h were also reported in other investigations [38,39]. It is likely that these differences are caused by the fact that the last sample in this research was extracted 72 h after drug intake, whereas a precise estimation of mirabegron t1/2 requires the collection of samples at least until 150 h after drug intake (i.e., three t1/2).
The lower exposure observed in females compared to males disappeared upon DW correction. Therefore, it is likely that these differences are caused by the lower weight of women compared to men, as shown in the demographic analysis, rather than by other sex differences. Actually, this higher exposure in women was previously documented in other studies with mirabegron [35,36,37]. If these differences are found to be clinically relevant in further research, the potential of adjusting dose to weight, rather than prescribing a unified dose, should be considered. Furthermore, a lower mirabegron exposure and a higher tmax were observed in fed volunteers (trial B) compared to those in fasting conditions (trials A and C). Consistently, it has been proposed that high-fat and low-fat meals reduce mirabegron AUC and Cmax and increase its tmax [7,40]. Since different transporters are involved in mirabegron influx into the lumen, such as the organic anion transporting polypeptide 2A1 (OATP2A1), the lower exposure and absorption observed may be attributable to the inhibition of these transporters by food [41,42]. The delayed mirabegron absorption in the presence of food might be explained by a delayed gastric emptying [36]. Despite these pharmacokinetic differences, safety and efficacy were achieved under both fed and fasting conditions in previous research [7]. Consequently, these differences are expected to be clinically irrelevant, and no dose adjustment ought to be required.
Mirabegron is extensively metabolized in the liver by CYP3A4 and CYP2D6 [43]. In this research, a higher t1/2 was observed in CYP2D6 IMs compared to NMs and UMs, and a progressive reduction in the AUC/DW, Cmax/DW and t1/2 and a progressive increase in CL/F when CYP2D6 activity augmented was observed. The lack of significance in exposure and CL might be caused by the reduced sample size available: only two UMs were included, with zero PMs. Indeed, to the best of our knowledge, differences in mirabegron exposure were previously only reported in CYP2D6 PMs, which showed higher drug exposure than NMs [18]. In accordance, individuals with a higher t1/2 (i.e., CYP2D6 IMs) presented a higher incidence of ADRs in this research. Thus, the augmented elimination rate in volunteers with higher CYP2D6 activity, along with their lower ADR incidence, suggests that CYP2D6 is involved in mirabegron metabolism and it might be a useful biomarker in the treatment with this drug. Nevertheless, a low contribution of CYP2D6 in mirabegron pharmacokinetics was proposed in previous research [44], which may also explain the lack of significance in CL and exposure. Consequently, further studies with increased sample size are required to study the impact of the CYP2D6 phenotype on the pharmacokinetics and safety of mirabegron.
On the other hand, no significant differences in mirabegron pharmacokinetics were observed based on the CYP3A4 phenotype, although it is proposed as its main metabolizing enzyme [45]. The exposure was higher in IMs than in NMs (see Supplementary Table S2), albeit not significantly, probably due to the reduced sample size; only two volunteers were CYP3A4 IMs and no PMs were found. Furthermore, the high inductility and inhibition of CYP3A4, and the different CYP3A4 expression levels between the sexes, may neutralize the impact of CYP3A4 genetic variation on its enzymatic activity [46,47]. Nevertheless, it should also be considered that other enzymes, such as CYP2D6 and UGTs, might play a more relevant role than CYP3A4 in mirabegron metabolism. Consequently, further research is needed to ascertain the relevance of CYP3A4 genetic variation in mirabegron pharmacokinetics.
Genetic variation in UGT1A3 and UGT1A4 was found to alter mirabegron pharmacokinetics in this study. A higher t1/2 was observed when the UGT1A4 rs2011425 T/G genotype was present, suggesting decreased elimination compared to the T/T genotype. This missense genetic variant was previously associated with both increased and decreased UGT1A4 glucuronidation activity in research with other drugs [48,49]. However, to the best of our knowledge, no previous associations between this genetic variant and differences in mirabegron pharmacokinetics or response are known. Nevertheless, other enzymes from the UGT1A gene cluster do participate in mirabegron metabolism, and it is known that enzymes from this family present substrate overlap [50]. Therefore, it could be proposed that the G allele in UGT1A4 rs2011425 causes a decreased mirabegron metabolization compared to the T allele, increasing the time needed for its elimination. However, a different study determined that the role of UGT1A4 in mirabegron metabolism is negligible [51]. Secondly, UGT1A3 rs2008584 was only associated with differences in tmax, which is a pharmacokinetic parameter more dependent on transporter activity, rather than that of hepatic metabolizing enzymes. In light of these considerations and the results obtained, it could be proposed that UGT1A4 rather than UGT1A3 might be involved in mirabegron metabolism, since genetic variation in UGT1A4 altered mirabegron t1/2. Nevertheless, further research on this association and on the impact of genetic variation in UGT1A4 activity should be performed. Lastly, UGT2B7 and UGT1A8 genetic variation was not found to impact mirabegron pharmacokinetics (see Supplementary Table S2). However, genetic variation in these genes and in their gene family is not well characterized (e.g., the impact of known genetic variants on enzyme activity is not clear and no haplotypes and phenotypes have been identified and defined). Therefore, a better structural and functional characterization of the genetic variation in these UGT genes is needed prior to evaluating their relevance as pharmacogenetic biomarkers.
The gene NAT2 codifies for the enzyme N-acetyltransferase 2, which is involved in the metabolism of diverse drugs [31,52,53]. A lower mirabegron Cmax/DW in individuals with lower NAT2 activity was observed. These results are contrary to previous results in the literature with other drugs, in which lower NAT2 activity was associated with reduced drug elimination and therefore higher exposure [54]. The observed differences might be caused by the different distribution of RAs, IAs and SAs in the three clinical trials. Actually, individuals with higher NAT2 activity were only present in clinical trials A and C, in which drug exposure was increased, compared to clinical trial B (Trial B: 0% RAs, 34.48% IAs and 65.52% SAs versus trials A + C: 14.29% RAs, 48.87% IAs and 36.73% SAs, p = 0.007). The lack of significant differences in the multivariate analysis supports this hypothesis. Nevertheless, the NAT2 role in mirabegron pharmacokinetics should be investigated in further research once the impact of each allele is clearly informed in the PharmVar database [55].
The solute carriers (SLC) transporters SLC6A2, SLC19A1 and SLCO1B1 were associated with differences in mirabegron pharmacokinetics and safety, respectively. Volunteers with the SLC6A2 rs12708954 C/C genotype presented lower exposure compared to those with the C/A genotype. SLC6A2 codifies for a noradrenaline transporter (NET) and NET transporters are involved in mirabegron uptake into the sympathetic nerve terminals [56]. Indeed, the SLC6A2 rs12708954 C/C genotype was previously associated with higher response and ADR incidence in other drugs that are also NET substrates, such as atomoxetine [57]. Despite the SLC6A2 expression in different regions of the gastrointestinal tract [58], it is unknown whether it has a role in drug pharmacokinetics. For these reasons and considering the reduced sample size available, further research is required to determine the relevance of this transporter to mirabegron pharmacokinetics. On the other side, SLC19A1 is a folate transporter with no previous associations with mirabegron pharmacokinetics [7,41,59]. Considering that the differences were observed in the heterozygous genotype, and their disappearance in the multivariate analysis, this result ought to be regarded as spurious.
Lastly, SLCO1B1 codifies for the OATP1B1, which is a liver specific transporter that contributes to the hepatic uptake of many drugs [60]. A higher ADR incidence in volunteers with lower SLCO1B1 activity was observed. To the best of our knowledge, no studies have proposed that OATP1B1 is involved in mirabegron transit, but previous research demonstrated that OATP1A2, a transporter of the same family, was capable of transporting mirabegron [41]. In research with other drugs, such as statins, a decreased SLCO1B1 activity was related to higher exposure, potentially reaching toxic concentrations that led to the appearance of ADRs [12]. However, differences in ADRs are expected to be accompanied by pharmacokinetic differences, which did not occur. Notwithstanding, when data from each clinical trial were analyzed individually, a trend towards a higher AUC/DW was observed in SLCO1B1 DFs (439.35 ± 125.50 kg·ng·h/mL·mg, n = 13) compared to NFs (364.01 ± 66.01 kg·ng·h/mL·mg, n = 15, puv = 0.080) in clinical trial B (fed conditions), but not in clinical trials A and C (fasting conditions). Thus, it could be proposed that mirabegron hepatic uptake is modified in the presence of food, conferring a more relevant role to SLCO1B1 under these conditions. Certain foods and food additives are known to inhibit different OATP isoforms [42], such as OATP2B1 [61], but, to our knowledge, no correlation between food consumption and SLCO1B1 inhibition has been proposed. Further research on the role of SLCO1B1 in mirabegron transport and its interaction with food is required.
This study presents some limitations. First, the sample size was arbitrary and limited, as it was determined by the amount of data available from the clinical trials. Increased sample size in further research would be advisable to ensure the not inconsiderable presence of low-prevalence genotypes and phenotypes (e.g., CYP2D6 UMs and PMS or CYP3A4 PMs), which would improve the statistical power of the study. Secondly, the genetic variants analyzed were selected based on previous knowledge and on their role in drug pharmacokinetics. Therefore, there is risk of missing unknown variants present in the genome that could alter mirabegron pharmacokinetics. Moreover, the subjects were healthy volunteers, and the majority of them received a single mirabegron dose, which reduced population variability, prevented efficacy assessment and hampered the safety analysis due to the low incidence of ADRs. Nevertheless, this research also offered advantages. It was performed in a controlled environment that allowed avoiding confounding factors, e.g., smoking, comorbidities or concomitant treatment, which are present in investigation with patients. Additionally, correct compliance with medication intake, homogeneous feeding conditions and controlled and precise sampling were ensured.

5. Conclusions

Having a reduced CYP2D6 activity was associated with lower elimination rate and higher incidence of ADRs, a finding which paves the way towards the use of this biomarker in mirabegron treatment. The metabolizing enzyme UGT1A4 and the transporters SLC6A2 and SLCO1B1 might be involved in mirabegron metabolism and absorption and distribution, respectively. The impact of CYP3A4 and UGT2B7 genetic variation in mirabegron pharmacokinetics requires further research.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics16081077/s1, Table S1. Analyzed genetic variants. Table S2. No statistically significant associations between genetic variation and pharmacokinetic parameters.

Author Contributions

Conceptualization, P.S.-C., P.C.-R. and F.A.-S.; Methodology, P.S.-C., P.C.-R. and A.R.-L.; Software, P.S.-C. and P.C.-R.; Validation, D.O., G.M.-A. and S.L.-B.; Formal Analysis, P.S.-C. and P.C.-R.; Investigation, P.S.-C. and P.C.-R.; Resources, D.O. and F.A.-S.; Data Curation, E.G.-I., M.R., S.M.-V., S.L.-B. and G.M.-A.; Writing—Original Draft, P.S.-C. and P.C.-R.; Writing—Review and Editing, A.R.-L., E.G.-I., M.R., S.M.-V., S.L.-B. and G.M.-A.; Supervision, F.A.-S. and D.O.; Project Administration, D.O. and S.M.-V. All authors have read and agreed to the published version of the manuscript.

Funding

P.S.-C. is financed by the FPI-UAM-2021 predoctoral fellowship (UAM). A.R.-L. contract is financed by Programa Investigo (NextGenerationEU funds of the Recovery and Resilience Facility), fellowship number 2022-C23.I01.P03.S0020-0000031. E.G.-I. is financed by the predoctoral fellowship PIPF-2022/SAL-GL-25946.

Institutional Review Board Statement

The data used in this research were obtained from three bioequivalence clinical trials conducted at the Clinical Trials Unit of the Hospital Universitario de La Princesa (UECHUP) (EUDRA-CT: 2020-004710-35, 2021-000283-29 and 2021-000285-15). The CT protocols were revised and approved by the Hospital’s Research Ethics Committee and by the Spanish Drugs Agency (AEMPS). CT and data handling were performed complying with the Good Clinical Practice guidelines [13], the Declaration of Helsinki [14] and the current Spanish legislation regarding clinical trials.

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the clinical trials and in the pharmacogenetic study.

Data Availability Statement

Data belong to the clinical trials’ sponsors and may be accessible upon reasonable request to the corresponding authors. The data are not publicly available due to the instruction of clinical trial’s sponsors.

Conflicts of Interest

F.A.-S. and D.O. have been consultants or investigators in clinical trials sponsored by the following pharmaceutical companies: Abbott, Alter, Aptatargets, Chemo, Cinfa, FAES, Farmalíder, Ferrer, GlaxoSmithKline, Galenicum, Gilead, Italfarmaco, Janssen-Cilag, Kern, Moderna, MSD, Normon, Novartis, Servier, Silver Pharma, Teva and Zambon. The remaining authors declare no conflicts of interest.

References

  1. Hutchinson, A.; Nesbitt, A.; Joshi, A.; Clubb, A.; Perera, M. Overactive Bladder Syndrome: Management and Treatment Options. Aust. J. Gen. Pract. 2020, 49, 593–598. [Google Scholar] [CrossRef] [PubMed]
  2. Pautz, A.; Michel, M.C. Sex and Gender Differences in the Pharmacology of the Overactive Urinary Bladder. In Sex and Gender Effects in Pharmacology; Tsirka, S.E., Acosta-Martinez, M., Eds.; Handbook of Experimental Pharmacology; Springer International Publishing: Cham, Switzerland, 2023; pp. 57–74. Volume 282, ISBN 978-3-031-42647-6. [Google Scholar]
  3. Yang, C.-F.; Huang, C.-Y.; Wang, S.-Y.; Chang, S.-R. Prevalence of and Associated Factors for Overactive Bladder Subtypes in Middle-Aged Women: A Cross-Sectional Study. Medicina 2022, 58, 383. [Google Scholar] [CrossRef] [PubMed]
  4. El-Zawahry, A. Combination Pharmacotherapy for Treatment of Overactive Bladder (OAB). Curr. Urol. Rep. 2019, 20, 33. [Google Scholar] [CrossRef]
  5. European Medicines Agency Betmiga. Last Update: 15th October 2015. Available online: https://www.ema.europa.eu/en/medicines/human/EPAR/betmiga (accessed on 12 April 2024).
  6. Peyronnet, B.; Mironska, E.; Chapple, C.; Cardozo, L.; Oelke, M.; Dmochowski, R.; Amarenco, G.; Gamé, X.; Kirby, R.; Van Der Aa, F.; et al. A Comprehensive Review of Overactive Bladder Pathophysiology: On the Way to Tailored Treatment. Eur. Urol. 2019, 75, 988–1000. [Google Scholar] [CrossRef] [PubMed]
  7. Spanish Drug Agency CIMA. FICHA TECNICA BETMIGA 50mg Comprimidos de Liberacion Prolongada. Available online: https://cima.aemps.es/cima/dochtml/ft/112809010/FT_112809010 (accessed on 12 April 2024).
  8. Konishi, K.; Tenmizu, D.; Takusagawa, S. Identification of Uridine 5′-Diphosphate-Glucuronosyltransferases Responsible for the Glucuronidation of Mirabegron, a Potent and Selective Β3-Adrenoceptor Agonist, in Human Liver Microsomes. Eur. J. Drug Metab. Pharmacokinet. 2018, 43, 301–309. [Google Scholar] [CrossRef] [PubMed]
  9. Cacabelos, R.; Naidoo, V.; Corzo, L.; Cacabelos, N.; Carril, J.C. Genophenotypic Factors and Pharmacogenomics in Adverse Drug Reactions. Int. J. Mol. Sci. 2021, 22, 13302. [Google Scholar] [CrossRef] [PubMed]
  10. Cecchin, E.; Stocco, G. Pharmacogenomics and Personalized Medicine. Genes 2020, 11, 679. [Google Scholar] [CrossRef] [PubMed]
  11. Lima, J.J.; Thomas, C.D.; Barbarino, J.; Desta, Z.; Van Driest, S.L.; El Rouby, N.; Johnson, J.A.; Cavallari, L.H.; Shakhnovich, V.; Thacker, D.L.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2C19 and Proton Pump Inhibitor Dosing. Clin. Pharmacol. Ther. 2021, 109, 1417–1423. [Google Scholar] [CrossRef] [PubMed]
  12. Cooper-DeHoff, R.M.; Niemi, M.; Ramsey, L.B.; Luzum, J.A.; Tarkiainen, E.K.; Straka, R.J.; Gong, L.; Tuteja, S.; Wilke, R.A.; Wadelius, M.; et al. The Clinical Pharmacogenetics Implementation Consortium Guideline for SLCO1B1, ABCG2, and CYP2C9 Genotypes and Statin-Associated Musculoskeletal Symptoms. Clin. Pharmacol. Ther. 2022, 111, 1007–1021. [Google Scholar] [CrossRef]
  13. Desta, Z.; Gammal, R.S.; Gong, L.; Whirl-Carrillo, M.; Gaur, A.H.; Sukasem, C.; Hockings, J.; Myers, A.; Swart, M.; Tyndale, R.F.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2B6 and Efavirenz-Containing Antiretroviral Therapy. Clin. Pharmacol. Ther. 2019, 106, 726–733. [Google Scholar] [CrossRef]
  14. Goetz, M.P.; Sangkuhl, K.; Guchelaar, H.-J.; Schwab, M.; Province, M.; Whirl-Carrillo, M.; Symmans, W.F.; McLeod, H.L.; Ratain, M.J.; Zembutsu, H.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2D6 and Tamoxifen Therapy. Clin. Pharmacol. Ther. 2018, 103, 770–777. [Google Scholar] [CrossRef]
  15. Birdwell, K.; Decker, B.; Barbarino, J.; Peterson, J.; Stein, C.; Sadee, W.; Wang, D.; Vinks, A.; He, Y.; Swen, J.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guidelines for CYP3A5 Genotype and Tacrolimus Dosing. Clin. Pharmacol. Ther. 2015, 98, 19–24. [Google Scholar] [CrossRef] [PubMed]
  16. Relling, M.V.; Schwab, M.; Whirl-Carrillo, M.; Suarez-Kurtz, G.; Pui, C.-H.; Stein, C.M.; Moyer, A.M.; Evans, W.E.; Klein, T.E.; Antillon-Klussmann, F.G.; et al. Clinical Pharmacogenetics Implementation Consortium Guideline for Thiopurine Dosing Based on TPMT and NUDT15 Genotypes: 2018 Update. Clin. Pharmacol. Ther. 2019, 105, 1095–1105. [Google Scholar] [CrossRef]
  17. Gammal, R.; Court, M.; Haidar, C.; Iwuchukwu, O.; Gaur, A.; Alvarellos, M.; Guillemette, C.; Lennox, J.; Whirl-Carrillo, M.; Brummel, S.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for UGT1A1 and Atazanavir Prescribing. Clin. Pharmacol. Ther. 2016, 99, 363–369. [Google Scholar] [CrossRef]
  18. PharmGKB Mirabegron. Available online: https://www.pharmgkb.org/chemical/PA166177513 (accessed on 22 May 2024).
  19. Zubiaur, P.; Mejía-Abril, G.; Navares-Gómez, M.; Villapalos-García, G.; Soria-Chacartegui, P.; Saiz-Rodríguez, M.; Ochoa, D.; Abad-Santos, F. PriME-PGx: La Princesa University Hospital Multidisciplinary Initiative for the Implementation of Pharmacogenetics. J. Clin. Med. 2021, 10, 3772. [Google Scholar] [CrossRef]
  20. Vijayananthan, A.; Nawawi, O. The Importance of Good Clinical Practice Guidelines and Its Role in Clinical Trials. Biomed. Imaging Interv. J. 2008, 4, e5. [Google Scholar] [CrossRef] [PubMed]
  21. World Medical Association Declaration of Helsinki Ethical Principles for Medical Research Involving Human Subjects. Last Update: 19th October 2013. Available online: https://jamanetwork.com/journals/jama/fullarticle/1760318 (accessed on 20 May 2024).
  22. Chin, R.Y.; Lee, B.Y. Principles and Practice of Clinical Trial Medicine; Elsevier: Amsterdam, The Netherlands; Academic Press: Cambridge, MA, USA, 2008; ISBN 978-0-12-373695-6. [Google Scholar]
  23. European Medicines Agency Guideline on the Investigation of Bioequivalence. Available online: https://www.ema.europa.eu/en/documents/scientific-guideline/guideline-investigation-bioequivalence-rev1_en.pdf (accessed on 20 May 2024).
  24. Aguirre, C.; García, M. Causality assessment in adverse drug reaction reports. Algorithm of the Spanish Pharmacovigilance System. Med. Clínica 2016, 147, 461–464. [Google Scholar] [CrossRef] [PubMed]
  25. National Institute of Health RefSeq: NCBI Reference Sequence Database. Available online: https://www.ncbi.nlm.nih.gov/refseq/ (accessed on 26 February 2024).
  26. Soria-Chacartegui, P.; Zubiaur, P.; Ochoa, D.; Villapalos-García, G.; Román, M.; Matas, M.; Figueiredo-Tor, L.; Mejía-Abril, G.; Calleja, S.; De Miguel, A.; et al. Genetic Variation in CYP2D6 and SLC22A1 Affects Amlodipine Pharmacokinetics and Safety. Pharmaceutics 2023, 15, 404. [Google Scholar] [CrossRef]
  27. Clinical Pharmacogenetics Implementation Consortium CPIC Website. Available online: https://cpicpgx.org/ (accessed on 26 February 2020).
  28. PharmGKB PGx Gene-Specific Information Tables. Available online: https://www.pharmgkb.org/page/pgxGeneRef (accessed on 26 February 2024).
  29. DPWG Annotation of DPWG Guideline for Quetiapine and CYP3A4. 2021. Available online: https://www.pharmgkb.org/guidelineAnnotation/PA166265421 (accessed on 26 February 2024).
  30. Campodónico, D.M.; Zubiaur, P.; Soria-Chacartegui, P.; Casajús, A.; Villapalos-García, G.; Navares-Gómez, M.; Gómez-Fernández, A.; Parra-Garcés, R.; Mejía-Abril, G.; Román, M.; et al. CYP2C8 *3 and *4 Define CYP2C8 Phenotype: An Approach with the Substrate Cinitapride. Clin. Transl. Sci. 2022, 15, 2613–2624. [Google Scholar] [CrossRef]
  31. Villapalos-García, G.; Zubiaur, P.; Ochoa, D.; Soria-Chacartegui, P.; Navares-Gómez, M.; Matas, M.; Mejía-Abril, G.; Casajús-Rey, A.; Campodónico, D.; Román, M.; et al. NAT2 Phenotype Alters Pharmacokinetics of Rivaroxaban in Healthy Volunteers. Biomed. Pharmacother. 2023, 165, 115058. [Google Scholar] [CrossRef]
  32. Humburg, J. Die Urininkontinenz der Frau: Einführung in die Diagnostik und Therapie. Ther. Umsch. 2018, 75, 535–540. [Google Scholar] [CrossRef] [PubMed]
  33. Scarneciu, I.; Lupu, S.; Bratu, O.; Teodorescu, A.; Maxim, L.; Brinza, A.; Laculiceanu, A.; Rotaru, R.; Lupu, A.-M.; Scarneciu, C. Overactive Bladder: A Review and Update. Exp. Ther. Med. 2021, 22, 1444. [Google Scholar] [CrossRef]
  34. Deeks, E.D. Mirabegron: A Review in Overactive Bladder Syndrome. Drugs 2018, 78, 833–844. [Google Scholar] [CrossRef]
  35. Krauwinkel, W.; Van Dijk, J.; Schaddelee, M.; Eltink, C.; Meijer, J.; Strabach, G.; Van Marle, S.; Kerbusch, V.; Van Gelderen, M. Pharmacokinetic Properties of Mirabegron, a Β3-Adrenoceptor Agonist: Results from Two Phase I, Randomized, Multiple-Dose Studies in Healthy Young and Elderly Men and Women. Clin. Ther. 2012, 34, 2144–2160. [Google Scholar] [CrossRef]
  36. Iitsuka, H.; Tokuno, T.; Amada, Y.; Matsushima, H.; Katashima, M.; Sawamoto, T.; Takusagawa, S.; Van Gelderen, M.; Tanaka, T.; Miyahara, H. Pharmacokinetics of Mirabegron, a Β3-Adrenoceptor Agonist for Treatment of Overactive Bladder, in Healthy Japanese Male Subjects: Results from Single- and Multiple-Dose Studies. Clin. Drug Investig. 2014, 34, 27–35. [Google Scholar] [CrossRef] [PubMed]
  37. Eltink, C.; Lee, J.; Schaddelee, M.; Zhang, W.; Kerbusch, V.; Meijer, J.; Marle, S.V.; Grunenberg, N.; Kowalski, D.; Drogendijk, T.; et al. Single Dose Pharmacokinetics and Absolute Bioavailability of Mirabegron, a β3-Adrenoceptor Agonist for Treatment of Overactive Bladder. Int. J. Clin. Pharmacol. Ther. 2012, 50, 838–850. [Google Scholar] [CrossRef] [PubMed]
  38. Warren, K.; Burden, H.; Abrams, P. Mirabegron in Overactive Bladder Patients: Efficacy Review and Update on Drug Safety. Ther. Adv. Drug Saf. 2016, 7, 204–216. [Google Scholar] [CrossRef]
  39. Tyagi, P.; Tyagi, V.; Chancellor, M. Mirabegron: A Safety Review. Expert. Opin. Drug Saf. 2011, 10, 287–294. [Google Scholar] [CrossRef] [PubMed]
  40. Paśko, P.; Rodacki, T.; Domagała-Rodacka, R.; Owczarek, D. A Short Review of Drug–Food Interactions of Medicines Treating Overactive Bladder Syndrome. Int. J. Clin. Pharm. 2016, 38, 1350–1356. [Google Scholar] [CrossRef]
  41. Takusagawa, S.; Ushigome, F.; Nemoto, H.; Takahashi, Y.; Li, Q.; Kerbusch, V.; Miyashita, A.; Iwatsubo, T.; Usui, T. Intestinal Absorption Mechanism of Mirabegron, a Potent and Selective β3-Adrenoceptor Agonist: Involvement of Human Efflux and/or Influx Transport Systems. Mol. Pharm. 2013, 10, 1783–1794. [Google Scholar] [CrossRef]
  42. Ali, Y.; Shams, T.; Wang, K.; Cheng, Z.; Li, Y.; Shu, W.; Bao, X.; Zhu, L.; Murray, M.; Zhou, F. The Involvement of Human Organic Anion Transporting Polypeptides (OATPs) in Drug-Herb/Food Interactions. Chin. Med. 2020, 15, 71. [Google Scholar] [CrossRef] [PubMed]
  43. Konishi, K.; Minematsu, T.; Nagasaka, Y.; Tabata, K. Physiologically-Based Pharmacokinetic Modeling for Mirabegron: A Multi-Elimination Pathway Mediated by Cytochrome P450 3A4, Uridine 5′-Diphosphate-Glucuronosyltransferase 2B7, and Butyrylcholinesterase. Xenobiotica 2019, 49, 912–921. [Google Scholar] [CrossRef] [PubMed]
  44. Lee, J.; Moy, S.; Meijer, J.; Krauwinkel, W.; Sawamoto, T.; Kerbusch, V.; Kowalski, D.; Roy, M.; Marion, A.; Takusagawa, S.; et al. Role of Cytochrome P450 Isoenzymes 3A and 2D6 in the In Vivo Metabolism of Mirabegron, a Β3-Adrenoceptor Agonist. Clin. Drug Investig. 2013, 33, 429–440. [Google Scholar] [CrossRef] [PubMed]
  45. Takusagawa, S.; Miyashita, A.; Iwatsubo, T.; Usui, T. In Vitro Inhibition and Induction of Human Cytochrome P450 Enzymes by Mirabegron, a Potent and Selective Β3-Adrenoceptor Agonist. Xenobiotica 2012, 42, 1187–1196. [Google Scholar] [CrossRef] [PubMed]
  46. Minegishi, G.; Kobayashi, Y.; Fujikura, M.; Sano, A.; Kazuki, Y.; Kobayashi, K. Induction of Hepatic CYP3A4 Expression by Cholesterol and Cholic Acid: Alterations of Gene Expression, Microsomal Activity, and Pharmacokinetics. Pharmacol. Res. Perspect. 2024, 12, e1197. [Google Scholar] [CrossRef] [PubMed]
  47. Chetty, M.; Mattison, D.; Rostami-Hodjegan, A. Sex Differences in the Clearance of CYP3A4 Substrates: Exploring Possible Reasons for the Substrate Dependency and Lack of Consensus. Curr. Drug Metab. 2012, 13, 778–786. [Google Scholar] [CrossRef] [PubMed]
  48. Pardiñas, A.F.; Nalmpanti, M.; Pocklington, A.J.; Legge, S.E.; Medway, C.; King, A.; Jansen, J.; Helthuis, M.; Zammit, S.; MacCabe, J.; et al. Pharmacogenomic Variants and Drug Interactions Identified through the Genetic Analysis of Clozapine Metabolism. Am. J. Psychiatry 2019, 176, 477–486. [Google Scholar] [CrossRef] [PubMed]
  49. Ghotbi, R.; Mannheimer, B.; Aklillu, E.; Suda, A.; Bertilsson, L.; Eliasson, E.; Ösby, U. Carriers of the UGT1A4 142T>G Gene Variant Are Predisposed to Reduced Olanzapine Exposure—An Impact Similar to Male Gender or Smoking in Schizophrenic Patients. Eur. J. Clin. Pharmacol. 2010, 66, 465–474. [Google Scholar] [CrossRef]
  50. Fujiwara, R.; Yokoi, T.; Nakajima, M. Structure and Protein–Protein Interactions of Human UDP-Glucuronosyltransferases. Front. Pharmacol. 2016, 7, 388. [Google Scholar] [CrossRef]
  51. Sodhi, J.K.; Benet, L.Z. The Necessity of Using Changes in Absorption Time to Implicate Intestinal Transporter Involvement in Oral Drug-Drug Interactions. AAPS J. 2020, 22, 111. [Google Scholar] [CrossRef]
  52. Zubiaur, P.; Figueiredo-Tor, L.; Villapalos-García, G.; Soria-Chacartegui, P.; Navares-Gómez, M.; Novalbos, J.; Matas, M.; Calleja, S.; Mejía-Abril, G.; Román, M.; et al. Association between CYP2C19 and CYP2B6 Phenotypes and the Pharmacokinetics and Safety of Diazepam. Biomed. Pharmacother. 2022, 155, 113747. [Google Scholar] [CrossRef] [PubMed]
  53. Thomas, L.; Raju, A.P.; Chaithra, M.S.S.; Varma, M.; Saravu, K.; Banerjee, M.; Sv, C.S.; Mallayasamy, S.; Rao, M. Influence of N-Acetyltransferase 2 (NAT2) Genotype/Single Nucleotide Polymorphisms on Clearance of Isoniazid in Tuberculosis Patients: A Systematic Review of Population Pharmacokinetic Models. Eur. J. Clin. Pharmacol. 2022, 78, 1535–1553. [Google Scholar] [CrossRef] [PubMed]
  54. PharmGKB NAT2. Available online: https://www.pharmgkb.org/gene/PA18 (accessed on 9 June 2024).
  55. PharmVar NAT2. Available online: https://www.pharmvar.org/gene/NAT2 (accessed on 9 June 2024).
  56. Mo, W.; Michel, M.C.; Lee, X.W.; Kaumann, A.J.; Molenaar, P. The β3-adrenoceptor Agonist Mirabegron Increases Human Atrial Force through β1-adrenoceptors: An Indirect Mechanism? Br. J. Pharmacol. 2017, 174, 2706–2715. [Google Scholar] [CrossRef]
  57. Gul, M.K.; Sener, E.F.; Onal, M.G.; Demirci, E. Role of the Norepinephrine Transporter Polymorphisms in Atomoxetine Treatment: From Response to Side Effects in Children with ADHD. J. Psychopharmacol. 2022, 36, 715–722. [Google Scholar] [CrossRef] [PubMed]
  58. Bgee database SLC6A2 Expression in Humans. Available online: https://www.bgee.org/gene/ENSG00000103546 (accessed on 31 May 2024).
  59. Groen-Wijnberg, M.; Van Dijk, J.; Krauwinkel, W.; Kerbusch, V.; Meijer, J.; Tretter, R.; Zhang, W.; Van Gelderen, M. Pharmacokinetic Interactions between Mirabegron and Metformin, Warfarin, Digoxin or Combined Oral Contraceptives. Eur. J. Drug Metab. Pharmacokinet. 2017, 42, 417–429. [Google Scholar] [CrossRef] [PubMed]
  60. Ramsey, L.B.; Gong, L.; Lee, S.; Wagner, J.B.; Zhou, X.; Sangkuhl, K.; Adams, S.M.; Straka, R.J.; Empey, P.E.; Boone, E.C.; et al. PHARMVAR GENEFOCUS: SLCO1B1. Clin. Pharmacol. Ther. 2022, 113, 782–793. [Google Scholar] [CrossRef]
  61. Tikkanen, A.; Pierrot, E.; Deng, F.; Sánchez, V.B.; Hagström, M.; Koenderink, J.B.; Kidron, H. Food Additives as Inhibitors of Intestinal Drug Transporter OATP2B1. Mol. Pharm. 2020, 17, 3748–3758. [Google Scholar] [CrossRef]
Table 1. Characteristics of the clinical trials included in the pharmacogenetic study.
Table 1. Characteristics of the clinical trials included in the pharmacogenetic study.
Clinical TrialEUDRA-CTConditionsDoseSample Size * (BE Study)Sample Size * (PG Study)
A2020-004710-35FastingSingle,
50 mg
36 (19/17)22 (12/10)
B2021-000283-29FedSingle,
50 mg
36 (18/18)30 (15/15)
C2021-000285-15FastingMultiple, 50 mg every 24 h36 (19/17)27 (11/16)
* Sample size is shown as total number (number of males/number of females). BE: bioequivalence, PG: pharmacogenetic. In each clinical trial, both formulations (T and R) were administered to every volunteer under the conditions and at the dose specified.
Table 2. Demographic characteristics according to sex, biogeographic origin and clinical trial.
Table 2. Demographic characteristics according to sex, biogeographic origin and clinical trial.
nAge
Years
Height
m
Weight
kg
BMI
kg/m2
Sex
Male3828.50 (24.00–34.25)1.76 (0.06)75.96 (11.64)24.85 (22.30–27.20)
Female4130.00 (24.00–33.50)1.62 (0.06) *63.92 (10.21) *24.33 (21.88–27.37)
Origin
European1826.00 (23.75–31.75)1.73 (0.09)76.81 (14.61)26.24 (22.40–27.82)
Latin American6131.00 (25.00–34.50)1.67 (0.09) *67.61 (10.99) *24.33 (21.84–26.55)
Trial
A2225.00 (23.00–31.50)1.70 (0.09)68.30 (11.71)24.49 (21.04–26.38)
B3028.50 (24.00–43.00)1.68 (0.09)68.63 (12.32)23.66 (21.01–27.53)
C2732.00 (26.00–33.00)1.68 (0.10)72.06 (13.21)25.64 (23.15–27.58)
Total7930.00 (24.00–34.00)1.68 (0.09)69.71 (12.43)24.58 (22.12–27.25)
Results are shown as mean (standard deviation) for normally distributed variables and as median (interquartile range) for not normally distributed variables. BMI: body mass index; *: p < 0.01 compared to males or Europeans.
Table 3. Pharmacokinetic data according to sex, biogeographic origin and clinical trial.
Table 3. Pharmacokinetic data according to sex, biogeographic origin and clinical trial.
nAUC/DW
kg·ng·h/mL·mg
Cmax/DW
kg·ng/mL·mg
tmax
h
t1/2
h
CL/F
L/kg·h
Sex
Female41577.02 (189.54)60.84 (29.19–74.58)4.50 (3.94–5.31)26.50 (8.21)1.81 (1.44–2.41)
Male38529.27 (177.04)48.31 (30.22–67.24)4.44 (3.76–5.13)26.34 (5.22)2.03 (1.64–2.61)
Origin
European18582.84 (202.00)54.80 (40.65–72.87)4.69 (4.19–4.91)25.27 (5.16)1.75 (1.55–2.33)
Latin American61545.56 (179.29)48.74 (27.06–71.30)4.41 (3.75–5.31)26.75 (7.33)1.92 (1.52–2.57)
Trial
A22633.34 (149.93)58.58 (50.19–67.74)4.31 (3.97–5.00)27.63 (3.89)1.77 (1.40–1.94)
B30397.84 (106.91) *24.85 (19.40–33.98) *5.13 (4.41–6.25) *30.62 (7.04)2.72 (2.21–3.35) *
C27663.03 (157.77)73.88 (63.14–93.71)4.25 (3.63–4.50)20.77 (4.51) $1.64 (1.26–1.81)
Total79554.06 (184.03)50.72 (29.50–72.54)4.50 (3.88–5.13)26.42 (6.90)1.91 (1.52–2.51)
Results are shown as mean (standard deviation) for normally distributed variables and as median (interquartile range) for not normally distributed variables. AUC/DW: dose/weight corrected area under the concentration-time curve; Cmax/DW: dose/weight corrected maximum drug plasma concentration; tmax: time to reach Cmax; t1/2: elimination half-life; CL/F: drug clearance adjusted by bioavailability. *: puv < 0.05 compared to clinical trials A and C; $: puv < 0.05 compared to clinical trials A and B; underlined: pmv < 0.05.
Table 4. Associations between genotypes or phenotypes and pharmacokinetic parameters.
Table 4. Associations between genotypes or phenotypes and pharmacokinetic parameters.
GeneGenotype/PhenotypenAUC/DWCmax/DWtmaxt1/2CL/F
kg·ng·h/mL·mgkg·ng/mL·mghhL/kg·h
CYP2D6PhenotypeUM2451.77 (231.42)45.37 [19.63–71.11]4.94 [4.75–5.13]20.42 (3.90)2.60 [1.64–3.55]
NM55547.97 (199.06)48.74 (28.59–72.57)4.38 (3.75–5.13)25.80 (6.01)1.94 (1.45–2.71)
IM18585.61 (188.13)50.24 (36.42–67.93)4.88 (3.97–5.78)29.97 (8.68) *1.73 (1.50–2.70)
UGT1A3rs2008584A/A16548.38 (162.14)58.17 (31.72–83.83)4.06 (3.38–4.48) *24.35 (6.62)1.84 (1.67–2.49)
A/G34564.37 (212.65)47.89 (25.47–72.54)4.69 (3.94–5.66)25.70 (5.44)1.84 (1.33–2.74)
G/G26537.28 (170.40)48.67 (28.80–67.93)4.69 (3.88–5.16)29.03 (8.39)1.92 (1.64–2.46)
UGT1A4rs2011425 T/T64541.38 (187.75)48.49 (28.39–73.09)4.50 (3.75–5.25)25.28 (6.67)1.93 (1.55–2.69)
T/G14587.99 (147.20)57.32 (30.52–66.27)4.31 (3.97–4.88)31.14 (5.96) *1.75 (1.52–1.98)
NAT2PhenotypeRA7602.63 (159.64)60.84 (48.74–65.11)4.00 (3.26–4.25)27.80 (7.77)1.75 (1.68–1.93)
IA34589.33 (196.52)61.37 (36.10–76.70)4.38 (3.75–5.16)24.61 (4.98)1.81 (1.33–2.38)
SA37516.32 (173.40)40.91 (24.85–68.97) *4.75 (4.06–5.56)27.74 (8.04)2.19 (1.66–2.67)
SLC6A2rs12708954C/C56522.63 (182.61) *48.31 (24.23–70.35)4.45 (3.75–5.09)26.56 (7.34)1.93 (1.65–2.81) *
C/A21648.45 (162.98)65.11 (41.05–78.25)4.63 (3.94–5.56)26.29 (6.01)1.67 (1.31–1.99)
A/A2442.73 (95.66)47.11 [36.51–57.72]5.06 [4.88–5.25] 23.94 (2.07)2.35 [1.98–2.71]
SLC19A1rs1051266A/A16498.20 (190.51)37.45 (20.53–53.46) *4.50 (4.03–6.25)27.92 (5.87)2.31 (1.50–3.26)
A/G34598.00 (181.87)60.74 (41.65–74.22)4.33 (3.75–5.03)25.72 (5.63)1.81 (1.43–2.21)
G/G27537.43 (182.46)48.18 (23.61–71.11)4.50 (3.76–5.13)26.12 (8.87)1.92 (1.61–2.74)
Results are shown as mean (standard deviation) for normally distributed variables and median (interquartile range) for not normally distributed variables. For categories with two volunteers, the two values in square brackets are shown instead of the interquartile range. AUC/DW: dose/weight corrected area under the concentration-time curve; Cmax/DW: dose/weight corrected maximum drug plasma concentration; tmax: time to reach Cmax; t1/2: elimination half-life; UM: ultrarapid metabolizers; NM: normal metabolizers; IM: intermediate metabolizers; RA: rapid acetilators, IA: intermediate acetilators, SA: slow acetilators. *: puv < 0.05 compared to CYP2D6 UMs + NMs, UGT1A3 rs2008584 A/G genotype, UGT1A4 rs2011425 T/T genotype, NAT2 RAs + IAs, SLC6A2 rs12708954 C/A genotype and SLC19A1 rs1051266 A/G genotype; underlined: pmv < 0.05.
Table 5. AUC and Cmax in volunteers with or without ADRs.
Table 5. AUC and Cmax in volunteers with or without ADRs.
Clinical TrialADRsnAUC
ng·h/mL
Cmax
ng/mL
Ano20470.43 (108.04)45.42 (34.53–51.28)
yes2449.80 (19.61)44.90 [43.71- 46.09]
Bno24282.60 (85.10)16.81 (12.28–24.08)
yes6359.96 (92.85)20.28 (12.15–24.60)
Cno25476.37 (161.56)52.75 (43.61–63.79)
yes2505.72 (22.20)54.95 [51.33–58.57]
Totalno69407.25 (152.77)26.21 (17.99–47.39)
yes10407.08 (94.52)40.52 (20.99–51.82)
Results are shown as mean (standard deviation) for normally distributed variables and median (interquartile range) for not normally distributed variables. For categories with two volunteers, the two values in square brackets are shown instead of the interquartile range. ADR: adverse reactions; AUC: area under the concentration-time curve, Cmax: maximum drug plasma concentration.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Soria-Chacartegui, P.; Cendoya-Ramiro, P.; González-Iglesias, E.; Martín-Vílchez, S.; Rodríguez-Lopez, A.; Mejía-Abril, G.; Román, M.; Luquero-Bueno, S.; Ochoa, D.; Abad-Santos, F. Genetic Variation in CYP2D6, UGT1A4, SLC6A2 and SLCO1B1 Alters the Pharmacokinetics and Safety of Mirabegron. Pharmaceutics 2024, 16, 1077. https://doi.org/10.3390/pharmaceutics16081077

AMA Style

Soria-Chacartegui P, Cendoya-Ramiro P, González-Iglesias E, Martín-Vílchez S, Rodríguez-Lopez A, Mejía-Abril G, Román M, Luquero-Bueno S, Ochoa D, Abad-Santos F. Genetic Variation in CYP2D6, UGT1A4, SLC6A2 and SLCO1B1 Alters the Pharmacokinetics and Safety of Mirabegron. Pharmaceutics. 2024; 16(8):1077. https://doi.org/10.3390/pharmaceutics16081077

Chicago/Turabian Style

Soria-Chacartegui, Paula, Patricia Cendoya-Ramiro, Eva González-Iglesias, Samuel Martín-Vílchez, Andrea Rodríguez-Lopez, Gina Mejía-Abril, Manuel Román, Sergio Luquero-Bueno, Dolores Ochoa, and Francisco Abad-Santos. 2024. "Genetic Variation in CYP2D6, UGT1A4, SLC6A2 and SLCO1B1 Alters the Pharmacokinetics and Safety of Mirabegron" Pharmaceutics 16, no. 8: 1077. https://doi.org/10.3390/pharmaceutics16081077

APA Style

Soria-Chacartegui, P., Cendoya-Ramiro, P., González-Iglesias, E., Martín-Vílchez, S., Rodríguez-Lopez, A., Mejía-Abril, G., Román, M., Luquero-Bueno, S., Ochoa, D., & Abad-Santos, F. (2024). Genetic Variation in CYP2D6, UGT1A4, SLC6A2 and SLCO1B1 Alters the Pharmacokinetics and Safety of Mirabegron. Pharmaceutics, 16(8), 1077. https://doi.org/10.3390/pharmaceutics16081077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop