Exploring the impact of CYP2D6 and UGT2B7 gene-drug interactions, and CYP-mediated DDI on oxycodone and oxymorphone pharmacokinetics using physiologically-based pharmacokinetic modeling and simulation.

Oxycodone is one of the most commonly used opioids to treat moderate to severe pain. It is metabolized mainly by CYP3A4 and CYP2D6, while only a small fraction of the dose is excreted unchanged into the urine. Oxymorphone, the metabolite primarily formed by CYP2D6, has a 40- to 60-fold higher mu-opioid receptor affinity than the parent compound. While CYP2D6-mediated gene-drug-interactions (GDIs) and drug-drug interactions (DDIs) are well-studied, they only account for a portion of the variability in oxycodone and oxymorphone exposure. The combined impact of CYP2D6-mediated GDIs and DDIs, CYP3A4-mediated DDIs, and UGT2B7 GDIs is not fully understood yet and hard to study in head-to-head clinical trials given the relatively large number of scenarios. Instead, we propose the use of a physiologically-based pharmacokinetic model that integrates available information on oxycodone's metabolism to characterize and predict the impact of DDIs and GDIs on the exposure of oxycodone and its major, pharmacologically-active metabolite oxymorphone. To this end, we first developed and verified a PBPK model for oxycodone and its metabolites using published clinical data. The verified model was then applied to determine the dose-exposure relationship of oxycodone and oxymorphone stratified by CYP2D6 and UGT2B7 phenotypes respectively, and administered perpetrators of CYP-based drug interactions. Our simulations demonstrate that the combination of CYP2D6 UM and a UGT2B7Y (268) mutation may lead to a 2.3-fold increase in oxymorphone exposure compared to individuals who are phenotyped as CYP2D6 NM / UGT2B7 NM. The extent of oxymorphone exposure increases up to 3.2-fold in individuals concurrently taking CYP3A4 inhibitors, such as ketoconazole. Inhibition of the CYP3A4 pathway results in a relative increase in the partial metabolic clearance of oxycodone to oxymorphone. Oxymorphone is impacted to a higher extent by GDIs and DDIs than oxycodone. We predict oxymorphone exposure to be highest in CYP2D6 UMs/UGT2B7 PMs in the presence of ketoconazole (strong CYP3A4 index inhibitor) and lowest in CYP2D6 PMs/UGT2B7 NMs in the presence of rifampicin (strong CYP3A4 index inducer) covering a 55-fold exposure range.

Identification of Novel and Early Biomarkers for Cisplatin-induced Nephrotoxicity and the Nephroprotective Role of Cimetidine using a Pharmacometabolomic-based Approach Coupled with In Vitro Toxicodynamic Modeling and Simulation.

Cisplatin is widely used for the treatment of various types of cancer. However, cisplatin-induced nephrotoxicity (CIN) is frequently observed in patients receiving cisplatin therapy which poses a challenge in its clinical utility. Currently used clinical biomarkers for CIN are not adequate for early detection of nephrotoxicity, hence there is a need to identify potential early biomarkers in predicting CIN. In the current study, a combination of in vitro toxicodynamic (TD) modeling and untargeted global metabolomics approach was used to identify novel potential metabolite biomarkers for early detection of CIN. In addition, we investigated the protective role of cimetidine (CIM), an inhibitor of the organic cation transporter 2 (OCT2), in suppressing CIN. We first characterized the time-course of nephrotoxic effects of cisplatin (CIS) and the protective effects of CIM in a human pseudo-immortalized renal proximal tubule epithelial cell line (RPTEC), SA7K cell line. Secondly, we used a mathematical cell-level, in vitro TD modeling approach to quantitatively characterize the time-course effects of CIS and CIM as single agents and combination in SA7K cells. Based on the experimental and modeling results, we selected relevant concentrations of CIS and CIM for our metabolomics study. With the help of PCA (Principal Component Analysis) and PLS-DA (Projection to Latent Structure - Discriminate Analysis) analyses, we confirmed global metabolome changes for different groups (CIS, CIM, CIS+CIM vs control) in SA7K cells. Based on the criterion of a p-value ≤ 0.05 and a fold change ≥ 2 or ≤ 0.5, we identified 20 top metabolites that were significantly changed during the early phase i.e. within first 12 h of CIS treatment. Finally, pathway analysis was conducted that revealed the key metabolic pathways that were most impacted in CIN.

Advancing drug safety science by integrating molecular knowledge with post-marketing adverse event reports.

Promising drug development efforts may frequently fail due to unintended adverse reactions. Several methods have been developed to analyze such data, aiming to improve pharmacovigilance and drug safety. In this work, we provide a brief review of key directions to quantitatively analyzing adverse events and explore the potential of augmenting these methods using additional molecular data descriptors. We argue that molecular expansion of adverse event data may provide a path to improving the insights gained through more traditional pharmacovigilance approaches. Examples include the ability to assess statistical relevance with respect to underlying biomolecular mechanisms, the ability to generate plausible causative hypotheses and/or confirmation where possible, the ability to computationally study potential clinical trial designs and/or results, as well as the further provision of advanced features incorporated in innovative methods, such as machine learning. In summary, molecular data expansion provides an elegant way to extend mechanistic modeling, systems pharmacology, and patient-centered approaches for the assessment of drug safety. We anticipate that such advances in real-world data informatics and outcome analytics will help to better inform public health, via the improved ability to prospectively understand and predict various types of drug-induced molecular perturbations and adverse events.

Application of a patient-centered reverse translational systems-based approach to understand mechanisms of an adverse drug reaction of immune checkpoint inhibitors.

Immunotherapy became a key pillar of cancer therapeutics with the approvals of ipilimumab, nivolumab, and pembrolizumab, which inhibit either cytotoxic T-lymphocyte antigen-4 (CTLA-4) or programmed death-1 (PD-1) that are negative regulators of T-cell activation. However, boosting T-cell activation is often accompanied by autoimmunity, leading to adverse drug reactions (ADRs), including high grade 3-4 colitis and its severe complications whose prevalence may reach 14% for combination checkpoint inhibitors. In this research, we investigated how mechanistic differences between anti-CTLA-4 (ipilimumab) and anti-PD-1 (nivolumab and pembrolizumab) affect colitis, a general class toxicity. The data analytical platform Molecular Health Effect was utilized to map population ADR data from the US Food and Drug Administration (FDA) Adverse Event Reporting System to chemical and biological databases for hypothesis generation regarding the underlying molecular mechanisms causing colitis. Disproportionality analysis was used to assess the statistical relevance between adverse events of interest and molecular causation. We verified that the anti-CTLA-4 drug is associated with an approximately three-fold higher proportional reporting ratio associated with colitis than those of the anti-PD-1 drugs. The signal of the molecular mechanisms, including signaling pathways of inflammatory cytokines, was statistically insignificant to test the hypothesis that the severer rate of colitis associated with ipilimumab would be due to a greater magnitude of T-cell activation as a result of earlier response of the anti-CTLA-4 drug in the immune response. This patient-centered systems-based approach provides an exploratory process to better understand drug pair adverse events at pathway and target levels through reverse translation from postmarket surveillance safety reports.

A case study of a patient-centered reverse translational systems-based approach to understand adverse event profiles in drug development.

Adverse drug reactions (ADRs) of targeted therapy drugs (TTDs) are frequently unexpected and long-term toxicities detract from exceptional efficacy of new TTDs. In this proof-of-concept study, we explored how molecular causation involved in trastuzumab-induced cardiotoxicity changes when trastuzumab was given in combination with doxorubicin, tamoxifen, paroxetine, or lapatinib. The data analytical platform Molecular Health Effect was utilized to map population ADR data from the US Food and Drug Administration (FDA) Adverse Event Reporting System to chemical and biological databases (such as UniProt and Reactome), for hypothesis generation regarding the underlying molecular mechanisms causing cardiotoxicity. Disproportionality analysis was used to assess the statistical relevance between adverse events of interest and molecular causation. Literature search was performed to compare the established hypotheses to published experimental findings. We found that the combination therapy of trastuzumab and doxorubicin may affect mitochondrial dysfunction in cardiomyocytes through different molecular pathways such as BCL-X and PGC-1α proteins, leading to a synergistic effect of cardiotoxicity. We found, on the other hand, that trastuzumab-induced cardiotoxicity would be diminished by concomitant use of tamoxifen, paroxetine, and/or lapatinib. Tamoxifen and paroxetine may cause less cardiotoxicity through an increase in antioxidant activities, such as glutathione conjugation. Lapatinib may decrease the apoptotic effects in cardiomyocytes by altering the effects of trastuzumab on BCL-X proteins. This patient-centered systems-based approach provides, based on the trastuzumab-induced ADR cardiotoxicity, an example of how to apply reverse translation to investigate ADRs at the molecular pathway and target level to understand the causality and prevalence during drug development of novel therapeutics.

Brexpiprazole Pharmacokinetics in CYP2D6 Poor Metabolizers: Using Physiologically Based Pharmacokinetic Modeling to Optimize Time to Effective Concentrations.

Brexpiprazole is an oral antipsychotic agent indicated for use in patients with schizophrenia, or as adjunctive treatment for major depressive disorder. As cytochrome P450 (CYP) 2D6 contributes significantly to brexpiprazole metabolism, there is a label-recommended 50% reduction in dose among patients with the CYP2D6 poor metabolizer phenotype. This study uses a whole-body physiologically based pharmacokinetic (PBPK) model to compare the pharmacokinetics of brexpiprazole in patients known to be extensive metabolizers (EMs) and poor metabolizers (PMs). A PBPK model was constructed, verified, and validated against brexpiprazole clinical data, and simulations of 500 subjects were performed to establish the median time to effective concentrations in EMs and PMs. The PBPK simulations captured brexpiprazole PK well and demonstrated significant differences in the time to effective concentrations between EMs and PMs according to the label-recommended titration. Additionally, these simulations suggest that CYP2D6 PMs consistently achieve lower minimum concentrations during the dosing interval than CYP2D6 EMs. Simulations using an alternative dosing strategy of twice-daily dosing (as opposed to once daily) in PMs during the first week of brexpiprazole dosing yielded more consistent plasma concentrations between EMs and PMs, without exceeding the area under the plasma concentration-time curve observed in the EMs. Taken together, the results of these PBPK simulations suggest that product labeling for brexpiprazole titration in CYP2D6 PMs likely overcompensates for the decreased clearance seen in this population. We propose an alternative dosing strategy that decreases the time to effective concentrations and recommend a reevaluation of steady-state PK in this population to potentially allow for higher daily doses in CYP2D6 PMs.

Dosing Regimen Prediction and Confirmation With Rivaroxaban for Thromboprophylaxis in Children After the Fontan Procedure: Insights From the Phase III UNIVERSE Study.

Thrombosis remains an important complication for children with single-ventricle physiology following the Fontan procedure, and effective thromboprophylaxis is an important unmet medical need. To obviate conventional dose-finding studies and expedite clinical development, a rivaroxaban dose regimen for this indication was determined using a model-informed drug development approach. A physiologically based pharmacokinetic rivaroxaban model was used to predict a pediatric dosing regimen that would produce drug exposures similar to that of 10 mg once daily in adults. This regimen was used in an open-label, multicenter phase III study, which investigated the use of rivaroxaban for thromboprophylaxis in post-Fontan patients 2 to 8 years of age. The pharmacokinetics (PK) of rivaroxaban was assessed in part A (n = 12) and in part B (n = 64) of the UNIVERSE study. The safety and efficacy in the rivaroxaban group were compared to those in the acetylsalicylic acid group for 12 months. Pharmacodynamic end points were assessed in both parts of the study. Rivaroxaban exposures achieved in parts A and B were similar to the adult reference exposures. Prothrombin time also showed similarity to the adult reference. Exposure-response analysis did not identify a quantitative relationship between rivaroxaban exposures and efficacy/safety outcomes within the observed exposure ranges. A body weight-based dose regimen selected by physiologically based pharmacokinetic modeling was shown in the UNIVERSE study to be appropriate for thromboprophylaxis in the post-Fontan pediatric population. Model-based dose selection can support pediatric drug development and bridge adult dose data to pediatrics, thereby obviating the need for dose-finding studies in pediatric programs.

Pharmacodynamic Modeling to Evaluate the Impact of Cimetidine, an OCT2 Inhibitor, on the Anticancer Effects of Cisplatin.

Despite potent anticancer activity, the clinical utilization of cisplatin is limited due to nephrotoxicity. As Organic Cation Transporter 2 (OCT2) has been shown to be one of the key transporters involved in the uptake of cisplatin into renal proximal tubules, OCT2 inhibitors such as cimetidine have been explored to suppress cisplatin-induced nephrotoxicity. Nonetheless, the impact of OCT2 inhibition or cimetidine on the anti-cancer effects of cisplatin has not been extensively examined. The main objective of the present study was to quantitatively characterize the anticancer effects of cisplatin and cimetidine and determine their nature of interactions in two cancer cell lines, OCT2-negative hepatocellular carcinoma (HCC) cell line, Huh7, and OCT2-positive breast cancer cell line, MDA-MB-468. First, we determined the static concentration-response curves of cisplatin and cimetidine as single agents. Next, with the help of three-dimensional (3D) response surface analyses and a competitive interaction model, we determined their nature of interactions at static concentrations to be modestly synergistic or additive in Huh7 and antagonistic in MDA-MB-468. These results were consistent with the cell-level pharmacodynamic (PD) modeling analysis which leveraged the time-course effects of drugs as single agents and drug combinations. Our developed PD model can be further used to design future preclinical studies to further investigate the cisplatin and cimetidine combinations in different in vitro and in vivo cancer models.

Perspective on model-informed drug development.

Model-informed drug development (MIDD) is a process intended to expedite drug development, enhance regulatory science, and produce benefits for patients. Quantitative modeling and simulation-principally by population pharmacokinetics (PK), exposure-response, and physiologically based pharmacokinetic (PBPK) analysis-is the technology that provides the capability to deploy MIDD across a range of applications. MIDD was codified in the 2017 Prescription Drug User Fee Act Reauthorization 1 (PDUFA VI, 2018-2022) and a performance goal was a MIDD pilot program to hold 2 to 4 industry-U.S. Food and Drug Administration (FDA) paired meetings quarterly through 2022.

Physiologically-based pharmacokinetics modeling to investigate formulation factors influencing the generic substitution of dabigatran etexilate.

The exposure-response relationship of direct acting oral anti-coagulants (DOACs) for bleeding risk is steep relative to ischemic stroke reduction. As a result, small changes in exposure may lead to bleeding events. The overall goal of this project was to determine the effect of critical formulation parameters on the pharmacokinetics (PKs) and thus safety and efficacy of generic DOACs. In this first installment of our overall finding, we developed and verified a physiologically-based PK (PBPK) model for dabigatran etexilate (DABE) and its metabolites. The model was developed following a middle out approach leveraging available in vitro and in vivo data. External validity of the model was confirmed by overlapping predicted and observed PK profiles for DABE as well as free and total dabigatran for a dataset not used during model development. The verified model was applied to interrogate the impact of modulating the microenvironment pH on DABE systemic exposure. The PBPK exploratory analyses highlighted the high sensitivity of DABE exposure to supersaturation ratio and precipitation kinetics.

Model-Informed Precision Dosing: Background, Requirements, Validation, Implementation, and Forward Trajectory of Individualizing Drug Therapy.

Model-informed precision dosing (MIPD) has become synonymous with modern approaches for individualizing drug therapy, in which the characteristics of each patient are considered as opposed to applying a one-size-fits-all alternative. This review provides a brief account of the current knowledge, practices, and opinions on MIPD while defining an achievable vision for MIPD in clinical care based on available evidence. We begin with a historical perspective on variability in dose requirements and then discuss technical aspects of MIPD, including the need for clinical decision support tools, practical validation, and implementation of MIPD in health care. We also discuss novel ways to characterize patient variability beyond the common perceptions of genetic control. Finally, we address current debates on MIPD from the perspectives of the new drug development, health economics, and drug regulations.

A Review on Drug-Induced Nephrotoxicity: Pathophysiological Mechanisms, Drug Classes, Clinical Management, and Recent Advances in Mathematical Modeling and Simulation Approaches.

A variety of marketed drugs belonging to various therapeutic classes are known to cause nephrotoxicity. Nephrotoxicity can manifest itself in several forms depending on the specific site involved as well as the underlying pathophysiological mechanisms. As they often coexist with other pathophysiological conditions, the steps that can be taken to treat them are often limited. Thus, drug-induced nephrotoxicity remains a major clinical challenge. Prior knowledge of risk factors associated with special patient populations and specific classes of drugs, combined with early diagnosis, therapeutic drug monitoring with dose adjustments, as well as timely prospective treatments are essential to prevent and manage them better. Most incident drug-induced renal toxicity is reversible only if diagnosed at an early stage and treated promptly. Hence, diagnosis at an early stage is the need of the hour to counter it. Significant recent advances in the identification of novel early biomarkers of nephrotoxicity are not beyond limitations. In such a scenario, mathematical modeling and simulation (M&S) approaches may help to better understand and predict toxicities in a clinical setting. This review summarizes pathophysiological mechanisms of drug-induced nephrotoxicity, classes of nephrotoxic drugs, management, prevention, and diagnosis in clinics. Finally, it also highlights some of the recent advancements in mathematical M&S approaches that could be used to better understand and predict drug-induced nephrotoxicity.

Quantitative Systems Pharmacology Model-Based Predictions of Clinical Endpoints to Optimize Warfarin and Rivaroxaban Anti-Thrombosis Therapy.

BACKGROUND: Tight monitoring of efficacy and safety of anticoagulants such as warfarin is imperative to optimize the benefit-risk ratio of anticoagulants in patients. The standard tests used are measurements of prothrombin time (PT), usually expressed as international normalized ratio (INR), and activated partial thromboplastin time (aPTT). OBJECTIVE: To leverage a previously developed quantitative systems pharmacology (QSP) model of the human coagulation network to predict INR and aPTT for warfarin and rivaroxaban, respectively. METHODS: A modeling and simulation approach was used to predict INR and aPTT measurements of patients receiving steady-state anticoagulation therapy. A previously developed QSP model was leveraged for the present analysis. The effect of genetic polymorphisms known to influence dose response of warfarin (CYP2C9, VKORC1) were implemented into the model by modifying warfarin clearance (CYP2C9 *1: 0.2 L/h; *2: 0.14 L/h, *3: 0.04 L/h) and the concentration of available vitamin K (VKORC1 GA: -22% from normal vitamin K concentration; AA: -44% from normal vitamin K concentration). Virtual patient populations were used to assess the ability of the model to accurately predict routine INR and aPTT measurements from patients under long-term anticoagulant therapy. RESULTS: The introduced model accurately described the observed INR of patients receiving long-term warfarin treatment. The model was able to demonstrate the influence of genetic polymorphisms of CYP2C9 and VKORC1 on the INR measurements. Additionally, the model was successfully used to predict aPTT measurements for patients receiving long-term rivaroxaban therapy. CONCLUSION: The QSP model accurately predicted INR and aPTT measurements observed during routine therapeutic drug monitoring. This is an exemplar of how a QSP model can be adapted and used as a model-based precision dosing tool during clinical practice and drug development to predict efficacy and safety of anticoagulants to ultimately help optimize anti-thrombotic therapy.

A clinical population pharmacokinetic/pharmacodynamic model for BIIB059, a monoclonal antibody for the treatment of systemic and cutaneous lupus erythematosus.

A population pharmacokinetic/pharmacodynamic (popPK/PD) model for BIIB059 (anti-blood dendritic cell antigen 2 [anti-BDCA2]), a humanized immunoglobulin G1 monoclonal antibody currently under development for the treatment of SLE and CLE, is presented. BIIB059 binds BDCA2, a plasmacytoid dendritic cell (pDC)-specific receptor that inhibits the production of IFN-I and other inflammatory mediators when ligated. Phase 1 PK and PD data of healthy adult volunteers (HV, n = 87) and SLE subjects (n = 22) were utilized for the development of the popPK/PD model. The data included single and multiple dosing of intravenous and subcutaneous BIIB059. BDCA2 internalization (PD marker) was measured for all subjects by monitoring reduction of BDCA2 on pDC cell surface and used for development of the popPD model. A two-compartment popPK model with linear plus non-linear elimination was found to best describe BIIB059 PK. BDCA2 levels were best captured using an indirect response model with stimulation of the elimination of BDCA2. Clearance in SLE subjects was 25% higher compared to HV (6.87 vs 5.52 mL/h). Bodyweight was identified as only other covariate on clearance and central volume. The estimates of EC50 and Emax were 0.35 µg/mL and 8.92, respectively. No difference in EC50 and Emax was observed between SLE and HV. The popPK/PD model described the data accurately, as evaluated by pcVPCs and bootstrap. The presented popPK/PD model for BIIB059 provides valuable insight into the dynamics and dose-response relationship of BIIB059 for the treatment of SLE and CLE and was used to guide dose selection for the Phase 2 clinical study (NCT02847598).

Semi-mechanistic modelling platform to assess cardiac contractility and haemodynamics in preclinical cardiovascular safety profiling of new molecular entities.

BACKGROUND AND PURPOSE: Cardiovascular safety is one of the most frequent causes of safety-related attrition both preclinically and clinically. Preclinical cardiovascular safety is routinely assessed using dog telemetry monitoring key cardiovascular functions. The present research was to develop a semi-mechanistic modelling platform to simultaneously assess changes in contractility (dPdtmax ), heart rate (HR) and mean arterial pressure (MAP) in preclinical studies. EXPERIMENTAL APPROACH: Data from dPdtmax , HR, preload (left ventricular end-diastolic pressure [LVEDP]) and MAP were available from dog telemetry studies after dosing with atenolol (n = 27), salbutamol (n = 5), L-NG -nitroarginine methyl ester (L-NAME; n = 4), milrinone (n = 4), verapamil (n = 12), dofetilide (n = 8), flecainide (n = 4) and AZ001 (n = 14). Literature model for rat CV function was used for the structural population pharmacodynamic model development. LVEDP was evaluated as covariate to account for the effect of preload on dPdtmax . KEY RESULTS: The model was able to describe drug-induced changes in dPdtmax , HR and MAP for all drugs included in the developed framework adequately, by incorporating appropriate drug effects on dPdtmax , HR and/or total peripheral resistance. Consistent with the Starling's law, incorporation of LVEDP as a covariate on dPdtmax to correct for the preload effect was found to be statistically significant. CONCLUSIONS AND IMPLICATIONS: The contractility and haemodynamics semi-mechanistic modelling platform accounts for diurnal variation, drug-induced changes and inter-animal variation. It can be used to hypothesize and evaluate pharmacological effects and provide a holistic cardiovascular safety profile for new drugs.

Association of Oral Anticoagulants and Verapamil or Diltiazem With Adverse Bleeding Events in Patients With Nonvalvular Atrial Fibrillation and Normal Kidney Function.

Importance: Direct oral anticoagulants (DOACs) are purported to have fewer drug-drug interactions than warfarin. However, potential interactions with coprescribed medications are still a safety concern. Verapamil hydrochloride and diltiazem hydrochloride are combined P-glycoprotein (P-gp) and CYP3A4 inhibitors and may be associated with increases in the risk of bleeding with DOACs. Objective: To evaluate the risk of bleeding with DOACs and verapamil or diltiazem using an active comparator design. Design, Setting, and Participants: A comparative effectiveness active comparator cohort study was conducted using US population-based data (2010-2015) analyzed between January 1 and July 15, 2019. Data were obtained on 48 442 patients with nonvalvular atrial fibrillation who had received an index prescription of dabigatran, rivaroxaban, or apixaban between October 19, 2010, through June 30, 2015, with final follow-up on October 1, 2015. Analysis was restricted to individuals with no history of kidney disease who were receiving standard doses of the DOACs. Exposures: Patients with initial prescriptions of DOACs who were receiving verapamil or diltiazem were compared with those receiving amlodipine or metoprolol. Main Outcomes and Measures: Overall and gastrointestinal major, moderate, and minor bleeding using primary or secondary diagnoses. Hazard ratios and 95% CIs were estimated using inverse probability of treatment weights in Cox proportional hazards regression models. Results: Of the 48 442 patients reviewed, analysis was conducted on 1764 patients receiving DOACs with verapamil or diltiazem compared with 3105 receiving amlodipine and 1793 patients receiving DOACs with verapamil or diltiazem compared with 3224 receiving metoprolol. Depending on the comparison, approximately 60% of the cohort were younger than 65 years and male, which differed by treatment group. Rivaroxaban and apixaban were not associated with increased rates of bleeding for patients receiving verapamil or diltiazem compared with those receiving amlodipine or metoprolol. Among patients receiving dabigatran etexilate, the overall bleeding rate was 52% higher (hazard ratio, 1.52; 95% CI, 1.05-2.20) with verapamil or diltiazem vs amlodipine and 43% higher (hazard ratio, 1.43; 95% CI, 1.02-2.00) vs metoprolol. Bleeding rates for dabigatran with verapamil or diltiazem were higher overall for other bleeding types (244.9 vs 158.4 per 1000 person-years; adjusted hazard ratios of overall GI bleeding: 2.16; 95% CI, 1.30-3.60; minor bleeding: 1.56; 95% CI, 1.07-2.27; and minor GI bleeding: 2.16; 95% CI, 1.29-3.63). Sensitivity analyses showed consistent results for dabigatran when used with verapamil and diltiazem, with magnitudes ranging from 50% to 100% increased hazard rates and no significant results for apixaban or rivaroxaban. Conclusions and Relevance: Current US prescribing information only recommends prescribing changes with dabigatran and P-gp inhibitors with lower kidney function. This study found increased bleeding risk associated with dabigatran when used concomitantly with the P-gp inhibitors verapamil and diltiazem in individuals with normal kidney function. Clinicians and patients may need to consider these drug-drug interactions when choosing oral anticoagulation.

Study of pharmacogenomic information in FDA-approved drug labeling to facilitate application of precision medicine.

Pharmacogenomics (PGx), studying the relationship between drug response and genetic makeup of an individual, is accelerating advances in precision medicine. The FDA includes PGx information in the labeling of approved drugs to better inform on their safety and effectiveness. We herein present a summary of PGx information found in 261 prescription drug labeling documents by querying the publicly available FDALabel database. A total of 362 drug-biomarker pairs (DBPs) were identified. We profiled DBPs using frequency of the biomarkers and their therapeutic classes. Four categories of applications (indication, safety, dosing and information) were discussed according to information in labeling. This analysis facilitates better understanding, utilization and translation of PGx information in drug labeling among researchers, healthcare professionals and the public.