1
|
Shukla A, Biswal AS, Chowdhury A, Halder R, Chatterjee S. Aggregation-Induced Modulation of Ground and Excited State Photophysics of 5-( tert-Butyl)-2-Hydroxy-1,3-Isophthalaldehyde (5- tBHI). J Phys Chem B 2024; 128:5437-5453. [PMID: 38662934 DOI: 10.1021/acs.jpcb.4c00258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2024]
Abstract
5-(tert-Butyl)-2-hydroxy-1,3-isophthalaldehyde (5-tBHI) is a photochromic material susceptible to either excited state proton transfer or excited state intramolecular proton transfer, depending upon the solvent. However, it has also been found to aggregate in the presence of sodium dodecyl sulfate. In this current study, based on the steady-state and time-resolved spectroscopy, supported by crystallography, quantum chemical density functional theory calculation, and molecular dynamics (MD) simulation, we report on the aggregation of this potential single benzene-based emitter (SBBE) in neat solvents as well as solid phase to modulate its photophysics. It has been found that 5-tBHI forms mixed aggregates of different orders, owing to the presence of both enolic and tautomeric forms, to yield tunable emission, although the emission intensity is quenched. These findings suggest that the intramolecular hydrogen bonding of 5-tBHI not only limits intermolecular interactions but also promotes nonradiative deactivation pathways. Hence, designing and structural engineering, with a focus to suppressing intramolecular hydrogen bonding as well as increasing through space conjugation by replacing the aldehydic moieties with bulky aliphatic or aromatic ketonic groups, can be a plausible approach to yielding improved probes with tunable emission and higher fluorescence quantum yields.
Collapse
Affiliation(s)
- Aparna Shukla
- Department of Chemistry and Chemical Biology, Indian Institute of Technology (Indian School of Mines) Dhanbad, Dhanbad 826004 Jharkhand, India
| | - Abhipsa Sekhar Biswal
- Department of Chemistry and Chemical Biology, Indian Institute of Technology (Indian School of Mines) Dhanbad, Dhanbad 826004 Jharkhand, India
| | - Arkaprava Chowdhury
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra 400076, India
| | - Ritaban Halder
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, United States
| | - Soumit Chatterjee
- Department of Chemistry and Chemical Biology, Indian Institute of Technology (Indian School of Mines) Dhanbad, Dhanbad 826004 Jharkhand, India
| |
Collapse
|
2
|
Meelua W, Wanjai T, Thinkumrob N, Oláh J, Cairns JRK, Hannongbua S, Ryde U, Jitonnom J. A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase. Phys Chem Chem Phys 2023; 25:8767-8778. [PMID: 36912034 DOI: 10.1039/d2cp05262h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
Abstract
Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of β-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-β-amino acids. This enzyme has the potential to be used as a tool in the production of β-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.
Collapse
Affiliation(s)
- Wijitra Meelua
- Demonstration School, University of Phayao, Phayao 56000, Thailand
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Tanchanok Wanjai
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Natechanok Thinkumrob
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Julianna Oláh
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Műegyetem rakpart 3, Budapest H-1111, Hungary
| | - James R Ketudat Cairns
- Center for Biomolecular Structure, Function and Application and School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
| | - Supa Hannongbua
- Department of Chemistry, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, Lund SE-221 00, Sweden
| | - Jitrayut Jitonnom
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| |
Collapse
|
3
|
Nakamura S, Akaki T, Nishiwaki K, Nakatani M, Kawase Y, Takahashi Y, Nakanishi I. System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study. J Comput Chem 2023; 44:824-831. [PMID: 36444861 DOI: 10.1002/jcc.27044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2022] [Revised: 11/01/2022] [Accepted: 11/03/2022] [Indexed: 11/30/2022]
Abstract
The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.
Collapse
Affiliation(s)
- Shinya Nakamura
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Tatsuo Akaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan.,Chemical Research Laboratories, Central Pharmaceutical Research Institute, Japan Tobacco Inc., Osaka, Japan
| | - Keiji Nishiwaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Midori Nakatani
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuji Kawase
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuki Takahashi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Isao Nakanishi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| |
Collapse
|
4
|
Mahdizadeh SJ, Pålsson E, Carlesso A, Chevet E, Eriksson LA. QM/MM Well-Tempered Metadynamics Study of the Mechanism of XBP1 mRNA Cleavage by Inositol Requiring Enzyme 1α RNase. J Chem Inf Model 2022; 62:4247-4260. [PMID: 35960929 PMCID: PMC9472280 DOI: 10.1021/acs.jcim.2c00735] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A range of in silico methodologies were herein employed to study the unconventional XBP1 mRNA cleavage mechanism performed by the unfolded protein response (UPR) mediator Inositol Requiring Enzyme 1α (IRE1). Using Protein-RNA molecular docking along with a series of extensive restrained/unrestrained atomistic molecular dynamics (MD) simulations, the dynamical behavior of the system was evaluated and a reliable model of the IRE1/XBP1 mRNA complex was constructed. From a series of well-converged quantum mechanics molecular mechanics well-tempered metadynamics (QM/MM WT-MetaD) simulations using the Grimme dispersion interaction corrected semiempirical parametrization method 6 level of theory (PM6-D3) and the AMBER14SB-OL3 force field, the free energy profile of the cleavage mechanism was determined, along with intermediates and transition state structures. The results show two distinct reaction paths based on general acid-general base type mechanisms, with different activation energies that perfectly match observations from experimental mutagenesis data. The study brings unique atomistic insights into the cleavage mechanism of XBP1 mRNA by IRE1 and clarifies the roles of the catalytic residues H910 and Y892. Increased understanding of the details in UPR signaling can assist in the development of new therapeutic agents for its modulation.
Collapse
Affiliation(s)
- Sayyed Jalil Mahdizadeh
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30 Göteborg, Sweden
| | - Emil Pålsson
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30 Göteborg, Sweden
| | - Antonio Carlesso
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30 Göteborg, Sweden.,Faculty of Biomedical Sciences, Euler Institute, Università della Svizzera Italiana (USI),, Lugano 6904, Switzerland.,Department of Pharmacology, Sahlgrenska Academy, University of Gothenburg, 405 30 Göteborg, Sweden
| | - Eric Chevet
- INSERM U1242, University of Rennes 1, 35000 Rennes, France.,Centre Eugène Marquis, 35000 Rennes, France
| | - Leif A Eriksson
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30 Göteborg, Sweden
| |
Collapse
|
5
|
Prasad S, Tantillo DJ. Roads Not Taken: Mechanism and Origins of Regio- and Chemoselectivity of Directed Co III-Catalyzed Alkenylation of N-Pyridyl 2-Pyridone. Organometallics 2022. [DOI: 10.1021/acs.organomet.1c00662] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Supreeth Prasad
- Department of Chemistry, University of California─Davis, Davis, California 95616, United States
| | - Dean J. Tantillo
- Department of Chemistry, University of California─Davis, Davis, California 95616, United States
| |
Collapse
|
6
|
QM and QM/MM study on inhibition mechanism of polyphenolic compounds as non-classical inhibitors of α-human carbonic anhydrase (II). Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02839-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
7
|
Smola M, Gutten O, Dejmek M, Kožíšek M, Evangelidis T, Tehrani ZA, Novotná B, Nencka R, Birkuš G, Rulíšek L, Boura E. Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein. Angew Chem Int Ed Engl 2021; 60:10172-10178. [PMID: 33616279 PMCID: PMC8251555 DOI: 10.1002/anie.202016805] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Indexed: 12/19/2022]
Abstract
STING (stimulator of interferon genes) is a key regulator of innate immunity that has recently been recognized as a promising drug target. STING is activated by cyclic dinucleotides (CDNs) which eventually leads to expression of type I interferons and other cytokines. Factors underlying the affinity of various CDN analogues are poorly understood. Herein, we correlate structural biology, isothermal calorimetry (ITC) and computational modeling to elucidate factors contributing to binding of six CDNs-three pairs of natural (ribo) and fluorinated (2'-fluororibo) 3',3'-CDNs. X-ray structural analyses of six {STING:CDN} complexes did not offer any explanation for the different affinities of the studied ligands. ITC showed entropy/enthalpy compensation up to 25 kcal mol-1 for this set of similar ligands. The higher affinities of fluorinated analogues are explained with help of computational methods by smaller loss of entropy upon binding and by smaller strain (free) energy.
Collapse
Affiliation(s)
- Miroslav Smola
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Milan Dejmek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Milan Kožíšek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Thomas Evangelidis
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Zahra Aliakbar Tehrani
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Barbora Novotná
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Radim Nencka
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Gabriel Birkuš
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Evzen Boura
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| |
Collapse
|
8
|
Smola M, Gutten O, Dejmek M, Kožíšek M, Evangelidis T, Tehrani ZA, Novotná B, Nencka R, Birkuš G, Rulíšek L, Boura E. Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202016805] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Miroslav Smola
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Milan Dejmek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Milan Kožíšek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Thomas Evangelidis
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Zahra Aliakbar Tehrani
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Barbora Novotná
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Radim Nencka
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Gabriel Birkuš
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Evzen Boura
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| |
Collapse
|
9
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|
10
|
Neese F, Wennmohs F, Becker U, Riplinger C. The ORCA quantum chemistry program package. J Chem Phys 2020; 152:224108. [DOI: 10.1063/5.0004608] [Citation(s) in RCA: 697] [Impact Index Per Article: 174.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
- FAccTs GmbH, Rolandstr. 67, 50677 Köln, Germany
| | - Frank Wennmohs
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Ute Becker
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | |
Collapse
|
11
|
Lang J, Antalík A, Veis L, Brandejs J, Brabec J, Legeza Ö, Pittner J. Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T). J Chem Theory Comput 2020; 16:3028-3040. [PMID: 32275424 DOI: 10.1021/acs.jctc.0c00065] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Sciences, Charles University, Albertov 6, 128 00 Prague 2, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
12
|
Liao S, Zhao M, Luo J, Luo K, Wu J, Liu R, Wang S, Jia P, Bai Y, Zheng X. The interaction mechanism between alkaloids and pepsin based on lum-AuNPs in the chemiluminescence analysis. RSC Adv 2019; 9:25569-25575. [PMID: 35530091 PMCID: PMC9070008 DOI: 10.1039/c9ra02978h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2019] [Accepted: 07/20/2019] [Indexed: 12/14/2022] Open
Abstract
Herein, novel luminol functional gold nanoparticles (lum-AuNPs) were quickly prepared in an alkaline luminol solution with HAuCl4, which had the unique characteristics of uniform size and excellent luminescence properties. A self-made flow injection-chemiluminescence (FI-CL) system was established to study the interaction between pepsin (Pep) and five alkaloids (anisodamine, berberine, reserpine, jatrorrhizine and matrine) using lum-AuNPs as the CL probe. Based on the abovementioned home-made CL system, the possible interaction mechanisms of Pep with five alkaloids have been comprehensively discussed by molecular docking simulation, chemical thermodynamics and kinetic studies. The results indicated that there were obvious CL enhancement and inhibition effects on the lum-AuNPs CL system for the Pep and the complex of Pep/alkaloids, respectively. The possible mechanism for the interaction of Pep–five alkaloids was mainly mediated by the hydrophobic force. The binding constant K and binding site n for the Pep–alkaloid interaction are consistent with the list of Ber > Res > Ani, Jat > Mat, which is relative to the potential of groups of alkaloids interacting with the active site of Pep. A flow injection-chemiluminescence system was established to study the interaction between pepsin and five alkaloids, combined with lum-AuNPs.![]()
Collapse
Affiliation(s)
- Sha Liao
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Meimei Zhao
- School of Pharmacy, Shaanxi Institute of International Trade & Commerce Xi'an 712046 China
| | - Jing Luo
- Shaanxi Traditional Chinese Medicine Hospital Xi'an 710004 China
| | - Kai Luo
- Department of Chemistry, Fudan University Shanghai 200438 China
| | - Jingni Wu
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Ruimin Liu
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Shixiang Wang
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Pu Jia
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Yajun Bai
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| | - Xiaohui Zheng
- Key Laboratory of Resource Biology and Biotechnology in Western China, Northwest University Xi'an 710069 China
| |
Collapse
|
13
|
Nutho B, Mulholland AJ, Rungrotmongkol T. Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate. J Phys Chem B 2019; 123:2889-2903. [PMID: 30845796 DOI: 10.1021/acs.jpcb.9b02157] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.
Collapse
Affiliation(s)
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol BS8 1TS , U.K
| | | |
Collapse
|
14
|
Lang J, Brabec J, Saitow M, Pittner J, Neese F, Demel O. Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method. Phys Chem Chem Phys 2019; 21:5022-5038. [PMID: 30762044 DOI: 10.1039/c8cp03577f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic.
| | | | | | | | | | | |
Collapse
|
15
|
Yingcharoen P, Kongtes C, Arayachukiat S, Suvarnapunya K, Vummaleti SVC, Wannakao S, Cavallo L, Poater A, D' Elia V. Assessing the pKa
-Dependent Activity of Hydroxyl Hydrogen Bond Donors in the Organocatalyzed Cycloaddition of Carbon Dioxide to Epoxides: Experimental and Theoretical Study. Adv Synth Catal 2018. [DOI: 10.1002/adsc.201801093] [Citation(s) in RCA: 73] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Prapussorn Yingcharoen
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
| | - Chutima Kongtes
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
| | - Sunatda Arayachukiat
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
| | - Kittipong Suvarnapunya
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
- Department of Immunology, Faculty of Medicine, Siriraj Hospital; Mahidol University; Bangkok 10700 Thailand
| | - Sai V. C. Vummaleti
- King Abdullah University of Science & Technology; KAUST Catalysis Center (KCC); 23955-6900 Thuwal Saudi Arabia
| | - Sippakorn Wannakao
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
| | - Luigi Cavallo
- King Abdullah University of Science & Technology; KAUST Catalysis Center (KCC); 23955-6900 Thuwal Saudi Arabia
| | - Albert Poater
- Institut de Química Computacional i Catàlisi and Departament de Química; Universitat de Girona, c/Maria Aurèlia Capmany 69; 17003 Girona, Catalonia Spain
| | - Valerio D' Elia
- Department of Materials Science and Engineering, School of Molecular Science and Engineering; Vidyasirimedhi Institute of Science and Technology, (VISTEC); Wang Chan Thailand
| |
Collapse
|
16
|
Thapa B, Beckett D, Erickson J, Raghavachari K. Theoretical Study of Protein–Ligand Interactions Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2018; 14:5143-5155. [DOI: 10.1021/acs.jctc.8b00531] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Jon Erickson
- Lilly Research Laboratories, Eli Lilly & Co., Indianapolis, Indiana 47285, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
17
|
Sárosi MB, Lybrand TP. Molecular Dynamics Simulation of Cyclooxygenase-2 Complexes with Indomethacin closo-Carborane Analogs. J Chem Inf Model 2018; 58:1990-1999. [PMID: 30067351 DOI: 10.1021/acs.jcim.8b00275] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Molecular dynamics simulation of carborane-containing ligands in complex with target enzymes is a challenging task due to the unique structure and properties of the carborane substituents and relative lack of appropriate experimental data to help assess the quality of carborane force field parameters. Here, we report results from energy minimization calculations for a series of carborane-amino acid complexes using carborane force field parameters published previously in the literature and adapted for use with the AMBER ff99SB and ff14SB potential functions. These molecular mechanics results agree well with quantum mechanical geometry optimization calculations obtained using dispersion-corrected density functional theory methods, suggesting that the carborane force field parameters should be suitable for more detailed calculations. We then performed molecular dynamics simulations for the 1,2-, 1,7-, and 1,12-dicarba- closo-dodecaborane(12) derivatives of indomethacin methyl ester bound with cyclooxygenase-2. The simulation results suggest that only the ortho-carborane derivative forms a stable complex, in agreement with experimental findings, and provide insight into the possible molecular basis for isomer binding selectivity.
Collapse
Affiliation(s)
- Menyhárt-Botond Sárosi
- Institute of Inorganic Chemistry, Faculty of Chemistry and Mineralogy , Leipzig University , Johannisallee 29 , D-04103 Leipzig , Germany
| | - Terry P Lybrand
- Departments of Chemistry and Pharmacology, Center for Structural Biology , Vanderbilt University , Nashville , Tennessee 37235-1822 , United States
| |
Collapse
|
18
|
Okimoto N, Otsuka T, Hirano Y, Taiji M. Use of the Multilayer Fragment Molecular Orbital Method to Predict the Rank Order of Protein-Ligand Binding Affinities: A Case Study Using Tankyrase 2 Inhibitors. ACS OMEGA 2018; 3:4475-4485. [PMID: 31458673 PMCID: PMC6641631 DOI: 10.1021/acsomega.8b00175] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/28/2018] [Accepted: 03/23/2018] [Indexed: 06/10/2023]
Abstract
In computational drug discovery, ranking a series of compound analogues in the order that is consistent with the experimental binding affinities remains a challenge. Many of the computational methods available for evaluating binding affinities have adopted molecular mechanics (MM)-based force fields, although they cannot completely describe protein-ligand interactions. By contrast, quantum mechanics (QM) calculations play an important role in understanding the protein-ligand interactions; however, their huge computational costs hinder their application in drug discovery. In this study, we have evaluated the ability to rank the binding affinities of tankyrase 2 ligands by combining both MM and QM calculations. Our computational approach uses the protein-ligand binding energies obtained from a cost-effective multilayer fragment molecular orbital (MFMO) method combined with the solvation energy obtained from the MM-Poisson-Boltzmann/surface area (MM-PB/SA) method to predict the binding affinity. This approach enabled us to rank tankyrase 2 inhibitor analogues, outperforming several MM-based methods, including rescoring by molecular docking and the MM-PB/SA method alone. Our results show that this computational approach using the MFMO method is a promising tool for predicting the rank order of the binding affinities of inhibitor analogues.
Collapse
|
19
|
Chandramouli B, Del Galdo S, Mancini G, Tasinato N, Barone V. Tailor-made computational protocols for precise characterization of small biological building blocks using QM and MM approaches. Biopolymers 2018. [DOI: 10.1002/bip.23109] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Balasubramanian Chandramouli
- Scuola Normale Superiore, Piazza dei Cavalieri 7; Pisa 56126 Italy
- Compunet, Istituto Italiano di Tecnologia, via Morego 30; Genova Italy
| | - Sara Del Galdo
- Scuola Normale Superiore, Piazza dei Cavalieri 7; Pisa 56126 Italy
| | - Giordano Mancini
- Scuola Normale Superiore, Piazza dei Cavalieri 7; Pisa 56126 Italy
- Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3; Pisa 56127 Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7; Pisa 56126 Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7; Pisa 56126 Italy
- Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3; Pisa 56127 Italy
| |
Collapse
|
20
|
Brabec J, Lang J, Saitow M, Pittner J, Neese F, Demel O. Domain-Based Local Pair Natural Orbital Version of Mukherjee’s State-Specific Coupled Cluster Method. J Chem Theory Comput 2018; 14:1370-1382. [DOI: 10.1021/acs.jctc.7b01184] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiri Brabec
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
- Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University in Prague, Hlavova 2030, 12840 Prague 2, Czech Republic
| | - Masaaki Saitow
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
21
|
Assemble-And-Match: A Novel Hybrid Tool for Enhancing Education and Research in Rational Structure Based Drug Design. Sci Rep 2018; 8:849. [PMID: 29339792 PMCID: PMC5770410 DOI: 10.1038/s41598-017-18151-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Accepted: 12/05/2017] [Indexed: 11/08/2022] Open
Abstract
Rational drug design is the process of finding new medication that can activate or inhibit the biofunction of a target molecule by binding to it and forming a molecular complex. Here, shape and charge complementarities between drug and target are key. To help find effective drug molecules out of a huge pool of possibilities, physical and computer aided tools have been developed. Former offers a tangible experience of the molecular interactions yet lacks measurement and evaluation capabilities. Latter enables accurate and fast evaluations, but does not deliver the interactive tangible experience of physical models. We introduce a novel hybrid model called "Assemble-And-Match" where, we enhance and combine the unique features of the two categories. Assemble-And-Match works based on fabrication of customized molecular fragments using our developed software and a 3D printer. Fragments are hinged to each other in different combinations and form flexible peptide chains, conformable to tertiary structures, to fit in the binding pocket of a (3D printed) target molecule. Through embedded measurement marks, the molecular model is reconstructed in silico and its properties are evaluated. We expect Assemble-And-Match tool can enable combination of visuospatial perception with in silico computational power to aid research and education in drug design.
Collapse
|
22
|
Lima Costa AH, Clemente WS, Bezerra KS, Lima Neto JX, Albuquerque EL, Fulco UL. Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists. NEW J CHEM 2018. [DOI: 10.1039/c8nj03521k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We present the energy profiles of estrogen receptor–agonist ligand interactions in atomic detail using a quantum biochemical approach.
Collapse
Affiliation(s)
- Aranthya H. Lima Costa
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - Washington S. Clemente
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - Katyanna S. Bezerra
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - José X. Lima Neto
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | | | - Umberto L. Fulco
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| |
Collapse
|
23
|
Ajani H, Pecina A, Eyrilmez SM, Fanfrlík J, Haldar S, Řezáč J, Hobza P, Lepšík M. Superior Performance of the SQM/COSMO Scoring Functions in Native Pose Recognition of Diverse Protein-Ligand Complexes in Cognate Docking. ACS OMEGA 2017; 2:4022-4029. [PMID: 30023710 PMCID: PMC6044937 DOI: 10.1021/acsomega.7b00503] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/24/2017] [Accepted: 07/18/2017] [Indexed: 06/08/2023]
Abstract
General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.
Collapse
Affiliation(s)
- Haresh Ajani
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
- Department
of Physical Chemistry, Palacký University, tř. 17. listopadu 1192/12, 77146 Olomouc, Czech Republic
| | - Adam Pecina
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
| | - Saltuk M. Eyrilmez
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
- Department
of Physical Chemistry, Palacký University, tř. 17. listopadu 1192/12, 77146 Olomouc, Czech Republic
| | - Jindřich Fanfrlík
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
| | - Susanta Haldar
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
| | - Jan Řezáč
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
| | - Pavel Hobza
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
- Department
of Physical Chemistry, Regional Centre of Advanced Technologies and
Materials, Palacký University, 77146 Olomouc, Czech Republic
| | - Martin Lepšík
- Department
of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, v.v.i., Flemingovo nam. 2, 16610 Praha 6, Czech Republic
| |
Collapse
|
24
|
Sárosi MB, Neumann W, Lybrand TP, Hey-Hawkins E. Molecular Modeling of the Interactions between Carborane-Containing Analogs of Indomethacin and Cyclooxygenase-2. J Chem Inf Model 2017. [DOI: 10.1021/acs.jcim.7b00113] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Menyhárt-Botond Sárosi
- Institute
of Inorganic Chemistry, Faculty of Chemistry and Mineralogy, Leipzig University, Johannisallee 29, D-04103 Leipzig, Germany
| | - Wilma Neumann
- Institute
of Inorganic Chemistry, Faculty of Chemistry and Mineralogy, Leipzig University, Johannisallee 29, D-04103 Leipzig, Germany
| | - Terry P. Lybrand
- Center
for Structural Biology, Departments of Chemistry and Pharmacology, Vanderbilt University, Nashville, Tennessee 37232-8725, United States
| | - Evamarie Hey-Hawkins
- Institute
of Inorganic Chemistry, Faculty of Chemistry and Mineralogy, Leipzig University, Johannisallee 29, D-04103 Leipzig, Germany
| |
Collapse
|
25
|
Grimme S, Bannwarth C, Shushkov P. A Robust and Accurate Tight-Binding Quantum Chemical Method for Structures, Vibrational Frequencies, and Noncovalent Interactions of Large Molecular Systems Parametrized for All spd-Block Elements (Z = 1–86). J Chem Theory Comput 2017; 13:1989-2009. [DOI: 10.1021/acs.jctc.7b00118] [Citation(s) in RCA: 702] [Impact Index Per Article: 100.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical
Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical
Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Philip Shushkov
- Mulliken Center for Theoretical
Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
26
|
Pecina A, Haldar S, Fanfrlík J, Meier R, Řezáč J, Lepšík M, Hobza P. SQM/COSMO Scoring Function at the DFTB3-D3H4 Level: Unique Identification of Native Protein–Ligand Poses. J Chem Inf Model 2017; 57:127-132. [DOI: 10.1021/acs.jcim.6b00513] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Adam Pecina
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Susanta Haldar
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Jindřich Fanfrlík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - René Meier
- Leibniz Institute of Plant Biochemistry, Weinberg 3, 06120 Halle, Germany
| | - Jan Řezáč
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Martin Lepšík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Pavel Hobza
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v.v.i., Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Palacký University, 77146 Olomouc, Czech Republic
| |
Collapse
|
27
|
Stasyuk OA, Jakubec D, Vondrášek J, Hobza P. Noncovalent Interactions in Specific Recognition Motifs of Protein-DNA Complexes. J Chem Theory Comput 2017; 13:877-885. [PMID: 27992205 DOI: 10.1021/acs.jctc.6b00775] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
In view of the importance of protein-DNA interactions in biological processes, we extracted from the Protein Data Bank several one-to-one complexes of amino acids with nucleotides that matched certain geometric and energetic specificity criteria and investigated them using quantum chemistry methods. The CCSD(T)/CBS interaction energies were used as a benchmark to compare the performance of the MP2.5, MP2-F12, DFT-D3, and PM6-D3H4 methods. All methods yielded good agreement with the reference values, with declining accuracy from MP2.5 to PM6-D3H4. Regardless of the site of interaction, the minima found after full optimization in implicit solvent with high dielectric constant were close to the structures experimentally detected in protein-DNA complexes. According to DFT-SAPT analysis, the nature of noncovalent interactions strongly depends on the type of amino acid. The negatively charged sugar-phosphate backbone of DNA heavily influences the strength of interactions and must be included in the computational model, especially in the case of interactions with charged amino acids.
Collapse
Affiliation(s)
- Olga A Stasyuk
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - David Jakubec
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University in Prague , Albertov 6, 128 43 Prague, Czech Republic
| | - Jiří Vondrášek
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University , 771 46 Olomouc, Czech Republic
| |
Collapse
|
28
|
Kokoschka M, Galgonek J, Vondrasek J, Hobza P. Computational methods for the description of pharmacologically relevant platinum complexes--molecular structure and bond dissociation. Phys Chem Chem Phys 2016; 18:4051-62. [PMID: 26777459 DOI: 10.1039/c5cp03391h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cancer is after cardiovascular disease the most frequent cause of death in Europe. In 28 of 53 countries considered in this area it is already the leading cause of death and expected to gain even more importance until the year 2020. Amongst the large arsenal of different anti-cancer drugs, platinum drugs belong to the first developed anticancer drugs and still have a large impact on cancer therapy. Nevertheless therapy with platinum-anticancer drugs is accompanied by severe adverse effects caused by frequent interactions with the amino acids of different human proteins. Computational chemistry offers methods to study such interactions and even those of not yet synthesized drugs in silico. For such studies a profound knowledge of the prediction quality of various computational methods towards platinum-drug-like complexes is necessary. By this article we are aiming on delivering important accuracy information of the frequently used computational methods. Most important findings are the high performance of the double hybrid functional B2PLYP for the calculation of geometries, even in small basis sets, followed by BP86 and PBE and the still acceptable performance of the semi-empirical Method PM6-D3H4X for extremely large systems. To follow absolute energies of the dissociation process, LPNO-CEPA and B3LYP-D3 can be suggested while SCS-MP2 shows an extremely narrow standard deviation and a low maximum error, which make it an ideal candidate for relative energy calculations in the exploration of reaction mechanisms.
Collapse
Affiliation(s)
- Malte Kokoschka
- Gilead Sciences & IOCB Research Center, Department of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic. and Department of Physical Chemistry, Palacky University, 771 46 Olomouc, Czech Republic
| | - Jakub Galgonek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic
| | - Jiri Vondrasek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic
| | - Pavel Hobza
- Gilead Sciences & IOCB Research Center, Department of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic. and Department of Physical Chemistry, Palacky University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
29
|
Kuechler ER, Giese TJ, York DM. Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations. J Chem Phys 2016; 143:234111. [PMID: 26696050 DOI: 10.1063/1.4937166] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom's local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the SN2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM/MM interactions.
Collapse
Affiliation(s)
- Erich R Kuechler
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, USA
| | - Timothy J Giese
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, USA
| | - Darrin M York
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, USA
| |
Collapse
|
30
|
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
|
31
|
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
|
32
|
Berg L, Mishra BK, Andersson CD, Ekström F, Linusson A. The Nature of Activated Non-classical Hydrogen Bonds: A Case Study on Acetylcholinesterase-Ligand Complexes. Chemistry 2016; 22:2672-81. [DOI: 10.1002/chem.201503973] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Indexed: 01/25/2023]
Affiliation(s)
- Lotta Berg
- Department of Chemistry; Umeå University; 901 87 Umeå Sweden
| | | | | | - Fredrik Ekström
- CBRN Defense and Security; Swedish Defense Research Agency; 906 21 Umeå Sweden
| | - Anna Linusson
- Department of Chemistry; Umeå University; 901 87 Umeå Sweden
| |
Collapse
|
33
|
Dixit M, Das S, Mhashal A, Eitan R, Major D. Practical Aspects of Multiscale Classical and Quantum Simulations of Enzyme Reactions. Methods Enzymol 2016; 577:251-86. [DOI: 10.1016/bs.mie.2016.05.046] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
34
|
Pecina A, Meier R, Fanfrlík J, Lepšík M, Řezáč J, Hobza P, Baldauf C. The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein–ligand interactions and implicit COSMO solvation. Chem Commun (Camb) 2016; 52:3312-5. [DOI: 10.1039/c5cc09499b] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Strictly uphill – in cognate docking experiments we show that a quantum mechanical description of interaction and solvation outperforms established scoring functions in sharply distinguishing the native state from decoy poses.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry and Biochemistry (IOCB) and Gilead Sciences and IOCB Research Center
- 16610 Prague 6
- Czech Republic
| | - René Meier
- Institut für Biochemie
- Fakultät für Biowissenschaften
- Pharmazie und Psychologie
- Universität Leipzig
- D-04109 Leipzig
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry (IOCB) and Gilead Sciences and IOCB Research Center
- 16610 Prague 6
- Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry and Biochemistry (IOCB) and Gilead Sciences and IOCB Research Center
- 16610 Prague 6
- Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry (IOCB) and Gilead Sciences and IOCB Research Center
- 16610 Prague 6
- Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry (IOCB) and Gilead Sciences and IOCB Research Center
- 16610 Prague 6
- Czech Republic
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft
- D-14195 Berlin
- Germany
| |
Collapse
|
35
|
Rodríguez JI, Matta CF, Uribe EA, Götz AW, Castillo-Alvarado F, Molina-Brito B. A QTAIM topological analysis of the P3HTPCBM dimer. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2015.11.052] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
36
|
Otsuka T, Okimoto N, Taiji M. Assessment and acceleration of binding energy calculations for protein-ligand complexes by the fragment molecular orbital method. J Comput Chem 2015; 36:2209-18. [DOI: 10.1002/jcc.24055] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2015] [Revised: 07/17/2015] [Accepted: 07/30/2015] [Indexed: 12/16/2022]
Affiliation(s)
- Takao Otsuka
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| | - Noriaki Okimoto
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| | - Makoto Taiji
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| |
Collapse
|
37
|
Genheden S, Ryde U, Söderhjelm P. Binding affinities by alchemical perturbation using QM/MM with a large QM system and polarizable MM model. J Comput Chem 2015; 36:2114-24. [PMID: 26280564 DOI: 10.1002/jcc.24048] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Revised: 06/13/2015] [Accepted: 07/17/2015] [Indexed: 12/19/2022]
Abstract
The most general way to improve the accuracy of binding-affinity calculations for protein-ligand systems is to use quantum-mechanical (QM) methods together with rigorous alchemical-perturbation (AP) methods. We explore this approach by calculating the relative binding free energy of two synthetic disaccharides binding to galectin-3 at a reasonably high QM level (dispersion-corrected density functional theory with a triple-zeta basis set) and with a sufficiently large QM system to include all short-range interactions with the ligand (744-748 atoms). The rest of the protein is treated as a collection of atomic multipoles (up to quadrupoles) and polarizabilities. Several methods for evaluating the binding free energy from the 3600 QM calculations are investigated in terms of stability and accuracy. In particular, methods using QM calculations only at the endpoints of the transformation are compared with the recently proposed non-Boltzmann Bennett acceptance ratio (NBB) method that uses QM calculations at several stages of the transformation. Unfortunately, none of the rigorous approaches give sufficient statistical precision. However, a novel approximate method, involving the direct use of QM energies in the Bennett acceptance ratio method, gives similar results as NBB but with better precision, ∼3 kJ/mol. The statistical error can be further reduced by performing a greater number of QM calculations.
Collapse
Affiliation(s)
- Samuel Genheden
- School of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, United Kingdom
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| |
Collapse
|
38
|
Demel O, Pittner J, Neese F. A Local Pair Natural Orbital-Based Multireference Mukherjee’s Coupled Cluster Method. J Chem Theory Comput 2015; 11:3104-14. [DOI: 10.1021/acs.jctc.5b00334] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ondřej Demel
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Pittner
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institute
of Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
39
|
Yilmazer ND, Heitel P, Schwabe T, Korth M. Benchmark of electronic structure methods for protein–ligand interactions based on high-level reference data. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s0219633615400015] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The accurate prediction of the strength of protein–ligand interactions is a very difficult problem despite impressive advances in the field of biomolecular modeling. There are good reasons to believe that quantum mechanical methods can help with this task, but the application of such methods in the context of scoring is still in its infancy. Here we benchmark several wave function theory (WFT), density functional theory (DFT) and semiempirical quantum mechanical (SQM) approaches against high-level theoretical references for realistic test cases. Based on our findings for systematically generated model systems of real protein/ligand complexes from the PDB-bind database, we can recommend SCS-MP2 and B2-PLYP-D3 as reference methods, TPSS-D3+Dabc/def-TZVPP as the best DFT approach and PM6-DH+ as a fast and accurate alternative to full ab initio treatments.
Collapse
Affiliation(s)
- Nusret Duygu Yilmazer
- Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany
| | - Pascal Heitel
- Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany
| | - Tobias Schwabe
- Center for Bioinformatics and Institute of Physical Chemistry, University of Hamburg, Bundesstraße 43, 20146 Hamburg, Germany
| | - Martin Korth
- Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany
| |
Collapse
|
40
|
Dantas DS, Oliveira JIN, Lima Neto JX, da Costa RF, Bezerra EM, Freire VN, Caetano EWS, Fulco UL, Albuquerque EL. Quantum molecular modelling of ibuprofen bound to human serum albumin. RSC Adv 2015. [DOI: 10.1039/c5ra04395f] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The total interaction energies of the ibuprofen complexed with FA3/FA4 and FA6 binding sites of human serum albumin are in agreement with the hypothesis that the Sudlow's site II is the main binding pocket for ibuprofen.
Collapse
Affiliation(s)
- Diego S. Dantas
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal
- Brazil
| | - Jonas I. N. Oliveira
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal
- Brazil
| | - José X. Lima Neto
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal
- Brazil
| | - Roner F. da Costa
- Departamento de Física
- Universidade Federal Rural do Semi-Árido
- Brazil
| | - Eveline M. Bezerra
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal
- Brazil
| | - Valder N. Freire
- Departamento de Física
- Universidade Federal do Ceará
- Fortaleza
- Brazil
| | | | - Umberto L. Fulco
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal
- Brazil
| | | |
Collapse
|
41
|
Kumar M, Busch DH, Subramaniam B, Thompson WH. Role of Tunable Acid Catalysis in Decomposition of α-Hydroxyalkyl Hydroperoxides and Mechanistic Implications for Tropospheric Chemistry. J Phys Chem A 2014; 118:9701-11. [DOI: 10.1021/jp505100x] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Manoj Kumar
- Department
of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| | - Daryle H. Busch
- Department
of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| | - Bala Subramaniam
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
- Department
of Chemical and Petroleum Engineering, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H. Thompson
- Department
of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| |
Collapse
|
42
|
Pandey KK. Shared electron versus donor–acceptor bonding description of Fe–ER bonds in [(η5-C5H5)(L)2Fe(ER)] (L = CO, PMe3; E = Si, Ge, Sn, Pb; R = Ph, Me). J Organomet Chem 2014. [DOI: 10.1016/j.jorganchem.2014.03.024] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
43
|
Kumar M, Busch DH, Subramaniam B, Thompson WH. Organic Acids Tunably Catalyze Carbonic Acid Decomposition. J Phys Chem A 2014; 118:5020-8. [DOI: 10.1021/jp5037469] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Manoj Kumar
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| | - Daryle H. Busch
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| | - Bala Subramaniam
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| | - Ward H. Thompson
- Center for Environmentally Beneficial Catalysis, 1501 Wakarusa Drive, Lawrence, Kansas 66047, United States
| |
Collapse
|
44
|
Kasende OE, Matondo A, Muzomwe M, Muya JT, Scheiner S. Interaction between temozolomide and water: Preferred binding sites. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.02.005] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
45
|
Michel J. Current and emerging opportunities for molecular simulations in structure-based drug design. Phys Chem Chem Phys 2014; 16:4465-77. [PMID: 24469595 PMCID: PMC4256725 DOI: 10.1039/c3cp54164a] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2013] [Accepted: 01/10/2014] [Indexed: 01/29/2023]
Abstract
An overview of the current capabilities and limitations of molecular simulation of biomolecular complexes in the context of computer-aided drug design is provided. Steady improvements in computer hardware coupled with more refined representations of energetics are leading to a new appreciation of the driving forces of molecular recognition. Molecular simulations are poised to more frequently guide the interpretation of biophysical measurements of biomolecular complexes. Ligand design strategies emerge from detailed analyses of computed structural ensembles. The feasibility of routine applications to ligand optimization problems hinges upon successful extensive large scale validation studies and the development of protocols to intelligently automate computations.
Collapse
Affiliation(s)
- Julien Michel
- EaStCHEM School of Chemistry, Joseph Black Building, The King's Buildings, Edinburgh, EH9 3JJ, UK.
| |
Collapse
|
46
|
Roos K, Viklund J, Meuller J, Kaspersson K, Svensson M. Potency prediction of β-secretase (BACE-1) inhibitors using density functional methods. J Chem Inf Model 2014; 54:818-25. [PMID: 24456077 DOI: 10.1021/ci400374z] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Scoring potency is a main challenge for structure based drug design. Inductive effects of subtle variations in the ligand are not possible to accurately predict by classical computational chemistry methods. In this study, the problem of predicting potency of ligands with electronic variations participating in key interactions with the protein was addressed. The potency was predicted for a large set of cyclic amidine and guanidine cores extracted from β-secretase (BACE-1) inhibitors. All cores were of similar size and had equal interaction motifs but were diverse with respect to electronic substitutions. A density functional theory approach, in combination with a representation of the active site of a protein using only key residues, was shown to be predictive. This computational approach was used to guide and support drug design, within the time frame of a normal drug discovery design cycle.
Collapse
Affiliation(s)
- Katarina Roos
- Department of Medicinal Chemistry and §Discovery Sciences, AstraZeneca R&D Mölndal , SE-431 83 Mölndal, Sweden
| | | | | | | | | |
Collapse
|
47
|
Moreira da Costa L, Stoyanov SR, Gusarov S, Seidl PR, Walkimar de M. Carneiro J, Kovalenko A. Computational Study of the Effect of Dispersion Interactions on the Thermochemistry of Aggregation of Fused Polycyclic Aromatic Hydrocarbons as Model Asphaltene Compounds in Solution. J Phys Chem A 2014; 118:896-908. [DOI: 10.1021/jp408005h] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Leonardo Moreira da Costa
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Organic Processes, School of Chemistry, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-901, Brazil
- Institute
of Chemistry, Universidade Federal Fluminense, Niteroi, RJ 24220-900, Brazil
- Centro Universitário da Zona Oeste, Avenida Manuel Caldeira de Alvarenga 1203, Campo Grande, RJ 23070-200, Brazil
| | - Stanislav R. Stoyanov
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 2V4, Canada
| | - Sergey Gusarov
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Peter R. Seidl
- Department
of Organic Processes, School of Chemistry, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-901, Brazil
| | | | - Andriy Kovalenko
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| |
Collapse
|
48
|
Sparta M, Neese F. Chemical applications carried out by local pair natural orbital based coupled-cluster methods. Chem Soc Rev 2014; 43:5032-41. [DOI: 10.1039/c4cs00050a] [Citation(s) in RCA: 106] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Local pair natural orbital coupled-cluster methods are suitable for accurate assessment of the energies of large molecular and supramolecular systems.
Collapse
Affiliation(s)
- Manuel Sparta
- Max Planck Institute for Chemical Energy Conversion
- D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion
- D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
49
|
Maurer SA, Beer M, Lambrecht DS, Ochsenfeld C. Linear-scaling symmetry-adapted perturbation theory with scaled dispersion. J Chem Phys 2013; 139:184104. [DOI: 10.1063/1.4827297] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Ehrlich S, Bettinger HF, Grimme S. Dispersion-Driven Conformational Isomerism in σ-Bonded Dimers of Larger Acenes. Angew Chem Int Ed Engl 2013; 52:10892-5. [DOI: 10.1002/anie.201304674] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Indexed: 11/06/2022]
|