Prevalence and Accuracy of Information on CYP2D6, CYP2C19, and CYP2C9 Related Substrate and Inhibitor Co-Prescriptions in the General Population: A Cross‐Sectional Descriptive Study as Part of the PharmLines Initiative

Background: Drug-drug interaction (DDI) is one of the main contributors to adverse drug reactions and therefore, it is important to study its frequency in the population. We aimed to investigate frequency and concordance on CYP2D6, CYP2C19, and CYP2C9 (CYP2D6/2C19/2C9)-mediated potential DDIs at the Lifelines cohort and linked data from the pharmacy database IADB.nl.

Methods: As part of the University of Groningen PharmLines Initiative, data were collected on CYP2D6/2C19/2C9-related substrate/inhibitors from entry questionnaires of Lifelines participants and linked information from the pharmacy database IADB.nl. CYP2D6/2C19/2C9 related co-prescriptions were divided based on the type of drugs i.e. chronically used medication (CM) or occasionally used medication (OM). This resulted in the combination of two chronically used drugs (CM-CM), chronically and occasionally used medication (CM-OM), and two occasionally used drugs (OM-OM). To measure the agreement level, cohen’s kappa statistics and test characteristics were used. Results were stratified by time window, gender, and age.

Results: Among 80,837 medicine users in the Lifelines, about 1–2 per hundred participants were exposed to a CYP2D6/2C19/2C9-mediated potential DDI. Overall, the overlapping time window of three months produced the highest mean kappa values between the databases i.e. 0.545 (95% CI:0.544–0.545), 0.512 (95% CI:0.511–0.512), and 0.374 (95% CI:0.373–0.375), respectively. CM-CM had a better level of agreement (good) than CM-OM (fair to moderate) and OM-OM combination (poor to moderate). The influence of gender on concordance values was different for different CYPs. Among older persons, agreement levels were higher than for the younger population.

Conclusions: CYP2D6/2C19/2C9-mediated potential DDIs were frequent and concordance of data varied by time window, type of combination, sex and age. Subsequent studies should rather use a combination of self-reported and pharmacy database information.

Introduction

Drug-drug interactions (DDIs) are an important contributor to adverse drug reaction leading to hospitalization or mortality (Doucet et al., 1996; Montane et al., 2018). CYP2D6, CYP2C19, and CYP2C9 (CYP2D6/2C19/2C9), subtypes of CYP450 drug metabolizing enzymes, are commonly involved in mediating partly inappropriate DDIs as these enzymes metabolize a wide variety of drugs in clinical practice (Flockhart and Oesterheld, 2000; Bahar et al., 2017b). CYP2D6/2C19/2C9 are highly polymorphic enzymes and the genetic polymorphisms produce inter-individual variabilities in drug metabolisms ranging from poor to accelerated metabolic activities (Zanger and Schwab, 2013). Consequently, the clinical impact of CYP2D6/2C19/2C9-mediated DDIs might be variable from person to person and depends on his/her genetic profile (Bahar et al., 2017a). The information on the clinical relevance and management of DDIs mediated by different CYP2D6/2C19/2C9 genotypes is therefore needed. In order to provide the information, the first step that needs to be done is to generate the data about the burden and type of potential DDIs mediated by these enzymes in the general population.

Estimation of the prevalence rate of a potential DDI is commonly performed using self-reporting methods in which patients are interviewed or filled out a self-administered questionnaire (Van den Brandt et al., 1991; Classen et al., 2007; Secoli et al., 2010). However, this kind of assessment is prone to information bias, because of inaccurate recall, which may influence the validity of results (Rockenbauer et al., 2001; West et al., 2005). Hence, it is important to validate drug information collected with self-reporting methods (Haapea et al., 2010; Hafferty et al., 2018).

The Lifelines cohort is a Dutch three-generation population cohort that provides a wide variety of medical and non-medical data, genomic information, and data on medication use (Stolk et al., 2008; Scholtens et al., 2014). The Lifelines cohort, as a prospective and long-term database, offers possibilities in pharmaco-epidemiological studies, such as assessing the impact of gene polymorphism on the magnitude of DDIs in the population. However, currently not much is known about the frequency, type and validity of potential DDIs in the open population.

This study has both a methodological and an epidemiological aim: we studied the frequency of potentially interacting substrates and inhibitors of the CYP2D6/2C19/2C9 and the concordance level of the information derived by self-reported drug use and an analysis of data from a drug-use database. For the latter aim, information as observed in the Lifelines cohort was compared with data from a prescription database, the University of Groningen prescription database IADB.nl, across type of medications, sex, and age (Visser et al., 2013; Sediq et al., 2018). A prescription database is regarded as an accurate database and not to be influenced by so-called recall bias (Monster et al., 2002; Schneeweiss and Avorn, 2005). Additionally, IADB.nl has been proven a reliable database in many pharmaco-epidemiological studies (Daud et al., 2017; Alfian et al., 2018; Bahar et al., 2018b).

Materials and Methods

The PharmLines Initiative is a university wide project in which the data of the Lifelines Cohort study have been linked to the University of Groningen prescription database IADB.nl. The project was started in 2017 by the Groningen Research Institute of Pharmacy, Departments of Epidemiology and Clinical Pharmacy, Department of Pharmacology of the University Medical Center Groningen and the Lifelines Cohort Study (https://www.lifelines.nl/researcher/cohort-and-biobank) (Sediq et al., 2018).

The Lifelines Cohort

The Lifelines cohort covers 167,729 participants from the Northern part of the Netherlands, aged 6 months until 93 years old, which were recruited from 2006 until 2013 (Stolk et al., 2008; Scholtens et al., 2014). It is an observational cohort study intended to facilitate research on the contribution and interaction between environmental, genetic, and phenotypic aspects in the development of chronic diseases and healthy aging (Stolk et al., 2008; Scholtens et al., 2014). The recruited participants will be followed for at least 30 years and are asked to complete a questionnaire every 1.5 years. In addition, once every 5 years, the participants have a comprehensive physical examination (Stolk et al., 2008; Scholtens et al., 2014). Baseline questionnaires included questions about general information, lifestyle and environment, psychosocial aspects, and health (including medication use) (Stolk et al., 2008; Scholtens et al., 2014; Klijs et al., 2015). The medication use information were collected in two ways i.e. a) patients filled out a questionnaire or b) patients carried the medication at the time of interview (Sediq et al., 2018). The medication data regarding their current prescription and dose were recorded and classified using the Anatomical Therapeutic Chemical (ATC) coding scheme (Stolk et al., 2008; Scholtens et al., 2014). The Lifelines population is multigenerational and generally representative of the Dutch population resided in the Northern part of the Netherlands (Klijs et al., 2015).

University of Groningen IADB.nl Database

The University of Groningen prescription database IADB.nl has recorded prescriptions from community pharmacies in the Netherlands since 1994, and is updated annually (Visser et al., 2013; Sediq et al., 2018). In 2017, it contained prescription data of approximately 700,000 individuals from around 72 pharmacies that are located in most of the area where the Lifelines cohort is also resident. The study population was reported to represent the general population in Netherlands (Visser et al., 2013; Sediq et al., 2018). In the IADB.nl, each patient has a unique and anonymous identifier. Each record contains information about patient’s sex, date of birth, and information about his/her prescribed medication such as ATC code, duration, daily dose, amount prescribed, and dispensing date (Visser et al., 2013; Sediq et al., 2018). The IADB.nl has no information about over-the-counter (OTC) drugs and prescriptions from the hospital.

Study Population and Linkage of Databases

The study population consists of all medicine users (≥18 years) in the Lifelines cohort. A Trusted Third Party, Statistics Netherlands (Dutch: Centraal Bureau voor de Statistiek; CBS), carried out the linkage of the Lifelines and the IADB.nl records at the patient level based on postal code in combination with sex and date of birth. The unique identifiers from both databases were removed, and once the linkage was completed, each patient was assigned a new unique code that cannot be traced back to their previous identifier. Using the new identifier, the data from both databases could be combined. The complete linking process was described in more detail by Sediq et al. (Sediq et al., 2018).

Exposures

Exposures were defined as substrates and inhibitors of CYP2D6/2C19/2C9. We defined a potential DDI as each combination of a substrate and inhibitor listed in the international standard and local guideline, Flockhart Table for CYP-mediated drug interactions and the Dutch Commentaren Medicatiebewaking book, respectively (Borgsteede, 2015; Flockhart, 2018). Based on the main indication according to the official product information, the exposures were classified as: 1) chronically used medication (CM) for example CYP2D6 substrates such as beta-blockers (metoprolol), and 2) occasionally used medication (OM) for example CYP2D6 substrates such as opioids (tramadol). The full list of medications including their classification can be found in Supplementary Material 1.

Outcomes

Outcome measures were defined as frequency of potential CYP2D6/2C19/2C9-mediated DDIs as well as the levels of agreement between the self-reported information from the Lifelines cohort and the IADB.nl prescription data on these potential DDIs across type of medications, age, and sex. If one patient was exposed to different types of CYP2C9/2D6/2C19 mediated DDIs, we calculated them as one participant with several incidences of potential DDIs. If the potential DDI was only found in the Lifelines cohort records, it was categorized as over-reporting (false positive). If the potential DDI was only found in the IADB.nl, it was categorized as under-reporting (false negative). We also provided data on sensitivity, specificity, negative predictive value/NPV, and positive predictive value/PPV for the top five potential DDIs detected in the lifelines database. Different overlapping time windows (i.e. 1 month, 3 months, 6 months, 9 months, and 1 year) between baseline date of self-reporting medication in the Lifelines cohort and dispensing date of prescription in the IADB.nl were applied to determine the optimum time window for assessing the agreement of both databases. Subgroup analyses by the type of medication (CM vs OM), age, and sex were performed to observe the potential influence of these factors on the agreement. Additionally, we also presented information about the clinical relevance of the potential DDIs based on the suggested management provided by Epocrates^® i.e. “contraindicated, avoid combination/use alternative, modify treatment/monitor and caution”. If Epocrates^® had no recommendation for the potential DDI, we checked whether Drugs.com, another online drug interactions screening software, provided suggestions for the potential DDI. Both of them were reported to have a high sensitivity for detection of potential DDIs (Perkins et al., 2006; Bossaer and Thomas, 2017).

Statistical Methods

Comparisons of the prevalence of potential CYP2D6/2C19/2C9-mediated DDIs [mean (SD)] and the frequency of participants with the potential DDIs [number (%)] across age groups (18–59 vs >=60 years old) and sex (men vs women) were performed by using independent sample t-test. A p-value which is less than 0.05 (< 0.05) is considered to indicate a statistically significant difference between comparison groups. Multivariate analysis of the influence of age and sex on the risk of having the potential DDIs was conducted by using a binary logistic regression method to obtain the crude and covariate-adjusted odds ratios as a measure of association. A p-value < 0.05 and 95% confidence interval (CI) not including 1 are considered as indicators for significant associations. To determine the agreement values between the databases on the potential DDIs, we used Cohen’s kappa statistics and 95% CI. Altman et al. provided some guidelines to define the Cohen’s kappa values i.e. poor (< 0.20), fair (0.20–0.40), moderate (0.41–0.60), good (0.61–0.80), and very good (0.81–1.00) (Altman, 1990).

Results

Among of 167,729 Lifelines participants, 80,837 adults were recorded with self-reported medicine use (mean age 46 years and 68.5% women) in the cohort at entry (Table 1). Among the subjects, there were 1,125 (1.4%) self-reported medicine users exposed to 1,199 potential CYP2D6/2C19/2C9-mediated DDIs (Figure 1). The prevalence of potential CYP2D6/2C19/2C9-mediated DDIs was 488, 513, and 198 respectively (Table 2). Older women had a significantly higher risk to be exposed to CYP2D6 (OR: 2.159, 95% CI: 1.386–3.363) and CYP2C19 (OR: 1.691, 95% CI: 1.184–2.416) mediated DDIs than older men but the comparable risks were observed among younger group. There was also a tendency that women had an increased risk to be exposed to CYP2C9 mediated DDIs than men (Table 3).

TABLE 1

Table 1 Characteristics of participants with self-reported medication use at entry in the Lifelines cohort database and overlap with IADB.nl prescription database.

FIGURE 1

Figure 1 Selection of the population study. DDI, drug-drug interaction.

TABLE 2

Table 2 Prevalence and participants with potential DDIs in the Lifelines cohort.

TABLE 3

Table 3 Multivariate analysis on the influence of age and sex on risk of having potential CYP2C9/2D6/2C19 mediated DDIs.

There were 24% and 47% of CYP2D6 and CYP2C19 mediated co-prescriptions, respectively, which were in the category of “avoid combination/use alternative”. Additionally, about 65%, 43%, and 93% of CYP2D6/2C19/2C9-mediated combinations were in the category of “modify treatment/monitor” according to the knowledgebase (Figure 2).

FIGURE 2

Figure 2 Proportion of potential DDIs based on the suggested managements provided by Epocrates® and Drugs.com.

Information from 45,160 Lifelines participants could be linked to the IADB.nl database. Among this linked population, there were 25,387 self-reported medicine users with comparable age and sex distribution (mean age 45.5 years and 68.6% women) as observed in the total medicine users in the Lifelines cohort (Table 1). Metoprolol-paroxetine (83 events), citalopram-omeprazole (173 events), and diclofenac-paroxetine (51 events) were the most prevalent potential DDIs mediated by CYP2D6/2C19/2C9, with good, moderate, and fair agreement of questionnaire and prescription data, respectively. Data on kappa, sensitivity, specificity, PPV and PPV values of the top five most frequent potential DDIs in the Lifelines database can be found in Table 4. Information on self-reported combinations of chronically used medications such as metoprolol-fluoxetine had very good agreement, high sensitivity and specificity as well as high PPV and NPV. Meanwhile, information on self-reported combinations with occasionally used medication such as ibuprofen-paroxetine tended to have fair kappa, sensitivity, and PPV but high specificity and NPV. The complete list of the potential DDIs in the Lifelines database and their kappa values can be found in the Supplementary Material 2 and 3, respectively.

TABLE 4

Table 4 Top five potential DDIs in the Lifelines cohort detected in the IADB.nl database with their kappa, sensitivity, specificity, PPV, and NPV, as well as 95% CI (time window: 3 months).

The application of different time windows resulted in different agreement levels of the potential DDIs (Figure 3). Overall, the time window of three months produced the highest mean kappa values among potential CYP2D6/2C19/2C9-mediated DDIs i.e. moderate [0.545 (95% CI: 0.544–0.545)], moderate [0.512 (95% CI: 0.511–0.512)], and fair [0.374 (95% CI: 0.373–0.375)], respectively. Extension of the time windows to 6, 9, and 12 months decreased the mean kappa values. The time window of 1 month also produced a low kappa value. For the time window of 3 months, subgroup analysis for the type of medication indicated the potential DDIs in CM-CM had better level of agreements (good) than CM-OM (fair to moderate) and OM-OM (poor to moderate). For the CYP2D6 and CYP2C9 mediated DDIs, CM-OM combination had better kappa values (fair agreement) than OM-OM combination (poor agreement). Meanwhile, for the CYP2C19 mediated DDIs, both the CM-OM and OM-OM combination had comparable agreement level (moderate). The summary of the results can be found in Supplementary Material 3.

FIGURE 3

Figure 3 The effect of different time windows on the agreement between the Lifelines cohort and the IADB.nl database.

Subgroup analysis of agreement by sex showed mixed results (Figure 4). In CYP2D6 mediated potential DDIs, females appeared to have a better level of agreement than males. The opposite result was observed in CYP2C19 and CYP2C9 mediated potential DDIs where males mostly had a better kappa value compared to females. Stratification by age indicated that people aged 60 years or older had a generally better kappa value than the younger population in CYP2D6/2C19/2C9 mediated potential DDIs (Figure 5).

FIGURE 4

Figure 4 The effect of sex on the agreement between the Lifelines cohort and the IADB.nl database.

FIGURE 5

Figure 5 The effect of age on the agreement between the Lifelines cohort and the IADB.nl database.

Discussion

In this cross-sectional study, CYP2D6/2C19/2C9-mediated potential DDIs were frequent and concordance of data varied by time window, type of medication, sex, and age. We found that one to two per hundred drug users in the Lifelines cohort were exposed to a potential CYP2D6/2C19/2C9-mediated DDI at a short moment in life time. Some of these potential DDIs are regarded as clinically relevant DDIs such as metoprolol and CYP2D6 inhibitors combinations. The DDIs may lead to bradycardia, hypotension, and atrioventricular block (Walley et al., 1993; Konig et al., 1996; Onalan et al., 2008; Bahar et al., 2018a). Other relevant DDIs were the combination of CYP2C9 inhibitors that consist of selective serotonin inhibitors (SSRIs), and nonsteroidal anti-inflammatory drugs (NSAIDs). The combination of SSRIs and NSAIDs was reported to increase risk of gastrointestinal bleedings (de Abajo et al., 1999; De Jong et al., 2003). Yet, the interaction between SSRIs and NSAIDs might be not solely a pharmacokinetic interaction but also involves a pharmacodynamic interaction (Moore et al., 2015). Our findings on the burden of DDIs might have potential clinical as well as economic implications. A DDI is one of the main contributors of an adverse drug reaction (ADR) which is one of the leading causes of hospitalisation and it can cost at minimum around €200 to €9,000 per hospitalisation (Formica et al., 2018).

Older women tended to have a higher risk to be exposed to potential CYP2D6/2C19/2C9-mediated DDIs than older men. A survey study from the United States about the pattern of drug use among adults in the outpatient setting indicated that elderly women (≥65 years old) had the highest burden of medication use in which 23% and 12% of them used at least five and 10 drugs, respectively (Kaufman et al., 2002). The risk of experiencing DDI is increased as the number of drugs taken also increased (Åstrand et al., 2006). Taking five to seven drugs and 10 to 14 drugs enhanced the risk of potentially relevant DDI by about 20%–30% and 40%–60%, respectively (Johnell and Klarin, 2007). Other studies also reported that being women and old are risk factors associated with DDIs (Grattagliano et al., 2010; Magro et al., 2012).

Based on this study, a three-month time window appeared to result in the best agreement level. This is consistent with the previous study by Sediq et al. about the validation of single drug used in the Lifelines (Sediq et al., 2018). Additionally, Lau et al. also found similar finding in their work on the validation of pharmacy records in Amsterdam, Netherlands (Lau et al., 1997). One of the possible reasons for this finding is that the Dutch reimbursement system only allows drugs to be prescribed for a maximum of 3 months period of supply (Lau et al., 1997). A comparable study to validate self-reported medication using a national prescription database in the Danish population also found a fixed 3 month time window was suitable for checking the agreement between the two sources of information on medicine use (Nielsen et al., 2008). Considering the time window is an important aspect, because a long time window may hamper the analysis of drugs used on as needed basis and a short time window may impair the analysis of drugs used chronically (Johnson and Vollmer, 1991; Nielsen et al., 2008).

Sub-analysis by type of medication indicates self-reported information on the CM-CM combination was more reliable than information on CM-OM and OM-OM combination. It offers the possibility to use co-prescription data of the CM-CM combination from the Lifelines cohort in research. This may because routinely used medication is more easily remembered by patients than drugs used occasionally. These results are consistent with previous studies (Nielsen et al., 2008; Sarangarm et al., 2012; Cohen et al., 2018).

Furthermore, for CYP2D6 and CYP2C9 related co-prescriptions, the kappa value of the OM-CM combination (fair agreement) was higher than OM-OM (poor agreement) except for CYP2C19. For the latest, the agreement level of CYP2C19 mediated OM-CM seems comparable with those of the OM-OM combination (moderate agreement). There are some possible explanations for this finding. One is the inclusion of proton pump inhibitors (PPIs) in the OM groups which are the main drugs in this group. PPIs have wide therapeutic indications and some of these indications need a chronic use of PPIs such as Zollinger-Ellisone syndrome, Barrett’s esophagus, and esophagitis (KNMP, 2018b). Another explanation is the inclusion of diazepam in the OM groups. Diazepam may be used chronically for treating patients with panic disorder and generalized anxiety disorder (KNMP, 2018a). Consequently, OM groups not solely consisted of drugs used ‘as needed’ but also may include chronically used drugs.

We found that the effect of sex on agreement is not consistent. Previous studies also reported mixed results. Some reports showed that men had a better recall accuracy than women (Linet et al., 1989; Haapea et al., 2010). Meanwhile, other studies indicated that sex had no influence on the agreement between self-reported medication use and prescription database (Van den Brandt et al., 1991; West et al., 1995). Therefore, more research is needed to determine the effect of sex on the recall accuracy, and concordance between self-reported medication use and information from a prescription database.

Our study found the agreement between the Lifelines and the IADB.nl database records is better in the older population (aged 60 years and older) than in younger adults. This result is in contrast to previous reports which found aging led to a low agreement between self-reported medication use and information from a drug database (Van den Brandt et al., 1991; West et al., 1995). A decrease in cognitive function and polypharmacy may cause poor recall information by old patients (Van den Brandt et al., 1991). However, other studies reported that age did not influence the agreement level (Sjahid et al., 1998; Lamiae et al., 2010). The method used to collect drug information may determine the influence of aging in recall bias. If the interviewers visit the patient’s house to ascertain the consumed drugs or if the patients are helped by their family in completing the questionnaire, the impact of self-reporting bias in the old participants (≥ 60 years old) can be reduced (Johnson and Vollmer, 1991; Lau et al., 1997; Richardson et al., 2013). In the Lifelines cohort, some participants filled out the questionnaire at home before visiting the premises. Therefore, the participants were potentially assisted by their relatives or may directly check their medication while completing the questionnaire. Meanwhile, some patients brought their medication at the time of interview so that interviewers could ascertain their medication list in the questionnaire.

Another possible culprit of conflicting reports is the type of medication. Most of the drugs related to CYP2D6/219/2C9 are mainly groups of drugs used by the old population chronically and were reported to be associated with a good recall such as cardiac therapy, antidiabetic agents, anti-thrombotic drugs, anticancer agents, antidepressant, and antipsychotic agents (Van den Brandt et al., 1991; Haukka et al., 2007; Gupta et al., 2011; Hafferty et al., 2018; Sediq et al., 2018). For the last two agents, Haukka et al. reported a good recall because patients brought their medication at the time of interview (Haukka et al., 2007). Lastly, the other possible explanation was the differential distribution of the population in each subgroup of age which may give a wrong impression about the influence of different age in the agreement (West et al., 1997). In our study, about 80% of the population is in the 18–59 years old subgroup. Therefore, a larger pharmaco-epidemiological study with a sub-group analysis is needed to elucidate the impact of age in the concordance of self-reported medication use and data from a prescription database.

Some strengths of our study are worth to be mentioned. Firstly, the linkage process between both databases is reliable since it was performed by CBS on individual level. Secondly, the population in our cohort is large and not limited to a certain group of population with diseases or using specific medications. Some other studies were conducted by using a limited sample and only in some particular groups of patients such as in elderly, pregnant women, patients with specific medical conditions or using certain drugs (West et al., 1995; Sjahid et al., 1998; Rockenbauer et al., 2001; Richardson et al., 2013). Thirdly, we included all types of drugs which may potentially trigger CYP2D6/219/2C9-mediated DDIs. However, there are also some limitations from our study. Firstly, we only checked the agreement of prescribed medication but not OTC drugs since the IADB.nl database has no information on OTC drugs. For example, ibuprofen is also available over the counter which may explain lower kappa values in potential DDI combinations. Next, we only had drug information from community pharmacies and, therefore, if the drugs recorded in the self-reported questionnaire were obtained from a hospital, it will not be detected in the IADB.nl and will be categorized as over-reporting information. Further, some potential DDIs included in our study were recommended to manage either by dose adjustments or monitoring of the possible potential side effects which have been possibly done by the responsible clinicians. Meanwhile, some other potential DDIs might not produce important side effects and only need caution on their use. However, we still kept them in our analysis because the influence of genetic polymorphisms on CYP2D6/2C19/2C9 may enhance the magnitude of the clinical impact of those interactions (Bahar et al., 2017a). Therefore, it would be valuable to research the interaction of CYP2D6/2C19/2C9 polymorphisms and CYP2D6/2C19/2C9-mediated DDIs in the next study. Additionally, we did not include the combination of substrates and inducers since the prevalence was too low to allow further analysis. Next, we only limited the focus of our study on the contribution of the three main phase I drug metabolizing enzymes (CYP2D6/2C19/2C9) since they cumulatively metabolize about 42% of drugs currently used in the clinical practice and mounting evidence has shown that clinical consequences of genetic polymorphisms are different among the CYP450 subfamily with CYP2D6/2C19/2C9 polymorphisms reported to have the most important clinically relevant implications (Zanger et al., 2008; Zanger and Schwab, 2013; Bahar et al., 2017a). Therefore, we assume that the CYP2D6/2C19/2C9 mediated DDIs will be the most frequent and relevant DDIs which will be found to be modified by genetic polymorphism in clinical practice. However, we would like to emphasize that drug interactions might also be facilitated by phase II drug metabolizing enzymes and drug transporters which are also subject to genetic polymorphisms and still are not optimally investigated (Board et al., 1998; Kerb, 2006). Furthermore, we only limited the analysis of the potential DDIs to the pairwise combination of medications since it reflects current guidelines and practice related to the management of DDIs in the Netherlands (van Roon et al., 2005; Heringa et al., 2016). However, the DDI might occur not only in the form of bimodal (involving two drugs) interaction but also in multimodal (involving more than two drugs) interactions especially in drugs metabolized by multiple metabolic pathways (Grönlund et al., 2011). Multimodal drug interactions were reported to produce more severe outcomes than bimodal interaction since all the metabolic pathways of the drugs are impaired (Niemi et al., 2003; Niemi et al., 2006). Lastly, our study had no clinical outcomes of the potential DDIs since in the current study our focus was limited to study the prevalence of the potential DDIs and the agreement of drug information between both databases. This study is pivotal in order to design valid follow-up studies with the aim to determine the clinical impact of the observed potential DDIs especially for chronically used medications.

Conclusion

In conclusion, CYP2D6/2C19/2C9-mediated potential DDIs were frequent and the agreement between the Lifelines cohort and the IADB.nl differed between time windows. The best concordance level was achieved at a 3-month time window. CM-CM co-prescription had a better agreement than CM-OM and OM-OM combinations. Sex had no consistent influence on the discordance between the databases. Meanwhile, the older population had a better kappa value than the younger population. For the next drug study, the self-reporting data should be complemented with the pharmacy data in order to achieve a better accuracy in capturing the real word information on medication use.

Data Availability Statement

The datasets for this manuscript are not publicly available, and used under license for the current study, in order to protect information that could potentially compromise the privacy of research participants. Requests to access the datasets should be directed to the Pharmlines Initiatives [email: cmVzZWFyY2hAbGlmZWxpbmVzLm5s]. Some data such as frequency, type and agreements (kappa values) of drug-drug-interactions are available as Supplementary Materials.

Ethics Statement

Ethics approval was not needed according to institutional guidelines and national legislation, since the data generated in this manuscript relied exclusively on the research database with pseudonymized information and informed consent was obtained at the time of original data collection

Author Contributions

Conceptualization and design (MB, JB, SB, BW, EH), data curation (MB, JB), investigation (MB, JB), resources (JB, AD, RA), statistical analysis (MB), interpretation of data (MB, SB, AD, RA, BW, EH), drafting manuscript (MB), critical evaluation of the manuscript, editing, and supervision (JB, SB, AD, RA, BW, EH). All authors read and agreed on the final version of the manuscript.

Funding

Lifelines is financially supported by several parties such as the Dutch Government, The Netherlands Organization of Scientific Research NWO (grant 175.010.2007.006), the European fund for regional development, Dutch Ministry of Economic Affairs, the Northern Netherlands Collaboration of Provinces (SNN), Provinces of Groningen and Drenthe, Pieken in de Delta, University Medical Center Groningen, and University of Groningen, Netherlands. Meanwhile, the prescription database IADB.nl and the PharmLines Initiative are financially supported by the Groningen Research Institute of Pharmacy, University of Groningen. MB obtained a DIKTI scholarship from the Ministry of Research, Technology and Higher Education of Indonesia. The funding organizations had no role and influence in the study design and results.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank all participants in the IADB.nl and Lifelines cohort for providing the data used in this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.00624/full#supplementary-material

References

Åstrand, B., Åstrand, E., Antonov, K., Petersson, G. (2006). Detection of potential drug interactions–a model for a national pharmacy register. Eur. J. Clin. Pharmacol. 62, 749–756. doi: 10.1007/s00228-006-0143-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfian, S. D., Worawutputtapong, P., Schuiling-Veninga, C. C., van der Schans, J., Bos, J. H., Hak, E., et al. (2018). Pharmacy-based predictors of non-persistence with and non-adherence to statin treatment among patients on oral diabetes medication in the Netherlands. Curr. Med. Res. Opin. 34, 1013–1019. doi: 10.1080/03007995.2017.1417242

PubMed Abstract | CrossRef Full Text | Google Scholar

Altman, D. G. (1990). Practical statistics for medical research (Chapman & Hall, London: CRC press).

Google Scholar

Bahar, M. A., Setiawan, D., Hak, E., Wilffert, B. (2017a). Pharmacogenetics of drug–drug interaction and drug–drug–gene interaction: a systematic review on CYP2C9, CYP2C19 and CYP2D6. Pharmacogenomics 18, 701–739. doi: 10.2217/pgs-2017-0194

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahar, M., Hak, E., Bos, J. H., Borgsteede, S. D., Wilffert, B. (2017b). The burden and management of cytochrome P450 2D6 (CYP2D6)-mediated drug–drug interaction (DDI): co-medication of metoprolol and paroxetine or fluoxetine in the elderly. Pharmacoepidemiol. Drug Saf. 26, 752–765. doi: 10.1002/pds.4200

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahar, M. A., Kamp, J., Borgsteede, S. D., Hak, E., Wilffert, B. (2018a). The Impact of CYP2D6 mediated Drug-Drug Interaction: A Systematic Review on a Combination of Metoprolol and Paroxetine/Fluoxetine. Br. J. Clin. Pharmacol. doi: 10.1111/bcp.13741

CrossRef Full Text | Google Scholar

Bahar, M. A., Wang, Y., Bos, J. H., Wilffert, B., Hak, E. (2018b). Discontinuation and dose adjustment of metoprolol after metoprolol-paroxetine/fluoxetine co-prescription in Dutch elderly. Pharmacoepidemiol. Drug Saf. 27, 621–629. doi: 10.1002/pds.4422

PubMed Abstract | CrossRef Full Text | Google Scholar

Board, P., Blackburn, A., Jermiin, L. S., Chelvanayagam, G. (1998). Polymorphism of phase II enzymes: identification of new enzymes and polymorphic variants by database analysis. Toxicol. Lett. 102, 149–154. doi: 10.1016/S0378-4274(98)00300-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Borgsteede, S. (2015). Commentaren medicatiebewaking.

Google Scholar

Bossaer, J. B., Thomas, C. M. (2017). Drug interaction database sensitivity with oral antineoplastics: An exploratory analysis. J. Oncol. Pract. 13, e217–e222. doi: 10.1200/JOP.2016.016212

PubMed Abstract | CrossRef Full Text | Google Scholar

Classen, S., Meuleman, J., Garvan, C., Ried, L. D., Mann, W., Asal, N. (2007). Review of prescription medications in home-based older adults with stroke: a pilot study. Res. Social. Adm. Pharm. 3, 104–122. doi: 10.1016/j.sapharm.2006.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, J. M., Wood, M. E., Hernandez-Diaz, S., Nordeng, H. (2018). Agreement between paternal self-reported medication use and records from a national prescription database. Pharmacoepidemiol. Drug Saf. 27, 413–421. doi: 10.1002/pds.4411

PubMed Abstract | CrossRef Full Text | Google Scholar

Daud, A. N., Bergman, J. E., Oktora, M. P., Kerstjens-Frederikse, W. S., Groen, H., Bos, J. H., et al. (2017). Maternal use of drug substrates of placental transporters and the effect of transporter-mediated drug interactions on the risk of congenital anomalies. PloS One 12, e0173530. doi: 10.1371/journal.pone.0173530

PubMed Abstract | CrossRef Full Text | Google Scholar

de Abajo, F. J., Rodriguez, L. A., Montero, D. (1999). Association between selective serotonin reuptake inhibitors and upper gastrointestinal bleeding: population based case-control study. BMJ 319, 1106–1109. doi: 10.1136/bmj.319.7217.1106

PubMed Abstract | CrossRef Full Text | Google Scholar

De Jong, J. C., Van Den, B., Paul, B., Tobi, H., De Jong, L. T., Van Den Berg, -. (2003). Combined use of SSRIs and NSAIDs increases the risk of gastrointestinal adverse effects. Br. J. Clin. Pharmacol. 55, 591–595. doi: 10.1046/j.0306-5251.2002.01770.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Doucet, J., Chassagne, P., Trivalle, C., Landrin, I., Pauty, M., Kadri, N., et al. (1996). Drug-drug interactions related to hospital admissions in older adults: A prospective study of 1000 patients. J. Am. Geriatr. Soc 44, 944–948. doi: 10.1111/j.1532-5415.1996.tb01865.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Flockhart, D. A., Oesterheld, J. R. (2000). Cytochrome P450-mediated drug interactions. Child Adolesc. Psychiatr. Clin. N. Am. 9, 43–76. doi: 10.1016/S1056-4993(18)30135-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Flockhart, D. A. (2018). Drug Interactions Flockhart Table™ . Indiana University School of Medicine. Available at: https://drug-interactions.medicine.iu.edu/MainTable.aspx

Google Scholar

Formica, D., Sultana, J., Cutroneo, P., Lucchesi, S., Angelica, R., Crisafulli, S., et al. (2018). The economic burden of preventable adverse drug reactions: a systematic review of observational studies. Expert. Opin. Drug Saf. 17, 681–695. doi: 10.1080/14740338.2018.1491547

PubMed Abstract | CrossRef Full Text | Google Scholar

Grönlund, J., Saari, T. I., Hagelberg, N. M., Neuvonen, P. J., Laine, K., Olkkola, K. T. (2011). Effect of inhibition of cytochrome P450 enzymes 2D6 and 3A4 on the pharmacokinetics of intravenous oxycodone. Clin. Drug Investig. 31, 143–153. doi: 10.2165/11539950-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Grattagliano, I., Portincasa, P., D’Ambrosio, G., Palmieri, V. O., Palasciano, G. (2010). Avoiding drug interactions: here’s help: you can’t count on an electronic prescribing system to catch all potential drug-drug interactions. These at-a-glance tables will help you minimize risk. J. Fam. Pract. 59, 322–329.

PubMed Abstract | Google Scholar

Gupta, V., Gu, K., Chen, Z., Lu, W., Shu, X. O., Zheng, Y. (2011). Concordance of self-reported and medical chart information on cancer diagnosis and treatment. BMC Med. Res. Methodol. 11, 72. doi: 10.1186/1471-2288-11-72

PubMed Abstract | CrossRef Full Text | Google Scholar

Haapea, M., Miettunen, J., Lindeman, S., Joukamaa, M., Koponen, H. (2010). Agreement between self-reported and pharmacy data on medication use in the Northern Finland 1966 Birth Cohort. Int. J. Methods Psychiatr. Res. 19, 88–96. doi: 10.1002/mpr.304

CrossRef Full Text | Google Scholar

Hafferty, J. D., Campbell, A. I., Navrady, L. B., Adams, M. J., MacIntyre, D., Lawrie, S. M., et al. (2018). Self-reported medication use validated through record linkage to national prescribing data. J. Clin. Epidemiol. 94, 132–142. doi: 10.1016/j.jclinepi.2017.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Haukka, J., Suvisaari, J., Tuulio-Henriksson, A., Lönnqvist, J. (2007). High concordance between self-reported medication and official prescription database information. Eur. J. Clin. Pharmacol. 63, 1069–1074. doi: 10.1007/s00228-007-0349-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Heringa, M., Floor-Schreudering, A., Tromp, P. C., de Smet, P. A., Bouvy, M. L. (2016). Nature and frequency of drug therapy alerts generated by clinical decision support in community pharmacy. Pharmacoepidemiol. Drug Saf. 25, 82–89. doi: 10.1002/pds.3915

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnell, K., Klarin, I. (2007). The relationship between number of drugs and potential drug-drug interactions in the elderly. Drug Saf. 30, 911–918. doi: 10.2165/00002018-200730100-00009

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, R. E., Vollmer, W. M. (1991). Comparing sources of drug data about the elderly. J. Am. Geriatr. Soc 39, 1079–1084. doi: 10.1111/j.1532-5415.1991.tb02872.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaufman, D. W., Kelly, J. P., Rosenberg, L., Anderson, T. E., Mitchell, A. A. (2002). Recent patterns of medication use in the ambulatory adult population of the United States: the Slone survey. JAMA 287, 337–344. doi: 10.1001/jama.287.3.337

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerb, R. (2006). Implications of genetic polymorphisms in drug transporters for pharmacotherapy. Cancer Lett. 234, 4–33. doi: 10.1016/j.canlet.2005.06.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Klijs, B., Scholtens, S., Mandemakers, J. J., Snieder, H., Stolk, R. P., Smidt, N. (2015). Representativeness of the LifeLines cohort study. PloS one 10, e0137203. doi: 10.1371/journal.pone.0137203

PubMed Abstract | CrossRef Full Text | Google Scholar

KNMP (2018a). Informatorium Medicamentorum: Diazepam. 2018.

Google Scholar

KNMP (2018b). MAAGDARMMIDDELEN. 2018.

Google Scholar

Konig, F., Hafele, M., Hauger, B., Loble, M., Wossner, S., Wolfersdorf, M. (1996). Bradycardia after beginning therapy with metoprolol and paroxetine. Psychiatr. Prax. 23, 244–245.

PubMed Abstract | Google Scholar

Lamiae, G., Michel, R., Elodie, A., Nabil, E. K., Jacques, B., Lucien, A. (2010). Agreement between patients’ self-report and physicians’ prescriptions on cardiovascular drug exposure: the PGRx database experience. Pharmacoepidemiol. Drug Saf. 19, 591–595. doi: 10.1002/pds.1952

PubMed Abstract | CrossRef Full Text | Google Scholar

Lau, H. S., de Boer, A., Beuning, K. S., Porsius, A. (1997). Validation of pharmacy records in drug exposure assessment. J. Clin. Epidemiol. 50, 619–625. doi: 10.1016/S0895-4356(97)00040-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Linet, M. S., Harlow, S. D., McLaughlin, J. K., McCaffrey, L. D. (1989). A comparison of interview data and medical records for previous medical conditions and surgery. J. Clin. Epidemiol. 42, 1207–1213. doi: 10.1016/0895-4356(89)90119-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Magro, L., Moretti, U., Leone, R. (2012). Epidemiology and characteristics of adverse drug reactions caused by drug–drug interactions. Expert. Opin. Drug Saf. 11, 83–94. doi: 10.1517/14740338.2012.631910

PubMed Abstract | CrossRef Full Text | Google Scholar

Monster, T. B., Janssen, W. M., de Jong, P. E., de Jong-van den Berg, L. T. W., PREVEND Study Group (2002). Pharmacy data in epidemiological studies: an easy to obtain and reliable tool. Pharmacoepidemiol. Drug Saf. 11, 379–384. doi: 10.1002/pds.722

PubMed Abstract | CrossRef Full Text | Google Scholar

Montane, E., Arellano, A. L., Sanz, Y., Roca, J., Farre, M. (2018). Drug-related deaths in hospital inpatients: A retrospective cohort study. Br. J. Clin. Pharmacol. 84, 542–552. doi: 10.1111/bcp.13471

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, N., Pollack, C., Butkerait, P. (2015). Adverse drug reactions and drug-drug interactions with over-the-counter NSAIDs. Ther. Clin. Risk Manage. 11, 1061–1075. doi: 10.2147/TCRM.S79135

CrossRef Full Text | Google Scholar

Nielsen, M. W., Søndergaard, B., Kjøller, M., Hansen, E. H. (2008). Agreement between self-reported data on medicine use and prescription records vary according to method of analysis and therapeutic group. J. Clin. Epidemiol. 61, 919–924. doi: 10.1016/j.jclinepi.2007.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Niemi, M., Backman, J. T., Neuvonen, M., Neuvonen, P. J. (2003). Effects of gemfibrozil, itraconazole, and their combination on the pharmacokinetics and pharmacodynamics of repaglinide: potentially hazardous interaction between gemfibrozil and repaglinide. Diabetologia 46, 347–351. doi: 10.1007/s00125-003-1034-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Niemi, M., Tornio, A., Pasanen, M. K., Fredrikson, H., Neuvonen, P. J., Backman, J. T. (2006). Itraconazole, gemfibrozil and their combination markedly raise the plasma concentrations of loperamide. Eur. J. Clin. Pharmacol. 62, 463–472. doi: 10.1007/s00228-006-0133-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Onalan, O., Cumurcu, B. E., Bekar, L. (2008). Complete atrioventricular block associated with concomitant use of metoprolol and paroxetine. Mayo Clin. Proc. 83, 595–599. doi: 10.4065/83.5.595

PubMed Abstract | CrossRef Full Text | Google Scholar

Perkins, N. A., Murphy, J. E., Malone, D. C., Armstrong, E. P. (2006). Performance of drug–drug interaction software for personal digital assistants. Ann. Pharmacother. 40, 850–855. doi: 10.1345/aph.1G603

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, K., Kenny, R. A., Peklar, J., Bennett, K. (2013). Agreement between patient interview data on prescription medication use and pharmacy records in those aged older than 50 years varied by therapeutic group and reporting of indicated health conditions. J. Clin. Epidemiol. 66, 1308–1316. doi: 10.1016/j.jclinepi.2013.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Rockenbauer, M., Olsen, J., Czeizel, A. E., Pedersen, L., Sørensen, H. T., The EuroMAP Group (2001). Recall bias in a case-control surveillance system on the use of medicine during pregnancy. Epidemiology 12 (4), 461–466. doi: 10.1097/00001648-200107000-00017

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarangarm, P., Young, B., Rayburn, W., Jaiswal, P., Dodd, M., Phelan, S., et al. (2012). Agreement between self-report and prescription data in medical records for pregnant women. Birth Defects Res. A Clin. Mol. Teratol. 94, 153–161. doi: 10.1002/bdra.22888

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneeweiss, S., Avorn, J. (2005). A review of uses of health care utilization databases for epidemiologic research on therapeutics. J. Clin. Epidemiol. 58, 323–337. doi: 10.1016/j.jclinepi.2004.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Scholtens, S., Smidt, N., Swertz, M. A., Bakker, S. J., Dotinga, A., Vonk, J. M., et al. (2014). Cohort Profile: LifeLines, a three-generation cohort study and biobank. Int. J. Epidemiol. 44, 1172–1180. doi: 10.1093/ije/dyu229

PubMed Abstract | CrossRef Full Text | Google Scholar

Secoli, S., Figueras, A., Lebrao, M. L., de Lima, F. D., Santos, J. L. F. (2010). Risk of potential drug-drug interactions among Brazilian elderly. Drugs Aging 27, 759–770. doi: 10.2165/11538460-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Sediq, R., van der Schans, J., Dotinga, A., Alingh, R. A., Wilffert, B., Bos, J. H., et al. (2018). Concordance Assessment of Self-Reported Medication Use in the Netherlands Three-Generation Lifelines Cohort Study With the Pharmacy Database iaDB.nl: The PharmLines Initiative. Clin. Epidemiol. 10, 981–989. doi: 10.2147/CLEP.S163037

PubMed Abstract | CrossRef Full Text | Google Scholar

Sjahid, S. I., Van Der, L., Paul, D., Stricker, B. H. C. (1998). Agreement between the pharmacy medication history and patient interview for cardiovascular drugs: the Rotterdam elderly study. Br. J. Clin. Pharmacol. 45, 591–595. doi: 10.1046/j.1365-2125.1998.00716.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Stolk, R. P., Rosmalen, J. G., Postma, D. S., de Boer, R. A., Navis, G., Slaets, J. P., et al. (2008). Universal risk factors for multifactorial diseases. Eur. J. Epidemiol. 23, 67–74. doi: 10.1007/s10654-007-9204-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Van den Brandt, P. A., Petri, H., Dorant, E., Goldbohm, R. A., Van de Crommert, S. (1991). Comparison of questionnaire information and pharmacy data on drug use. Pharm. Weekbl. 13, 91–96. doi: 10.1007/BF01974987

CrossRef Full Text | Google Scholar

van Roon, E. N., Flikweert, S., le Comte, M., Langendijk, P. N., Kwee-Zuiderwijk, W. J., Smits, P., et al. (2005). Clinical relevance of drug-drug interactions. Drug Saf. 28, 1131–1139. doi: 10.2165/00002018-200528120-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Visser, S. T., Schuiling-Veninga, C. C., Bos, J. H., den Berg, de J.-v., TW, L., Postma, M. J. (2013). The population-based prescription database IADB.nl: its development, usefulness in outcomes research and challenges. Expert. Rev. Pharmacoecon. Outcomes Res. 13, 285–292. doi: 10.1586/erp.13.20

PubMed Abstract | CrossRef Full Text | Google Scholar

Walley, T., Pirmohamed, M., Proudlove, C., Maxwell, D. (1993). Interaction of metoprolol and fluoxetine. Lancet 341, 967–968. doi: 10.1016/0140-6736(93)91265-N

CrossRef Full Text | Google Scholar

West, S. L., Savitz, D. A., Koch, G., Strom, B. L., Guess, H. A., Hartzema, A. (1995). Recall accuracy for prescription medications: self-report compared with database information. Am. J. Epidemiol. 142, 1103–1112. doi: 10.1093/oxfordjournals.aje.a117563

PubMed Abstract | CrossRef Full Text | Google Scholar

West, S. L., Savitz, D. A., Koch, G., Sheff, K. L., Strom, B. L., Guess, H. A., et al. (1997). Demographics, health behaviors, and past drug use as predictors of recall accuracy for previous prescription medication use. J. Clin. Epidemiol. 50, 975–980. doi: 10.1016/S0895-4356(97)00026-7

PubMed Abstract | CrossRef Full Text | Google Scholar

West, S. L., Strom, B. L., Poole, C. (2005). "Validity of pharmacoepidemiology drug and diagnosis data," in Pharmacoepidemiology Ed. Strom, B. L. (Chichester, England: John Wiley & Sons Ltd), 709–765.4.

Google Scholar

Zanger, U. M., Schwab, M. (2013). Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol. Ther. 138, 103–141. doi: 10.1016/j.pharmthera.2012.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanger, U. M., Turpeinen, M., Klein, K., Schwab, M. (2008). Functional pharmacogenetics/genomics of human cytochromes P450 involved in drug biotransformation. Anal. Bioanal. Chem. 392, 1093–1108. doi: 10.1007/s00216-008-2291-6

PubMed Abstract | CrossRef Full Text | Google Scholar