1
|
Ertl P, Gerebtzoff G, Lewis RA, Muenkler H, Schneider N, Sirockin F, Stiefl N, Tosco P. Chemical reactivity prediction: current methods and different application areas. Mol Inform 2021; 41:e2100277. [PMID: 34964302 DOI: 10.1002/minf.202100277] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 12/28/2021] [Indexed: 11/10/2022]
Abstract
The ability to predict chemical reactivity of a molecule is highly desirable in drug discovery, both ex vivo (synthetic route planning, formulation, stability) and in vivo: metabolic reactions determine pharmacodynamics, pharmacokinetics and potential toxic effects, and early assessment of liabilities is vital to reduce attrition rates in later stages of development. Quantum mechanics offer a precise description of the interactions between electrons and orbitals in the breaking and forming of new bonds. Modern algorithms and faster computers have allowed the study of more complex systems in a punctual and accurate fashion, and answers for chemical questions around stability and reactivity can now be provided. Through machine learning, predictive models can be built out of descriptors derived from quantum mechanics and cheminformatics, even in the absence of experimental data to train on. In this article, current progress on computational reactivity prediction is reviewed: applications to problems in drug design, such as modelling of metabolism and covalent inhibition, are highlighted and unmet challenges are posed.
Collapse
Affiliation(s)
| | | | - Richard A Lewis
- Computer-Aided Drug Design, Eli Lilly and Company Limited, Windlesham, SWITZERLAND
| | - Hagen Muenkler
- Novartis Institutes for BioMedical Research Inc, SWITZERLAND
| | | | | | | | - Paolo Tosco
- Novartis Institutes for BioMedical Research Inc, SWITZERLAND
| |
Collapse
|
2
|
Holmer M, de Bruyn Kops C, Stork C, Kirchmair J. CYPstrate: A Set of Machine Learning Models for the Accurate Classification of Cytochrome P450 Enzyme Substrates and Non-Substrates. Molecules 2021; 26:molecules26154678. [PMID: 34361831 PMCID: PMC8347321 DOI: 10.3390/molecules26154678] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Revised: 07/23/2021] [Accepted: 07/26/2021] [Indexed: 11/16/2022] Open
Abstract
The interaction of small organic molecules such as drugs, agrochemicals, and cosmetics with cytochrome P450 enzymes (CYPs) can lead to substantial changes in the bioavailability of active substances and hence consequences with respect to pharmacological efficacy and toxicity. Therefore, efficient means of predicting the interactions of small organic molecules with CYPs are of high importance to a host of different industries. In this work, we present a new set of machine learning models for the classification of xenobiotics into substrates and non-substrates of nine human CYP isozymes: CYPs 1A2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D6, 2E1, and 3A4. The models are trained on an extended, high-quality collection of known substrates and non-substrates and have been subjected to thorough validation. Our results show that the models yield competitive performance and are favorable for the detection of CYP substrates. In particular, a new consensus model reached high performance, with Matthews correlation coefficients (MCCs) between 0.45 (CYP2C8) and 0.85 (CYP3A4), although at the cost of coverage. The best models presented in this work are accessible free of charge via the "CYPstrate" module of the New E-Resource for Drug Discovery (NERDD).
Collapse
Affiliation(s)
- Malte Holmer
- Center for Bioinformatics (ZBH), Department of Informatics, Universität Hamburg, 20146 Hamburg, Germany; (M.H.); (C.d.B.K.); (C.S.)
| | - Christina de Bruyn Kops
- Center for Bioinformatics (ZBH), Department of Informatics, Universität Hamburg, 20146 Hamburg, Germany; (M.H.); (C.d.B.K.); (C.S.)
| | - Conrad Stork
- Center for Bioinformatics (ZBH), Department of Informatics, Universität Hamburg, 20146 Hamburg, Germany; (M.H.); (C.d.B.K.); (C.S.)
| | - Johannes Kirchmair
- Center for Bioinformatics (ZBH), Department of Informatics, Universität Hamburg, 20146 Hamburg, Germany; (M.H.); (C.d.B.K.); (C.S.)
- Division of Pharmaceutical Chemistry, Department of Pharmaceutical Sciences, University of Vienna, 1090 Vienna, Austria
- Correspondence:
| |
Collapse
|
3
|
Šícho M, Stork C, Mazzolari A, de Bruyn Kops C, Pedretti A, Testa B, Vistoli G, Svozil D, Kirchmair J. FAME 3: Predicting the Sites of Metabolism in Synthetic Compounds and Natural Products for Phase 1 and Phase 2 Metabolic Enzymes. J Chem Inf Model 2019; 59:3400-3412. [PMID: 31361490 DOI: 10.1021/acs.jcim.9b00376] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In this work we present the third generation of FAst MEtabolizer (FAME 3), a collection of extra trees classifiers for the prediction of sites of metabolism (SoMs) in small molecules such as drugs, druglike compounds, natural products, agrochemicals, and cosmetics. FAME 3 was derived from the MetaQSAR database ( Pedretti et al. J. Med. Chem. 2018 , 61 , 1019 ), a recently published data resource on xenobiotic metabolism that contains more than 2100 substrates annotated with more than 6300 experimentally confirmed SoMs related to redox reactions, hydrolysis and other nonredox reactions, and conjugation reactions. In tests with holdout data, FAME 3 models reached competitive performance, with Matthews correlation coefficients (MCCs) ranging from 0.50 for a global model covering phase 1 and phase 2 metabolism, to 0.75 for a focused model for phase 2 metabolism. A model focused on cytochrome P450 metabolism yielded an MCC of 0.57. Results from case studies with several synthetic compounds, natural products, and natural product derivatives demonstrate the agreement between model predictions and literature data even for molecules with structural patterns clearly distinct from those present in the training data. The applicability domains of the individual models were estimated by a new, atom-based distance measure (FAMEscore) that is based on a nearest-neighbor search in the space of atom environments. FAME 3 is available via a public web service at https://nerdd.zbh.uni-hamburg.de/ and as a self-contained Java software package, free for academic and noncommercial research.
Collapse
Affiliation(s)
- Martin Šícho
- Faculty of Mathematics, Informatics and Natural Sciences, Department of Informatics, Center for Bioinformatics , Universität Hamburg , 20146 Hamburg , Germany.,Faculty of Chemical Technology, Department of Informatics and Chemistry, CZ-OPENSCREEN: National Infrastructure for Chemical Biology , University of Chemistry and Technology Prague , 166 28 Prague 6 , Czech Republic
| | - Conrad Stork
- Faculty of Mathematics, Informatics and Natural Sciences, Department of Informatics, Center for Bioinformatics , Universität Hamburg , 20146 Hamburg , Germany
| | - Angelica Mazzolari
- Facoltà di Scienze del Farmaco, Dipartimento di Scienze Farmaceutiche "Pietro Pratesi" , Università degli Studi di Milano , I-20133 Milan , Italy
| | - Christina de Bruyn Kops
- Faculty of Mathematics, Informatics and Natural Sciences, Department of Informatics, Center for Bioinformatics , Universität Hamburg , 20146 Hamburg , Germany
| | - Alessandro Pedretti
- Facoltà di Scienze del Farmaco, Dipartimento di Scienze Farmaceutiche "Pietro Pratesi" , Università degli Studi di Milano , I-20133 Milan , Italy
| | - Bernard Testa
- University of Lausanne , 1015 Lausanne , Switzerland
| | - Giulio Vistoli
- Facoltà di Scienze del Farmaco, Dipartimento di Scienze Farmaceutiche "Pietro Pratesi" , Università degli Studi di Milano , I-20133 Milan , Italy
| | - Daniel Svozil
- Faculty of Chemical Technology, Department of Informatics and Chemistry, CZ-OPENSCREEN: National Infrastructure for Chemical Biology , University of Chemistry and Technology Prague , 166 28 Prague 6 , Czech Republic
| | - Johannes Kirchmair
- Faculty of Mathematics, Informatics and Natural Sciences, Department of Informatics, Center for Bioinformatics , Universität Hamburg , 20146 Hamburg , Germany
| |
Collapse
|
4
|
Kaitoh K, Kotera M, Funatsu K. Novel Electrotopological Atomic Descriptors for the Prediction of Xenobiotic Cytochrome P450 Reactions. Mol Inform 2019; 38:e1900010. [PMID: 31187601 DOI: 10.1002/minf.201900010] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2019] [Accepted: 04/28/2019] [Indexed: 01/06/2023]
Abstract
Cytochrome P450 (CYP) is an enzyme family that plays a crucial role in metabolism, mainly metabolizing xenobiotics to produce non-toxic structures, however, some metabolized products can cause hepatotoxicity. Hence, predicting the structures of CYP products is an important task in designing non-hepatotoxic drugs. Here, we have developed novel atomic descriptors to predict the sites of metabolism (SoM) in CYP substrates. We proposed descriptors that describe topological and electrostatic characteristics of CYP substrates using Gasteiger charge. The proposed descriptors were applied to CYP3A4 data analysis as a case study. As a result of the descriptor selection, we obtained a gradient boosting decision tree-based SoM classification model that used 139 existing descriptors and the proposed 45 descriptors, and the model performed well in terms of the Matthews correlation coefficient. We also developed a structure converter to predict CYP products. This converter correctly generated 51 structural formulas of experimentally observed CYP3A4 products according to a manual evaluation.
Collapse
Affiliation(s)
- Kazuma Kaitoh
- Department of Chemical System Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
| | - Masaaki Kotera
- Department of Chemical System Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
| | - Kimito Funatsu
- Department of Chemical System Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
| |
Collapse
|
5
|
Nair PC, McKinnon RA, Miners JO. Computational Prediction of the Site(s) of Metabolism and Binding Modes of Protein Kinase Inhibitors Metabolized by CYP3A4. Drug Metab Dispos 2019; 47:616-631. [DOI: 10.1124/dmd.118.085167] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Accepted: 03/18/2019] [Indexed: 01/13/2023] Open
|
6
|
Tyzack JD, Kirchmair J. Computational methods and tools to predict cytochrome P450 metabolism for drug discovery. Chem Biol Drug Des 2019; 93:377-386. [PMID: 30471192 PMCID: PMC6590657 DOI: 10.1111/cbdd.13445] [Citation(s) in RCA: 95] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2018] [Revised: 11/05/2018] [Accepted: 11/11/2018] [Indexed: 01/08/2023]
Abstract
In this review, we present important, recent developments in the computational prediction of cytochrome P450 (CYP) metabolism in the context of drug discovery. We discuss in silico models for the various aspects of CYP metabolism prediction, including CYP substrate and inhibitor predictors, site of metabolism predictors (i.e., metabolically labile sites within potential substrates) and metabolite structure predictors. We summarize the different approaches taken by these models, such as rule‐based methods, machine learning, data mining, quantum chemical methods, molecular interaction fields, and docking. We highlight the scope and limitations of each method and discuss future implications for the field of metabolism prediction in drug discovery.
Collapse
Affiliation(s)
| | - Johannes Kirchmair
- Department of Chemistry, University of Bergen, Bergen, Norway.,Computational Biology Unit (CBU), University of Bergen, Bergen, Norway.,Center for Bioinformatics, Universität Hamburg, Hamburg, Germany
| |
Collapse
|
7
|
Cavasin AT, Hillisch A, Uellendahl F, Schneckener S, Göller AH. Reliable and Performant Identification of Low-Energy Conformers in the Gas Phase and Water. J Chem Inf Model 2018; 58:1005-1020. [PMID: 29717870 DOI: 10.1021/acs.jcim.8b00151] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Prediction of compound properties from structure via quantitative structure-activity relationship and machine-learning approaches is an important computational chemistry task in small-molecule drug research. Though many such properties are dependent on three-dimensional structures or even conformer ensembles, the majority of models are based on descriptors derived from two-dimensional structures. Here we present results from a thorough benchmark study of force field, semiempirical, and density functional methods for the calculation of conformer energies in the gas phase and water solvation as a foundation for the correct identification of relevant low-energy conformers. We find that the tight-binding ansatz GFN-xTB shows the lowest error metrics and highest correlation to the benchmark PBE0-D3(BJ)/def2-TZVP in the gas phase for the computationally fast methods and that in solvent OPLS3 becomes comparable in performance. MMFF94, AM1, and DFTB+ perform worse, whereas the performance-optimized but far more expensive functional PBEh-3c yields energies almost perfectly correlated to the benchmark and should be used whenever affordable. On the basis of our findings, we have implemented a reliable and fast protocol for the identification of low-energy conformers of drug-like molecules in water that can be used for the quantification of strain energy and entropy contributions to target binding as well as for the derivation of conformer-ensemble-dependent molecular descriptors.
Collapse
Affiliation(s)
| | - Alexander Hillisch
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| | - Felix Uellendahl
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| | - Sebastian Schneckener
- Bayer AG , Engineering & Technology, Applied Mathematics , 51368 Leverkusen , Germany
| | - Andreas H Göller
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| |
Collapse
|
8
|
Zheng M, Zhao J, Cui C, Fu Z, Li X, Liu X, Ding X, Tan X, Li F, Luo X, Chen K, Jiang H. Computational chemical biology and drug design: Facilitating protein structure, function, and modulation studies. Med Res Rev 2018; 38:914-950. [DOI: 10.1002/med.21483] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2017] [Revised: 12/13/2017] [Accepted: 12/15/2017] [Indexed: 12/12/2022]
Affiliation(s)
- Mingyue Zheng
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Jihui Zhao
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Chen Cui
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Zunyun Fu
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Xutong Li
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Xiaohong Liu
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
- School of Life Science and Technology; ShanghaiTech University; Shanghai China
| | - Xiaoyu Ding
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Xiaoqin Tan
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Fei Li
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
- Department of Chemistry, College of Sciences; Shanghai University; Shanghai China
| | - Xiaomin Luo
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Kaixian Chen
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
- School of Life Science and Technology; ShanghaiTech University; Shanghai China
| | - Hualiang Jiang
- State Key Laboratory of Drug Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| |
Collapse
|
9
|
Šícho M, de Bruyn Kops C, Stork C, Svozil D, Kirchmair J. FAME 2: Simple and Effective Machine Learning Model of Cytochrome P450 Regioselectivity. J Chem Inf Model 2017; 57:1832-1846. [DOI: 10.1021/acs.jcim.7b00250] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Martin Šícho
- Faculty
of Mathematics, Informatics and Natural Sciences, Department of Computer
Science, Center for Bioinformatics, Universität Hamburg, Hamburg, 20146, Germany
- CZ-OPENSCREEN:
National Infrastructure for Chemical Biology, Laboratory of Informatics
and Chemistry, Faculty of Chemical Technology, University of Chemistry and Technology Prague, 166 28 Prague 6, Czech Republic
| | - Christina de Bruyn Kops
- Faculty
of Mathematics, Informatics and Natural Sciences, Department of Computer
Science, Center for Bioinformatics, Universität Hamburg, Hamburg, 20146, Germany
| | - Conrad Stork
- Faculty
of Mathematics, Informatics and Natural Sciences, Department of Computer
Science, Center for Bioinformatics, Universität Hamburg, Hamburg, 20146, Germany
| | - Daniel Svozil
- CZ-OPENSCREEN:
National Infrastructure for Chemical Biology, Laboratory of Informatics
and Chemistry, Faculty of Chemical Technology, University of Chemistry and Technology Prague, 166 28 Prague 6, Czech Republic
| | - Johannes Kirchmair
- Faculty
of Mathematics, Informatics and Natural Sciences, Department of Computer
Science, Center for Bioinformatics, Universität Hamburg, Hamburg, 20146, Germany
| |
Collapse
|