1
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
2
|
Akkus E, Tayfuroglu O, Yildiz M, Kocak A. Accurate Binding Free Energy Method from End-State MD Simulations. J Chem Inf Model 2022; 62:4095-4106. [PMID: 35972783 PMCID: PMC9472276 DOI: 10.1021/acs.jcim.2c00601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
Herein, we introduce a new strategy to estimate binding
free energies
using end-state molecular dynamics simulation trajectories. The method
is adopted from linear interaction energy (LIE) and ANI-2x neural
network potentials (machine learning) for the atomic simulation environment
(ASE). It predicts the single-point interaction energies between ligand–protein
and ligand–solvent pairs at the accuracy of the wb97x/6-31G*
level for the conformational space that is sampled by molecular dynamics
(MD) simulations. Our results on 54 protein–ligand complexes
show that the method can be accurate and have a correlation of R = 0.87–0.88 to the experimental binding free energies,
outperforming current end-state methods with reduced computational
cost. The method also allows us to compare BFEs of ligands with different
scaffolds. The code is available free of charge (documentation and
test files) at https://github.com/otayfuroglu/deepQM.
Collapse
Affiliation(s)
- Ebru Akkus
- Department of Bioengineering, Gebze Technical University, 41400 Gebze, Kocaeli, Turkey
| | - Omer Tayfuroglu
- Department of Chemistry, Gebze Technical University, 41400 Gebze, Kocaeli, Turkey
| | - Muslum Yildiz
- Department of Molecular Biology and Genetics, Gebze Technical University, 41400 Gebze, Kocaeli, Turkey
| | - Abdulkadir Kocak
- Department of Chemistry, Gebze Technical University, 41400 Gebze, Kocaeli, Turkey
| |
Collapse
|
3
|
Nunes RS, Vila-Viçosa D, Costa PJ. Halogen Bonding: An Underestimated Player in Membrane–Ligand Interactions. J Am Chem Soc 2021; 143:4253-4267. [DOI: 10.1021/jacs.0c12470] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Rafael Santana Nunes
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
- Centro de Química Estrutural, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
| | - Diogo Vila-Viçosa
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| | - Paulo J. Costa
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| |
Collapse
|
4
|
Gundelach L, Fox T, Tautermann CS, Skylaris CK. Protein–ligand free energies of binding from full-protein DFT calculations: convergence and choice of exchange–correlation functional. Phys Chem Chem Phys 2021; 23:9381-9393. [DOI: 10.1039/d1cp00206f] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum mechanical binding free energies based on thousands of full-protein DFT calculations are tractable, reproducible and converge well.
Collapse
Affiliation(s)
- Lennart Gundelach
- University of Southampton Faculty of Engineering Science and Mathematics, Chemistry
- University Road
- Southampton
- UK
| | - Thomas Fox
- Boehringer Ingelheim Pharma GmbH & Co KG
- Medicinal Chemistry
- 88397 Biberach an der Riss
- Germany
| | | | - Chris-Kriton Skylaris
- University of Southampton Faculty of Engineering Science and Mathematics, Chemistry
- University Road
- Southampton
- UK
| |
Collapse
|
5
|
Costa PJ, Nunes R, Vila-Viçosa D. Halogen bonding in halocarbon-protein complexes and computational tools for rational drug design. Expert Opin Drug Discov 2019; 14:805-820. [PMID: 31131651 DOI: 10.1080/17460441.2019.1619692] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Introduction: Halogens have a prominent role in drug design. Often used as a mean to improve ADME properties, they are also becoming a tool in protein-ligand recognition given their ability to form a non-covalent interaction, termed halogen bond, where halogens act as electrophilic species interacting with electron-rich partners. Rational drug design of halogen-bonding lead molecules requires an accurate description of halocarbon-protein complexes by computational tools though not all methods are able to tackle this non-covalent interaction. Areas covered: The authors present a review of computational methodologies that can be used to properly describe halogen bonds in the context of protein-ligand complexes, providing also insights on how these methods can be used in the context of computer-aided drug design. Expert opinion: Although in the last few years many computational tools, ranging from fast screening methods to the more expensive QM calculations, have been developed to tackle the halogen bonding phenomenon, they are not yet standard in the literature. This will eventually change as official software distributions are including support for halogen bonding in their methods. Tackling desolvation of halogenated species seems to be a good strategy to improve the accuracy of computational methods, that will be more commonly used prior to laboratory work in the future.
Collapse
Affiliation(s)
- Paulo J Costa
- a Centro de Quı́mica e Bioquı́mica, Departamento de Quı́mica e Bioquı́mica , Faculdade de Ciências da Universidade de Lisboa, Campo Grande , Lisboa , Portugal.,b University of Lisboa, Faculty of Sciences , BioISI - Biosystems & Integrative Sciences Institute , Lisboa , Portugal
| | - Rafael Nunes
- a Centro de Quı́mica e Bioquı́mica, Departamento de Quı́mica e Bioquı́mica , Faculdade de Ciências da Universidade de Lisboa, Campo Grande , Lisboa , Portugal.,b University of Lisboa, Faculty of Sciences , BioISI - Biosystems & Integrative Sciences Institute , Lisboa , Portugal
| | - Diogo Vila-Viçosa
- a Centro de Quı́mica e Bioquı́mica, Departamento de Quı́mica e Bioquı́mica , Faculdade de Ciências da Universidade de Lisboa, Campo Grande , Lisboa , Portugal.,b University of Lisboa, Faculty of Sciences , BioISI - Biosystems & Integrative Sciences Institute , Lisboa , Portugal
| |
Collapse
|
6
|
Glycan binding and specificity of viral influenza neuraminidases by classical molecular dynamics and replica exchange molecular dynamics simulations. J Biomol Struct Dyn 2018; 37:3354-3365. [DOI: 10.1080/07391102.2018.1514326] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
7
|
Looking Back, Looking Forward at Halogen Bonding in Drug Discovery. Molecules 2017; 22:molecules22091397. [PMID: 28837116 PMCID: PMC6151711 DOI: 10.3390/molecules22091397] [Citation(s) in RCA: 102] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2017] [Accepted: 08/18/2017] [Indexed: 11/25/2022] Open
Abstract
Halogen bonding has emerged at the forefront of advances in improving ligand: receptor interactions. In particular the newfound ability of this extant non-covalent-bonding phenomena has revolutionized computational approaches to drug discovery while simultaneously reenergizing synthetic approaches to the field. Here we survey, via examples of classical applications involving halogen atoms in pharmaceutical compounds and their biological hosts, the unique advantages that halogen atoms offer as both Lewis acids and Lewis bases.
Collapse
|
8
|
Rowe RK, Ho PS. Relationships between hydrogen bonds and halogen bonds in biological systems. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2017; 73:255-264. [PMID: 28362290 DOI: 10.1107/s2052520617003109] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Accepted: 02/24/2017] [Indexed: 06/07/2023]
Abstract
The recent recognition that halogen bonding (XB) plays important roles in the recognition and assembly of biological molecules has led to new approaches in medicinal chemistry and biomolecular engineering. When designing XBs into strategies for rational drug design or into a biomolecule to affect its structure and function, we must consider the relationship between this interaction and the more ubiquitous hydrogen bond (HB). In this review, we explore these relationships by asking whether and how XBs can replace, compete against or behave independently of HBs in various biological systems. The complex relationships between the two interactions inform us of the challenges we face in fully utilizing XBs to control the affinity and recognition of inhibitors against their therapeutic targets, and to control the structure and function of proteins, nucleic acids and other biomolecular scaffolds.
Collapse
Affiliation(s)
- Rhianon K Rowe
- Department of Biochemistry and Molecular Biology, Colorado State University, Fort Collins, Colorado 80523-1870, USA
| | - P Shing Ho
- Department of Biochemistry and Molecular Biology, Colorado State University, Fort Collins, Colorado 80523-1870, USA
| |
Collapse
|
9
|
González-Vera JA, Medina RA, Martín-Fontecha M, Gonzalez A, de la Fuente T, Vázquez-Villa H, García-Cárceles J, Botta J, McCormick PJ, Benhamú B, Pardo L, López-Rodríguez ML. A new serotonin 5-HT 6 receptor antagonist with procognitive activity - Importance of a halogen bond interaction to stabilize the binding. Sci Rep 2017; 7:41293. [PMID: 28117458 PMCID: PMC5259792 DOI: 10.1038/srep41293] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Accepted: 12/16/2016] [Indexed: 12/13/2022] Open
Abstract
Serotonin 5-HT6 receptor has been proposed as a promising therapeutic target for cognition enhancement though the development of new antagonists is still needed to validate these molecules as a drug class for the treatment of Alzheimer's disease and other pathologies associated with memory deficiency. As part of our efforts to target the 5-HT6 receptor, new benzimidazole-based compounds have been designed and synthesized. Site-directed mutagenesis and homology models show the importance of a halogen bond interaction between a chlorine atom of the new class of 5-HT6 receptor antagonists identified herein and a backbone carbonyl group in transmembrane domain 4. In vitro pharmacological characterization of 5-HT6 receptor antagonist 7 indicates high affinity and selectivity over a panel of receptors including 5-HT2B subtype and hERG channel, which suggests no major cardiac issues. Compound 7 exhibited in vivo procognitive activity (1 mg/kg, ip) in the novel object recognition task as a model of memory deficit.
Collapse
Affiliation(s)
- Juan A. González-Vera
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Rocío A. Medina
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Mar Martín-Fontecha
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Angel Gonzalez
- Laboratori de Medicina Computacional, Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Spain
| | - Tania de la Fuente
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Henar Vázquez-Villa
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Javier García-Cárceles
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Joaquín Botta
- School of Pharmacy, University of East Anglia, NR4 7TJ Norwich, UK
| | | | - Bellinda Benhamú
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| | - Leonardo Pardo
- Laboratori de Medicina Computacional, Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Spain
| | - María L. López-Rodríguez
- Departamento de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
| |
Collapse
|
10
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
11
|
Chen J, Wang J, Zhang Q, Chen K, Zhu W. Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation. Sci Rep 2015; 5:17421. [PMID: 26616018 PMCID: PMC4663504 DOI: 10.1038/srep17421] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2015] [Accepted: 10/29/2015] [Indexed: 12/13/2022] Open
Abstract
Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250014, China
| | - Jinan Wang
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Kaixian Chen
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Weiliang Zhu
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| |
Collapse
|
12
|
Affiliation(s)
- Melissa Coates Ford
- Department of Biochemistry & Molecular Biology, Colorado State University, Fort Collins, Colorado 80523-1870, United States
| | - P. Shing Ho
- Department of Biochemistry & Molecular Biology, Colorado State University, Fort Collins, Colorado 80523-1870, United States
| |
Collapse
|
13
|
Yang Z, Liu Y, Chen Z, Xu Z, Shi J, Chen K, Zhu W. A quantum mechanics-based halogen bonding scoring function for protein-ligand interactions. J Mol Model 2015; 21:138. [DOI: 10.1007/s00894-015-2681-6] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2015] [Accepted: 04/13/2015] [Indexed: 03/08/2023]
|
14
|
Yilmazer ND, Korth M. Enhanced semiempirical QM methods for biomolecular interactions. Comput Struct Biotechnol J 2015; 13:169-75. [PMID: 25848495 PMCID: PMC4372622 DOI: 10.1016/j.csbj.2015.02.004] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2015] [Revised: 02/17/2015] [Accepted: 02/19/2015] [Indexed: 12/21/2022] Open
Abstract
Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.
Collapse
Affiliation(s)
| | - Martin Korth
- Institute of Theoretical Chemistry, Ulm University, D-89069 Ulm, Germany
| |
Collapse
|
15
|
Chen J, Wang J, Zhang Q, Chen K, Zhu W. A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation. J Biomol Struct Dyn 2015; 33:2606-18. [PMID: 25562613 DOI: 10.1080/07391102.2014.1003146] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.
Collapse
Affiliation(s)
- Jianzhong Chen
- a School of Science , Shandong Jiaotong University , Jinan , 250014 , China
| | - Jinan Wang
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| | - Qinggang Zhang
- c College of Physics and Electronics , Shandong Normal University , Jinan , 250014 , China
| | - Kaixian Chen
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| | - Weiliang Zhu
- b Discovery and Design Center , CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences , 555 Zuchongzhi Road, Shanghai , 201203 , China
| |
Collapse
|
16
|
Luo Y, Yin S, Lai W, Wang Y. Effects of global orbital cutoff value and numerical basis set size on accuracies of theoretical atomization energies. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1580-8] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
17
|
Wang L, Gao J, Bi F, Song B, Liu C. Toward the Development of the Potential with Angular Distortion for Halogen Bond: A Comparison of Potential Energy Surfaces between Halogen Bond and Hydrogen Bond. J Phys Chem A 2014; 118:9140-7. [DOI: 10.1021/jp502739c] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
| | | | | | - Bo Song
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Laboratory of Physical Biology, Shanghai, 201800, P. R. China
| | | |
Collapse
|
18
|
Wang YT, Chen YC. Insights from QM/MM Modeling the 3D Structure of the 2009 H1N1 Influenza A Virus Neuraminidase and Its Binding Interactions with Antiviral Drugs. Mol Inform 2014; 33:240-9. [PMID: 27485692 DOI: 10.1002/minf.201300117] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Accepted: 11/15/2014] [Indexed: 11/07/2022]
Abstract
Control of drug release through the inhibition of neuraminidase-1 has been identified as a potential target for the treatment of H1N1 influenza; however, the drug binding mode of neuraminidase is not yet completely understood. In this work, we propose a model for a neuraminidase-1 complex based on four known X-ray structures of drug/neuraminidase-1 complexes. Specifically, H1N1 neuraminidase-1 complexed with 4 drugs (zanamivir, laninamivir, laninamivir octanoate and oseltamivir) was first investigated using a combined quantum mechanical and molecular mechanical (QM/MM) approach. Based on these structures, a model for the H1N1 neuraminidase-1 complex was proposed and simulated using the same computational protocol. Implications to drug/H1N1 neuraminidase-1 binding modes are discussed. From our simulations, the predicted binding free energies of the four drugs are in good agreement with the experimental results, with the correlation coefficient being 0.84.
Collapse
Affiliation(s)
- Yeng-Tseng Wang
- Department of Biochemistry, College of Medicine, Kaohsiung Medical University, No. 100, Shih-Chuan 1st Road, Kaohsiung 80708, Taiwan, R.O.C.
| | - Yu-Ching Chen
- Department of Biomedical Informatics of Asia University, No. 500, Lioufeng Road, Wufeng, Taichung 41354, Taiwan R.O.C
| |
Collapse
|
19
|
Bi F, Gao J, Wang L, Du L, Song B, Liu C. Polarization-enhanced bonding process of halogen bond, a theoretical study on F–H/F–X (X=F,Cl,Br,I) and ammonia. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2013.09.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
20
|
Liu Y, Xu Z, Yang Z, Chen K, Zhu W. A knowledge-based halogen bonding scoring function for predicting protein-ligand interactions. J Mol Model 2013; 19:5015-30. [PMID: 24072554 DOI: 10.1007/s00894-013-2005-7] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2013] [Accepted: 09/08/2013] [Indexed: 11/30/2022]
Abstract
Halogen bonding, a non-covalent interaction between the halogen σ-hole and Lewis bases, could not be properly characterized by majority of current scoring functions. In this study, a knowledge-based halogen bonding scoring function, termed XBPMF, was developed by an iterative method for predicting protein-ligand interactions. Three sets of pairwise potentials were derived from two training sets of protein-ligand complexes from the Protein Data Bank. It was found that two-dimensional pairwise potentials could characterize appropriately the distance and angle profiles of halogen bonding, which is superior to one-dimensional pairwise potentials. With comparison to six widely used scoring functions, XBPMF was evaluated to have moderate power for predicting protein-ligand interactions in terms of "docking power", "ranking power" and "scoring power". Especially, it has a rather satisfactory performance for the systems with typical halogen bonds. To the best of our knowledge, XBPMF is the first halogen bonding scoring function that is not dependent on any dummy atom, and is practical for high-throughput virtual screening. Therefore, this scoring function should be useful for the study and application of halogen bonding interactions like molecular docking and lead optimization.
Collapse
Affiliation(s)
- Yingtao Liu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Structure and Function, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China,
| | | | | | | | | |
Collapse
|
21
|
Kubař T, Bodrog Z, Gaus M, Köhler C, Aradi B, Frauenheim T, Elstner M. Parametrization of the SCC-DFTB Method for Halogens. J Chem Theory Comput 2013; 9:2939-49. [PMID: 26583977 DOI: 10.1021/ct4001922] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Parametrization of the approximative DFT method SCC-DFTB for halogen elements is presented. The new parameter set is intended to describe halogenated organic as well as inorganic molecules, and it is compatible with the established parametrization of SCC-DFTB for carbon, hydrogen, oxygen, and nitrogen. The performance of the parameter set is tested on a representative set of molecules and discussed.
Collapse
Affiliation(s)
- Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Zoltán Bodrog
- Bremen Center for Computational Materials Science, University of Bremen , Am Fallturm 1, 28359 Bremen, Germany
| | - Michael Gaus
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Christof Köhler
- Bremen Center for Computational Materials Science, University of Bremen , Am Fallturm 1, 28359 Bremen, Germany
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen , Am Fallturm 1, 28359 Bremen, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen , Am Fallturm 1, 28359 Bremen, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , Kaiserstr. 12, 76131 Karlsruhe, Germany
| |
Collapse
|
22
|
Vikas, Reenu, Chayawan. Does electron-correlation has any role in the quantitative structure-activity relationships? J Mol Graph Model 2013; 42:7-16. [PMID: 23501159 DOI: 10.1016/j.jmgm.2013.02.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2012] [Revised: 02/05/2013] [Accepted: 02/12/2013] [Indexed: 01/31/2023]
Abstract
For developing quantitative structure-activity relationships (QSARs), quantum-mechanical molecular descriptors based on the state-of-the-art quantum-mechanical methods such as Hartree-Fock (HF) method and density-functional theory (DFT), are now routinely employed. The validity of these quantum-mechanical methods, however, rests on the accurate estimation of electron-correlation energy. This work analyses the role of electron-correlation, using correlation energy as a molecular descriptor, in the QSARs. In particular, QSAR models, for the mutagenic activity of a set of nitrated polycyclic aromatic hydrocarbons (nitro-PAHs), are examined for the role of electron-correlation through state-of-the-art external validation parameters such as concordance correlation coefficient and recently proposed predictive squared correlation coefficients, namely, QF1(2), QF2(2), and QF3(2) etc. The electron-correlation contribution to the highest occupied and lowest unoccupied molecular orbital (HOMO/LUMO) energies is also analyzed. QSAR models based on the semi-empirical quantum-mechanical methods like PM6 and RM1 are also compared. It is found that the models, developed using electron-correlation contribution of the quantum-mechanical descriptors, are not only robust but also relatively more predictive than those developed with the HF and DFT descriptors. The latter are found to be even less reliable than PM6 and RM1 descriptors based models, which show comparable robustness and predictivity with those developed using electron correlation based descriptors. The external predictivity of model based on semi-empirical descriptors can be improved if electron-correlation contribution of the quantum-mechanical descriptors is explicitly included in the model. This work reports the first-ever use of electron-correlation energy and its contribution to the HOMO/LUMO energies as molecular descriptors.
Collapse
Affiliation(s)
- Vikas
- Quantum Chemistry Group, Department of Chemistry & Centre of Advanced Studies in Chemistry, Panjab University, Chandigarh 160 014, India.
| | | | | |
Collapse
|
23
|
Poznański J, Shugar D. Halogen bonding at the ATP binding site of protein kinases: preferred geometry and topology of ligand binding. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:1381-6. [PMID: 23376187 DOI: 10.1016/j.bbapap.2013.01.026] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2012] [Accepted: 01/23/2013] [Indexed: 11/25/2022]
Abstract
Halogenated ligands have been widely developed as potent, and frequently selective, inhibitors of protein kinases (PK). Herein, all structures of protein kinases complexed with a halogenated ligand, identified in the PDB, were analyzed in the context of eventual contribution of halogen bonding to protein-ligand interactions. Global inspection shows that two carbonyl groups of residues located in the hinge region are the most abundant halogen bond acceptors. In contrast to solution data, well-defined water molecules, located at sites conserved across most PK structures, are also involved in halogen bonding. Analysis of cumulative distributions of halogen-acceptor distances shows that structures displaying short contacts involving a halogen atom are overpopulated, contributing together to clearly defined maxima of 2.82, 2.91 and 2.94Å for chlorine, bromine and iodine, respectively. The angular preference of a halogen bond favors ideal topology (180°, 120°) for iodine. For bromine the distribution is much more dispersed, and no such preference was found for chlorine. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).
Collapse
Affiliation(s)
- Jarosław Poznański
- Institute of Biochemistry and Biophysics PAS, Pawińskiego 5a, 02-106 Warszawa, Poland.
| | | |
Collapse
|
24
|
Wilcken R, Zimmermann MO, Lange A, Joerger AC, Boeckler FM. Principles and Applications of Halogen Bonding in Medicinal Chemistry and Chemical Biology. J Med Chem 2013; 56:1363-88. [DOI: 10.1021/jm3012068] [Citation(s) in RCA: 839] [Impact Index Per Article: 76.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Rainer Wilcken
- Laboratory for Molecular Design
and Pharmaceutical Biophysics, Department of Pharmaceutical and Medicinal
Chemistry, Institute of Pharmacy, Eberhard Karls University, Tuebingen, Auf der Morgenstelle 8, 72076 Tuebingen,
Germany
- MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 0QH,
United Kingdom
| | - Markus O. Zimmermann
- Laboratory for Molecular Design
and Pharmaceutical Biophysics, Department of Pharmaceutical and Medicinal
Chemistry, Institute of Pharmacy, Eberhard Karls University, Tuebingen, Auf der Morgenstelle 8, 72076 Tuebingen,
Germany
| | - Andreas Lange
- Laboratory for Molecular Design
and Pharmaceutical Biophysics, Department of Pharmaceutical and Medicinal
Chemistry, Institute of Pharmacy, Eberhard Karls University, Tuebingen, Auf der Morgenstelle 8, 72076 Tuebingen,
Germany
| | - Andreas C. Joerger
- MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 0QH,
United Kingdom
| | - Frank M. Boeckler
- Laboratory for Molecular Design
and Pharmaceutical Biophysics, Department of Pharmaceutical and Medicinal
Chemistry, Institute of Pharmacy, Eberhard Karls University, Tuebingen, Auf der Morgenstelle 8, 72076 Tuebingen,
Germany
| |
Collapse
|
25
|
Carter M, Rappé AK, Ho PS. Scalable Anisotropic Shape and Electrostatic Models for Biological Bromine Halogen Bonds. J Chem Theory Comput 2012; 8:2461-73. [DOI: 10.1021/ct3001969] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Megan Carter
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| | - Anthony K. Rappé
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| | - P. Shing Ho
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| |
Collapse
|
26
|
Molecular mechanical perspective on halogen bonding. J Mol Model 2012; 18:4625-38. [DOI: 10.1007/s00894-012-1454-8] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2012] [Accepted: 04/30/2012] [Indexed: 11/30/2022]
|
27
|
Lu Y, Liu Y, Xu Z, Li H, Liu H, Zhu W. Halogen bonding for rational drug design and new drug discovery. Expert Opin Drug Discov 2012; 7:375-83. [DOI: 10.1517/17460441.2012.678829] [Citation(s) in RCA: 218] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
28
|
Ibrahim MAA. AMBER Empirical Potential Describes the Geometry and Energy of Noncovalent Halogen Interactions Better than Advanced Semiempirical Quantum Mechanical Method PM6-DH2X. J Phys Chem B 2012; 116:3659-69. [DOI: 10.1021/jp3003905] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mahmoud A. A. Ibrahim
- School of
Chemistry, University of Manchester, Oxford
Road, Manchester M139PL, United Kingdom
- Chemistry Department, Faculty
of Science, Minia University, Minia 61519, Egypt
| |
Collapse
|