1
|
Chang J, Fan X, Tian B. DeepP450: Predicting Human P450 Activities of Small Molecules by Integrating Pretrained Protein Language Model and Molecular Representation. J Chem Inf Model 2024; 64:3149-3160. [PMID: 38587937 DOI: 10.1021/acs.jcim.4c00115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
Cytochrome P450 enzymes (CYPs) play a crucial role in Phase I drug metabolism in the human body, and CYP activity toward compounds can significantly affect druggability, making early prediction of CYP activity and substrate identification essential for therapeutic development. Here, we established a deep learning model for assessing potential CYP substrates, DeepP450, by fine-tuning protein and molecule pretrained models through feature integration with cross-attention and self-attention layers. This model exhibited high prediction accuracy (0.92) on the test set, with area under the receiver operating characteristic curve (AUROC) values ranging from 0.89 to 0.98 in substrate/nonsubstrate predictions across the nine major human CYPs, surpassing current benchmarks for CYP activity prediction. Notably, DeepP450 uses only one model to predict substrates/nonsubstrates for any of the nine CYPs and exhibits certain generalizability on novel compounds and different categories of human CYPs, which could greatly facilitate early stage drug design by avoiding CYP-reactive compounds.
Collapse
Affiliation(s)
- Jiamin Chang
- MOE Key Laboratory of Bioinformatics, State Key Laboratory of Molecular Oncology, School of Pharmaceutical Sciences, Tsinghua University, Beijing 100084, China
| | - Xiaoyu Fan
- MOE Key Laboratory of Bioinformatics, State Key Laboratory of Molecular Oncology, School of Pharmaceutical Sciences, Tsinghua University, Beijing 100084, China
| | - Boxue Tian
- MOE Key Laboratory of Bioinformatics, State Key Laboratory of Molecular Oncology, School of Pharmaceutical Sciences, Tsinghua University, Beijing 100084, China
| |
Collapse
|
2
|
Feng Y, Gong C, Zhu J, Liu G, Tang Y, Li W. Prediction of Sites of Metabolism of CYP3A4 Substrates Utilizing Docking-Derived Geometric Features. J Chem Inf Model 2023. [PMID: 37336765 DOI: 10.1021/acs.jcim.3c00549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
Cytochrome P450 3A4 (CYP3A4) is one of the major drug-metabolizing enzymes in the human body and is responsible for the metabolism of ∼50% of clinically used drugs. Therefore, the identification of the compound's sites of metabolism (SOMs) mediated by CYP3A4 is of utmost importance in the early stage of drug discovery and development. Herein, docking-based approaches incorporating geometric features were used for SOMs prediction of CYP3A4 substrates. The cross-docking poses of a relatively large data set containing 474 substrates were analyzed in depth, and a widely observed geometric pattern called the close proximity of SOMs was derived from the poses. On the basis of the close proximity, several structure-based models have been constructed, which demonstrated better performance than those structure-based models using the criterion of Fe-SOM distance. For further improving the prediction performance, the structure-based models were also combined with the well-known ligand-based model SMARTCyp. One combined model exhibited good performance on the SOMs prediction of an external substrate set containing kinase inhibitors, PROTACs, approved drugs, and some lead compounds.
Collapse
Affiliation(s)
- Yanjun Feng
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Changda Gong
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Jieyu Zhu
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Guixia Liu
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Yun Tang
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Weihua Li
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| |
Collapse
|
3
|
Guengerich FP. Roles of cytochrome P450 enzymes in pharmacology and toxicology: Past, present, and future. ADVANCES IN PHARMACOLOGY (SAN DIEGO, CALIF.) 2022; 95:1-47. [PMID: 35953152 PMCID: PMC9869358 DOI: 10.1016/bs.apha.2021.12.001] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The development of the cytochrome P450 (P450) field has been remarkable in the areas of pharmacology and toxicology, particularly in drug development. Today it is possible to use the knowledge base and relatively straightforward assays to make intelligent predictions about drug disposition prior to human dosing. Much is known about the structures, regulation, chemistry of catalysis, and the substrate and inhibitor specificity of human P450s. Many aspects of drug-drug interactions and side effects can be understood in terms of P450s. This knowledge has also been useful in pharmacy practice, as well as in the pharmaceutical industry and medical practice. However, there are still basic and practical questions to address regarding P450s and their roles in pharmacology and toxicology. Another aspect is the discovery of drugs that inhibit P450 to treat diseases.
Collapse
|
4
|
Conan M, Théret N, Langouet S, Siegel A. Constructing xenobiotic maps of metabolism to predict enzymes catalyzing metabolites capable of binding to DNA. BMC Bioinformatics 2021; 22:450. [PMID: 34548010 PMCID: PMC8454073 DOI: 10.1186/s12859-021-04363-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 08/28/2021] [Indexed: 12/22/2022] Open
Abstract
Background The liver plays a major role in the metabolic activation of xenobiotics (drugs, chemicals such as pollutants, pesticides, food additives...). Among environmental contaminants of concern, heterocyclic aromatic amines (HAA) are xenobiotics classified by IARC as possible or probable carcinogens (2A or 2B). There exist little information about the effect of these HAA in humans. While HAA is a family of more than thirty identified chemicals, the metabolic activation and possible DNA adduct formation have been fully characterized in human liver for only a few of them (MeIQx, PhIP, A\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha$$\end{document}αC). Results We have developed a modeling approach in order to predict all the possible metabolites of a xenobiotic and enzymatic profiles that are linked to the production of metabolites able to bind DNA. Our prediction of metabolites approach relies on the construction of an enriched and annotated map of metabolites from an input metabolite.The pipeline assembles reaction prediction tools (SyGMa), sites of metabolism prediction tools (Way2Drug, SOMP and Fame 3), a tool to estimate the ability of a xenobotics to form DNA adducts (XenoSite Reactivity V1), and a filtering procedure based on Bayesian framework. This prediction pipeline was evaluated using caffeine and then applied to HAA. The method was applied to determine enzymes profiles associated with the maximization of metabolites derived from each HAA which are able to bind to DNA. The classification of HAA according to enzymatic profiles was consistent with their chemical structures. Conclusions Overall, a predictive toxicological model based on an in silico systems biology approach opens perspectives to estimate the genotoxicity of various chemical classes of environmental contaminants. Moreover, our approach based on enzymes profile determination opens the possibility of predicting various xenobiotics metabolites susceptible to bind to DNA in both normal and physiopathological situations. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04363-6.
Collapse
Affiliation(s)
- Mael Conan
- Institut de Recherche en Santé, Environnement et Travail, Univ Rennes, Inserm, EHESP, IRSET, Rennes, France.,Institut de Recherche en Informatique et Systèmes Aléatoires, Univ Rennes, Inria, CNRS, IRISA, Rennes, France
| | - Nathalie Théret
- Institut de Recherche en Santé, Environnement et Travail, Univ Rennes, Inserm, EHESP, IRSET, Rennes, France.,Institut de Recherche en Informatique et Systèmes Aléatoires, Univ Rennes, Inria, CNRS, IRISA, Rennes, France
| | - Sophie Langouet
- Institut de Recherche en Santé, Environnement et Travail, Univ Rennes, Inserm, EHESP, IRSET, Rennes, France.
| | - Anne Siegel
- Institut de Recherche en Informatique et Systèmes Aléatoires, Univ Rennes, Inria, CNRS, IRISA, Rennes, France.
| |
Collapse
|
5
|
He Z, Wang Z, Gao B, Liu S, Zhao X, Shi H, Wang M. Stereostructure-activity mechanism of cyproconazole by cytochrome P450 in rat liver microsomes: A combined experimental and computational study. JOURNAL OF HAZARDOUS MATERIALS 2021; 416:125764. [PMID: 33827004 DOI: 10.1016/j.jhazmat.2021.125764] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 03/22/2021] [Accepted: 03/23/2021] [Indexed: 06/12/2023]
Abstract
Cyproconazole (CPZ), representing the chiral triazole fungicides, is widely used in the pharmaceutical and agricultural fields. To clarify its potential adverse effects on the generalized CYP-mediated processes within mammalian, a comparative experimental and computational approach was employed to investigate the CYP-mediated metabolism processes of CPZ stereoisomers in rat liver microsomes (RLMs). The depletion rate of CPZ stereoisomers in vitro incubation system with RLMs followed the order RR-> SS-> SR-> RS-CPZ. The results of kinetic assays were in line with the depletion rate results. Further inhibition assay confirmed the stereoselective metabolism of CPZ stereoisomers by different CYP isoforms. Molecular dynamics (MD) simulation revealed the stereoselective metabolism mechanism. Several hydrogen bonds and π-stacking restrict the position of CPZ isomers in the active cavity of CYPs so that the 4'-nitrogen on the triazole ring can bind closely to the heme of CYP, which results in the metabolism of CPZ isomers. By combining the computational and experimental approaches, the structure-activity relationship of CPZ and CYP was elucidated, and this method can be further applied to predict the degree of uncertainty in the process of xenobiotic biotransformation of triazole fungicides and serve as a basis for risk assessment.
Collapse
Affiliation(s)
- Zongzhe He
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China
| | - Zhen Wang
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China
| | - Beibei Gao
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China; Toxicological Center, University of Antwerp, Wilrijk, Belgium
| | - Shiling Liu
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China
| | - Xuejun Zhao
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China
| | - Haiyan Shi
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China
| | - Minghua Wang
- Department of Pesticide Science, College of Plant Protection, Nanjing Agricultural University, State & Local Joint Engineering Research Center of Green Pesticide Invention and Application, Nanjing 210095, China.
| |
Collapse
|
6
|
Bruno A, Costantino G, Sartori L, Radi M. The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization. Curr Med Chem 2019; 26:3838-3873. [PMID: 29110597 DOI: 10.2174/0929867324666171107101035] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2017] [Revised: 09/27/2017] [Accepted: 09/28/2017] [Indexed: 01/04/2023]
Abstract
BACKGROUND Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug. METHODS In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan. RESULTS A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given. CONCLUSION The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.
Collapse
Affiliation(s)
- Agostino Bruno
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Gabriele Costantino
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| | - Luca Sartori
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Marco Radi
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| |
Collapse
|
7
|
Diao X, Huestis MA. New Synthetic Cannabinoids Metabolism and Strategies to Best Identify Optimal Marker Metabolites. Front Chem 2019; 7:109. [PMID: 30886845 PMCID: PMC6409358 DOI: 10.3389/fchem.2019.00109] [Citation(s) in RCA: 76] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Accepted: 02/11/2019] [Indexed: 11/21/2022] Open
Abstract
Synthetic cannabinoids (SCs) were initially developed as pharmacological tools to probe the endocannabinoid system and as novel pharmacotherapies, but are now highly abused. This is a serious public health and social problem throughout the world and it is highly challenging to identify which SC was consumed by the drug abusers, a necessary step to tie adverse health effects to the new drug's toxicity. Two intrinsic properties complicate SC identification, their often rapid and extensive metabolism, and their generally high potency relative to the natural psychoactive Δ9-tetrahydrocannabinol in cannabis. Additional challenges are the lack of reference standards for the major urinary metabolites needed for forensic verification, and the sometimes differing illicit and licit status and, in some cases, identical metabolites produced by closely related SC pairs, i.e., JWH-018/AM-2201, THJ-018/THJ-2201, and BB-22/MDMB-CHMICA/ADB-CHMICA. We review current SC prevalence, establish the necessity for SC metabolism investigation and contrast the advantages and disadvantages of multiple metabolic approaches. The human hepatocyte incubation model for determining a new SC's metabolism is highly recommended after comparison to human liver microsomes incubation, in silico prediction, rat in vivo, zebrafish, and fungus Cunninghamella elegans models. We evaluate SC metabolic patterns, and devise a practical strategy to select optimal urinary marker metabolites for SCs. New SCs are incubated first with human hepatocytes and major metabolites are then identified by high-resolution mass spectrometry. Although initially difficult to obtain, authentic human urine samples following the specified SC exposure are hydrolyzed and analyzed by high-resolution mass spectrometry to verify identified major metabolites. Since some SCs produce the same major urinary metabolites, documentation of the specific SC consumed may require identification of the SC parent itself in either blood or oral fluid. An encouraging trend is the recent reduction in the number of new SC introduced per year. With global collaboration and communication, we can improve education of the public about the toxicity of new SC and our response to their introduction.
Collapse
Affiliation(s)
- Xingxing Diao
- Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Marilyn A. Huestis
- The Lambert Center for the Study of Medicinal Cannabis and Hemp, Institute for Emerging Health Professions, Thomas Jefferson University, Philadelphia, PA, United States
| |
Collapse
|
8
|
Saikia S, Bordoloi M, Sarmah R, Kolita B. Antiviral compound screening, peptide designing, and protein network construction of influenza a virus (strain a/Puerto Rico/8/1934 H1N1). Drug Dev Res 2018; 80:106-124. [PMID: 30276835 DOI: 10.1002/ddr.21475] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 09/03/2018] [Accepted: 09/05/2018] [Indexed: 11/11/2022]
Abstract
Plant-based antiviral therapy is the current need for holistic health care management, which can be achieved through screening of phytochemicals and designing of antiviral peptides. There exist certain host's factors which are directly involved for rapid viral replication causing worldwide pandemic. A total of 177 phytochemicals from Ocimum sanctum (L.), Tinospora cordifolia (Thunb.) Miers, Cinnamomum camphora (L.) J. Presl., Allium sativum (L.), Curcuma longa (L.), and Aloe vera (L.) Burm. f. were evaluated for their affinity to all viral proteins of H1N1. Applying drug filters and keeping the threshold of such filters relative to the standards, 82 compounds were found suitable for further analysis. Consensus scoring system was used for screening top ligands from 82 compounds, which screened the top 12 compounds. Highly conserved regions (>80%) which were hydrophilic, flexible, antigenic, and also charged were screened out as potent antiviral peptides. The viral proteins were taken as the targets for the modeled peptides for protein-protein docking. Further, host-pathogen interacting network was constructed to unveil host factors involved in viral replication, from which unique protein clusters representing their involvement in viral reproduction were selected through mapping with pathway databases. Twelve compounds and five peptides were found to be highly effective against all the proteins of H1N1. Based on the uniqueness, 13 clusters of proteins were obtained which are engaged in cellular process, namely, viral reproduction, fructose-6-phosphate metabolism, nitrogen compound metabolism, biosynthesis, cellular process, oligodendrocyte development, localization, multiorganism process, primary metabolism, response to unfolded protein, metabolism, and response to protein and catabolism.
Collapse
Affiliation(s)
- Surovi Saikia
- Chemical Sciences & Technology (Natural Products Chemistry) Division, CSIR North East Institute of Science & Technology, Jorhat, Assam, India
| | - Manobjyoti Bordoloi
- Chemical Sciences & Technology (Natural Products Chemistry) Division, CSIR North East Institute of Science & Technology, Jorhat, Assam, India
| | - Rajeev Sarmah
- Allied Health Sciences, Assam Down Town University, Panikhaiti, Guwahati, Assam, India
| | - Bhaskor Kolita
- Chemical Sciences & Technology (Natural Products Chemistry) Division, CSIR North East Institute of Science & Technology, Jorhat, Assam, India
| |
Collapse
|
9
|
Huang X, Groves JT. Oxygen Activation and Radical Transformations in Heme Proteins and Metalloporphyrins. Chem Rev 2018; 118:2491-2553. [PMID: 29286645 PMCID: PMC5855008 DOI: 10.1021/acs.chemrev.7b00373] [Citation(s) in RCA: 591] [Impact Index Per Article: 98.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Indexed: 12/20/2022]
Abstract
As a result of the adaptation of life to an aerobic environment, nature has evolved a panoply of metalloproteins for oxidative metabolism and protection against reactive oxygen species. Despite the diverse structures and functions of these proteins, they share common mechanistic grounds. An open-shell transition metal like iron or copper is employed to interact with O2 and its derived intermediates such as hydrogen peroxide to afford a variety of metal-oxygen intermediates. These reactive intermediates, including metal-superoxo, -(hydro)peroxo, and high-valent metal-oxo species, are the basis for the various biological functions of O2-utilizing metalloproteins. Collectively, these processes are called oxygen activation. Much of our understanding of the reactivity of these reactive intermediates has come from the study of heme-containing proteins and related metalloporphyrin compounds. These studies not only have deepened our understanding of various functions of heme proteins, such as O2 storage and transport, degradation of reactive oxygen species, redox signaling, and biological oxygenation, etc., but also have driven the development of bioinorganic chemistry and biomimetic catalysis. In this review, we survey the range of O2 activation processes mediated by heme proteins and model compounds with a focus on recent progress in the characterization and reactivity of important iron-oxygen intermediates. Representative reactions initiated by these reactive intermediates as well as some context from prior decades will also be presented. We will discuss the fundamental mechanistic features of these transformations and delineate the underlying structural and electronic factors that contribute to the spectrum of reactivities that has been observed in nature as well as those that have been invented using these paradigms. Given the recent developments in biocatalysis for non-natural chemistries and the renaissance of radical chemistry in organic synthesis, we envision that new enzymatic and synthetic transformations will emerge based on the radical processes mediated by metalloproteins and their synthetic analogs.
Collapse
Affiliation(s)
- Xiongyi Huang
- Department
of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
- Department
of Chemistry, California Institute of Technology, Pasadena, California 91125, United States
| | - John T. Groves
- Department
of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
10
|
Meng J, Li S, Liu X, Zheng M, Li H. RD-Metabolizer: an integrated and reaction types extensive approach to predict metabolic sites and metabolites of drug-like molecules. Chem Cent J 2017; 11:65. [PMID: 29086838 PMCID: PMC5515729 DOI: 10.1186/s13065-017-0290-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2017] [Accepted: 07/03/2017] [Indexed: 11/10/2022] Open
Abstract
Background Experimental approaches for determining the metabolic properties of the drug candidates are usually expensive, time-consuming and labor intensive. There is a great deal of interest in developing computational methods to accurately and efficiently predict the metabolic decomposition of drug-like molecules, which can provide decisive support and guidance for experimentalists. Results Here, we developed an integrated, low false positive and reaction types extensive metabolism prediction approach called RD-Metabolizer (Reaction Database-based Metabolizer). RD-Metabolizer firstly employed the detailed reaction SMARTS patterns to encode different metabolism reaction types with the aim of covering larger chemical reaction space. 2D fingerprint similarity calculation model was built to calculate the metabolic probability of each site in a molecule. RDKit was utilized to act on pre-written reaction SMARTS patterns to correct the metabolic ranking of each site in a molecule generated by the 2D fingerprint similarity calculation model as well as generate corresponding structures of metabolites, thus helping to reduce the false positive metabolites. Two test sets were adopted to evaluate the performance of RD-Metabolizer in predicting SOMs and structures of metabolites. The results indicated that RD-Metabolizer was better than or at least as good as several widely used SOMs prediction methods. Besides, the number of false positive metabolites was obviously reduced compared with MetaPrint2D-React. Conclusions The accuracy and efficiency of RD-Metabolizer was further illustrated by a metabolism prediction case of AZD9291, which is a mutant-selective EGFR inhibitor. RD-Metabolizer will serve as a useful toolkit for the early metabolic properties assessment of drug-like molecules at the preclinical stage of drug discovery.A visual example of the metabolic site and the corresponding metabolite of Chloroquine predicted by RD-Metabolizer ![]()
Collapse
Affiliation(s)
- Jiajia Meng
- State Key Laboratory of Bioreactor Engineering, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 200237, China
| | - Shiliang Li
- State Key Laboratory of Bioreactor Engineering, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 200237, China.,Shanghai Key Laboratory of Chemical Biology, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai, 200237, China
| | - Xiaofeng Liu
- State Key Laboratory of Bioreactor Engineering, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 200237, China
| | - Mingyue Zheng
- Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Honglin Li
- State Key Laboratory of Bioreactor Engineering, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 200237, China.
| |
Collapse
|
11
|
Dixit VA, Deshpande S. Advances in Computational Prediction of Regioselective and Isoform-Specific Drug Metabolism Catalyzed by CYP450s. ChemistrySelect 2016. [DOI: 10.1002/slct.201601051] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Vaibhav A. Dixit
- Department of Pharmaceutical Chemistry; School of Pharmacy and Technology Management (SPTM), Shri Vile Parle Kelavani Mandal's (SVKM's) Narsee Monjee Institute of Management Studies (NMIMS), Mukesh Patel Technology Park, Babulde, Bank of Tapi River; Mumbai-Agra Road Shirpur, Dist. Dhule−425405 India
| | - Shirish Deshpande
- Department of Pharmaceutical Chemistry; School of Pharmacy and Technology Management (SPTM), Shri Vile Parle Kelavani Mandal's (SVKM's) Narsee Monjee Institute of Management Studies (NMIMS), Mukesh Patel Technology Park, Babulde, Bank of Tapi River; Mumbai-Agra Road Shirpur, Dist. Dhule−425405 India
| |
Collapse
|
12
|
Metabolite identification of bentysrepinine (Y101), a novel anti-HBV agent in rats using a five-step strategy based on a combined workflow with two different platforms of liquid chromatography-tandem mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci 2016; 1040:118-128. [PMID: 27978466 DOI: 10.1016/j.jchromb.2016.12.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2016] [Revised: 11/27/2016] [Accepted: 12/05/2016] [Indexed: 01/05/2023]
Abstract
Bentysrepinine (Y101), a derivative of repensine (a compound isolated from Dichondrarepens Forst), is a novel phenyalanine dipeptide inhibiting DNA-HBV and cccDNA activities and is currently under development for the treatment of hepatitis B virus (HBV)-infected hepatitis. Our previous study implied that there might be an existence of extensive metabolism of Y101 in rats. Therefore, it is necessary to perform metabolic profiling study to further evaluate its safety and drug-like properties. In this study, the metabolism of Y101 in rats was investigated by a convincible five-step strategy to characterize metabolites in plasma and that excreted into urine, bile and feces. The five-step strategy was realized by using an combined workflow on two different MS platforms, including various scan modes of liquid chromatography with hybrid quadruple-linear ion trap mass spectrometry (LC-QTRAP-MS/MS) and various post-acquiring data mining tools of liquid chromatography-quadrupole time-of-flight tandem mass spectrometry (LC-QTOF-MS/MS). QTOF MS/MS was employed as a powerful complementary tool to enable high confidence of metabolites identification using its functions of accurate MS and MS/MS fragmentation. As a result, a total of 30 metabolites were detected, including 25 phase I and 5 phase II metabolites. Among them, four primary metabolites (M6-M9) were further identified by comparing with the authentic standards chemically synthesized. The possible metabolic pathways of Y101 in rats were proposed to be amide hydrolysis, monohydroxylation, dihydroxylation, N-oxidation, demethylation, methylation, glucosidation and glucuronidation. This is the first study of the metabolism of Y101 in rats. The five-step strategy was successfully used to systematically characterize metabolites of Y101 in rats, and it would be generally applied for metabolite identification of new drug candidate.
Collapse
|
13
|
Diao X, Huestis MA. Approaches, Challenges, and Advances in Metabolism of New Synthetic Cannabinoids and Identification of Optimal Urinary Marker Metabolites. Clin Pharmacol Ther 2016; 101:239-253. [DOI: 10.1002/cpt.534] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2016] [Revised: 09/16/2016] [Accepted: 10/07/2016] [Indexed: 12/15/2022]
Affiliation(s)
- X Diao
- Department of Chemistry and Drug Metabolism, IRP, National Institute on Drug Abuse; National Institutes of Health; Baltimore Maryland USA
| | - MA Huestis
- Department of Chemistry and Drug Metabolism, IRP, National Institute on Drug Abuse; National Institutes of Health; Baltimore Maryland USA
- University of Maryland School of Medicine; Baltimore Maryland USA
| |
Collapse
|
14
|
Mukherjee G, Lal Gupta P, Jayaram B. Predicting the binding modes and sites of metabolism of xenobiotics. MOLECULAR BIOSYSTEMS 2016; 11:1914-24. [PMID: 25913019 DOI: 10.1039/c5mb00118h] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Metabolism studies are an essential integral part of ADMET profiling of drug candidates to evaluate their safety and efficacy. Cytochrome P-450 (CYP) metabolizes a wide variety of xenobiotics/drugs. The binding modes of these compounds with CYP and their intrinsic reactivities decide the metabolic products. We report here a novel computational protocol, which comprises docking of ligands to heme-containing CYPs and prediction of binding energies through a newly developed scoring function, followed by analyses of the docked structures and molecular orbitals of the ligand molecules, for predicting the sites of metabolism (SOM) of ligands. The calculated binding free energies of 121 heme-containing protein-ligand docked complexes yielded a correlation coefficient of 0.84 against experiment. Molecular orbital analyses of the resultant top three unique poses of the docked complexes provided a success rate of 87% in identifying the experimentally known sites of metabolism of the xenobiotics. The SOM prediction methodology is freely accessible at .
Collapse
Affiliation(s)
- Goutam Mukherjee
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi-110016, India.
| | | | | |
Collapse
|
15
|
Perryman AL, Stratton TP, Ekins S, Freundlich JS. Predicting Mouse Liver Microsomal Stability with "Pruned" Machine Learning Models and Public Data. Pharm Res 2016; 33:433-49. [PMID: 26415647 PMCID: PMC4712113 DOI: 10.1007/s11095-015-1800-5] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2015] [Accepted: 09/22/2015] [Indexed: 02/07/2023]
Abstract
PURPOSE Mouse efficacy studies are a critical hurdle to advance translational research of potential therapeutic compounds for many diseases. Although mouse liver microsomal (MLM) stability studies are not a perfect surrogate for in vivo studies of metabolic clearance, they are the initial model system used to assess metabolic stability. Consequently, we explored the development of machine learning models that can enhance the probability of identifying compounds possessing MLM stability. METHODS Published assays on MLM half-life values were identified in PubChem, reformatted, and curated to create a training set with 894 unique small molecules. These data were used to construct machine learning models assessed with internal cross-validation, external tests with a published set of antitubercular compounds, and independent validation with an additional diverse set of 571 compounds (PubChem data on percent metabolism). RESULTS "Pruning" out the moderately unstable / moderately stable compounds from the training set produced models with superior predictive power. Bayesian models displayed the best predictive power for identifying compounds with a half-life ≥1 h. CONCLUSIONS Our results suggest the pruning strategy may be of general benefit to improve test set enrichment and provide machine learning models with enhanced predictive value for the MLM stability of small organic molecules. This study represents the most exhaustive study to date of using machine learning approaches with MLM data from public sources.
Collapse
Affiliation(s)
- Alexander L Perryman
- Division of Infectious Disease, Department of Medicine, and the Ruy V. Lourenço Center for the Study of Emerging and Re-emerging Pathogens, Rutgers University-New Jersey Medical School, Newark, New Jersey, 07103, USA
| | - Thomas P Stratton
- Department of Pharmacology & Physiology, Rutgers University-New Jersey Medical School, Medical Sciences Building, I-503, 185 South Orange Ave., Newark, New Jersey, 07103, USA
| | - Sean Ekins
- Collaborations in Chemistry, 5616 Hilltop Needmore Road, Fuquay-Varina, NC, 27526, USA
- Collaborative Drug Discovery, 1633 Bayshore Highway, Suite 342, Burlingame, CA, 94010, USA
| | - Joel S Freundlich
- Division of Infectious Disease, Department of Medicine, and the Ruy V. Lourenço Center for the Study of Emerging and Re-emerging Pathogens, Rutgers University-New Jersey Medical School, Newark, New Jersey, 07103, USA.
- Department of Pharmacology & Physiology, Rutgers University-New Jersey Medical School, Medical Sciences Building, I-503, 185 South Orange Ave., Newark, New Jersey, 07103, USA.
| |
Collapse
|
16
|
Li J, Cai J, Su H, Du H, Zhang J, Ding S, Liu G, Tang Y, Li W. Effects of protein flexibility and active site water molecules on the prediction of sites of metabolism for cytochrome P450 2C19 substrates. MOLECULAR BIOSYSTEMS 2016; 12:868-78. [DOI: 10.1039/c5mb00784d] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Structure-based prediction of sites of metabolism (SOMs) mediated by cytochrome P450s (CYPs) is of great interest in drug discovery and development.
Collapse
Affiliation(s)
- Junhao Li
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Jinya Cai
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Haixia Su
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Hanwen Du
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Juan Zhang
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Shihui Ding
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Guixia Liu
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Yun Tang
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Weihua Li
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| |
Collapse
|
17
|
Jurva U, Weidolf L. Electrochemical generation of drug metabolites with applications in drug discovery and development. Trends Analyt Chem 2015. [DOI: 10.1016/j.trac.2015.04.010] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
18
|
Olsen L, Oostenbrink C, Jørgensen FS. Prediction of cytochrome P450 mediated metabolism. Adv Drug Deliv Rev 2015; 86:61-71. [PMID: 25958010 DOI: 10.1016/j.addr.2015.04.020] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2014] [Revised: 03/30/2015] [Accepted: 04/27/2015] [Indexed: 10/23/2022]
Abstract
Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.
Collapse
|
19
|
Scharkoi O, Becker R, Esslinger S, Weber M, Nehls I. Predicting sites of cytochrome P450-mediated hydroxylation applied to CYP3A4 and hexabromocyclododecane. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2014.898845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
20
|
Shah MB, Wilderman PR, Liu J, Jang HH, Zhang Q, Stout CD, Halpert JR. Structural and biophysical characterization of human cytochromes P450 2B6 and 2A6 bound to volatile hydrocarbons: analysis and comparison. Mol Pharmacol 2015; 87:649-59. [PMID: 25585967 PMCID: PMC4366795 DOI: 10.1124/mol.114.097014] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2014] [Accepted: 01/12/2015] [Indexed: 11/22/2022] Open
Abstract
X-ray crystal structures of complexes of cytochromes CYP2B6 and CYP2A6 with the monoterpene sabinene revealed two distinct binding modes in the active sites. In CYP2B6, sabinene positioned itself with the putative oxidation site located closer to the heme iron. In contrast, sabinene was found in an alternate conformation in the more compact CYP2A6, where the larger hydrophobic side chains resulted in a significantly reduced active-site cavity. Furthermore, results from isothermal titration calorimetry indicated a much more substantial contribution of favorable enthalpy to sabinene binding to CYP2B6 as opposed to CYP2A6, consistent with the previous observations with (+)-α-pinene. Structural analysis of CYP2B6 complexes with sabinene and the structurally similar (3)-carene and comparison with previously solved structures revealed how the movement of the F206 side chain influences the volume of the binding pocket. In addition, retrospective analysis of prior structures revealed that ligands containing -Cl and -NH functional groups adopted a distinct orientation in the CYP2B active site compared with other ligands. This binding mode may reflect the formation of Cl-π or NH-π bonds with aromatic rings in the active site, which serve as important contributors to protein-ligand binding affinity and specificity. Overall, the findings from multiple techniques illustrate how drugs metabolizing CYP2B6 and CYP2A6 handle a common hydrocarbon found in the environment. The study also provides insight into the role of specific functional groups of the ligand that may influence the binding to CYP2B6.
Collapse
Affiliation(s)
- Manish B Shah
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - P Ross Wilderman
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - Jingbao Liu
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - Hyun-Hee Jang
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - Qinghai Zhang
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - C David Stout
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| | - James R Halpert
- Department of Pharmaceutical Sciences, The University of Connecticut, Storrs, Connecticut (M.B.S., P.R.W., J.L., J.R.H.); School of Biological Sciences and Technology, Chonnam National University, Gwangju, Republic of Korea (H.-H.J.); and Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California (Q.Z., C.D.S.)
| |
Collapse
|
21
|
Oostenbrink C. Structure‐Based Methods for Predicting the Sites and Products of Metabolism. ACTA ACUST UNITED AC 2014. [DOI: 10.1002/9783527673261.ch10] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
22
|
Rydberg P. Reactivity‐Based Approaches and Machine Learning Methods for Predicting the Sites of Cytochrome P450‐Mediated Metabolism. ACTA ACUST UNITED AC 2014. [DOI: 10.1002/9783527673261.ch11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
23
|
Pragyan P, Kesharwani SS, Nandekar PP, Rathod V, Sangamwar AT. Predicting drug metabolism by CYP1A1, CYP1A2, and CYP1B1: insights from MetaSite, molecular docking and quantum chemical calculations. Mol Divers 2014; 18:865-78. [DOI: 10.1007/s11030-014-9534-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2014] [Accepted: 06/19/2014] [Indexed: 12/13/2022]
|
24
|
Cruciani G, Baroni M, Benedetti P, Goracci L, Fortuna CG. Exposition and reactivity optimization to predict sites of metabolism in chemicals. DRUG DISCOVERY TODAY. TECHNOLOGIES 2014; 10:e155-65. [PMID: 24050245 DOI: 10.1016/j.ddtec.2012.11.001] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Chemical modifications of drugs induced by phase I biotransformations significantly affect their pharmacokinetic properties. Because the metabolites produced can themselves have a pharmacological effect and an intrinsic toxicity, medicinal chemists need to accurately predict the sites of metabolism (SoM) of drugs as early as possible. However, site of metabolism prediction is rarely accompanied by a prediction of the relative abundance of the various metabolites. Such a prediction would be a great help in the study of drug– drug interactions and in the process of reducing the toxicity of potential drug candidates. The aim of this paper is to present recent developments in the prediction of xenobiotic metabolism and to use concrete examples to explain the computational mechanism employed.
Collapse
|
25
|
Brändén G, Sjögren T, Schnecke V, Xue Y. Structure-based ligand design to overcome CYP inhibition in drug discovery projects. Drug Discov Today 2014; 19:905-11. [PMID: 24642031 DOI: 10.1016/j.drudis.2014.03.012] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2013] [Revised: 02/26/2014] [Accepted: 03/11/2014] [Indexed: 01/01/2023]
Abstract
Cytochrome P450 (CYP) enzymes are key players in xenobiotic metabolism, and inhibition of CYPs can therefore result in unwanted drug-drug interactions. Within drug discovery, CYP inhibition can cause delays in the progression of candidate drugs, or even premature closure of projects. During the past decade, a massive effort in the pharmaceutical industry and academic research has produced a wealth of structural information in the CYP field. In this short review, we will describe how structure-based approaches can be used in the pharmaceutical industry to work away from CYP inhibition, with a focus on the opportunities and challenges. We will show two examples from our own work where structural information on CYP2C9 and CYP3A4 inhibitor complexes have been successfully exploited in ongoing drug discovery projects.
Collapse
Affiliation(s)
- Gisela Brändén
- Department of Chemistry and Molecular Biology, University of Gothenburg, Göteborg S-405 30, Sweden.
| | - Tove Sjögren
- Discovery Sciences, AstraZeneca R&D Mölndal, Mölndal S-431 83, Sweden
| | - Volker Schnecke
- CVMD iMed, AstraZeneca R&D Mölndal, Mölndal S-431 83, Sweden
| | - Yafeng Xue
- Discovery Sciences, AstraZeneca R&D Mölndal, Mölndal S-431 83, Sweden
| |
Collapse
|
26
|
Ekins S. Progress in computational toxicology. J Pharmacol Toxicol Methods 2013; 69:115-40. [PMID: 24361690 DOI: 10.1016/j.vascn.2013.12.003] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Accepted: 12/08/2013] [Indexed: 01/02/2023]
Abstract
INTRODUCTION Computational methods have been widely applied to toxicology across pharmaceutical, consumer product and environmental fields over the past decade. Progress in computational toxicology is now reviewed. METHODS A literature review was performed on computational models for hepatotoxicity (e.g. for drug-induced liver injury (DILI)), cardiotoxicity, renal toxicity and genotoxicity. In addition various publications have been highlighted that use machine learning methods. Several computational toxicology model datasets from past publications were used to compare Bayesian and Support Vector Machine (SVM) learning methods. RESULTS The increasing amounts of data for defined toxicology endpoints have enabled machine learning models that have been increasingly used for predictions. It is shown that across many different models Bayesian and SVM perform similarly based on cross validation data. DISCUSSION Considerable progress has been made in computational toxicology in a decade in both model development and availability of larger scale or 'big data' models. The future efforts in toxicology data generation will likely provide us with hundreds of thousands of compounds that are readily accessible for machine learning models. These models will cover relevant chemistry space for pharmaceutical, consumer product and environmental applications.
Collapse
Affiliation(s)
- Sean Ekins
- Collaborations in Chemistry, 5616 Hilltop Needmore Road, Fuquay Varina, NC 27526, USA; Department of Pharmaceutical Sciences, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA; Department of Pharmacology, Rutgers University-Robert Wood Johnson Medical School, 675 Hoes Lane, Piscataway, NJ 08854, USA; Division of Chemical Biology and Medicinal Chemistry, UNC Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, NC 27599-7355, USA.
| |
Collapse
|
27
|
Rydberg P, Jørgensen FS, Olsen L. Use of density functional theory in drug metabolism studies. Expert Opin Drug Metab Toxicol 2013; 10:215-27. [PMID: 24295134 DOI: 10.1517/17425255.2014.864278] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
INTRODUCTION The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes. AREAS COVERED We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom. EXPERT OPINION DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.
Collapse
Affiliation(s)
- Patrik Rydberg
- University of Copenhagen, Department of Drug Design and Pharmacology , Denmark
| | | | | |
Collapse
|
28
|
Tyzack JD, Williamson MJ, Torella R, Glen RC. Prediction of cytochrome P450 xenobiotic metabolism: tethered docking and reactivity derived from ligand molecular orbital analysis. J Chem Inf Model 2013; 53:1294-305. [PMID: 23701380 DOI: 10.1021/ci400058s] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Metabolism of xenobiotic and endogenous compounds is frequently complex, not completely elucidated, and therefore often ambiguous. The prediction of sites of metabolism (SoM) can be particularly helpful as a first step toward the identification of metabolites, a process especially relevant to drug discovery. This paper describes a reactivity approach for predicting SoM whereby reactivity is derived directly from the ground state ligand molecular orbital analysis, calculated using Density Functional Theory, using a novel implementation of the average local ionization energy. Thus each potential SoM is sampled in the context of the whole ligand, in contrast to other popular approaches where activation energies are calculated for a predefined database of molecular fragments and assigned to matching moieties in a query ligand. In addition, one of the first descriptions of molecular dynamics of cytochrome P450 (CYP) isoforms 3A4, 2D6, and 2C9 in their Compound I state is reported, and, from the representative protein structures obtained, an analysis and evaluation of various docking approaches using GOLD is performed. In particular, a covalent docking approach is described coupled with the modeling of important electrostatic interactions between CYP and ligand using spherical constraints. Combining the docking and reactivity results, obtained using standard functionality from common docking and quantum chemical applications, enables a SoM to be identified in the top 2 predictions for 75%, 80%, and 78% of the data sets for 3A4, 2D6, and 2C9, respectively, results that are accessible and competitive with other recently published prediction tools.
Collapse
Affiliation(s)
- Jonathan D Tyzack
- Unilever Centre for Molecular Science Informatics, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, United Kingdom
| | | | | | | |
Collapse
|
29
|
Artese A, Cross S, Costa G, Distinto S, Parrotta L, Alcaro S, Ortuso F, Cruciani G. Molecular interaction fields in drug discovery: recent advances and future perspectives. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1150] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Anna Artese
- Dipartimento di Scienze della Salute; Università degli Studi “Magna Graecia” di Catanzaro; Campus “S. Venuta”; Viale Europa Catanzaro Italy
| | - Simon Cross
- Molecular Discovery Ltd, Pinner; Middlesex London United Kingdom
| | - Giosuè Costa
- Dipartimento di Scienze della Salute; Università degli Studi “Magna Graecia” di Catanzaro; Campus “S. Venuta”; Viale Europa Catanzaro Italy
| | - Simona Distinto
- Dipartimento di Scienze della Vita e dell'Ambiente; Università di Cagliari; Cagliari Italy
| | - Lucia Parrotta
- Dipartimento di Scienze della Salute; Università degli Studi “Magna Graecia” di Catanzaro; Campus “S. Venuta”; Viale Europa Catanzaro Italy
| | - Stefano Alcaro
- Dipartimento di Scienze della Salute; Università degli Studi “Magna Graecia” di Catanzaro; Campus “S. Venuta”; Viale Europa Catanzaro Italy
| | - Francesco Ortuso
- Dipartimento di Scienze della Salute; Università degli Studi “Magna Graecia” di Catanzaro; Campus “S. Venuta”; Viale Europa Catanzaro Italy
| | - Gabriele Cruciani
- Laboratory for Chemometrics and Cheminformatics; Chemistry Department; University of Perugia; Perugia Italy
| |
Collapse
|
30
|
Leach AG. Tactics to Avoid Inhibition of Cytochrome P450s. TOPICS IN MEDICINAL CHEMISTRY 2013. [DOI: 10.1007/7355_2013_25] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
|
31
|
Hutzler JM, Obach RS, Dalvie D, Zientek MA. Strategies for a comprehensive understanding of metabolism by aldehyde oxidase. Expert Opin Drug Metab Toxicol 2012; 9:153-68. [DOI: 10.1517/17425255.2013.738668] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
32
|
Campagna-Slater V, Pottel J, Therrien E, Cantin LD, Moitessier N. Development of a computational tool to rival experts in the prediction of sites of metabolism of xenobiotics by p450s. J Chem Inf Model 2012; 52:2471-83. [PMID: 22916680 DOI: 10.1021/ci3003073] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The metabolism of xenobiotics--and more specifically drugs--in the liver is a critical process controlling their half-life. Although there exist experimental methods, which measure the metabolic stability of xenobiotics and identify their metabolites, developing higher throughput predictive methods is an avenue of research. It is expected that predicting the chemical nature of the metabolites would be an asset for designing safer drugs and/or drugs with modulated half-lives. We have developed IMPACTS (In-silico Metabolism Prediction by Activated Cytochromes and Transition States), a computational tool combining docking to metabolic enzymes, transition state modeling, and rule-based substrate reactivity prediction to predict the site of metabolism (SoM) of xenobiotics. Its application to sets of CYP1A2, CYP2C9, CYP2D6, and CYP3A4 substrates and comparison to experts' predictions demonstrates its accuracy and significance. IMPACTS identified an experimentally observed SoM in the top 2 predicted sites for 77% of the substrates, while the accuracy of biotransformation experts' prediction was 65%. Application of IMPACTS to external sets and comparison of its accuracy to those of eleven other methods further validated the method implemented in IMPACTS.
Collapse
Affiliation(s)
- Valérie Campagna-Slater
- Department of Chemistry, McGill University, 801 Sherbrooke St W, Montreal, QC H3A 0B8, Canada
| | | | | | | | | |
Collapse
|
33
|
Ballard P, Brassil P, Bui KH, Dolgos H, Petersson C, Tunek A, Webborn PJH. The right compound in the right assay at the right time: an integrated discovery DMPK strategy. Drug Metab Rev 2012; 44:224-52. [DOI: 10.3109/03602532.2012.691099] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
34
|
Shahrokh K, Cheatham TE, Yost GS. Conformational dynamics of CYP3A4 demonstrate the important role of Arg212 coupled with the opening of ingress, egress and solvent channels to dehydrogenation of 4-hydroxy-tamoxifen. Biochim Biophys Acta Gen Subj 2012; 1820:1605-17. [PMID: 22677141 DOI: 10.1016/j.bbagen.2012.05.011] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2012] [Revised: 05/22/2012] [Accepted: 05/23/2012] [Indexed: 12/12/2022]
Abstract
BACKGROUND Structure-based methods for P450 substrates are commonly used during drug development to identify sites of metabolism. However, docking studies using available X-ray structures for the major drug-metabolizing P450, CYP3A4, do not always identify binding modes supportive of the production of high-energy toxic metabolites. Minor pathways such as P450-catalyzed dehydrogenation have been experimentally shown to produce reactive products capable of forming biomolecular adducts which can lead to increased risk toxicities. 4-Hydroxy-tamoxifen (4OHT) is metabolized by CYP3A4 via competing hydroxylation and dehydrogenation reactions. METHODS Ab initio gas-phase electronic structural characterization of 4OHT was used to develop a docking scoring scheme. Conformational sampling of CYP3A4 with molecular dynamics simulations along multiple trajectories were used to generate representative structures for docking studies using recently published heme parameters. A key predicted binding mode was tested experimentally using site-directed mutagenesis of CYP3A4 and liquid chromatography-mass spectroscopy analysis. RESULTS Docking with MD-refined CYP3A4 structures incorporating hexa-coordinate heme parameters identifies a unique binding mode involving ARG212 and channel 4, unobserved in the starting PDB ID: 1TQN X-ray structure. The models supporting dehydrogenation are consistent with results from in vitro incubations. GENERAL SIGNIFICANCE Our models indicate that coupled structural contributions of the ingress, egress and solvent channels to the CYP3A4 active site geometries play key roles in the observed 4OHT binding modes. Thus adequate sampling of the conformational space of these drug-metabolizing promiscuous enzymes is important for substrates that may bind in malleable regions of the enzyme active-site.
Collapse
Affiliation(s)
- Kiumars Shahrokh
- Department of Pharmacology and Toxicology, College of Pharmacy, Skaggs Hall 201, University of Utah, Salt Lake City, UT 84112, USA
| | | | | |
Collapse
|
35
|
Liu R, Liu J, Tawa G, Wallqvist A. 2D SMARTCyp Reactivity-Based Site of Metabolism Prediction for Major Drug-Metabolizing Cytochrome P450 Enzymes. J Chem Inf Model 2012; 52:1698-712. [DOI: 10.1021/ci3001524] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ruifeng Liu
- DoD Biotechnology High Performance
Computing Software Applications Institute, Telemedicine and Advanced
Technology Research Center, U.S. Army Medical Research and Material Command, Fort Detrick, Maryland 21702,
United States
| | - Jin Liu
- DoD Biotechnology High Performance
Computing Software Applications Institute, Telemedicine and Advanced
Technology Research Center, U.S. Army Medical Research and Material Command, Fort Detrick, Maryland 21702,
United States
| | - Greg Tawa
- DoD Biotechnology High Performance
Computing Software Applications Institute, Telemedicine and Advanced
Technology Research Center, U.S. Army Medical Research and Material Command, Fort Detrick, Maryland 21702,
United States
| | - Anders Wallqvist
- DoD Biotechnology High Performance
Computing Software Applications Institute, Telemedicine and Advanced
Technology Research Center, U.S. Army Medical Research and Material Command, Fort Detrick, Maryland 21702,
United States
| |
Collapse
|
36
|
Rydberg P, Olsen L. Predicting drug metabolism by cytochrome P450 2C9: comparison with the 2D6 and 3A4 isoforms. ChemMedChem 2012; 7:1202-9. [PMID: 22593031 DOI: 10.1002/cmdc.201200160] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2012] [Revised: 04/23/2012] [Indexed: 11/10/2022]
Abstract
By the use of knowledge gained through modeling of drug metabolism mediated by the cytochrome P450 2D6 and 3A4 isoforms, we constructed a 2D-based model for site-of-metabolism prediction for the cytochrome P450 2C9 isoform. The similarities and differences between the models for the 2C9 and 2D6 isoforms are discussed through structural knowledge from the X-ray crystal structures and trends in experimental data. The final model was validated on an independent test set, resulting in an area under the curve value of 0.92, and a site of metabolism was found among the top two ranked atoms for 77% of the compounds.
Collapse
Affiliation(s)
- Patrik Rydberg
- Department of Drug Design and Pharmacology, University of Copenhagen, Universitetsparken 2, Denmark.
| | | |
Collapse
|
37
|
Dalvie D, Sun H, Xiang C, Hu Q, Jiang Y, Kang P. Effect of Structural Variation on Aldehyde Oxidase-Catalyzed Oxidation of Zoniporide. Drug Metab Dispos 2012; 40:1575-87. [DOI: 10.1124/dmd.112.045823] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
38
|
McPhail B, Tie Y, Hong H, Pearce BA, Schnackenberg LK, Ge W, Fuscoe JC, Tong W, Buzatu DA, Wilkes JG, Fowler BA, Demchuk E, Beger RD. Modeling chemical interaction profiles: I. Spectral data-activity relationship and structure-activity relationship models for inhibitors and non-inhibitors of cytochrome P450 CYP3A4 and CYP2D6 isozymes. Molecules 2012; 17:3383-406. [PMID: 22421792 PMCID: PMC6268752 DOI: 10.3390/molecules17033383] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2012] [Revised: 02/27/2012] [Accepted: 02/28/2012] [Indexed: 02/07/2023] Open
Abstract
An interagency collaboration was established to model chemical interactions that may cause adverse health effects when an exposure to a mixture of chemicals occurs. Many of these chemicals—drugs, pesticides, and environmental pollutant—interact at the level of metabolic biotransformations mediated by cytochrome P450 (CYP) enzymes. In the present work, spectral data-activity relationship (SDAR) and structure-activity relationship (SAR) approaches were used to develop machine-learning classifiers of inhibitors and non-inhibitors of the CYP3A4 and CYP2D6 isozymes. The models were built upon 602 reference pharmaceutical compounds whose interactions have been deduced from clinical data, and 100 additional chemicals that were used to evaluate model performance in an external validation (EV) test. SDAR is an innovative modeling approach that relies on discriminant analysis applied to binned nuclear magnetic resonance (NMR) spectral descriptors. In the present work, both 1D 13C and 1D 15N-NMR spectra were used together in a novel implementation of the SDAR technique. It was found that increasing the binning size of 1D 13C-NMR and 15N-NMR spectra caused an increase in the tenfold cross-validation (CV) performance in terms of both the rate of correct classification and sensitivity. The results of SDAR modeling were verified using SAR. For SAR modeling, a decision forest approach involving from 6 to 17 Mold2 descriptors in a tree was used. Average rates of correct classification of SDAR and SAR models in a hundred CV tests were 60% and 61% for CYP3A4, and 62% and 70% for CYP2D6, respectively. The rates of correct classification of SDAR and SAR models in the EV test were 73% and 86% for CYP3A4, and 76% and 90% for CYP2D6, respectively. Thus, both SDAR and SAR methods demonstrated a comparable performance in modeling a large set of structurally diverse data. Based on unique NMR structural descriptors, the new SDAR modeling method complements the existing SAR techniques, providing an independent estimator that can increase confidence in a structure-activity assessment. When modeling was applied to hazardous environmental chemicals, it was found that up to 20% of them may be substrates and up to 10% of them may be inhibitors of the CYP3A4 and CYP2D6 isoforms. The developed models provide a rare opportunity for the environmental health branch of the public health service to extrapolate to hazardous chemicals directly from human clinical data. Therefore, the pharmacological and environmental health branches are both expected to benefit from these reported models.
Collapse
Affiliation(s)
- Brooks McPhail
- Division of Toxicology and Environmental Medicine, Agency for Toxic Substances and Disease Registry, Atlanta, GA 30333, USA; (B.M.); (Y.T.); (B.A.F.)
| | - Yunfeng Tie
- Division of Toxicology and Environmental Medicine, Agency for Toxic Substances and Disease Registry, Atlanta, GA 30333, USA; (B.M.); (Y.T.); (B.A.F.)
| | - Huixiao Hong
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Bruce A. Pearce
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Laura K. Schnackenberg
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Weigong Ge
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - James C. Fuscoe
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Weida Tong
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Dan A. Buzatu
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Jon G. Wilkes
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| | - Bruce A. Fowler
- Division of Toxicology and Environmental Medicine, Agency for Toxic Substances and Disease Registry, Atlanta, GA 30333, USA; (B.M.); (Y.T.); (B.A.F.)
| | - Eugene Demchuk
- Division of Toxicology and Environmental Medicine, Agency for Toxic Substances and Disease Registry, Atlanta, GA 30333, USA; (B.M.); (Y.T.); (B.A.F.)
- Department of Basic Pharmaceutical Sciences, West Virginia University, Morgantown, WV 26506-9530, USA
- Author to whom correspondence should be addressed; ; Tel.: +1-770-488-3327; Fax: +1-404-248-4142
| | - Richard D. Beger
- Division of Systems Biology, National Center for Toxicological Research, U.S. Food and Drug Administration, Jefferson, AR 72079, USA; (H.H.); (B.A.P.); (L.K.S.); (W.G.); (J.C.F.); (W.T.); (D.A.B.); (J.G.W.); (R.D.B.)
| |
Collapse
|
39
|
Kirchmair J, Williamson MJ, Tyzack JD, Tan L, Bond PJ, Bender A, Glen RC. Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms. J Chem Inf Model 2012; 52:617-48. [PMID: 22339582 PMCID: PMC3317594 DOI: 10.1021/ci200542m] [Citation(s) in RCA: 187] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Metabolism of xenobiotics remains a central challenge
for the discovery
and development of drugs, cosmetics, nutritional supplements, and
agrochemicals. Metabolic transformations are frequently related to
the incidence of toxic effects that may result from the emergence
of reactive species, the systemic accumulation of metabolites, or
by induction of metabolic pathways. Experimental investigation of
the metabolism of small organic molecules is particularly resource
demanding; hence, computational methods are of considerable interest
to complement experimental approaches. This review provides a broad
overview of structure- and ligand-based computational methods for
the prediction of xenobiotic metabolism. Current computational approaches
to address xenobiotic metabolism are discussed from three major perspectives:
(i) prediction of sites of metabolism (SOMs), (ii) elucidation of
potential metabolites and their chemical structures, and (iii) prediction
of direct and indirect effects of xenobiotics on metabolizing enzymes,
where the focus is on the cytochrome P450 (CYP) superfamily of enzymes,
the cardinal xenobiotics metabolizing enzymes. For each of these domains,
a variety of approaches and their applications are systematically
reviewed, including expert systems, data mining approaches, quantitative
structure–activity relationships (QSARs), and machine learning-based
methods, pharmacophore-based algorithms, shape-focused techniques,
molecular interaction fields (MIFs), reactivity-focused techniques,
protein–ligand docking, molecular dynamics (MD) simulations,
and combinations of methods. Predictive metabolism is a developing
area, and there is still enormous potential for improvement. However,
it is clear that the combination of rapidly increasing amounts of
available ligand- and structure-related experimental data (in particular,
quantitative data) with novel and diverse simulation and modeling
approaches is accelerating the development of effective tools for
prediction of in vivo metabolism, which is reflected by the diverse
and comprehensive data sources and methods for metabolism prediction
reviewed here. This review attempts to survey the range and scope
of computational methods applied to metabolism prediction and also
to compare and contrast their applicability and performance.
Collapse
Affiliation(s)
- Johannes Kirchmair
- Unilever Centre for Molecular Science Informatics, Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW, Cambridge, United Kingdom
| | | | | | | | | | | | | |
Collapse
|
40
|
Accessing, using, and creating chemical property databases for computational toxicology modeling. Methods Mol Biol 2012; 929:221-41. [PMID: 23007432 DOI: 10.1007/978-1-62703-050-2_10] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Toxicity data is expensive to generate, is increasingly seen as precompetitive, and is frequently used for the generation of computational models in a discipline known as computational toxicology. Repositories of chemical property data are valuable for supporting computational toxicologists by providing access to data regarding potential toxicity issues with compounds as well as for the purpose of building structure-toxicity relationships and associated prediction models. These relationships use mathematical, statistical, and modeling computational approaches and can be used to understand the mechanisms by which chemicals cause harm and, ultimately, enable prediction of adverse effects of these chemicals to human health and/or the environment. Such approaches are of value as they offer an opportunity to prioritize chemicals for testing. An increasing amount of data used by computational toxicologists is being published into the public domain and, in parallel, there is a greater availability of Open Source software for the generation of computational models. This chapter provides an overview of the types of data and software available and how these may be used to produce predictive toxicology models for the community.
Collapse
|
41
|
Li J, Schneebeli ST, Bylund J, Farid R, Friesner RA. IDSite: An accurate approach to predict P450-mediated drug metabolism. J Chem Theory Comput 2011; 7:3829-3845. [PMID: 22247702 PMCID: PMC3254112 DOI: 10.1021/ct200462q] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate prediction of drug metabolism is crucial for drug design. Since a large majority of drugs metabolism involves P450 enzymes, we herein describe a computational approach, IDSite, to predict P450-mediated drug metabolism. To model induced-fit effects, IDSite samples the conformational space with flexible docking in Glide followed by two refinement stages using the Protein Local Optimization Program (PLOP). Sites of metabolism (SOMs) are predicted according to a physical-based score that evaluates the potential of atoms to react with the catalytic iron center. As a preliminary test, we present in this paper the prediction of hydroxylation and O-dealkylation sites mediated by CYP2D6 using two different models: a physical-based simulation model, and a modification of this model in which a small number of parameters are fit to a training set. Without fitting any parameters to experimental data, the Physical IDSite scoring recovers 83% of the experimental observations for 56 compounds with a very low false positive rate. With only 4 fitted parameters, the Fitted IDSite was trained with the subset of 36 compounds and successfully applied to the other 20 compounds, recovering 94% of the experimental observations with high sensitivity and specificity for both sets.
Collapse
Affiliation(s)
- Jianing Li
- Department of Chemistry, Columbia University, New York, NY
| | | | - Joseph Bylund
- Department of Chemistry, Columbia University, New York, NY
| | - Ramy Farid
- Schrödinger, Inc., 120 W. 45 St., New York, NY
| | | |
Collapse
|
42
|
Moroy G, Martiny VY, Vayer P, Villoutreix BO, Miteva MA. Toward in silico structure-based ADMET prediction in drug discovery. Drug Discov Today 2011; 17:44-55. [PMID: 22056716 DOI: 10.1016/j.drudis.2011.10.023] [Citation(s) in RCA: 167] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2011] [Revised: 10/07/2011] [Accepted: 10/21/2011] [Indexed: 12/12/2022]
Abstract
Quantitative structure-activity relationship (QSAR) methods and related approaches have been used to investigate the molecular features that influence the absorption, distribution, metabolism, excretion and toxicity (ADMET) of drugs. As the three-dimensional structures of several major ADMET proteins become available, structure-based (docking-scoring) computations can be carried out to complement or to go beyond QSAR studies. Applying docking-scoring methods to ADMET proteins is a challenging process because they usually have a large and flexible binding cavity; however, promising results relating to metabolizing enzymes have been reported. After reviewing current trends in the field we applied structure-based methods in the context of receptor flexibility in a case study involving the phase II metabolizing sulfotransferases. Overall, the explored concepts and results suggested that structure-based ADMET profiling will probably join the mainstream during the coming years.
Collapse
Affiliation(s)
- Gautier Moroy
- Inserm UMR-S 973, Molécules Thérapeutiques In Silico, Université Paris Diderot, Sorbonne Paris Cité, 35 Rue Helene Brion, 75013 Paris, France
| | | | | | | | | |
Collapse
|
43
|
Pelkonen O, Turpeinen M, Raunio H. In vivo-in vitro-in silico pharmacokinetic modelling in drug development: current status and future directions. Clin Pharmacokinet 2011; 50:483-91. [PMID: 21740072 DOI: 10.2165/11592400-000000000-00000] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Although clinical drug trials are indispensable in providing an appropriate background for dosage recommendations, they can provide mechanistic pharmacokinetic information only indirectly with the help of certain biomarkers for pathological, physiological and pharmacological determinants. Thus, to provide such mechanistic information of clinical value, various in vitro and in silico tests and approaches are increasingly employed in drug discovery and development. Integration of the results of these primarily preclinical studies has been made possible by various computational models, such as in vitro-in vivo extrapolation of hepatic clearance or physiologically based pharmacokinetic modelling. In this article, the current status of these modelling approaches is surveyed and some examples are given, highlighting advantages and disadvantages in applying them at various phases of drug development. A new paradigm of model-based drug development is briefly described, and the importance of the approach of integrating all of the information coming from different investigations at all levels--be it in vivo, in vitro or in silico--is emphasized.
Collapse
Affiliation(s)
- Olavi Pelkonen
- Department of Pharmacology and Toxicology, Institute of Biomedicine, University of Oulu, Oulu, Finland.
| | | | | |
Collapse
|
44
|
Shahrokh K, Orendt A, Yost GS, Cheatham TE. Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle. J Comput Chem 2011; 33:119-33. [PMID: 21997754 DOI: 10.1002/jcc.21922] [Citation(s) in RCA: 199] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2011] [Revised: 07/28/2011] [Accepted: 07/30/2011] [Indexed: 01/31/2023]
Abstract
Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.
Collapse
Affiliation(s)
- Kiumars Shahrokh
- Department of Pharmacology and Toxicology, University of Utah, Salt Lake City, Utah 84112, USA
| | | | | | | |
Collapse
|
45
|
Mishra NK. Computational modeling of P450s for toxicity prediction. Expert Opin Drug Metab Toxicol 2011; 7:1211-31. [DOI: 10.1517/17425255.2011.611501] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
46
|
T'jollyn H, Boussery K, Mortishire-Smith RJ, Coe K, De Boeck B, Van Bocxlaer JF, Mannens G. Evaluation of Three State-of-the-Art Metabolite Prediction Software Packages (Meteor, MetaSite, and StarDrop) through Independent and Synergistic Use. Drug Metab Dispos 2011; 39:2066-75. [DOI: 10.1124/dmd.111.039982] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
47
|
Danielson ML, Desai PV, Mohutsky MA, Wrighton SA, Lill MA. Potentially increasing the metabolic stability of drug candidates via computational site of metabolism prediction by CYP2C9: The utility of incorporating protein flexibility via an ensemble of structures. Eur J Med Chem 2011; 46:3953-63. [PMID: 21703735 DOI: 10.1016/j.ejmech.2011.05.067] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2011] [Revised: 05/24/2011] [Accepted: 05/26/2011] [Indexed: 10/18/2022]
Abstract
Cytochrome P450 enzymes are responsible for metabolizing many endogenous and xenobiotic molecules encountered by the human body. It has been estimated that 75% of all drugs are metabolized by cytochrome P450 enzymes. Thus, predicting a compound's potential sites of metabolism (SOM) is highly advantageous early in the drug development process. We have combined molecular dynamics, AutoDock Vina docking, the neighboring atom type (NAT) reactivity model, and a solvent-accessible surface-area term to form a reactivity-accessibility model capable of predicting SOM for cytochrome P450 2C9 substrates. To investigate the importance of protein flexibility during the ligand-binding process, the results of SOM prediction using a static protein structure for docking were compared to SOM prediction using multiple protein structures in ensemble docking. The results reported here indicate that ensemble docking increases the number of ligands that can be docked in a bioactive conformation (ensemble: 96%, static: 85%) but only leads to a slight improvement (49% vs. 44%) in predicting an experimentally known SOM in the top-1 position for a ligand library of 75 CYP2C9 substrates. Using ensemble docking, the reactivity-accessibility model accurately predicts SOM in the top-1 ranked position for 49% of the ligand library and considering the top-3 predicted sites increases the prediction success rate to approximately 70% of the ligand library. Further classifying the substrate library according to K(m) values leads to an improvement in SOM prediction for substrates with low K(m) values (57% at top-1). While the current predictive power of the reactivity-accessibility model still leaves significant room for improvement, the results illustrate the usefulness of this method to identify key protein-ligand interactions and guide structural modifications of the ligand to increase its metabolic stability.
Collapse
Affiliation(s)
- Matthew L Danielson
- Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, 575 Stadium Mall Drive, West Lafayette, IN 47907, USA
| | | | | | | | | |
Collapse
|
48
|
Dixit VA, Bharatam PV. Toxic Metabolite Formation from Troglitazone (TGZ): New Insights from a DFT Study. Chem Res Toxicol 2011; 24:1113-22. [DOI: 10.1021/tx200110h] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Vaibhav A. Dixit
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), S. A. S. Nagar, Punjab-160062, India
| | - Prasad V. Bharatam
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), S. A. S. Nagar, Punjab-160062, India
| |
Collapse
|
49
|
Abstract
'It is better to be useful than perfect'. This review attempts to critically cover and assess the currently available approaches and tools to answer the crucial question: Is it possible (and if it is, to what extent is it possible) to predict in vivo metabolites and their abundances on the basis of in vitro and preclinical animal studies? In preclinical drug development, it is possible to produce metabolite patterns from a candidate drug by virtual means (i.e., in silico models), but these are not yet validated. However, they may be useful to cover the potential range of metabolites. In vitro metabolite patterns and apparent relative abundances are produced by various in vitro systems employing tissue preparations (mainly liver) and in most cases using liquid chromatography-mass spectrometry analytical techniques for tentative identification. The pattern of the metabolites produced depends on the enzyme source; the most comprehensive source of drug-metabolizing enzymes is cultured human hepatocytes, followed by liver homogenate fortified with appropriate cofactors. For specific purposes, such as the identification of metabolizing enzyme(s), recombinant enzymes can be used. Metabolite data from animal in vitro and in vivo experiments, despite known species differences, may help pinpoint metabolites that are not apparently produced in in vitro human systems, or suggest alternative experimental approaches. The range of metabolites detected provides clues regarding the enzymes attacking the molecule under study. We also discuss established approaches to identify the major enzymes. The last question, regarding reliability and robustness of metabolite extrapolations from in vitro to in vivo, both qualitatively and quantitatively, cannot be easily answered. There are a number of examples in the literature suggesting that extrapolations are generally useful, but there are only a few systematic and comprehensive studies to validate in vitro-in vivo extrapolations. In conclusion, extrapolation from preclinical metabolite data to the in vivo situation is certainly useful, but it is not known to what extent.
Collapse
|
50
|
Andrade CH, Freitas LMD, Oliveira VD. Twenty-six years of HIV science: an overview of anti-HIV drugs metabolism. BRAZ J PHARM SCI 2011. [DOI: 10.1590/s1984-82502011000200003] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
From the identification of HIV as the agent causing AIDS, to the development of effective antiretroviral drugs, the scientific achievements in HIV research over the past twenty-six years have been formidable. Currently, there are twenty-five anti-HIV compounds which have been formally approved for clinical use in the treatment of AIDS. These compounds fall into six categories: nucleoside reverse transcriptase inhibitors (NRTIs), nucleotide reverse transcriptase inhibitors (NtRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs), cell entry inhibitors or fusion inhibitors (FIs), co-receptor inhibitors (CRIs), and integrase inhibitors (INIs). Metabolism by the host organism is one of the most important determinants of the pharmacokinetic profile of a drug. Formation of active or toxic metabolites will also have an impact on the pharmacological and toxicological outcomes. Therefore, it is widely recognized that metabolism studies of a new chemical entity need to be addressed early in the drug discovery process. This paper describes an overview of the metabolism of currently available anti-HIV drugs.
Collapse
|