Revisiting Aldehyde Oxidase Mediated Metabolism in Drug-like Molecules: An Improved Computational Model

Transfer learning enhanced graph neural network for aldehyde oxidase metabolism prediction and its experimental application

Aldehyde oxidase (AOX) is a molybdoenzyme that is primarily expressed in the liver and is involved in the metabolism of drugs and other xenobiotics. AOX-mediated metabolism can result in unexpected outcomes, such as the production of toxic metabolites and high metabolic clearance, which can lead to the clinical failure of novel therapeutic agents. Computational models can assist medicinal chemists in rapidly evaluating the AOX metabolic risk of compounds during the early phases of drug discovery and provide valuable clues for manipulating AOX-mediated metabolism liability. In this study, we developed a novel graph neural network called AOMP for predicting AOX-mediated metabolism. AOMP integrated the tasks of metabolic substrate/non-substrate classification and metabolic site prediction, while utilizing transfer learning from 13C nuclear magnetic resonance data to enhance its performance on both tasks. AOMP significantly outperformed the benchmark methods in both cross-validation and external testing. Using AOMP, we systematically assessed the AOX-mediated metabolism of common fragments in kinase inhibitors and successfully identified four new scaffolds with AOX metabolism liability, which were validated through in vitro experiments. Furthermore, for the convenience of the community, we established the first online service for AOX metabolism prediction based on AOMP, which is freely available at https://aomp.alphama.com.cn.

INTEDE 2.0: the metabolic roadmap of drugs

The metabolic roadmap of drugs (MRD) is a comprehensive atlas for understanding the stepwise and sequential metabolism of certain drug in living organisms. It plays a vital role in lead optimization, personalized medication, and ADMET research. The MRD consists of three main components: (i) the sequential catalyses of drug and its metabolites by different drug-metabolizing enzymes (DMEs), (ii) a comprehensive collection of metabolic reactions along the entire MRD and (iii) a systematic description on efficacy & toxicity for all metabolites of a studied drug. However, there is no database available for describing the comprehensive metabolic roadmaps of drugs. Therefore, in this study, a major update of INTEDE was conducted, which provided the stepwise & sequential metabolic roadmaps for a total of 4701 drugs, and a total of 22 165 metabolic reactions containing 1088 DMEs and 18 882 drug metabolites. Additionally, the INTEDE 2.0 labeled the pharmacological properties (pharmacological activity or toxicity) of metabolites and provided their structural information. Furthermore, 3717 drug metabolism relationships were supplemented (from 7338 to 11 055). All in all, INTEDE 2.0 is highly expected to attract broad interests from related research community and serve as an essential supplement to existing pharmaceutical/biological/chemical databases. INTEDE 2.0 can now be accessible freely without any login requirement at: http://idrblab.org/intede/.

Large Computational Survey of Intrinsic Reactivity of Aromatic Carbon Atoms with Respect to a Model Aldehyde Oxidase

Aldehyde oxidase (AOX) and other related molybdenum-containing enzymes are known to oxidize the C-H bonds of aromatic rings. This process contributes to the metabolism of pharmaceutical compounds and, therefore, is of vital importance to drug pharmacokinetics. The present work describes an automated computational workflow and its use for the prediction of intrinsic reactivity of small aromatic molecules toward a minimal model of the active site of AOX. The workflow is based on quantum chemical transition state searches for the underlying single-step oxidation reaction, where the automated protocol includes identification of unique aromatic C-H bonds, creation of three-dimensional reactant and product complex geometries via a templating approach, search for a transition state, and validation of reaction end points. Conformational search on the reactants, products, and the transition states is performed. The automated procedure has been validated on previously reported transition state barriers and was used to evaluate the intrinsic reactivity of nearly three hundred heterocycles commonly found in approved drug molecules. The intrinsic reactivity of more than 1000 individual aromatic carbon sites is reported. Stereochemical and conformational aspects of the oxidation reaction, which have not been discussed in previous studies, are shown to play important roles in accurate modeling of the oxidation reaction. Observations on structural trends that determine the reactivity are provided and rationalized.

In Silico Prediction of Metabolic Reaction Catalyzed by Human Aldehyde Oxidase

Aldehyde oxidase (AOX) plays an important role in drug metabolism. Human AOX (hAOX) is widely distributed in the body, and there are some differences between species. Currently, animal models cannot accurately predict the metabolism of hAOX. Therefore, more and more in silico models have been constructed for the prediction of the hAOX metabolism. These models are based on molecular docking and quantum chemistry theory, which are time-consuming and difficult to automate. Therefore, in this study, we compared traditional machine learning methods, graph convolutional neural network methods, and sequence-based methods with limited data, and proposed a ligand-based model for the metabolism prediction catalyzed by hAOX. Compared with the published models, our model achieved better performance (ACC = 0.91, F1 = 0.77). What's more, we built a web server to predict the sites of metabolism (SOMs) for hAOX. In summary, this study provides a convenient and automatable model and builds a web server named Meta-hAOX for accelerating the drug design and optimization stage.

Consensus Prediction of Human Gut Microbiota-Mediated Metabolism Susceptibility for Small Molecules by Machine Learning, Structural Alerts, and Dietary Compounds-Based Average Similarity Methods

The human gut microbiota (HGM) colonizing human gastrointestinal tract (HGT) confers a repertoire of dynamic and unique metabolic capacities that are not possessed by the host and therefore is tentatively perceived as an alternative metabolic ″organ″ besides the liver in the host. Nevertheless, the significant contribution of HGM to the overall human metabolism is often overlooked in the modern drug discovery pipeline. Hence, a systematic evaluation of HGM-mediated drug metabolism is gradually important, and its computational prediction becomes increasingly necessary. In this work, a new data set containing both the HGM-mediated metabolism susceptible (HGMMS) and insusceptible (HGMMI) compounds (329 vs 320) was manually curated. Based on this data set, the first machine learning (ML) model, a new structural alerts (SA) model, and the K-nearest neighboring dietary compounds-based average similarity (AS) model were proposed to directly predict the HGM-mediated metabolism susceptibility for small molecules, and exhibit promising performance on three independent test sets. Finally, consensus prediction (ML/SA/AS) for DrugBank molecules revealed an intriguing phenomenon that a typical Michael acceptor ″α,β-unsaturated carbonyl group″ is a very common warhead for the design of covalent inhibitors and inclined to be metabolized by HGM in anaerobic HGT to generate the reduced metabolite without the reactive warhead, which could be a new concern to medicinal chemists. To the best of our knowledge, we gleaned the first HGMMS/HGMMI data set, developed the first HGMMS/HGMMI classification model, implemented a relatively comprehensive program based on ML/SA/AS approaches, and found a new phenomenon on the HGM-mediated deactivation of an extensively used warhead for covalent inhibitors.

Interrogating the Inhibition Mechanisms of Human Aldehyde Oxidase by X-ray Crystallography and NMR Spectroscopy: The Raloxifene Case

Human aldehyde oxidase (hAOX1) is mainly present in the liver and has an emerging role in drug metabolism, since it accepts a wide range of molecules as substrates and inhibitors. Herein, we employed an integrative approach by combining NMR, X-ray crystallography, and enzyme inhibition kinetics to understand the inhibition modes of three hAOX1 inhibitors-thioridazine, benzamidine, and raloxifene. These integrative data indicate that thioridazine is a noncompetitive inhibitor, while benzamidine presents a mixed type of inhibition. Additionally, we describe the first crystal structure of hAOX1 in complex with raloxifene. Raloxifene binds tightly at the entrance of the substrate tunnel, stabilizing the flexible entrance gates and elucidating an unusual substrate-dependent mechanism of inhibition with potential impact on drug-drug interactions. This study can be considered as a proof-of-concept for an efficient experimental screening of prospective substrates and inhibitors of hAOX1 relevant in drug discovery.

Recent developments in predicting CYP-independent metabolism

As lead optimization efforts have successfully reduced metabolic liabilities due to cytochrome P450 (CYP)-mediated metabolism, there has been an increase in the frequency of involvement of non-CYP enzymes in the metabolism of investigational compounds. Although there have been numerous notable advancements in the characterization of non-CYP enzymes with respect to their localization, reaction mechanisms, species differences and identification of typical substrates, accurate prediction of non-CYP-mediated clearance, with a particular emphasis with the difficulties in accounting for any extrahepatic contributions, remains a challenge. The current manuscript comprehensively summarizes the recent advancements in the prediction of drug metabolism and the in vitro to in vitro extrapolation of clearance for substrates of non-CYP drug metabolizing enzymes.