Search

Tyrosine kinase inhibitors (TKIs) are increasingly popular drugs used to treat more than a dozen different diseases including some forms of cancer. Despite having fewer adverse effects than traditional chemotherapies, they are not without risks. Liver injury is a particular concern. Of the FDA-approved TKIs, approximately 40% cause hepatotoxicity. However, little is known about the underlying pathophysiology. The leading hypothesis is that TKIs are converted by cytochrome P450 3A4 (CYP3A4) to reactive metabolites that damage proteins. Indeed, there is strong evidence for this bioactivation of TKIs in in vitro reactions. However, the actual toxic effects are underexplored. Here, we measured the cytotoxicity of several TKIs in primary mouse hepatocytes, HepaRG cells and HepG2 cells with and without CYP3A4 modulation. To our surprise, the data indicate that CYP3A4 increases resistance to sorafenib and lapatinib hepatotoxicity. The results have implications for the mechanism of toxicity of these drugs in patients and underline the importance of selecting an appropriate experimental model.

Meclofenamate is a nonsteroidal anti-inflammatory drug used in the treatment of mild-to-moderate pain yet poses a rare risk of hepatotoxicity through an unknown mechanism. Nonsteroidal anti-inflammatory drug (NSAID) bioactivation is a common molecular initiating event for hepatotoxicity. Thus, we hypothesized a similar mechanism for meclofenamate and leveraged computational and experimental approaches to identify and characterize its bioactivation. Analyses employing our XenoNet model indicated possible pathways to meclofenamate bioactivation into 19 reactive metabolites subsequently trapped into glutathione adducts. We describe the first reported bioactivation kinetics for meclofenamate and relative importance of those pathways using human liver microsomes. The findings validated only four of the many bioactivation pathways predicted by modeling. For experimental studies, dansyl glutathione was a critical trap for reactive quinone metabolites and provided a way to characterize adduct structures by mass spectrometry and quantitate yields during reactions. Of the four quinone adducts, we were able to characterize structures for three of them. Based on kinetics, the most efficient bioactivation pathway led to the monohydroxy para-quinone-imine followed by the dechloro-ortho-quinone-imine. Two very inefficient pathways led to the dihydroxy ortho-quinone and a likely multiply adducted quinone. When taken together, bioactivation pathways for meclofenamate accounted for approximately 13% of total metabolism. In sum, XenoNet facilitated prediction of reactive metabolite structures, whereas quantitative experimental studies provided a tractable approach to validate actual bioactivation pathways for meclofenamate. Our results provide a foundation for assessing reactive metabolite load more accurately for future comparative studies with other NSAIDs and drugs in general. Significance Statement Meclofenamate bioactivation may initiate hepatotoxicity, yet common risk assessment approaches are often cumbersome and inefficient and yield qualitative insights that do not scale relative bioactivation risks. We developed and applied innovative computational modeling and quantitative kinetics to identify and validate meclofenamate bioactivation pathways and relevance as a function of time and concentration. This strategy yielded novel insights on meclofenamate bioactivation and provides a tractable approach to more accurately and efficiently assess other drug bioactivations and correlate risks to toxicological outcomes.

The 12(th) International Society for the Study of Xenobiotics (ISSX) meeting, held in Portland, OR, USA from July 28–31, 2019, was attended by diverse members of the pharmaceutical sciences community. The ISSX New Investigators Group provides learning and professional growth opportunities for student and early career members of ISSX. To share meeting content with those who were unable to attend, the ISSX New Investigators herein elected to highlight the “Advances in the Study of Drug Metabolism” symposium, as it engaged attendees with diverse backgrounds. This session covered a wide range of current topics in drug metabolism research including predicting sites and routes of metabolism, metabolite identification, ligand docking, and medicinal and natural products chemistry, and highlighted approaches complemented by computational modeling. In silico tools have been increasingly applied in both academic and industrial settings, alongside traditional and evolving in vitro techniques, to strengthen and streamline pharmaceutical research. Approaches such as quantum mechanics simulations facilitate understanding of reaction energetics towards prediction of routes and sites of drug metabolism. Furthermore, in tandem with crystallographic and orthogonal wet lab techniques for structural validation of drug metabolizing enzymes, in silico models can aid understanding of substrate recognition by particular enzymes, identify metabolic soft spots and predict toxic metabolites for improved molecular design. Of note, integration of chemical synthesis and biosynthesis using natural products remains an important approach for identifying new chemical scaffolds in drug discovery. These subjects, compiled by the symposium organizers, presenters, and the ISSX New Investigators Group, are discussed in this review.

Diphenylamine NSAIDs are taken frequently for chronic pain conditions, yet their use may potentiate hepatotoxicity risks through poorly characterized metabolic mechanisms. Our previous work revealed that seven marketed or withdrawn diphenylamine NSAIDs undergo bioactivation into quinone-species metabolites, whose reaction specificities depended on halogenation and the type of acidic group on the diphenylamine. Herein, we identified cytochromes P450 responsible for those bioactivations, determined reaction specificities, and estimated relative contributions of enzymes to overall hepatic bioactivations and detoxifications. A qualitative activity screen revealed CYP2C8, 2C9, 2C19, and 3A4 played roles in drug bioactivation. Subsequent steady-state studies with recombinant CYPs recapitulated the importance of halogenation and acidic group type on bioactivations but importantly, showed patterns unique to each CYP. CYP2C9, 2C19 and 3A4 bioactivated all NSAIDs with CYP2C9 dominating all possible bioactivation pathways and their specificities. For each CYP, specificities for overall oxidative metabolism were not impacted significantly by differences in NSAID structures but the values themselves differed among the enzymes such that CYP2C9 and 3A4 were more efficient than others. When considering hepatic CYP abundance, CYP2C9 almost exclusively accounted for diphenylamine NSAID bioactivations, whereas CYP3A4 provided a critical counterbalance favoring their overall detoxification. Preference for either outcome would depend on molecular structures favoring metabolism by the CYPs as well as the influence of clinical factors altering their expression and/or activity. While focused on NSAIDs, these findings have broader implications on bioactivation risks given the expansion of the diphenylamine scaffold to other drug classes such as targeted cancer therapeutics.

The bromodomain and extra terminal (BET) protein family recognizes acetylated lysines within histones and transcription factors using two N-terminal bromodomains, D1 and D2. The protein–protein interactions between BET bromodomains, acetylated histones, and transcription factors are therapeutic targets for BET-related diseases, including inflammatory disease and cancer. Prior work demonstrated that methylated-1,2,3-triazoles are suitable N-acetyl lysine mimetics for BET inhibition. Here we describe a structure–activity relationship study of triazole-based inhibitors that improve affinity, D1 selectivity, and microsomal stability. These outcomes were accomplished by targeting a nonconserved residue, Asp144 and a conserved residue, Met149, on BRD4 D1. The lead inhibitors DW34 and 26 have a BRD4 D1 K(d) of 12 and 6.4 nM, respectively. Cellular activity was demonstrated through suppression of c-Myc expression in MM.1S cells and downregulation of IL-8 in TNF-α-stimulated A549 cells. These data indicate that DW34 and 26 are new leads to investigate the anticancer and anti-inflammatory activity of BET proteins.

In 2020, nearly one-third of new drugs on the global market were synthetic cannabinoids including the drug of abuse N-(1-adamantyl)-1-(5-pentyl)-1H-indazole-3-carboxamide (5F-APINACA, 5F-AKB48). Knowledge of 5F-APINACA metabolism provides a critical mechanistic basis to interpret and predict abuser outcomes. Prior qualitative studies identified which metabolic processes occur but not the order and extent of them and often relied on problematic &ldquo, semi-quantitative&rdquo, mass spectroscopic (MS) approaches. We capitalized on 5F-APINACA absorbance for quantitation while leveraging MS to characterize metabolite structures for measuring 5F-APINACA steady-state kinetics. We demonstrated the reliability of absorbance and not MS for inferring metabolite levels. Human liver microsomal reactions yielded eight metabolites by MS but only five by absorbance. Subsequent kinetic studies on primary and secondary metabolites revealed highly efficient mono- and dihydroxylation of the adamantyl group and much less efficient oxidative defluorination at the N-pentyl terminus. Based on regiospecificity and kinetics, we constructed pathways for competing and intersecting steps in 5F-APINACA metabolism. Overall efficiency for adamantyl oxidation was 17-fold higher than that for oxidative defluorination, showing significant bias in metabolic flux and subsequent metabolite profile compositions. Lastly, our analytical approach provides a powerful new strategy to more accurately assess metabolic kinetics for other understudied synthetic cannabinoids possessing the indazole chromophore.

Meloxicam is a thiazole-containing NSAID that was approved for marketing with favorable clinical outcomes despite being structurally similar to the hepatotoxic sudoxicam. Introduction of a single methyl group on the thiazole results in an overall lower toxic risk, yet the group’s impact on P450 isozyme bioactivation is unclear. Through analytical methods, we used inhibitor phenotyping and recombinant P450s to identify contributing P450s, and then measured steady-state kinetics for bioactivation of sudoxicam and meloxicam by the recombinant P450s to determine relative efficiencies. Experiments showed that CYP2C8, 2C19, and 3A4 catalyze sudoxicam bioactivation, and CYP1A2 catalyzes meloxicam bioactivation, indicating that the methyl group not only impacts enzyme affinity for the drugs, but also alters which isozymes catalyze the metabolic pathways. Scaling of relative P450 efficiencies based on average liver concentration revealed that CYP2C8 dominates the sudoxicam bioactivation pathway and CYP2C9 dominates meloxicam detoxification. Dominant P450s were applied for an informatics assessment of electronic health records to identify potential correlations between meloxicam drug-drug interactions and drug-induced liver injury. Overall, our findings provide a cautionary tale on assumed impacts of even simple structural modifications on drug bioactivation while also revealing specific targets for clinical investigations of predictive factors that determine meloxicam-induced idiosyncratic liver injury.

The 3,5-dimethylisoxazole motif has become a useful and popular acetyl-lysine mimic employed in isoxazole-containing bromodomain and extra-terminal (BET) inhibitors but may introduce the potential for bioactivations into toxic reactive metabolites. As a test, we coupled deep neural models for quinone formation, metabolite structures, and biomolecule reactivity to predict bioactivation pathways for 32 BET inhibitors and validate the bioactivation of select inhibitors experimentally. Based on model predictions, inhibitors were more likely to undergo bioactivation than reported non-bioactivated molecules containing isoxazoles. The model outputs varied with substituents indicating the ability to scale their impact on bioactivation. We selected OXFBD02, OXFBD04, and I-BET151 for more in-depth analysis. OXFBD’s bioactivations were evenly split between traditional quinones and novel extended quinone-methides involving the isoxazole yet strongly favored the latter quinones. Subsequent experimental studies confirmed the formation of both types of quinones for OXFBD molecules, yet traditional quinones were the dominant reactive metabolites. Modeled I-BET151 bioactivations led to extended quinone-methides, which were not verified experimentally. The differences in observed and predicted bioactivations reflected the need to improve overall bioactivation scaling. Nevertheless, our coupled modeling approach predicted BET inhibitor bioactivations including novel extended quinone methides, and we experimentally verified those pathways highlighting potential concerns for toxicity in the development of these new drug leads.

Thiazoles are biologically active aromatic heterocyclic rings occurring frequently in natural products and drugs. These molecules undergo typically harmless elimination; however, a hepatotoxic response can occur due to multistep bioactivation of the thiazole to generate a reactive thioamide. A basis for those differences in outcomes remains unknown. A textbook example is the high hepatotoxicity observed for sudoxicam in contrast to the relative safe use and marketability of meloxicam, which differs in structure from sudoxicam by the addition of a single methyl group. Both drugs undergo bioactivation, but meloxicam exhibits an additional detoxification pathway due to hydroxylation of the methyl group. We hypothesized that thiazole bioactivation efficiency is similar between sudoxicam and meloxicam due to the methyl group being a weak electron donator, and thus, the relevance of bioactivation depends on the competing detoxification pathway. For a rapid analysis, we modeled epoxidation of sudoxicam derivatives to investigate the impact of substituents on thiazole bioactivation. As expected, electron donating groups increased the likelihood for epoxidation with a minimal effect for the methyl group, but model predictions did not extrapolate well among all types of substituents. Through analytical methods, we measured steady-state kinetics for metabolic bioactivation of sudoxicam and meloxicam by human liver microsomes. Sudoxicam bioactivation was 6-fold more efficient than that for meloxicam, yet meloxicam showed a 6-fold higher efficiency of detoxification than bioactivation. Overall, sudoxicam bioactivation was 15-fold more likely than meloxicam considering all metabolic clearance pathways. Kinetic differences likely arise from different enzymes catalyzing respective metabolic pathways based on phenotyping studies. Rather than simply providing an alternative detoxification pathway, the meloxicam methyl group suppressed the bioactivation reaction. These findings indicate the impact of thiazole substituents on bioactivation is more complex than previously thought and likely contributes to the unpredictability of their toxic potential.