1
|
Öeren M, Hunt PA, Wharrick CE, Tabatabaei Ghomi H, Segall MD. Predicting routes of phase I and II metabolism based on quantum mechanics and machine learning. Xenobiotica 2024; 54:379-393. [PMID: 37966132 DOI: 10.1080/00498254.2023.2284251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 11/13/2023] [Indexed: 11/16/2023]
Abstract
Unexpected metabolism could lead to the failure of many late-stage drug candidates or even the withdrawal of approved drugs. Thus, it is critical to predict and study the dominant routes of metabolism in the early stages of research.We describe the development and validation of a 'WhichEnzyme' model that accurately predicts the enzyme families most likely to be responsible for a drug-like molecule's metabolism. Furthermore, we combine this model with our previously published regioselectivity models for Cytochromes P450, Aldehyde Oxidases, Flavin-containing Monooxygenases, UDP-glucuronosyltransferases and Sulfotransferases - the most important Phase I and Phase II drug metabolising enzymes - and a 'WhichP450' model that predicts the Cytochrome P450 isoform(s) responsible for a compound's metabolism.The regioselectivity models are based on a mechanistic understanding of these enzymes' actions and use quantum mechanical simulations with machine learning methods to accurately predict sites of metabolism and the resulting metabolites. We train heuristics based on the outputs of the 'WhichEnzyme', 'WhichP450', and regioselectivity models to determine the most likely routes of metabolism and metabolites to be observed experimentally.Finally, we demonstrate that this combination delivers high sensitivity in identifying experimentally reported metabolites and higher precision than other methods for predicting in vivo metabolite profiles.
Collapse
Affiliation(s)
- Mario Öeren
- Optibrium Limited, Cambridge Innovation Park, Cambridge, UK
| | - Peter A Hunt
- Optibrium Limited, Cambridge Innovation Park, Cambridge, UK
| | | | | | | |
Collapse
|
2
|
Pereira AMG, de Oliveira VM, da Rocha MN, Roberto CHA, Cajazeiras FFM, Guedes JM, Marinho MM, Teixeira AMR, Marinho ES, de Lima-Neto P, Dos Santos HS. Structure and Ligand Based Virtual Screening and MPO Topological Analysis of Triazolo Thiadiazepine-fused Coumarin Derivatives as Anti-Parkinson Drug Candidates. Mol Biotechnol 2024:10.1007/s12033-024-01200-y. [PMID: 38834896 DOI: 10.1007/s12033-024-01200-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 05/10/2024] [Indexed: 06/06/2024]
Abstract
Parkinson's disease (PD) is a debilitating condition that can cause locomotor problems in affected patients, such as tremors and body rigidity. PD therapy often includes the use of monoamine oxidase B (MAOB) inhibitors, particularly phenylhalogen compounds and coumarin-based semi-synthetic compounds. The objective of this study was to analyze the structural, pharmacokinetic, and pharmacodynamic profile of a series of Triazolo Thiadiazepine-fused Coumarin Derivatives (TDCDs) against MAOB, in comparison with the inhibitor safinamide. To achieve this goal, we utilized structure-based virtual screening techniques, including target prediction and absorption, distribution, metabolism, and excretion (ADME) prediction based on multi-parameter optimization (MPO) topological analysis, as well as ligand-based virtual screening techniques, such as docking and molecular dynamics. The findings indicate that the TDCDs exhibit structural similarity to other bioactive compounds containing coumarin and MAOB-binding azoles, which are present in the ChEMBL database. The topological analyses suggest that TDCD3 has the best ADME profile, particularly due to the alignment between low lipophilicity and high polarity. The coumarin and triazole portions make a strong contribution to this profile, resulting in a permeability with Papp estimated at 2.15 × 10-5 cm/s, indicating high cell viability. The substance is predicted to be metabolically stable. It is important to note that this is an objective evaluation based on the available data. Molecular docking simulations showed that the ligand has an affinity energy of - 8.075 kcal/mol with MAOB and interacts with biological substrate residues such as Pro102 and Phe103. The results suggest that the compound has a safe profile in relation to the MAOB model, making it a promising active ingredient for the treatment of PD.
Collapse
Affiliation(s)
- Antônio Mateus Gomes Pereira
- Doctoral Program in Biotechnology, Northeast Biotechnology Network, State University of Ceará, Fortaleza, CE, Brazil
- Center of Molecular Bioprospecting and Applied Experimentation, University Center INTA - UNINTA, Sobral, CE, Brazil
| | | | - Matheus Nunes da Rocha
- Postgraduate Program in Natural Sciences, State University of Ceará, Fortaleza, CE, Brazil
| | | | | | - Jesyka Macêdo Guedes
- Center of Exact Sciences and Technology, State University Vale Do Acaraú, Sobral, CE, Brazil
| | - Márcia Machado Marinho
- Center of Exact Sciences and Technology, State University Vale Do Acaraú, Sobral, CE, Brazil
| | | | - Emmanuel Silva Marinho
- Postgraduate Program in Natural Sciences, State University of Ceará, Fortaleza, CE, Brazil
| | - Pedro de Lima-Neto
- Department of Analytical Chemistry and Phisicochemistry, Federal University of Ceará, Campus Do Pici, Fortaleza, CE, Brazil
| | - Hélcio Silva Dos Santos
- Center of Exact Sciences and Technology, State University Vale Do Acaraú, Sobral, CE, Brazil.
| |
Collapse
|
3
|
Gabor-Worwa E, Kowal-Chwast A, Gaud N, Gogola D, Littlewood P, Smoluch M, Brzózka K, Kus K. Uridine 5'-Diphospho-glucuronosyltransferase 1A3 (UGT1A3) Prediction of Hepatic Clearance of Organic Anion Transporting Polypeptide 1B3 (OATP1B3) Substrate Telmisartan by Glucuronidation Using In Vitro-In Vivo Extrapolation (IVIVE). Eur J Drug Metab Pharmacokinet 2024; 49:393-403. [PMID: 38642299 DOI: 10.1007/s13318-024-00895-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/27/2024] [Indexed: 04/22/2024]
Abstract
BACKGROUND AND OBJECTIVE The prediction of pharmacokinetic parameters for drugs metabolised by cytochrome P450 enzymes has been the subject of active research for many years, while the application of in vitro-in vivo extrapolation (IVIVE) techniques for non-cytochrome P450 enzymes has not been thoroughly evaluated. There is still no established quantitative method for predicting hepatic clearance of drugs metabolised by uridine 5'-diphospho-glucuronosyltransferases (UGTs), not to mention those which undergo hepatic uptake. The objective of the study was to predict the human hepatic clearance for telmisartan based on in vitro metabolic stability and hepatic uptake results. METHODS Telmisartan was examined in liver systems, allowing to estimate intrinsic clearance (CLint, in vitro) based on the substrate disappearance rate with the use of liquid chromatography tandem mass spectrometry (LC-MS/MS) technique. Obtained CLint, in vitro values were corrected for corresponding unbound fractions. Prediction of human hepatic clearance was made from scaled unbound CLint, in vitro data with the use of the well-stirred model, and finally referenced to the literature value of observed clearance in humans, allowing determination of the essential scaling factors. RESULTS The in vitro scaled CLint, in vitro by UGT1A3 was assessed using three systems, human hepatocytes, liver microsomes, and recombinant enzymes. Obtained values were scaled and hepatic metabolism clearance was predicted, resulting in significant clearance underprediction. Utilization of the extended clearance concept (ECC) and hepatic uptake improved prediction of hepatic metabolism clearance. The scaling factors for hepatocytes, assessing the in vitro-in vivo difference, changed from sixfold difference to only twofold difference with the application of the ECC. CONCLUSIONS The study showed that taking into consideration hepatic uptake of a drug allows us to obtain satisfactory scaling factors, hence enabling the prediction of in vivo hepatic glucuronidation from in vitro data.
Collapse
Affiliation(s)
- Ewelina Gabor-Worwa
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland.
- Department of Analytical Chemistry and Biochemistry, Faculty of Materials Science and Ceramics, AGH University of Krakow, Mickiewicza 30 Street, 30-059, Krakow, Poland.
| | - Anna Kowal-Chwast
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
- Department of Analytical Chemistry and Biochemistry, Faculty of Materials Science and Ceramics, AGH University of Krakow, Mickiewicza 30 Street, 30-059, Krakow, Poland
| | - Nilesh Gaud
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
| | - Dawid Gogola
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
| | - Peter Littlewood
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
| | - Marek Smoluch
- Department of Analytical Chemistry and Biochemistry, Faculty of Materials Science and Ceramics, AGH University of Krakow, Mickiewicza 30 Street, 30-059, Krakow, Poland
| | - Krzysztof Brzózka
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
| | - Kamil Kus
- Department of Drug Metabolism and Pharmacokinetics, Ryvu Therapeutics S.A., Sternbacha 2 Street, 30-394, Krakow, Poland
| |
Collapse
|
4
|
Kumar N, He J, Rusling JF. Electrochemical transformations catalyzed by cytochrome P450s and peroxidases. Chem Soc Rev 2023; 52:5135-5171. [PMID: 37458261 DOI: 10.1039/d3cs00461a] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/01/2023]
Abstract
Cytochrome P450s (Cyt P450s) and peroxidases are enzymes featuring iron heme cofactors that have wide applicability as biocatalysts in chemical syntheses. Cyt P450s are a family of monooxygenases that oxidize fatty acids, steroids, and xenobiotics, synthesize hormones, and convert drugs and other chemicals to metabolites. Peroxidases are involved in breaking down hydrogen peroxide and can oxidize organic compounds during this process. Both heme-containing enzymes utilize active FeIVO intermediates to oxidize reactants. By incorporating these enzymes in stable thin films on electrodes, Cyt P450s and peroxidases can accept electrons from an electrode, albeit by different mechanisms, and catalyze organic transformations in a feasible and cost-effective way. This is an advantageous approach, often called bioelectrocatalysis, compared to their biological pathways in solution that require expensive biochemical reductants such as NADPH or additional enzymes to recycle NADPH for Cyt P450s. Bioelectrocatalysis also serves as an ex situ platform to investigate metabolism of drugs and bio-relevant chemicals. In this paper we review biocatalytic electrochemical reactions using Cyt P450s including C-H activation, S-oxidation, epoxidation, N-hydroxylation, and oxidative N-, and O-dealkylation; as well as reactions catalyzed by peroxidases including synthetically important oxidations of organic compounds. Design aspects of these bioelectrocatalytic reactions are presented and discussed, including enzyme film formation on electrodes, temperature, pH, solvents, and activation of the enzymes. Finally, we discuss challenges and future perspective of these two important bioelectrocatalytic systems.
Collapse
Affiliation(s)
- Neeraj Kumar
- Department of Chemistry, University of Connecticut, Storrs, CT 06269-3136, USA.
| | - Jie He
- Department of Chemistry, University of Connecticut, Storrs, CT 06269-3136, USA.
- Institute of Materials Science, University of Connecticut, Storrs, CT 06269-3136, USA
| | - James F Rusling
- Department of Chemistry, University of Connecticut, Storrs, CT 06269-3136, USA.
- Institute of Materials Science, University of Connecticut, Storrs, CT 06269-3136, USA
- Department of Surgery and Neag Cancer Center, Uconn Health, Farmington, CT 06030, USA
- School of Chemistry, National University of Ireland at Galway, Galway, Ireland
| |
Collapse
|
5
|
Öeren M, Kaempf SC, Ponting DJ, Hunt PA, Segall MD. Predicting Regioselectivity of Cytosolic Sulfotransferase Metabolism for Drugs. J Chem Inf Model 2023. [PMID: 37229540 DOI: 10.1021/acs.jcim.3c00275] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Cytosolic sulfotransferases (SULTs) are a family of enzymes responsible for the sulfation of small endogenous and exogenous compounds. SULTs contribute to the conjugation phase of metabolism and share substrates with the uridine 5'-diphospho-glucuronosyltransferase (UGT) family of enzymes. UGTs are considered to be the most important enzymes in the conjugation phase, and SULTs are an auxiliary enzyme system to them. Understanding how the regioselectivity of SULTs differs from that of UGTs is essential from the perspective of developing novel drug candidates. We present a general ligand-based SULT model trained and tested using high-quality experimental regioselectivity data. The current study suggests that, unlike other metabolic enzymes in the modification and conjugation phases, the SULT regioselectivity is not strongly influenced by the activation energy of the rate-limiting step of the catalysis. Instead, the prominent role is played by the substrate binding site of SULT. Thus, the model is trained only on steric and orientation descriptors, which mimic the binding pocket of SULT. The resulting classification model, which predicts whether a site is metabolized, achieved a Cohen's kappa of 0.71.
Collapse
Affiliation(s)
- Mario Öeren
- Cambridge Innovation Park, Optibrium Limited, Denny End Road, Cambridge CB25 9GL, U.K
| | - Sylvia C Kaempf
- Cambridge Innovation Park, Optibrium Limited, Denny End Road, Cambridge CB25 9GL, U.K
- School of Chemistry, North Haugh, University of St Andrews, St Andrews KY16 9ST, U.K
| | - David J Ponting
- Lhasa Limited, Granary Wharf House, 2 Canal Wharf, Leeds LS11 5PS, U.K
| | - Peter A Hunt
- Cambridge Innovation Park, Optibrium Limited, Denny End Road, Cambridge CB25 9GL, U.K
| | - Matthew D Segall
- Cambridge Innovation Park, Optibrium Limited, Denny End Road, Cambridge CB25 9GL, U.K
| |
Collapse
|
6
|
Hirlekar BU, Nuthi A, Singh KD, Murty US, Dixit VA. An overview of compound properties, multiparameter optimization, and computational drug design methods for PARP-1 inhibitor drugs. Eur J Med Chem 2023; 252:115300. [PMID: 36989813 DOI: 10.1016/j.ejmech.2023.115300] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 03/14/2023] [Accepted: 03/17/2023] [Indexed: 03/29/2023]
Abstract
Breast cancer treatment with PARP-1 inhibitors remains challenging due to emerging toxicities, drug resistance, and unaffordable costs of treatment options. How do we invent strategies to design better anti-cancer drugs? A part of the answer is in optimized compound properties, desirability functions, and modern computational drug design methods that drive selectivity and toxicity and have not been reviewed for PARP-1 inhibitors. Nonetheless, comparisons of these compound properties for PARP-1 inhibitors are not available in the literature. In this review, we analyze the physchem, PKPD space to identify inherent desirability functions characteristic of approved drugs that can be valuable for the design of better candidates. Recent literature utilizing ligand, structure-based drug design strategies and matched molecular pair analysis (MMPA) for the discovery of novel PARP-1 inhibitors are also reviewed. Thus, this perspective provides valuable insights into the medchem and multiparameter optimization of PARP-1 inhibitors that might be useful to other medicinal chemists.
Collapse
|
7
|
Dixit VA, Kulkarni A. Applications of Bond Energy‐Based Thermodynamic Analysis to the Feasibility of Unfunctionalized C−C Cross‐Coupling Reactions. ChemistrySelect 2022. [DOI: 10.1002/slct.202203111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Vaibhav A. Dixit
- Department of Medicinal Chemistry National Institute of Pharmaceutical Education and Research Guwahati NIPER Guwahati) Department of Pharmaceuticals Ministry of Chemicals & Fertilizers, Govt. of India, Sila Katamur (Halu-gurisuk) Changsari Kamrup 781101 Guwahati Assam India
| | - Aniket Kulkarni
- Department of Pharmacy Birla Institute of Technology and Sciences Pilani (BITS Pilani) Vidya Vihar Campus, 41 Pilani 333031 Rajasthan India
| |
Collapse
|
8
|
Dixit VA, Murty US, Bajaj P, Blumberger J, de Visser SP. Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3. J Phys Chem B 2022; 126:9737-9747. [PMID: 36384294 DOI: 10.1021/acs.jpcb.2c03967] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; ΔG°), and log ET rates (log kET) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate kET via ET ΔG°. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative ΔG°) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.
Collapse
Affiliation(s)
- Vaibhav A Dixit
- Department of Medicinal Chemistry, Department of Pharmaceuticals, Ministry of Chemicals & Fertilizers,, National Institute of Pharmaceutical Education and Research, Guwahati, (NIPER Guwahati) Govt. of India, Sila Katamur (Halugurisuk), P.O.: Changsari, Dist: Kamrup, 781101Guwahati, Assam, India
| | - Upadhyayula Suryanarayana Murty
- Department of Medicinal Chemistry, Department of Pharmaceuticals, Ministry of Chemicals & Fertilizers,, National Institute of Pharmaceutical Education and Research, Guwahati, (NIPER Guwahati) Govt. of India, Sila Katamur (Halugurisuk), P.O.: Changsari, Dist: Kamrup, 781101Guwahati, Assam, India
| | - Priyanka Bajaj
- National Institute of Pharmaceutical Education and Research, Hyderabad (NIPER Hyderabad), NH-9, Balanagar Main Road, Kukatpally Industrial Estate, Balanagar, Hyderabad500037, Telangana, India
| | - Jochen Blumberger
- Department of Physics and Astronomy, and Thomas Young Centre, University College London, Gower Street, LondonWC1E 6BT, U.K
| | - Sam P de Visser
- Manchester Institute of Biotechnology and Department of Chemical Engineering, The University of Manchester, 131 Princess Street, ManchesterM17DN, U.K
| |
Collapse
|
9
|
Öeren M, Walton PJ, Suri J, Ponting DJ, Hunt PA, Segall MD. Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning. J Med Chem 2022; 65:14066-14081. [PMID: 36239985 DOI: 10.1021/acs.jmedchem.2c01303] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.
Collapse
Affiliation(s)
- Mario Öeren
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge CB25 9GL, U.K
| | - Peter J Walton
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge CB25 9GL, U.K.,School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| | - James Suri
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge CB25 9GL, U.K.,School of Chemistry, University of St Andrews, North Haugh, St Andrews KY16 9ST, U.K
| | - David J Ponting
- Lhasa Limited, Granary Wharf House, 2 Canal Wharf, Leeds LS11 5PS, U.K
| | - Peter A Hunt
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge CB25 9GL, U.K
| | - Matthew D Segall
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge CB25 9GL, U.K
| |
Collapse
|
10
|
Recent advance of fluorescent probes for detection of drug-induced liver injury markers. CHINESE CHEM LETT 2022. [DOI: 10.1016/j.cclet.2021.12.043] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
11
|
Chai L, Zhang H, Song R, Yang H, Yu H, Paneth P, Kepp KP, Akamatsu M, Ji L. Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2021; 55:14037-14050. [PMID: 34663070 DOI: 10.1021/acs.est.1c03036] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcoholβ-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and βC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the βC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcoholβ-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehydeβ-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.
Collapse
Affiliation(s)
- Lihong Chai
- School of Environment Science and Spatial Informatics, China University of Mining and Technology, Daxue Road 1, Xuzhou 221116, China
- College of Environmental and Resource Sciences, Zhejiang University, Yuhangtang Road 866, Hangzhou 310058, China
| | - Huanni Zhang
- College of Environmental and Resource Sciences, Zhejiang University, Yuhangtang Road 866, Hangzhou 310058, China
| | - Runqian Song
- College of Environmental and Resource Sciences, Zhejiang University, Yuhangtang Road 866, Hangzhou 310058, China
| | - Haohan Yang
- School of Environment Science and Spatial Informatics, China University of Mining and Technology, Daxue Road 1, Xuzhou 221116, China
| | - Haiying Yu
- College of Geography and Environmental Sciences, Zhejiang Normal University, Jinhua 321004, China
| | - Piotr Paneth
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| | - Kasper P Kepp
- DTU Chemistry, Technical University of Denmark, Building 206, Kgs. Lyngby DK-2800, Denmark
| | - Miki Akamatsu
- Graduate School of Agriculture, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Li Ji
- School of Environment Science and Spatial Informatics, China University of Mining and Technology, Daxue Road 1, Xuzhou 221116, China
- College of Environmental and Resource Sciences, Zhejiang University, Yuhangtang Road 866, Hangzhou 310058, China
- Graduate School of Agriculture, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
12
|
Roach S, Faponle AS, Satpathy JK, Sastri CV, de Visser SP. Substrate sulfoxidation by a biomimetic cytochrome P450 Compound I mimic: How do porphyrin and phthalocyanine equatorial ligands compare? J CHEM SCI 2021. [DOI: 10.1007/s12039-021-01917-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
13
|
Öeren M, Walton PJ, Hunt PA, Ponting DJ, Segall MD. Predicting reactivity to drug metabolism: beyond P450s-modelling FMOs and UGTs. J Comput Aided Mol Des 2020; 35:541-555. [PMID: 32533369 DOI: 10.1007/s10822-020-00321-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2020] [Accepted: 06/07/2020] [Indexed: 11/28/2022]
Abstract
We present a study based on density functional theory calculations to explore the rate limiting steps of product formation for oxidation by Flavin-containing Monooxygenase (FMO) and glucuronidation by the UDP-glucuronosyltransferase (UGT) family of enzymes. FMOs are responsible for the modification phase of metabolism of a wide diversity of drugs, working in conjunction with Cytochrome P450 (CYP) family of enzymes, and UGTs are the most important class of drug conjugation enzymes. Reactivity calculations are important for prediction of metabolism by CYPs and reactivity alone explains around 70-85% of the experimentally observed sites of metabolism within CYP substrates. In the current work we extend this approach to propose model systems which can be used to calculate the activation energies, i.e. reactivity, for the rate-limiting steps for both FMO oxidation and glucuronidation of potential sites of metabolism. These results are validated by comparison with the experimentally observed reaction rates and sites of metabolism, indicating that the presented models are suitable to provide the basis of a reactivity component within generalizable models to predict either FMO or UGT metabolism.
Collapse
Affiliation(s)
- Mario Öeren
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge, CB25 9PB, UK.
| | - Peter J Walton
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge, CB25 9PB, UK.,School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, UK
| | - Peter A Hunt
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge, CB25 9PB, UK
| | - David J Ponting
- Lhasa Limited, Granary Wharf House, 2 Canal Wharf, Leeds, LS11 5PS, UK
| | - Matthew D Segall
- Optibrium Limited, Cambridge Innovation Park, Denny End Road, Cambridge, CB25 9PB, UK
| |
Collapse
|
14
|
Ji L. Synergy between Experiments and Computations: A Green Channel for Revealing Metabolic Mechanism of Xenobiotics in Chemical Toxicology. Chem Res Toxicol 2020; 33:1539-1550. [DOI: 10.1021/acs.chemrestox.9b00448] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Li Ji
- College of Environmental and Resource Sciences, Zhejiang University, Yuhangtang Road 866, Hangzhou 310058, China
- Graduate School of Agriculture, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
15
|
Ší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
|
16
|
Yang X, Wang Y, Byrne R, Schneider G, Yang S. Concepts of Artificial Intelligence for Computer-Assisted Drug Discovery. Chem Rev 2019; 119:10520-10594. [PMID: 31294972 DOI: 10.1021/acs.chemrev.8b00728] [Citation(s) in RCA: 350] [Impact Index Per Article: 70.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Artificial intelligence (AI), and, in particular, deep learning as a subcategory of AI, provides opportunities for the discovery and development of innovative drugs. Various machine learning approaches have recently (re)emerged, some of which may be considered instances of domain-specific AI which have been successfully employed for drug discovery and design. This review provides a comprehensive portrayal of these machine learning techniques and of their applications in medicinal chemistry. After introducing the basic principles, alongside some application notes, of the various machine learning algorithms, the current state-of-the art of AI-assisted pharmaceutical discovery is discussed, including applications in structure- and ligand-based virtual screening, de novo drug design, physicochemical and pharmacokinetic property prediction, drug repurposing, and related aspects. Finally, several challenges and limitations of the current methods are summarized, with a view to potential future directions for AI-assisted drug discovery and design.
Collapse
Affiliation(s)
- Xin Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital , Sichuan University , Chengdu , Sichuan 610041 , China
| | - Yifei Wang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital , Sichuan University , Chengdu , Sichuan 610041 , China
| | - Ryan Byrne
- ETH Zurich , Department of Chemistry and Applied Biosciences , Vladimir-Prelog-Weg 4 , CH-8093 Zurich , Switzerland
| | - Gisbert Schneider
- ETH Zurich , Department of Chemistry and Applied Biosciences , Vladimir-Prelog-Weg 4 , CH-8093 Zurich , Switzerland
| | - Shengyong Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital , Sichuan University , Chengdu , Sichuan 610041 , China
| |
Collapse
|
17
|
Lin HY, Ho Y, Liu HL. Structure-Based Pharmacophore Modeling to Discover Novel CCR5 Inhibitors for HIV-1/Cancers Therapy. ACTA ACUST UNITED AC 2019. [DOI: 10.4236/jbise.2019.121002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
18
|
Dixit VA. A simple model to solve a complex drug toxicity problem. Toxicol Res (Camb) 2018; 8:157-171. [PMID: 30997019 DOI: 10.1039/c8tx00261d] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2018] [Accepted: 11/28/2018] [Indexed: 02/06/2023] Open
Abstract
Linear drug toxicity models like therapeutic index (TI), physicochemical rules (rule of five, 3/75), ligand efficiency indices (LEI), ideal pharmacokinetic (PK) and pharmacodynamic (PD) profiles are widely used in drug discovery and development. In spite of this, predicting drug toxicity at various stages remains challenging and the overall productivity (<20%) and ultimate benefit to the patients remain low. A simple drug toxicity model, "Drug Toxicity Index" (DTI), is developed here using 711 oral drugs. DTI redefines drug toxicity as scaled biphasic and exponential functions of PD, PK and physicochemical parameters. PD, PK and physicochemical toxicity contributions were estimated from the on and off target IC50, maximum unbound plasma drug concentration (free C max), and log D values, respectively. These contributions are then scaled by molar dose and oral bioavailability and the logarithm of the sum of scaled contributions is DTI. Drugs with DTI above the WHO ATC drug category specific average values consistently have toxic profiles, while drugs with DTI below this average are relatively safe. DTI performs better than standard rules for lead optimization, LEI and exposure based TIs in identifying safe and toxic drugs. DTI classifies 392 drugs reported in the US-FDA's Liver Toxicity Knowledge Base (LTKB) with an AUC for ROC curves of 0.91-0.64 for different WHO ATC categories. DTI has been used to predict network meta-analysis results on relative toxicity within/across eight different therapeutic areas. It is useful in understanding PD, PK and physicochemical toxicity contributions and identifying potentially toxic drugs and the toxicity of recently approved drugs. Decision trees are proposed for applying the DTI concept in preclinical drug discovery and clinical trial settings. DTI can potentially reduce failure in drug discovery and might be useful in therapeutic drug monitoring and in xenobiotic and environmental toxicity studies.
Collapse
Affiliation(s)
- Vaibhav A Dixit
- Department of Pharmacy , Birla Institute of Technology and Sciences Pilani (BITS Pilani) , Vidya Vihar Campus , Street number 41 , Pilani , 333031 , Rajasthan , India . ; ; Tel: +91 1596 255652 ; Tel: +91-7709129400
| |
Collapse
|
19
|
Finkelmann AR, Goldmann D, Schneider G, Göller AH. MetScore: Site of Metabolism Prediction Beyond Cytochrome P450 Enzymes. ChemMedChem 2018; 13:2281-2289. [PMID: 30184341 DOI: 10.1002/cmdc.201800309] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Revised: 08/31/2018] [Indexed: 12/20/2022]
Abstract
The metabolism of xenobiotics by humans and other organisms is a complex process involving numerous enzymes that catalyze phase I (functionalization) and phase II (conjugation) reactions. Herein we introduce MetScore, a machine learning model that can predict both phase I and phase II reaction sites of drugs in a single prediction run. We developed cheminformatics workflows to filter and process reactions to obtain suitable phase I and phase II data sets for model training. Employing a recently developed molecular representation based on quantum chemical partial charges, we constructed random forest machine learning models for phase I and phase II reactions. After combining these models with our previous cytochrome P450 model and calibrating the combination against Bayer in-house data, we obtained the MetScore model that shows good performance, with Matthews correlation coefficients of 0.61 and 0.76 for diverse phase I and phase II reaction types, respectively. We validated its potential applicability to lead optimization campaigns for a new and independent data set compiled from recent publications. The results of this study demonstrate the usefulness of quantum-chemistry-derived molecular representations for reactivity prediction.
Collapse
Affiliation(s)
- Arndt R Finkelmann
- ETH Zurich, Department of Chemistry and Applied Biosciences, Zurich, Switzerland
| | - Daria Goldmann
- KNIME GmbH, Reichenaustrasse 11, 78467, Konstanz, Germany
| | - Gisbert Schneider
- ETH Zurich, Department of Chemistry and Applied Biosciences, Zurich, Switzerland
| | - Andreas H Göller
- Bayer AG, Pharmaceuticals, Research & Development, 42096, Wuppertal, Germany
| |
Collapse
|