1
|
Panecka-Hofman J, Poehner I. Structure and dynamics of pteridine reductase 1: the key phenomena relevant to enzyme function and drug design. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2023; 52:521-532. [PMID: 37608196 PMCID: PMC10618315 DOI: 10.1007/s00249-023-01677-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2023] [Revised: 07/08/2023] [Accepted: 07/24/2023] [Indexed: 08/24/2023]
Abstract
Pteridine reductase 1 (PTR1) is a folate and pterin pathway enzyme unique for pathogenic trypanosomatids. As a validated drug target, PTR1 has been the focus of recent research efforts aimed at finding more effective treatments against human parasitic diseases such as leishmaniasis or sleeping sickness. Previous PTR1-centered structural studies highlighted the enzyme characteristics, such as flexible regions around the active site, highly conserved structural waters, and species-specific differences in pocket properties and dynamics, which likely impacts the binding of natural substrates and inhibitors. Furthermore, several aspects of the PTR1 function, such as the substrate inhibition phenomenon and the level of ligand binding cooperativity in the enzyme homotetramer, likely related to the global enzyme dynamics, are poorly known at the molecular level. We postulate that future drug design efforts could greatly benefit from a better understanding of these phenomena through studying both the local and global PTR1 dynamics. This review highlights the key aspects of the PTR1 structure and dynamics relevant to structure-based drug design that could be effectively investigated by modeling approaches. Particular emphasis is given to the perspective of molecular dynamics, what has been accomplished in this area to date, and how modeling could impact the PTR1-targeted drug design in the future.
Collapse
Affiliation(s)
- Joanna Panecka-Hofman
- Division of Biophysics, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02-093, Warsaw, Poland.
| | - Ina Poehner
- School of Pharmacy, University of Eastern Finland, Yliopistonranta 1 C, 70211, Kuopio, Finland
| |
Collapse
|
2
|
Panecka-Hofman J, Poehner I, Wade R. Anti-trypanosomatid structure-based drug design - lessons learned from targeting the folate pathway. Expert Opin Drug Discov 2022; 17:1029-1045. [PMID: 36073204 DOI: 10.1080/17460441.2022.2113776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
INTRODUCTION Trypanosomatidic parasitic infections of humans and animals caused by Trypanosoma brucei, Trypanosoma cruzi, and Leishmania species pose a significant health and economic burden in developing countries. There are few effective and accessible treatments for these diseases, and the existing therapies suffer from problems such as parasite resistance and side effects. Structure-based drug design (SBDD) is one of the strategies that has been applied to discover new compounds targeting trypanosomatid-borne diseases. AREAS COVERED We review the current literature (mostly over the last 5 years, searched in PubMed database on Nov 11th 2021) on the application of structure-based drug design approaches to identify new anti-trypanosomatidic compounds that interfere with a validated target biochemical pathway, the trypanosomatid folate pathway. EXPERT OPINION The application of structure-based drug design approaches to perturb the trypanosomatid folate pathway has successfully provided many new inhibitors with good selectivity profiles, most of which are natural products or their derivatives or have scaffolds of known drugs. However, the inhibitory effect against the target protein(s) often does not translate to anti-parasitic activity. Further progress is hampered by our incomplete understanding of parasite biology and biochemistry, which is necessary to complement SBDD in a multiparameter optimization approach to discovering selective anti-parasitic drugs.
Collapse
Affiliation(s)
- Joanna Panecka-Hofman
- Division of Biophysics, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5a, 02-097 Warsaw, Poland
| | - Ina Poehner
- School of Pharmacy, University of Eastern Finland, Kuopio, Yliopistonranta 1C, PO Box 1627, FI-70211 Kuopio, Finland
| | - Rebecca Wade
- Center for Molecular Biology (ZMBH), Heidelberg University, Im Neuenheimer Feld 282, Heidelberg 69120, Germany.,Heidelberg Institute for Theoretical Studies (HITS), Schloß-Wolfsbrunnenweg 35, Heidelberg 69118, Germany.,DKFZ-ZMBH Alliance and Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| |
Collapse
|
3
|
Jedwabny W, Dyguda-Kazimierowicz E, Pernal K, Szalewicz K, Patkowski K. Extension of an Atom-Atom Dispersion Function to Halogen Bonds and Its Use for Rational Design of Drugs and Biocatalysts. J Phys Chem A 2021; 125:1787-1799. [PMID: 33620223 PMCID: PMC8028329 DOI: 10.1021/acs.jpca.0c11347] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Revised: 02/10/2021] [Indexed: 12/17/2022]
Abstract
A dispersion function Das in the form of a damped atom-atom asymptotic expansion fitted to ab initio dispersion energies from symmetry-adapted perturbation theory was improved and extended to systems containing heavier halogen atoms. To illustrate its performance, the revised Das function was implemented in the multipole first-order electrostatic and second-order dispersion (MED) scoring model. The extension has allowed applications to a much larger set of biocomplexes than it was possible with the original Das. A reasonable correlation between MED and experimentally determined inhibitory activities was achieved in a number of test cases, including structures featuring nonphysically shortened intermonomer distances, which constitute a particular challenge for binding strength predictions. Since the MED model is also computationally efficient, it can be used for reliable and rapid assessment of the ligand affinity or multidimensional scanning of amino acid side-chain conformations in the process of rational design of novel drugs or biocatalysts.
Collapse
Affiliation(s)
- Wiktoria Jedwabny
- Department
of Chemistry, Wrocław University of
Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Edyta Dyguda-Kazimierowicz
- Department
of Chemistry, Wrocław University of
Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Katarzyna Pernal
- Institute
of Physics, Łódź University
of Technology, Wólczańska
219, 90-924 Łódź, Poland
| | - Krzysztof Szalewicz
- Department
of Physics and Astronomy, University of
Delaware, Newark, Delaware 19716, United
States
| | - Konrad Patkowski
- Department
of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
4
|
Jedwabny W, Dyguda-Kazimierowicz E. Revisiting the halogen bonding between phosphodiesterase type 5 and its inhibitors. J Mol Model 2019; 25:29. [PMID: 30613843 PMCID: PMC6321839 DOI: 10.1007/s00894-018-3897-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Accepted: 12/04/2018] [Indexed: 01/12/2023]
Abstract
Halogenated ligands are nowadays commonly designed in order to increase their potency against protein targets. Although novel computational methods of evaluating the affinity of such halogenated inhibitors have emerged, they still lack the sufficient accuracy, which is especially noticeable in the case of empirical scoring functions, being the method of choice in the drug design process. Here, we evaluated a series of halogenated inhibitors of phosphodiesterase type 5 with ab initio methods, revealing the physical nature of ligand binding and determining the components of interaction energy that are essential for proper inhibitor ranking. In particular, a nonempirical scoring model combining long-range contributions to the interaction energy provided a significant correlation with experimental binding potency, outperforming a number of commonly used empirical scoring functions. Considering the low computational cost associated with remarkable predictive abilities of the aforementioned model, it could be used for rapid assessment of the ligand affinity in the process of rational design of novel halogenated compounds.
Collapse
Affiliation(s)
- Wiktoria Jedwabny
- Department of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | | |
Collapse
|
5
|
Haldar S, Comitani F, Saladino G, Woods C, van der Kamp MW, Mulholland AJ, Gervasio FL. A Multiscale Simulation Approach to Modeling Drug-Protein Binding Kinetics. J Chem Theory Comput 2018; 14:6093-6101. [PMID: 30208708 DOI: 10.1021/acs.jctc.8b00687] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Drug-target binding kinetics has recently emerged as a sometimes critical determinant of in vivo efficacy and toxicity. Its rational optimization to improve potency or reduce side effects of drugs is, however, extremely difficult. Molecular simulations can play a crucial role in identifying features and properties of small ligands and their protein targets affecting the binding kinetics, but significant challenges include the long time scales involved in (un)binding events and the limited accuracy of empirical atomistic force fields (lacking, e.g., changes in electronic polarization). In an effort to overcome these hurdles, we propose a method that combines state-of-the-art enhanced sampling simulations and quantum mechanics/molecular mechanics (QM/MM) calculations at the BLYP/VDZ level to compute association free energy profiles and characterize the binding kinetics in terms of structure and dynamics of the transition state ensemble. We test our combined approach on the binding of the anticancer drug Imatinib to Src kinase, a well-characterized target for cancer therapy with a complex binding mechanism involving significant conformational changes. The results indicate significant changes in polarization along the binding pathways, which affect the predicted binding kinetics. This is likely to be of widespread importance in binding of ligands to protein targets.
Collapse
Affiliation(s)
- Susanta Haldar
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , BS8 1TS , United Kingdom
| | | | | | - Christopher Woods
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , BS8 1TS , United Kingdom
| | - Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , BS8 1TS , United Kingdom
- School of Biochemistry , University of Bristol , Bristol , BS8 1TD , United Kingdom
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , BS8 1TS , United Kingdom
| | | |
Collapse
|
6
|
Jedwabny W, Cierpicki T, Grembecka J, Dyguda-Kazimierowicz E. Validation of approximate nonempirical scoring model for menin-mixed lineage leukemia inhibitors. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2350-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
7
|
Jedwabny W, Lodola A, Dyguda-Kazimierowicz E. Theoretical Model of EphA2-Ephrin A1 Inhibition. Molecules 2018; 23:molecules23071688. [PMID: 29997324 PMCID: PMC6099714 DOI: 10.3390/molecules23071688] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Revised: 07/05/2018] [Accepted: 07/06/2018] [Indexed: 02/03/2023] Open
Abstract
This work aims at the theoretical description of EphA2-ephrin A1 inhibition by small molecules. Recently proposed ab initio-based scoring models, comprising long-range components of interaction energy, is tested on lithocholic acid class inhibitors of this protein–protein interaction (PPI) against common empirical descriptors. We show that, although limited to compounds with similar solvation energy, the ab initio model is able to rank the set of selected inhibitors more effectively than empirical scoring functions, aiding the design of novel compounds.
Collapse
Affiliation(s)
- Wiktoria Jedwabny
- Department of Chemistry, Wrocław University of Science and Technology, 50370 Wrocław, Poland.
| | - Alessio Lodola
- Department of Food and Drug, University of Parma, 43100 Parma, Italy.
| | | |
Collapse
|
8
|
The influence of the dispersion corrections on the performance of DFT method in modeling HNgY noble gas molecules and their complexes. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2017.11.040] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|