1
|
Szél V, Zsidó BZ, Jeszenői N, Hetényi C. Target-ligand binding affinity from single point enthalpy calculation and elemental composition. Phys Chem Chem Phys 2023; 25:31714-31725. [PMID: 37964670 DOI: 10.1039/d3cp04483a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Reliable target-ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target-ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target-ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target-ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol-1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2-3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.
Collapse
Affiliation(s)
- Viktor Szél
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Balázs Zoltán Zsidó
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Norbert Jeszenői
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Csaba Hetényi
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| |
Collapse
|
2
|
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
|
3
|
Sabe VT, Ntombela T, Jhamba LA, Maguire GEM, Govender T, Naicker T, Kruger HG. Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review. Eur J Med Chem 2021; 224:113705. [PMID: 34303871 DOI: 10.1016/j.ejmech.2021.113705] [Citation(s) in RCA: 188] [Impact Index Per Article: 62.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 07/12/2021] [Accepted: 07/12/2021] [Indexed: 12/30/2022]
Abstract
Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.
Collapse
Affiliation(s)
- Victor T Sabe
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| | - Thandokuhle Ntombela
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| | - Lindiwe A Jhamba
- HIV Pathogenesis Program, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa; School of Chemistry and Physics, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Thavendran Govender
- Faculty of Science and Agriculture, Department of Chemistry, University of Zululand, KwaDlangezwa, 3886, South Africa
| | - Tricia Naicker
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| |
Collapse
|
4
|
Gupta B, Ambekar RS, Tromer RM, Ghosal PS, Sinha R, Majumder A, Kumbhakar P, Ajayan PM, Galvao DS, Gupta AK, Tiwary CS. Development of a schwarzite-based moving bed 3D printed water treatment system for nanoplastic remediation. RSC Adv 2021; 11:19788-19796. [PMID: 35479224 PMCID: PMC9033670 DOI: 10.1039/d1ra03097c] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Accepted: 05/20/2021] [Indexed: 01/22/2023] Open
Abstract
The impact of micro and nanoplastic debris on our aquatic ecosystem is among the most prominent environmental challenges we face today. In addition, nanoplastics create significant concern for environmentalists because of their toxicity and difficulty in separation and removal. Here we report the development of a 3D printed moving bed water filter (M-3DPWF), which can perform as an efficient nanoplastic scavenger. The enhanced separation of the nanoplastics happens due to the creation of a charged filter material that traps the more surface charged nanoparticles selectively. Synthetic contaminated water from polycarbonate waste has been tested with the filter, and enhanced nanoplastic removal has been achieved. The proposed filtration mechanism of surface-charge based water cleaning is further validated using density function theory (semi-empirical) based simulation. The filter has also shown good structural and mechanical stability in both static and dynamic water conditions. The field suitability of the novel treatment system has also been confirmed using water from various sources, such as sea, river, and pond. Our results suggest that the newly developed water filter can be used for the removal of floating nanoparticles in water as a robust advanced treatment system.
Collapse
Affiliation(s)
- Bramha Gupta
- School of Water Resources, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Rushikesh S Ambekar
- Metallurgical and Materials Engineering, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Raphael M Tromer
- Applied Physics Department, State University of Campinas - UNICAMP 13083-859-Campinas SP Brazil
| | - Partha Sarathi Ghosal
- School of Water Resources, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Rupal Sinha
- School of Water Resources, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Abhradeep Majumder
- School of Environmental Science and Engineering, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Partha Kumbhakar
- Metallurgical and Materials Engineering, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - P M Ajayan
- Department of Materials Science and Nanoengineering, Smalley-Curl Institute, Rice University Houston Texas 77005 USA
| | - Douglas S Galvao
- Applied Physics Department, State University of Campinas - UNICAMP 13083-859-Campinas SP Brazil
| | - Ashok Kumar Gupta
- Environmental Engineering Division, Department of Civil Engineering, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| | - Chandra Sekhar Tiwary
- Metallurgical and Materials Engineering, Indian Institute of Technology Kharagpur Kharagpur 721302 India
| |
Collapse
|
5
|
Enantioseparation, in vitro testing, and structural characterization of triple-binding reactivators of organophosphate-inhibited cholinesterases. Biochem J 2020; 477:2771-2790. [PMID: 32639532 DOI: 10.1042/bcj20200192] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Revised: 07/06/2020] [Accepted: 07/08/2020] [Indexed: 11/17/2022]
Abstract
The enantiomers of racemic 2-hydroxyimino-N-(azidophenylpropyl)acetamide-derived triple-binding oxime reactivators were separated, and tested for inhibition and reactivation of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibited with tabun (GA), cyclosarin (GF), sarin (GB), and VX. Both enzymes showed the greatest affinity toward the methylimidazole derivative (III) of 2-hydroxyimino-N-(azidophenylpropyl)acetamide (I). The crystal structure was determined for the complex of oxime III within human BChE, confirming that all three binding groups interacted with active site residues. In the case of BChE inhibited by GF, oximes I (kr = 207 M-1 min-1) and III (kr = 213 M-1 min-1) showed better reactivation efficiency than the reference oxime 2-PAM. Finally, the key mechanistic steps in the reactivation of GF-inhibited BChE with oxime III were modeled using the PM7R6 method, stressing the importance of proton transfer from Nε of His438 to Oγ of Ser203 for achieving successful reactivation.
Collapse
|
6
|
Bagheri S, Behnejad H, Firouzi R, Karimi-Jafari MH. Using the Semiempirical Quantum Mechanics in Improving the Molecular Docking: A Case Study with CDK2. Mol Inform 2020; 39:e2000036. [PMID: 32485047 DOI: 10.1002/minf.202000036] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2020] [Accepted: 05/28/2020] [Indexed: 11/12/2022]
Abstract
In this study, we use some modified semiempirical quantum mechanics (SQM) methods for improving the molecular docking process. To this end, the three popular SQM Hamiltonians, PM6, PM6-D3H4X, and PM7 are employed for geometry optimization of some binding modes of ligands docked into the human cyclin-dependent kinase 2 (CDK2) by two widely used docking tools, AutoDock and AutoDock Vina. The results were analyzed with two different evaluation metrics: the symmetry-corrected heavy-atom RMSD and the fraction of recovered ligand-protein contacts. It is shown that the evaluation of the fraction of recovered contacts is more useful to measure the similarity between two structures when interacting with a protein. It was also found that AutoDock is more successful than AutoDock Vina in producing the correct ligand poses (RMSD≤2.0 Å) and ranking of the poses. It is also demonstrated that the ligand optimization at the SQM level improves the docking results and the SQM structures have a significantly better fit to the observed crystal structures. Finally, the SQM optimizations reduce the number of close contacts in the docking poses and successfully remove most of the clash or bad contacts between ligand and protein.
Collapse
Affiliation(s)
- Saleh Bagheri
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran
| | - Hassan Behnejad
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran
| | - Rohoullah Firouzi
- Department of Physical Chemistry, Chemistry and Chemical Engineering Research Center of Iran, Tehran, Iran
| | | |
Collapse
|
7
|
Rocha SFLS, Sant'Anna CMR. A procedure combining molecular docking and semiempirical method PM7 for identification of selective Shp2 inhibitors. Biopolymers 2019; 110:e23320. [DOI: 10.1002/bip.23320] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Revised: 05/24/2019] [Accepted: 06/17/2019] [Indexed: 12/12/2022]
Affiliation(s)
- Sheisi F. L. S. Rocha
- Programa de Pós‐Graduação em Química, Instituto de QuímicaUniversidade Federal Rural do Rio de Janeiro Seropédica Brazil
| | - Carlos M. R. Sant'Anna
- Programa de Pós‐Graduação em Química, Instituto de QuímicaUniversidade Federal Rural do Rio de Janeiro Seropédica Brazil
- Departamento de Química Fundamental, Instituto de QuímicaUniversidade Federal Rural do Rio de Janeiro Seropédica Brazil
| |
Collapse
|
8
|
Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, Hou T. End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem Rev 2019; 119:9478-9508. [DOI: 10.1021/acs.chemrev.9b00055] [Citation(s) in RCA: 578] [Impact Index Per Article: 115.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ercheng Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Huiyong Sun
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Zhe Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Hui Liu
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200122, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Tingjun Hou
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
9
|
Interface Interactions of the Bowman-Birk Inhibitor BTCI in a Ternary Complex with Trypsin and Chymotrypsin Evaluated by Semiempirical Quantum Mechanical Calculations. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800754] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
10
|
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
|
11
|
Tutone M, Almerico AM. Recent advances on CDK inhibitors: An insight by means of in silico methods. Eur J Med Chem 2017; 142:300-315. [PMID: 28802482 DOI: 10.1016/j.ejmech.2017.07.067] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Revised: 07/19/2017] [Accepted: 07/28/2017] [Indexed: 02/06/2023]
Abstract
The cyclin dependent kinases (CDKs) are a small family of serine/threonine protein kinases that can act as a potential therapeutic target in several proliferative diseases, including cancer. This short review is a survey on the more recent research progresses in the field achieved by using in silico methods. All the "armamentarium" available to the medicinal chemists (docking protocols and molecular dynamics, fragment-based, de novo design, virtual screening, and QSAR) has been employed to the discovery of new, potent, and selective inhibitors of cyclin dependent kinases. The results cited herein can be useful to understand the nature of the inhibitor-target interactions, and furnish an insight on the structural/molecular requirements necessary to achieve the required selectivity against cyclin dependent kinases over other types of kinases.
Collapse
Affiliation(s)
- Marco Tutone
- Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche (STEBICEF), Università di Palermo, Via Archirafi 32, 90123 Palermo, Italy
| | - Anna Maria Almerico
- Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche (STEBICEF), Università di Palermo, Via Archirafi 32, 90123 Palermo, Italy.
| |
Collapse
|
12
|
Stewart JJP. An investigation into the applicability of the semiempirical method PM7 for modeling the catalytic mechanism in the enzyme chymotrypsin. J Mol Model 2017; 23:154. [PMID: 28378242 PMCID: PMC5380709 DOI: 10.1007/s00894-017-3326-8] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2017] [Accepted: 03/17/2017] [Indexed: 02/01/2023]
Abstract
The catalytic cycle for the serine protease α-chymotrypsin was investigated in an attempt to determine the suitability of using the semiempirical method PM7 in the program MOPAC for investigating enzyme-catalyzed reactions. All six classical intermediates were modeled using standard methods, and were characterized as stable minima on the potential energy surface. Using a modified saddle point optimization method, five transition states were located and verified both by vibrational and by intrinsic reaction coordinate analysis. Some individual features, such as the hydrogen bonds in the oxyanion hole, the nature of various electrostatic interactions, and the role of Met192, were examined. This involved designing and running computational experiments to model mutations that would allow features of interest, in particular the energies involved, to be isolated. Three features within the enzyme were examined in detail: the reaction site itself, where covalent bonds were made and broken, the electrostatic effects of the buried aspartate anion, a passive but essential component of the catalytic triad, and the oxyanion hole, where hydrogen bonds help stabilize charged intermediates. With one minor exception, all phenomena investigated agreed with previously-reported descriptions. This result, along with the fact that all the techniques used were relatively straightforward, leads to the recommendation that PM7 and related methods, such as PM6-D3H4, are appropriate for modeling similar enzyme-catalyzed reactions. Graphical abstract Fifth of six transition states, showing water splitting into hydroxyl anion and a proton, to form the second tetrahedral intermediate and histidinium ion. Atoms of the water molecule involved in the hydrolysis are indicated by halos.
Collapse
Affiliation(s)
- James J P Stewart
- Stewart Computational Chemistry, 15210 Paddington Circle, Colorado Springs, CO, 80921, USA.
| |
Collapse
|
13
|
Nekardová M, Vymětalová L, Khirsariya P, Kováčová S, Hylsová M, Jorda R, Kryštof V, Fanfrlík J, Hobza P, Paruch K. Structural Basis of the Interaction of Cyclin-Dependent Kinase 2 with Roscovitine and Its Analogues Having Bioisosteric Central Heterocycles. Chemphyschem 2017; 18:785-795. [DOI: 10.1002/cphc.201601319] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Revised: 01/24/2017] [Indexed: 11/11/2022]
Affiliation(s)
- Michaela Nekardová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences; Flemingovo nám. 2 166 10 Prague 6 Czech Republic
- Faculty of Mathematics and Physics; Charles University in Prague; Ke Karlovu 3 Prague 2 121 16 Czech Republic
| | - Ladislava Vymětalová
- Laboratory of Growth Regulators, Faculty of Science of; Palacky University and Institute of Experimental Botany of the Czech Academy of Sciences; Šlechtitelů 27 78371 Olomouc Czech Republic
| | - Prashant Khirsariya
- Department of Chemistry, CZ Openscreen; Masaryk University; Kamenice 5 62500 Brno Czech Republic
- International Center for Clinical Research; St. Anne's University Hospital Brno; Pekařská 53 656 91 Brno Czech Republic
| | - Silvia Kováčová
- Department of Chemistry, CZ Openscreen; Masaryk University; Kamenice 5 62500 Brno Czech Republic
| | - Michaela Hylsová
- Department of Chemistry, CZ Openscreen; Masaryk University; Kamenice 5 62500 Brno Czech Republic
| | - Radek Jorda
- Laboratory of Growth Regulators, Faculty of Science of; Palacky University and Institute of Experimental Botany of the Czech Academy of Sciences; Šlechtitelů 27 78371 Olomouc Czech Republic
| | - Vladimír Kryštof
- Laboratory of Growth Regulators, Faculty of Science of; Palacky University and Institute of Experimental Botany of the Czech Academy of Sciences; Šlechtitelů 27 78371 Olomouc Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences; Flemingovo nám. 2 166 10 Prague 6 Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences; Flemingovo nám. 2 166 10 Prague 6 Czech Republic
| | - Kamil Paruch
- Department of Chemistry, CZ Openscreen; Masaryk University; Kamenice 5 62500 Brno Czech Republic
- International Center for Clinical Research; St. Anne's University Hospital Brno; Pekařská 53 656 91 Brno Czech Republic
| |
Collapse
|
14
|
Hylsová M, Carbain B, Fanfrlík J, Musilová L, Haldar S, Köprülüoğlu C, Ajani H, Brahmkshatriya PS, Jorda R, Kryštof V, Hobza P, Echalier A, Paruch K, Lepšík M. Explicit treatment of active-site waters enhances quantum mechanical/implicit solvent scoring: Inhibition of CDK2 by new pyrazolo[1,5-a]pyrimidines. Eur J Med Chem 2016; 126:1118-1128. [PMID: 28039837 DOI: 10.1016/j.ejmech.2016.12.023] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2016] [Revised: 12/07/2016] [Accepted: 12/09/2016] [Indexed: 12/17/2022]
Abstract
We present comprehensive testing of solvent representation in quantum mechanics (QM)-based scoring of protein-ligand affinities. To this aim, we prepared 21 new inhibitors of cyclin-dependent kinase 2 (CDK2) with the pyrazolo[1,5-a]pyrimidine core, whose activities spanned three orders of magnitude. The crystal structure of a potent inhibitor bound to the active CDK2/cyclin A complex revealed that the biphenyl substituent at position 5 of the pyrazolo[1,5-a]pyrimidine scaffold was located in a previously unexplored pocket and that six water molecules resided in the active site. Using molecular dynamics, protein-ligand interactions and active-site water H-bond networks as well as thermodynamics were probed. Thereafter, all the inhibitors were scored by the QM approach utilizing the COSMO implicit solvent model. Such a standard treatment failed to produce a correlation with the experiment (R2 = 0.49). However, the addition of the active-site waters resulted in significant improvement (R2 = 0.68). The activities of the compounds could thus be interpreted by taking into account their specific noncovalent interactions with CDK2 and the active-site waters. In summary, using a combination of several experimental and theoretical approaches we demonstrate that the inclusion of explicit solvent effects enhance QM/COSMO scoring to produce a reliable structure-activity relationship with physical insights. More generally, this approach is envisioned to contribute to increased accuracy of the computational design of novel inhibitors.
Collapse
Affiliation(s)
- Michaela Hylsová
- Department of Chemistry, CZ Openscreen, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Benoit Carbain
- Department of Chemistry, CZ Openscreen, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekařská 53, 656 91 Brno, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Lenka Musilová
- Department of Chemistry, CZ Openscreen, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Susanta Haldar
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic; Regional Center of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic; Regional Center of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| | - Haresh Ajani
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic; Regional Center of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| | - Pathik S Brahmkshatriya
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Radek Jorda
- Laboratory of Growth Regulators, Faculty of Science, Palacký University, Institute of Experimental Botany, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Vladimír Kryštof
- Laboratory of Growth Regulators, Faculty of Science, Palacký University, Institute of Experimental Botany, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic; Regional Center of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| | - Aude Echalier
- Centre de Biochimie Structurale, CNRS UMR 5048 - UM - INSERM U 1054, 29 rue de Navacelles, 34090 Montpellier, France
| | - Kamil Paruch
- Department of Chemistry, CZ Openscreen, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekařská 53, 656 91 Brno, Czech Republic.
| | - Martin Lepšík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nam. 2, 166 10 Prague 6, Czech Republic.
| |
Collapse
|
15
|
Stewart JJP. A method for predicting individual residue contributions to enzyme specificity and binding-site energies, and its application to MTH1. J Mol Model 2016; 22:259. [PMID: 27714533 PMCID: PMC5054044 DOI: 10.1007/s00894-016-3119-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2016] [Accepted: 09/08/2016] [Indexed: 12/20/2022]
Abstract
A new method for predicting the energy contributions to substrate binding and to specificity has been developed. Conventional global optimization methods do not permit the subtle effects responsible for these properties to be modeled with sufficient precision to allow confidence to be placed in the results, but by making simple alterations to the model, the precisions of the various energies involved can be improved from about ±2 kcal mol−1 to ±0.1 kcal mol−1. This technique was applied to the oxidized nucleotide pyrophosphohydrolase enzyme MTH1. MTH1 is unusual in that the binding and reaction sites are well separated—an advantage from a computational chemistry perspective, as it allows the energetics involved in docking to be modeled without the need to consider any issues relating to reaction mechanisms. In this study, two types of energy terms were investigated: the noncovalent interactions between the binding site and the substrate, and those responsible for discriminating between the oxidized nucleotide 8-oxo-dGTP and the normal dGTP. Both of these were investigated using the semiempirical method PM7 in the program MOPAC. The contributions of the individual residues to both the binding energy and the specificity of MTH1 were calculated by simulating the effect of mutations. Where comparisons were possible, all calculated results were in agreement with experimental observations. This technique provides fresh insight into the binding mechanism that enzymes use for discriminating between possible substrates.
Collapse
Affiliation(s)
- James J P Stewart
- Stewart Computational Chemistry, 15210 Paddington Circle, Colorado Springs, CO, 80921, USA.
| |
Collapse
|
16
|
Bazgier V, Berka K, Otyepka M, Banáš P. Exponential repulsion improves structural predictability of molecular docking. J Comput Chem 2016; 37:2485-94. [PMID: 27620738 DOI: 10.1002/jcc.24473] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2016] [Revised: 06/21/2016] [Accepted: 07/16/2016] [Indexed: 01/13/2023]
Abstract
Molecular docking is a powerful tool for theoretical prediction of the preferred conformation and orientation of small molecules within protein active sites. The obtained poses can be used for estimation of binding energies, which indicate the inhibition effect of designed inhibitors, and therefore might be used for in silico drug design. However, the evaluation of ligand binding affinity critically depends on successful prediction of the native binding mode. Contemporary docking methods are often based on scoring functions derived from molecular mechanical potentials. In such potentials, nonbonded interactions are typically represented by electrostatic interactions between atom-centered partial charges and standard 6-12 Lennard-Jones potential. Here, we present implementation and testing of a scoring function based on more physically justified exponential repulsion instead of the standard Lennard-Jones potential. We found that this scoring function significantly improved prediction of the native binding modes in proteins bearing narrow active sites such as serine proteases and kinases. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Václav Bazgier
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research, Institute of Experimental Botany ASCR & Palacky University, Šlechtitelů 11, Olomouc, 783 71, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, Olomouc, 77146, Czech Republic
| | - Karel Berka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, Olomouc, 77146, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, Olomouc, 77146, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, Olomouc, 77146, Czech Republic.
| |
Collapse
|
17
|
Yilmazer ND, Korth M. Recent Progress in Treating Protein-Ligand Interactions with Quantum-Mechanical Methods. Int J Mol Sci 2016; 17:ijms17050742. [PMID: 27196893 PMCID: PMC4881564 DOI: 10.3390/ijms17050742] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2016] [Revised: 04/18/2016] [Accepted: 05/03/2016] [Indexed: 11/16/2022] Open
Abstract
We review the first successes and failures of a “new wave” of quantum chemistry-based approaches to the treatment of protein/ligand interactions. These approaches share the use of “enhanced”, dispersion (D), and/or hydrogen-bond (H) corrected density functional theory (DFT) or semi-empirical quantum mechanical (SQM) methods, in combination with ensemble weighting techniques of some form to capture entropic effects. Benchmark and model system calculations in comparison to high-level theoretical as well as experimental references have shown that both DFT-D (dispersion-corrected density functional theory) and SQM-DH (dispersion and hydrogen bond-corrected semi-empirical quantum mechanical) perform much more accurately than older DFT and SQM approaches and also standard docking methods. In addition, DFT-D might soon become and SQM-DH already is fast enough to compute a large number of binding modes of comparably large protein/ligand complexes, thus allowing for a more accurate assessment of entropic effects.
Collapse
Affiliation(s)
- Nusret Duygu Yilmazer
- Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany.
| | - Martin Korth
- Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany.
| |
Collapse
|
18
|
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
|
19
|
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
|
20
|
Spyrakis F, Cavasotto CN. Open challenges in structure-based virtual screening: Receptor modeling, target flexibility consideration and active site water molecules description. Arch Biochem Biophys 2015; 583:105-19. [DOI: 10.1016/j.abb.2015.08.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 08/03/2015] [Accepted: 08/03/2015] [Indexed: 01/05/2023]
|
21
|
Martin BP, Brandon CJ, Stewart JJP, Braun‐Sand SB. Accuracy issues involved in modeling in vivo protein structures using PM7. Proteins 2015; 83:1427-35. [PMID: 25973843 PMCID: PMC4744657 DOI: 10.1002/prot.24826] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 04/20/2015] [Accepted: 05/02/2015] [Indexed: 11/07/2022]
Abstract
Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X-ray structures. Three important contributory factors are the experimental limitations of X-ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop in heat of formation was calculated. Analysis of the changes showed that about 10% of this decrease in heat of formation was caused by faults in PM7, the balance being attributable to the X-ray structure and the difference between the crystal and solution environments. A previously unknown fault in PM7 was revealed during tests to validate the geometries generated using PM7. Clashscores generated by the Molprobity molecular mechanics structure validation program showed that PM7 was predicting unrealistically close contacts between nonbonding atoms in regions where the local geometry is dominated by very weak noncovalent interactions. The origin of this fault was traced to an underestimation of the core-core repulsion between atoms at distances smaller than the equilibrium distance.
Collapse
Affiliation(s)
- Benjamin P. Martin
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
| | | | - James J. P. Stewart
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
- Stewart Computational Chemistry15210 Paddington CircleColorado SpringsCO80921
| | - Sonja B. Braun‐Sand
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
- Stewart Computational Chemistry15210 Paddington CircleColorado SpringsCO80921
| |
Collapse
|
22
|
Nefedova IV, Martynov AG, Averin AA, Kirakosyan GA, Tsivadze AY, Gorbunova YG. New Octopus-like Phthalocyanines as Fullerene Receptors: Synthesis and Photophysical Investigation. Isr J Chem 2015. [DOI: 10.1002/ijch.201500024] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
23
|
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
|
24
|
Ma S, Tan S, Fang D, Zhang R, Zhou S, Wu W, Zheng K. Probing the binding mechanism of novel dual NF-κB/AP-1 inhibitors by 3D-QSAR, docking and molecular dynamics simulations. RSC Adv 2015. [DOI: 10.1039/c5ra10831d] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Potent dual NF-κB/AP-1 inhibitors could effectively treat immunoinflammatory diseases. An integrated computational study was carried out to identify the most favourable binding sites, the structural features and the interaction mechanisms.
Collapse
Affiliation(s)
- Shaojie Ma
- Department of Physical Chemistry
- College of Pharmacy
- Guangdong Pharmaceutical University
- Guangzhou 510006
- PR China
| | - Shepei Tan
- Department of Physical Chemistry
- College of Pharmacy
- Guangdong Pharmaceutical University
- Guangzhou 510006
- PR China
| | - Danqing Fang
- Department of Cardiothoracic Surgery
- Affiliated Second Hospital of Guangzhou Medical University
- Guangzhou 510260
- PR China
| | - Rong Zhang
- Department of Physical Chemistry
- College of Pharmacy
- Guangdong Pharmaceutical University
- Guangzhou 510006
- PR China
| | - Shengfu Zhou
- Department of Physical Chemistry
- College of Pharmacy
- Guangdong Pharmaceutical University
- Guangzhou 510006
- PR China
| | - Wenjuan Wu
- Department of Physical Chemistry
- College of Pharmacy
- Guangdong Pharmaceutical University
- Guangzhou 510006
- PR China
| | - Kangcheng Zheng
- School of Chemistry and Chemical Engineering
- SunYat-Sen University
- Guangzhou 510275
- PR China
| |
Collapse
|
25
|
Ma S, Zeng G, Fang D, Wang J, Wu W, Xie W, Tan S, Zheng K. Studies of N(9)-arenthenyl purines as novel DFG-in and DFG-out dual Src/Abl inhibitors using 3D-QSAR, docking and molecular dynamics simulations. MOLECULAR BIOSYSTEMS 2014; 11:394-406. [PMID: 25406390 DOI: 10.1039/c4mb00350k] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recently, the development of Src/Abl (c-Src/Bcr-Abl tyrosine kinases) dual inhibitors has attracted considerable attention from the research community for treatment of malignancies. In order to explore the different structural features impacting the Src and Abl inhibitory activities of N(9)-arenethenyl purines and to investigate the molecular mechanisms of ligand-receptor interactions, a molecular modeling study combining the three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations was performed. The obtained CoMFA (comparative molecular field analysis) models exhibited satisfactory internal and external predictability. The plots of the CoMFA fields could be used to investigate the structural differences between DFG-in (targeting the active enzyme conformation) and DFG-out (targeting the inactive enzyme conformation) inhibitors. The key amino acid residues were identified by docking studies, and the detailed binding modes of the compounds with different activities were determined by MD simulations. The binding free energies gave a good correlation with the experimental determined activities. In an energetic analysis, the MM-PBSA (molecular mechanics Poisson-Boltzmann surface) energy decomposition revealed that the van der Waals interactions were the major driving force for the binding of the DFG-in and DFG-out compounds to Src and Abl, especially the hydrophobic interactions between ligands and residues Ala403/380, Asp404/381, and Phe405/382 in DFG-out Src and Abl complexes. They also help to stabilize the DFG-out conformations. These results can offer useful references for designing novel potential DFG-in and DFG-out dual Src/Abl inhibitors.
Collapse
Affiliation(s)
- Shaojie Ma
- Department of Physical Chemistry, College of Pharmacy, Guangdong Pharmaceutical University, Guangzhou 510006, China.
| | | | | | | | | | | | | | | |
Collapse
|
26
|
Adasme-Carreño F, Muñoz-Gutierrez C, Caballero J, Alzate-Morales JH. Performance of the MM/GBSA scoring using a binding site hydrogen bond network-based frame selection: the protein kinase case. Phys Chem Chem Phys 2014; 16:14047-58. [DOI: 10.1039/c4cp01378f] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Conformational clustering using hydrogen bond network analysis improved the MM/GBSA scoring for some protein-kinase–ligand systems used as case studies.
Collapse
Affiliation(s)
- Francisco Adasme-Carreño
- Centro de Bioinformática y Simulación Molecular (CBSM)
- Escuela de Ingeniería en Bioinformática
- Facultad de Ingeniería
- Universidad de Talca
- Talca, Chile
| | - Camila Muñoz-Gutierrez
- Centro de Bioinformática y Simulación Molecular (CBSM)
- Escuela de Ingeniería en Bioinformática
- Facultad de Ingeniería
- Universidad de Talca
- Talca, Chile
| | - Julio Caballero
- Centro de Bioinformática y Simulación Molecular (CBSM)
- Escuela de Ingeniería en Bioinformática
- Facultad de Ingeniería
- Universidad de Talca
- Talca, Chile
| | - Jans H. Alzate-Morales
- Centro de Bioinformática y Simulación Molecular (CBSM)
- Escuela de Ingeniería en Bioinformática
- Facultad de Ingeniería
- Universidad de Talca
- Talca, Chile
| |
Collapse
|
27
|
Fanfrlík J, Brahmkshatriya PS, Řezáč J, Jílková A, Horn M, Mareš M, Hobza P, Lepšík M. Quantum mechanics-based scoring rationalizes the irreversible inactivation of parasitic Schistosoma mansoni cysteine peptidase by vinyl sulfone inhibitors. J Phys Chem B 2013; 117:14973-82. [PMID: 24195769 DOI: 10.1021/jp409604n] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The quantum mechanics (QM)-based scoring function that we previously developed for the description of noncovalent binding in protein-ligand complexes has been modified and extended to treat covalent binding of inhibitory ligands. The enhancements are (i) the description of the covalent bond breakage and formation using hybrid QM/semiempirical QM (QM/SQM) restrained optimizations and (ii) the addition of the new ΔG(cov)' term to the noncovalent score, describing the "free" energy difference between the covalent and noncovalent complexes. This enhanced QM-based scoring function is applied to a series of 20 vinyl sulfone-based inhibitory compounds inactivating the cysteine peptidase cathepsin B1 of the Schistosoma mansoni parasite (SmCB1). The available X-ray structure of the SmCB1 in complex with a potent vinyl sulfone inhibitor K11017 is used as a template to build the other covalently bound complexes and to model the derived noncovalent complexes. We present the correlation of the covalent score and its constituents with the experimental binding data. Four outliers are identified. They contain bulky R1' substituents structurally divergent from the template, which might induce larger protein rearrangements than could be accurately modeled. In summary, we propose a new computational approach and an optimal protocol for the rapid evaluation and prospective design of covalent inhibitors with a conserved binding mode.
Collapse
Affiliation(s)
- Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry, v.v.i., and Gilead Sciences and IOCB Research Center, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| | | | | | | | | | | | | | | |
Collapse
|
28
|
Fanfrlík J, Kolář M, Kamlar M, Hurný D, Ruiz FX, Cousido-Siah A, Mitschler A, Řezáč J, Munusamy E, Lepšík M, Matějíček P, Veselý J, Podjarny A, Hobza P. Modulation of aldose reductase inhibition by halogen bond tuning. ACS Chem Biol 2013; 8:2484-92. [PMID: 23988122 DOI: 10.1021/cb400526n] [Citation(s) in RCA: 74] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
In this paper, we studied a designed series of aldose reductase (AR) inhibitors. The series was derived from a known AR binder, which had previously been shown to form a halogen bond between its bromine atom and the oxygen atom of the Thr-113 side chain of AR. In the series, the strength of the halogen bond was modulated by two factors, namely bromine-iodine substitution and the fluorination of the aromatic ring in several positions. The role of the single halogen bond in AR-ligand binding was elucidated by advanced binding free energy calculations involving the semiempirical quantum chemical Hamiltonian. The results were complemented with ultrahigh-resolution X-ray crystallography and IC50 measurements. All of the AR inhibitors studied were shown by X-ray crystallography to bind in an identical manner. Further, it was demonstrated that it was possible to decrease the IC50 value by about 1 order of magnitude by tuning the strength of the halogen bond by a monoatomic substitution. The calculations revealed that the protein-ligand interaction energy increased upon the substitution of iodine for bromine or upon the addition of electron-withdrawing fluorine atoms to the ring. However, the effect on the binding affinity was found to be more complex due to the change of the solvation/desolvation properties within the ligand series. The study shows that it is possible to modulate the strength of a halogen bond in a protein-ligand complex as was designed based on the previous studies of low-molecular-weight complexes.
Collapse
Affiliation(s)
- Jindřich Fanfrlík
- Institute of Organic
Chemistry and Biochemistry and Gilead Science Research Center, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Michal Kolář
- Institute of Organic
Chemistry and Biochemistry and Gilead Science Research Center, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Martin Kamlar
- Department
of Organic Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - David Hurný
- Department
of Organic Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - Francesc X. Ruiz
- Department
of Integrative Biology, IGBMC, CNRS, INSERM, UdS, 1 rue Laurent Fries, 67404 Illkirch CEDEX, France
| | - Alexandra Cousido-Siah
- Department
of Integrative Biology, IGBMC, CNRS, INSERM, UdS, 1 rue Laurent Fries, 67404 Illkirch CEDEX, France
| | - André Mitschler
- Department
of Integrative Biology, IGBMC, CNRS, INSERM, UdS, 1 rue Laurent Fries, 67404 Illkirch CEDEX, France
| | - Jan Řezáč
- Institute of Organic
Chemistry and Biochemistry and Gilead Science Research Center, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Elango Munusamy
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Martin Lepšík
- Institute of Organic
Chemistry and Biochemistry and Gilead Science Research Center, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Pavel Matějíček
- Department
of Physical and Macromolecular Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - Jan Veselý
- Department
of Organic Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - Alberto Podjarny
- Department
of Integrative Biology, IGBMC, CNRS, INSERM, UdS, 1 rue Laurent Fries, 67404 Illkirch CEDEX, France
| | - Pavel Hobza
- Institute of Organic
Chemistry and Biochemistry and Gilead Science Research Center, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Palacký University, Olomouc, 771 46 Olomouc, Czech Republic
| |
Collapse
|
29
|
The Semiempirical Quantum Mechanical Scoring Function for In Silico Drug Design. Chempluschem 2013; 78:921-931. [DOI: 10.1002/cplu.201300199] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2013] [Indexed: 12/19/2022]
|
30
|
Ilatovskiy AV, Abagyan R, Kufareva I. Quantum Mechanics Approaches to Drug Research in the Era of Structural Chemogenomics. INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY 2013; 113:1669-1675. [PMID: 25414519 PMCID: PMC4235788 DOI: 10.1002/qua.24400] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
The rapid growth of the available crystallographic information about proteins and binding pockets creates remarkable opportunities for enriching the drug research pipelines with computational prediction of novel protein-ligand interactions. While ab initio quantum mechanical approaches are known to provide unprecedented accuracy in structure-based binding energy calculations, they are limited to only small systems of dozens of atoms. In the structural chemogenomics era, it is critical that new approaches are developed that enable application of QM methodologies to non-covalent interactions in systems as large as protein-ligand complexes and conformational ensembles. This perspective highlights recent advances towards bridging the gap between high accuracy and high volume computations in drug research.
Collapse
Affiliation(s)
- Andrey V. Ilatovskiy
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, USA, 92093
- Division of Molecular and Radiation Biophysics, Konstantinov Petersburg Nuclear Physics Institute, NRC Kurchatov Institute, Gatchina, Russia, 188300
| | - Ruben Abagyan
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, USA, 92093
| | - Irina Kufareva
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, USA, 92093
| |
Collapse
|
31
|
Kolář M, Fanfrlík J, Lepšík M, Forti F, Luque FJ, Hobza P. Assessing the accuracy and performance of implicit solvent models for drug molecules: conformational ensemble approaches. J Phys Chem B 2013; 117:5950-62. [PMID: 23600402 DOI: 10.1021/jp402117c] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accuracy and performance of implicit solvent methods for solvation free energy calculations were assessed on a set of 20 neutral drug molecules. Molecular dynamics (MD) provided ensembles of conformations in water and water-saturated octanol. The solvation free energies were calculated by popular implicit solvent models based on quantum mechanical (QM) electronic densities (COSMO-RS, MST, SMD) as well as on molecular mechanical (MM) point-charge models (GB, PB). The performance of the implicit models was tested by a comparison with experimental water-octanol transfer free energies (ΔG(ow)) by using single- and multiconformation approaches. MD simulations revealed difficulties in a priori estimation of the flexibility features of the solutes from simple structural descriptors, such as the number of rotatable bonds. An increasing accuracy of the calculated ΔG(ow) was observed in the following order: GB1 ~ PB < GB7 ≪ MST < SMD ~ COSMO-RS with a clear distinction identified between MM- and QM-based models, although for the set excluding three largest molecules, the differences among COSMO-RS, MST, and SMD were negligible. It was shown that the single-conformation approach applied to crystal geometries provides a rather accurate estimate of ΔG(ow) for rigid molecules yet fails completely for the flexible ones. The multiconformation approaches improved the performance, but only when the deformation contribution was ignored. It was revealed that for large-scale calculations on small molecules a recent GB model, GB7, provided a reasonable accuracy/speed ratio. In conclusion, the study contributes to the understanding of solvation free energy calculations for physical and medicinal chemistry applications.
Collapse
Affiliation(s)
- Michal Kolář
- Institute of Organic Chemistry and Biochemistry and Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, v.v.i., 166 10 Prague 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
32
|
Hostaš J, Řezáč J, Hobza P. On the performance of the semiempirical quantum mechanical PM6 and PM7 methods for noncovalent interactions. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.02.069] [Citation(s) in RCA: 94] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
33
|
Kiewisch K, Jacob CR, Visscher L. Quantum-Chemical Electron Densities of Proteins and of Selected Protein Sites from Subsystem Density Functional Theory. J Chem Theory Comput 2013; 9:2425-40. [DOI: 10.1021/ct3008759] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Karin Kiewisch
- Amsterdam Center for Multiscale
Modeling, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| | - Christoph R. Jacob
- Center for Functional Nanostructures
and Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe,
Germany
| | - Lucas Visscher
- Amsterdam Center for Multiscale
Modeling, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
34
|
Mucs D, Bryce RA. The application of quantum mechanics in structure-based drug design. Expert Opin Drug Discov 2013; 8:263-76. [DOI: 10.1517/17460441.2013.752812] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
35
|
Muddana HS, Gilson MK. Calculation of Host-Guest Binding Affinities Using a Quantum-Mechanical Energy Model. J Chem Theory Comput 2012; 8:2023-2033. [PMID: 22737045 PMCID: PMC3378313 DOI: 10.1021/ct3002738] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The prediction of protein-ligand binding affinities is of central interest in computer-aided drug discovery, but it is still difficult to achieve a high degree of accuracy. Recent studies suggesting that available force fields may be a key source of error motivate the present study, which reports the first mining minima (M2) binding affinity calculations based on a quantum mechanical energy model, rather than an empirical force field. We apply a semi-empirical quantum-mechanical energy function, PM6-DH+, coupled with the COSMO solvation model, to 29 host-guest systems with a wide range of measured binding affinities. After correction for a systematic error, which appears to derive from the treatment of polar solvation, the computed absolute binding affinities agree well with experimental measurements, with a mean error 1.6 kcal/mol and a correlation coefficient of 0.91. These calculations also delineate the contributions of various energy components, including solute energy, configurational entropy, and solvation free energy, to the binding free energies of these host-guest complexes. Comparison with our previous calculations, which used empirical force fields, point to significant differences in both the energetic and entropic components of the binding free energy. The present study demonstrates successful combination of a quantum mechanical Hamiltonian with the M2 affinity method.
Collapse
Affiliation(s)
- Hari S Muddana
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA 92093-0736
| | | |
Collapse
|
36
|
Antony J, Grimme S. Fully ab initio protein-ligand interaction energies with dispersion corrected density functional theory. J Comput Chem 2012; 33:1730-9. [PMID: 22570225 DOI: 10.1002/jcc.23004] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2012] [Revised: 04/10/2012] [Accepted: 04/11/2012] [Indexed: 01/31/2023]
Abstract
Dispersion corrected density functional theory (DFT-D3) is used for fully ab initio protein-ligand (PL) interaction energy calculation via molecular fractionation with conjugated caps (MFCC) and applied to PL complexes from the PDB comprising 3680, 1798, and 1060 atoms. Molecular fragments with n amino acids instead of one in the original MFCC approach are considered, thereby allowing for estimating the three-body and higher many-body terms. n > 1 is recommended both in terms of accuracy and efficiency of MFCC. For neutral protein side-chains, the computed PL interaction energy is visibly independent of the fragment length n. The MFCC fractionation error is determined by comparison to a full-system calculation for the 1060 atoms containing PL complex. For charged amino acid side-chains, the variation of the MFCC result with n is increased. For these systems, using a continuum solvation model with a dielectricity constant typical for protein environments (ϵ = 4) reduces both the variation with n and improves the stability of the DFT calculations considerably. The PL interaction energies for two typical complexes obtained ab initio for the first time are found to be rather large (-30 and -54 kcal/mol).
Collapse
Affiliation(s)
- Jens Antony
- Universität Münster, Organisch-Chemisches Institut, Germany
| | | |
Collapse
|
37
|
Hobza P. Calculations on noncovalent interactions and databases of benchmark interaction energies. Acc Chem Res 2012; 45:663-72. [PMID: 22225511 DOI: 10.1021/ar200255p] [Citation(s) in RCA: 144] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Although covalent interactions determine the primary structure of a molecule, the noncovalent interactions are responsible for the tertiary and quaternary structure of a molecule and create the fascinating world of the 3D architectures of biomacromolecules. For example, the double helical structure of DNA is of fundamental importance for the function of DNA: it allows it to store and transfer genetic information. To fulfill this role, the structure is rigid to maintain the double helix with a proper positioning of the complementary base, and floppy to allow for its opening. Very strong covalent interactions cannot fulfill both of these criteria, but noncovalent interactions, which are about 2 orders of magnitude weaker, can. This Account highlights the recent advances in the field of the design of novel wave function theory (WFT) methods applicable to noncovalent complexes ranging in size from less than 100 atoms, for which highly accurate ab initio methods are available, up to extended ones (several thousands atoms), which are the domain of semiempirical QM (SQM) methods. Accurate interaction energies for noncovalent complexes are generated by the coupled-cluster technique, taking single- and double-electron excitations iteratively and triple-electron excitation perturbatively with a complete basis set description (CCSD(T)/CBS). The procedure provides interaction energies with high accuracy (error less than 1 kcal/mol). Because the method is computationally demanding, its application is limited to complexes smaller than 30 atoms. But researchers would also like to use computational methods to determine these interaction energies accurately for larger biological and nanoscale structures. Standard QM methods such as MP2, MP3, CCSD, or DFT fail to describe various types of noncovalent systems (H-bonded, stacked, dispersion-controlled, etc.) with comparable accuracy. Therefore, novel methods are needed that have been parametrized toward noncovalent interactions, and existing benchmark data sets represent an important tool for the development of new methods providing reliable characteristics of noncovalent clusters. Our laboratory developed the first suitable data set of CCSD(T)/CBS interaction energies and geometries of various noncovalent complexes, called S22. Since its publication in 2006, it has frequently been applied in parametrization and/or verification of various wave function and density functional techniques. During the intense use of this data set, several inconsistencies emerged, such as the insufficient accuracy of the CCSD(T) correction term or its unbalanced character, which has triggered the introduction of a new, broader, and more accurate data set called the S66 data set. It contains not only 66 CCSD(T)/CBS interaction energies determined in the equilibrium geometries but also 1056 interaction energies calculated at the same level for nonequilibrium geometries. The S22 and S66 data sets have been used for the verification of various WFT methods, and the lowest RMSE (S66, in kcal/mol) was found for the recently introduced SCS-MI-CCSD/CBS (0.08), MP2.5/CBS (0.16), MP2.X/6-31G* (0.27), and SCS-MI-MP2/CBS (0.38) methods. Because of their computational economy, the MP2.5 and MP2.X/6-31G* methods can be recommended for highly accurate calculations of large complexes with up to 100 atoms. The evaluation of SQM methods was based only on the S22 data set, and because some of these methods have been parametrized toward the same data set, the respective results should be taken with caution. For really extended complexes such as protein-ligand systems, only the SMQ methods are applicable. After adding the corrections to the dispersion energy and H-bonding, several methods exhibit surprisingly low RMSE (even below 0.5 kcal/mol). Among the various SMQ methods, the PM6-DH2 can be recommended because of its computational efficiency and it can be used for optimization (which is not the case for other SQM methods). The PM6-DH2 is the base of our novel scoring function used in in silico drug design.
Collapse
Affiliation(s)
- Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 771 46 Olomouc, Czech Republic
- Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| |
Collapse
|
38
|
Mikulskis P, Genheden S, Wichmann K, Ryde U. A semiempirical approach to ligand-binding affinities: Dependence on the Hamiltonian and corrections. J Comput Chem 2012; 33:1179-89. [DOI: 10.1002/jcc.22949] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2011] [Revised: 01/10/2012] [Accepted: 01/16/2012] [Indexed: 12/30/2022]
|
39
|
Seddon G, Lounnas V, McGuire R, van den Bergh T, Bywater RP, Oliveira L, Vriend G. Drug design for ever, from hype to hope. J Comput Aided Mol Des 2012; 26:137-50. [PMID: 22252446 PMCID: PMC3268973 DOI: 10.1007/s10822-011-9519-9] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2011] [Accepted: 12/05/2011] [Indexed: 01/28/2023]
Abstract
In its first 25 years JCAMD has been disseminating a large number of techniques aimed at finding better medicines faster. These include genetic algorithms, COMFA, QSAR, structure based techniques, homology modelling, high throughput screening, combichem, and dozens more that were a hype in their time and that now are just a useful addition to the drug-designers toolbox. Despite massive efforts throughout academic and industrial drug design research departments, the number of FDA-approved new molecular entities per year stagnates, and the pharmaceutical industry is reorganising accordingly. The recent spate of industrial consolidations and the concomitant move towards outsourcing of research activities requires better integration of all activities along the chain from bench to bedside. The next 25 years will undoubtedly show a series of translational science activities that are aimed at a better communication between all parties involved, from quantum chemistry to bedside and from academia to industry. This will above all include understanding the underlying biological problem and optimal use of all available data.
Collapse
Affiliation(s)
| | - V. Lounnas
- CMBI, Radboud University Nijmegen Medical Centre, Geert Grooteplein 26–28, 6525 GA Nijmegen, The Netherlands
| | - R. McGuire
- BioAxis Research, Bergse Heihoek 56, Berghem, 5351 SL The Netherlands
| | - T. van den Bergh
- Bio-Prodict, Dreijenplein 10, 6703 HB Wageningen, The Netherlands
| | | | - L. Oliveira
- Sao Paulo Federal University (UNIFESP), Sao Paulo, Brazil
| | - G. Vriend
- CMBI, Radboud University Nijmegen Medical Centre, Geert Grooteplein 26–28, 6525 GA Nijmegen, The Netherlands
| |
Collapse
|
40
|
Molecular simulation methods in drug discovery: a prospective outlook. J Comput Aided Mol Des 2011; 26:81-6. [DOI: 10.1007/s10822-011-9506-1] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2011] [Accepted: 11/30/2011] [Indexed: 11/27/2022]
|
41
|
Caballero J, Alzate-Morales JH, Vergara-Jaque A. Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds. J Chem Inf Model 2011; 51:2920-31. [PMID: 22011048 DOI: 10.1021/ci200306w] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.
Collapse
Affiliation(s)
- Julio Caballero
- Centro de Bioinformática y Simulación Molecular, Universidad de Talca, 2 Norte 685, Casilla 721, Talca, Chile.
| | | | | |
Collapse
|
42
|
Dobeš P, Řezáč J, Fanfrlík J, Otyepka M, Hobza P. Semiempirical Quantum Mechanical Method PM6-DH2X Describes the Geometry and Energetics of CK2-Inhibitor Complexes Involving Halogen Bonds Well, While the Empirical Potential Fails. J Phys Chem B 2011; 115:8581-9. [DOI: 10.1021/jp202149z] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Petr Dobeš
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
- Center of Molecular Biology and Gene Therapy, Department of Internal Medicine − Hematooncology, University Hospital Brno, 625 00 Brno, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 771 46 Olomouc, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
43
|
Pecina A, Přenosil O, Fanfrlík J, Řezáč J, Granatier J, Hobza P, Lepšík M. On the reliability of the corrected semiempirical quantum chemical method (PM6-DH2) for assigning the protonation states in HIV-1 protease/inhibitor complexes. ACTA ACUST UNITED AC 2011. [DOI: 10.1135/cccc2011035] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
A novel computational protocol for determining the most probable protonation states in protein/ligand complexes is presented. The method consists in treating large parts of the enzyme using the corrected semiempirical quantum chemical (QM) method – PM6-D2 for optimization and PM6-DH2 for single-point energies – while the rest is calculated using molecular mechanics (MM) within a hybrid QM/MM fashion. The surrounding solvent is approximated by an implicit model. This approach is applied to two model systems, two different carboxylate pairs in one general and one unique HIV-1 protease/inhibitor complex. The effect of the size of the movable QM part is investigated in a series of several sizes, 3-, 6-, 8- and 10-Å regions surrounding the inhibitor. For the smallest region (< 450 atoms) the computationally more costly DFT QM/MM optimizations are performed as a check of the correctness. Proton transfer (PT) phenomena occur at both the PM6-D2 and DFT levels, which underlines the requirement for a QM approach. The barriers of PT are checked in model carboxylic acid pairs using the highly accurate MP2 and CCSD(T) values. An important result of this study is the fine-tuning of the protocol which can be used in further applications; its limitations are also shown, pointing to future developments. The calculations reveal which protonation variants of the active site are the most stable. In conclusion, the presented protocol can also be utilized for defining probable isomers in biomolecular systems. It can also serve as a preparatory step for further interaction-energy and binding-score calculations.
Collapse
|
44
|
Kolár M, Fanfrlík J, Hobza P. Ligand conformational and solvation/desolvation free energy in protein-ligand complex formation. J Phys Chem B 2011; 115:4718-24. [PMID: 21466174 DOI: 10.1021/jp2010265] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this study, an extensive sampling of the conformational space of nine HIV-1 protease inhibitors was performed to estimate the uncertainty with which a single-conformation scoring scheme approximates the ligand-protein binding free energy. The SMD implicit solvation/desolvation energy and gas-phase PM6-DH2 energy were calculated for a set of 1600 conformations of each ligand. The probability density functions of the energies were compared with the values obtained from the single-conformation approach and from a short ab initio molecular dynamics simulation. The relative uncertainty in the score within the set of nine inhibitors was calculated to be 3.5 kcal·mol(-1) and 2.7 kcal·mol(-1) for the single-conformation and short dynamics, respectively. These results, though limited to the consideration of flexible ligands, provide a valuable insight into the precision of rigid models in the current computer-aided drug design.
Collapse
Affiliation(s)
- Michal Kolár
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
| | | | | |
Collapse
|