1
|
Sharma H, Raju B, Narendra G, Motiwale M, Sharma B, Verma H, Silakari O. QM/MM Studies on Enzyme Catalysis and Insight into Designing of New Inhibitors by ONIOM Approach: Recent Update. ChemistrySelect 2023. [DOI: 10.1002/slct.202203319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Himani Sharma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Baddipadige Raju
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Gera Narendra
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Mohit Motiwale
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Bhavna Sharma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Himanshu Verma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Om Silakari
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| |
Collapse
|
2
|
Rozza AM, Papp M, McFarlane NR, Harvey JN, Oláh J. The Mechanism of Biochemical NO‐Sensing: Insights from Computational Chemistry. Chemistry 2022; 28:e202200930. [PMID: 35670519 PMCID: PMC9542423 DOI: 10.1002/chem.202200930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Indexed: 11/22/2022]
Abstract
The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO‐receptor in humans is soluble guanylyl cyclase (sGC) – a histidine‐ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO‐signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high‐level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H‐NOX protein from Shewanella oneidensis (So H‐NOX) homologous to the NO‐binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six‐coordinate intermediate heme‐NO complex, followed by relatively facile His decoordination yielding a five‐coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.
Collapse
Affiliation(s)
- Ahmed M. Rozza
- Department of Inorganic and Analytical Chemistry Budapest University of Technology and Economics 1111 Budapest Műegyetem rakpart 3. Hungary
- Department of Biotechnology Faculty of Agriculture Al-Azhar University Cairo 11651 Egypt
| | - Marcell Papp
- Department of Inorganic and Analytical Chemistry Budapest University of Technology and Economics 1111 Budapest Műegyetem rakpart 3. Hungary
| | - Neil R. McFarlane
- Department of Chemistry KU Leuven 3001 Leuven Celestijnenlaan 200 f- box 2404 Belgium
| | - Jeremy N. Harvey
- Department of Chemistry KU Leuven 3001 Leuven Celestijnenlaan 200 f- box 2404 Belgium
| | - Julianna Oláh
- Department of Inorganic and Analytical Chemistry Budapest University of Technology and Economics 1111 Budapest Műegyetem rakpart 3. Hungary
| |
Collapse
|
3
|
Zhang S, Li X, Wang Y, Wei J, Zhang X, Liu Y. Computational Study of the Peroxygenase Mechanism Catalyzed by Hemoglobin Dehaloperoxidase Involved in the Degradation of Chlorophenols. Inorg Chem 2022; 61:2628-2639. [PMID: 35080380 DOI: 10.1021/acs.inorgchem.1c03721] [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/29/2022]
Abstract
The biochemical evidence showed that hemoglobin dehaloperoxidase (DHP B) from Amphitrite Ornata is a multifunctional hemoprotein that catalyzes both dehalogenation and hydroxylation of halophenols via the peroxidase and peroxygenase mechanism, respectively, which sets the basis for the degradation of halophenols. In the peroxygenase mechanism, the reaction was previously suggested to be triggered either by the hydrogen atom abstraction by the Fe═O center or by the proton abstraction by His55. To illuminate the peroxygenase mechanism of DHP B at the atomistic level, on the basis of the high-resolution crystal structure, computational models were constructed, and a series of quantum mechanical/molecular mechanical calculations have been performed. According to the calculation results, the pathway (Path a) initiated by the H-abstraction by the Fe═O center is feasible. In another pathway (Path b), His55 can abstract the proton from the hydroxyl group of the substrate (4-Cl-o-cresol) to initiate the reaction; however, its feasibility depends on the prior electron transfer from the substrate to the porphyrin group. The rate-limiting step of Path a is the OH-rebound, which corresponds to an energy barrier of 14.7 kcal/mol at the quartet state. His55 acts as an acid-base catalyst and directly involves in the catalysis. Our mutant study indicates that His55 can be replaced by other titratable residues. These findings may provide useful information for further understanding of the catalytic reaction of DHP B and for the design of enzymes in the degradation of pollutants, in particular, halophenols.
Collapse
Affiliation(s)
- Shiqing Zhang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| | - Xinyi Li
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| | - Yijing Wang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| | - Jingjing Wei
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| | - Xue Zhang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| | - Yongjun Liu
- School of Chemistry and Chemical Engineering, Shandong University, Jinan, Shandong 250100, China
| |
Collapse
|
4
|
Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences. Top Catal 2021; 65:165-186. [DOI: 10.1007/s11244-021-01521-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
5
|
Zalloum WA, Zalloum N. Comparative QM/MM Molecular Dynamics and Umbrella Sampling Simulations: Interaction of the Zinc-Bound Intermediate Gem-Diolate Trapoxin A Inhibitor and Acetyl-l-lysine Substrate with Histone Deacetylase 8. J Phys Chem B 2021; 125:5321-5337. [PMID: 33998791 DOI: 10.1021/acs.jpcb.1c01696] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Targeting the genetic material without destruction is a priority to develop safe anticancer drugs. Histone deacetylase 8 (HDAC8), which is proved to be involved in carcinogenesis, is an enzyme associated with the chromatin for post-translational deacetylation of acetylated lysine. In this study, HDAC8 co-crystallized with the intermediate state tetrapeptide Trapoxin A (TA) inhibitor and the holoenzyme are utilized to find their conformational ensembles. Furthermore, the co-crystallized intermediate gem-diolate TA was used to find optimum interactions with the active site residues by conventional molecular dynamics (MD) simulation and QM/MM umbrella sampling. Finally, the intermediate state of the acetyl-l-lysine substrate was explored by QM/MM steered MD and compared to the binding of the intermediate state of the inhibitor. This research showed that HDAC8 is flexible and exists in conformational ensembles in its holoenzyme state. Binding of the intermediate state TA stabilizes its conformation. The optimum binding to the active site of HDAC8 for structures of gem-diolate TA (intermediate state) and acetyl-l-lysine (intermediate state) was determined according to the corresponding energy profiles. The use of these models will aid in the design of potentially reversible, potent, and selective inhibitors of HDAC8 for cancer treatment.
Collapse
Affiliation(s)
- Waleed A Zalloum
- Department of Pharmacy, Faculty of Health Science, American University of Madaba, P.O. Box 2882, Amman 11821, Jordan
| | - Needa Zalloum
- Department of Biopharmaceutics and Clinical Pharmacy, Faculty of Pharmacy, University of Jordan, Amman 11942, Jordan
| |
Collapse
|
6
|
Piccini G, Parrinello M. Accurate Quantum Chemical Free Energies at Affordable Cost. J Phys Chem Lett 2019; 10:3727-3731. [PMID: 31244270 DOI: 10.1021/acs.jpclett.9b01301] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Free energy sampling methods allow studying the full dynamics of activated processes. Unfortunately, the affordable accuracy of the potential describing the energy and forces of the system is usually rather low. Here we introduce a new method that by combining metadynamics and free energy perturbation allows calculating accurate quantum chemical free energies for chemical reactions. To prove the effectiveness of this new approach we study the SN2 reaction of CH3F + Cl- → CH3Cl + F- in vacuo and solvated by water. Comparisons are made with harmonic transition-state theory to show how this method could provide accurate equilibrium and rate constants for complex systems.
Collapse
Affiliation(s)
- GiovanniMaria Piccini
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della SvizzeraItaliana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della SvizzeraItaliana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genova , Italy
| |
Collapse
|
7
|
Böckers M, Franke R, Staemmler V. A heuristic estimate of molecular correlation energies using pair correlation energies of localized molecular orbitals. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2422-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
8
|
Limb MAL, Suardíaz R, Grant IM, Mulholland AJ. Quantum Mechanics/Molecular Mechanics Simulations Show Saccharide Distortion is Required for Reaction in Hen Egg-White Lysozyme. Chemistry 2019; 25:764-768. [PMID: 30347479 DOI: 10.1002/chem.201805250] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Indexed: 01/07/2023]
Abstract
Hybrid quantum mechanics/molecular mechanics (QM/MM) calculations on lysozyme show significant distortion of the bound saccharide is required to facilitate the catalytic reaction.
Collapse
Affiliation(s)
- Michael A L Limb
- School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS81TS, UK
| | - Reynier Suardíaz
- School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS81TS, UK
| | - Ian M Grant
- Present address: Institute of Psychiatry, Psychology & Neuroscience, King's College London, Denmark Hill, London, SE5 8AF, UK
| | - Adrian J Mulholland
- School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS81TS, UK
| |
Collapse
|
9
|
Proppe J, Husch T, Simm GN, Reiher M. Uncertainty quantification for quantum chemical models of complex reaction networks. Faraday Discuss 2018; 195:497-520. [PMID: 27730243 DOI: 10.1039/c6fd00144k] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
For the quantitative understanding of complex chemical reaction mechanisms, it is, in general, necessary to accurately determine the corresponding free energy surface and to solve the resulting continuous-time reaction rate equations for a continuous state space. For a general (complex) reaction network, it is computationally hard to fulfill these two requirements. However, it is possible to approximately address these challenges in a physically consistent way. On the one hand, it may be sufficient to consider approximate free energies if a reliable uncertainty measure can be provided. On the other hand, a highly resolved time evolution may not be necessary to still determine quantitative fluxes in a reaction network if one is interested in specific time scales. In this paper, we present discrete-time kinetic simulations in discrete state space taking free energy uncertainties into account. The method builds upon thermo-chemical data obtained from electronic structure calculations in a condensed-phase model. Our kinetic approach supports the analysis of general reaction networks spanning multiple time scales, which is here demonstrated for the example of the formose reaction. An important application of our approach is the detection of regions in a reaction network which require further investigation, given the uncertainties introduced by both approximate electronic structure methods and kinetic models. Such cases can then be studied in greater detail with more sophisticated first-principles calculations and kinetic simulations.
Collapse
Affiliation(s)
- Jonny Proppe
- Laboratory of Physical Chemistry, ETH Zürich, Zürich, Switzerland.
| | - Tamara Husch
- Laboratory of Physical Chemistry, ETH Zürich, Zürich, Switzerland.
| | - Gregor N Simm
- Laboratory of Physical Chemistry, ETH Zürich, Zürich, Switzerland.
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Zürich, Switzerland.
| |
Collapse
|
10
|
Rinaldi S, Van der Kamp MW, Ranaghan KE, Mulholland AJ, Colombo G. Understanding Complex Mechanisms of Enzyme Reactivity: The Case of Limonene-1,2-Epoxide Hydrolases. ACS Catal 2018. [DOI: 10.1021/acscatal.8b00863] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Silvia Rinaldi
- Istituto di Chimica del Riconoscimento Molecolare, C.N.R., Via Mario Bianco 9, 20131 Milano, Italy
| | - Marc W. Van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
- School of Biochemistry, University of Bristol, University Walk, Bristol BS8 1TD, United Kingdom
| | - Kara E. Ranaghan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Giorgio Colombo
- Istituto di Chimica del Riconoscimento Molecolare, C.N.R., Via Mario Bianco 9, 20131 Milano, Italy
- Dipartimento di Chimica, Università degli Studi di Pavia, Via Taramelli 12, 27100 Pavia, Italy
| |
Collapse
|
11
|
Wang XX, Yuan L, Jia CX, Qu HJ, Li BJ, Chi YJ, Yu HT. A combined density functional theory and numerical simulation investigation of levels of chirality transfer and regioselectivity for the radical cyclizations of N-methyl-, N-ethyl- and N-isopropyl-substituted ortho-halo- N-acryloylanilides. NEW J CHEM 2018. [DOI: 10.1039/c8nj01102h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
There is a high chirality transfer ratio when the N-alkyl substituent is larger in volume than isopropyl.
Collapse
Affiliation(s)
- Xiao-xu Wang
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Lang Yuan
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Cai-xin Jia
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Hong-jie Qu
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Bai-jian Li
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Yu-juan Chi
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| | - Hai-tao Yu
- Key Laboratory of Functional Inorganic Material Chemistry (Ministry of Education of China)
- School of Chemistry and Materials Science
- Heilongjiang University
- Harbin 150080
- China
| |
Collapse
|
12
|
Alavi FS, Zahedi M, Safari N, Ryde U. QM/MM Study of the Conversion of Oxophlorin into Verdoheme by Heme Oxygenase. J Phys Chem B 2017; 121:11427-11436. [PMID: 29090581 DOI: 10.1021/acs.jpcb.7b08332] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Heme oxygenase is an enzyme that degrades heme, thereby recycling iron in most organisms, including humans. Pervious density functional theory (DFT) calculations have suggested that iron(III) hydroxyheme, an intermediate generated in the first step of heme degradation by heme oxygenase, is converted to iron(III) superoxo oxophlorin in the presence of dioxygen. In this article, we have studied the detailed mechanism of conversion of iron(III) superoxo oxophlorin to verdoheme by using combined quantum mechanics and molecular mechanics (QM/MM) calculations. The calculations employed the B3LYP method and the def2-QZVP basis set, considering dispersion effects with the DFT-D3 approach, obtaining accurate energies with large QM regions of almost 1000 atoms. The reaction was found to be exothermic by -35 kcal/mol, with a rate-determining barrier of 19 kcal/mol in the doublet state. The protein environment and especially water in the enzyme pocket significantly affects the reaction by decreasing the reaction activation energies and changing the structures by providing strategic hydrogen bonds.
Collapse
Affiliation(s)
- Fatemeh Sadat Alavi
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Mansour Zahedi
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Nasser Safari
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University , Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
13
|
Pruitt SR, Steinmann C. Mapping Interaction Energies in Chorismate Mutase with the Fragment Molecular Orbital Method. J Phys Chem A 2017; 121:1797-1807. [DOI: 10.1021/acs.jpca.6b12830] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Spencer R. Pruitt
- Academic & Research Computing, Worcester Polytechnic Institute, Worcester, Massachusetts 01602, United States
| | - Casper Steinmann
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
14
|
Wang X, Zhu W, Liu Y. Tryptophan lyase (NosL): mechanistic insights into amine dehydrogenation and carboxyl fragment migration by QM/MM calculations. Catal Sci Technol 2017. [DOI: 10.1039/c7cy00573c] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
QM/MM calculations suggest two feasible pathways for the breaking of the C–C bond of the substrate. The breaking of the Cα–Cβ bond leads to the final product, whereas the cleavage of the Cα–C bond will terminate in the EPR-trapped radical intermediate.
Collapse
Affiliation(s)
- Xiya Wang
- School of Chemistry and Chemical Engineering
- Shandong University
- Jinan
- China
| | - Wenyou Zhu
- College of Chemistry and Chemical Engineering
- Xuzhou Institute of Technology
- Xuzhou
- China
| | - Yongjun Liu
- School of Chemistry and Chemical Engineering
- Shandong University
- Jinan
- China
| |
Collapse
|
15
|
Dhoke GV, Ensari Y, Davari MD, Ruff AJ, Schwaneberg U, Bocola M. What's My Substrate? Computational Function Assignment of Candida parapsilosis ADH5 by Genome Database Search, Virtual Screening, and QM/MM Calculations. J Chem Inf Model 2016; 56:1313-23. [PMID: 27387009 DOI: 10.1021/acs.jcim.6b00076] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Zinc-dependent medium chain reductase from Candida parapsilosis can be used in the reduction of carbonyl compounds to pharmacologically important chiral secondary alcohols. To date, the nomenclature of cpADH5 is differing (CPCR2/RCR/SADH) in the literature, and its natural substrate is not known. In this study, we utilized a substrate docking based virtual screening method combined with KEGG, MetaCyc pathway, and Candida genome databases search for the discovery of natural substrates of cpADH5. The virtual screening of 7834 carbonyl compounds from the ZINC database provided 94 aldehydes or methyl/ethyl ketones as putative carbonyl substrates. Out of which, 52 carbonyl substrates of cpADH5 with catalytically active docking pose were identified by employing mechanism based substrate docking protocol. Comparison of the virtual screening results with KEGG, MetaCyc database search, and Candida genome pathway analysis suggest that cpADH5 might be involved in the Ehrlich pathway (reduction of fusel aldehydes in leucine, isoleucine, and valine degradation). Our QM/MM calculations and experimental activity measurements affirmed that butyraldehyde substrates are the potential natural substrates of cpADH5, suggesting a carbonyl reductase role for this enzyme in butyraldehyde reduction in aliphatic amino acid degradation pathways. Phylogenetic tree analysis of known ADHs from Candida albicans shows that cpADH5 is close to caADH5. We therefore propose, according to the experimental substrate identification and sequence similarity, the common name butyraldehyde dehydrogenase cpADH5 for Candida parapsilosis CPCR2/RCR/SADH.
Collapse
Affiliation(s)
- Gaurao V Dhoke
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany
| | - Yunus Ensari
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany
| | - Mehdi D Davari
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany
| | - Anna Joëlle Ruff
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany
| | - Ulrich Schwaneberg
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany.,DWI-Leibniz Institut für Interaktive Materialien , Forckenbeckstraße 50, 52056 Aachen, Germany
| | - Marco Bocola
- Lehrstuhl für Biotechnologie, RWTH Aachen University , Worringerweg 3, 52074 Aachen, Germany
| |
Collapse
|
16
|
Kats D. Speeding up local correlation methods: System-inherent domains. J Chem Phys 2016; 145:014103. [DOI: 10.1063/1.4954963] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
17
|
Piccini G, Alessio M, Sauer J. Ab Initio Calculation of Rate Constants for Molecule-Surface Reactions with Chemical Accuracy. Angew Chem Int Ed Engl 2016. [DOI: 10.1002/ange.201601534] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- GiovanniMaria Piccini
- Institut für Chemie; Humboldt-Universität; Unter den Linden 6 10099 Berlin Germany
- Department of Chemistry and Applied Biosciences; ETH Zurich; c/o USI Campus; Via Giuseppe Buffi 13 6900 Lugano, Ticino Switzerland
| | - Maristella Alessio
- Institut für Chemie; Humboldt-Universität; Unter den Linden 6 10099 Berlin Germany
| | - Joachim Sauer
- Institut für Chemie; Humboldt-Universität; Unter den Linden 6 10099 Berlin Germany
| |
Collapse
|
18
|
Piccini G, Alessio M, Sauer J. Ab Initio Calculation of Rate Constants for Molecule-Surface Reactions with Chemical Accuracy. Angew Chem Int Ed Engl 2016; 55:5235-7. [PMID: 27008460 PMCID: PMC4834608 DOI: 10.1002/anie.201601534] [Citation(s) in RCA: 85] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Indexed: 11/10/2022]
Abstract
The ab initio prediction of reaction rate constants for systems with hundreds of atoms with an accuracy that is comparable to experiment is a challenge for computational quantum chemistry. We present a divide-and-conquer strategy that departs from the potential energy surfaces obtained by standard density functional theory with inclusion of dispersion. The energies of the reactant and transition structures are refined by wavefunction-type calculations for the reaction site. Thermal effects and entropies are calculated from vibrational partition functions, and the anharmonic frequencies are calculated separately for each vibrational mode. This method is applied to a key reaction of an industrially relevant catalytic process, the methylation of small alkenes over zeolites. The calculated reaction rate constants (free energies), pre-exponential factors (entropies), and enthalpy barriers show that our computational strategy yields results that agree with experiment within chemical accuracy limits (less than one order of magnitude).
Collapse
Affiliation(s)
- GiovanniMaria Piccini
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10099, Berlin, Germany.,Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, 6900, Lugano, Ticino, Switzerland
| | - Maristella Alessio
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10099, Berlin, Germany
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10099, Berlin, Germany.
| |
Collapse
|
19
|
Murugan V, Parasuraman P, Selvin JFA, Priyadarzini TRK, Gromiha MM, Fukui K, Veluraja K. Geometry Optimization of Carbohydrate Binding Sites of Influenza: A Quantum Mechanical Approach. J Carbohydr Chem 2015. [DOI: 10.1080/07328303.2015.1093135] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
20
|
Dhoke GV, Davari MD, Schwaneberg U, Bocola M. QM/MM Calculations Revealing the Resting and Catalytic States in Zinc-Dependent Medium-Chain Dehydrogenases/Reductases. ACS Catal 2015. [DOI: 10.1021/cs501524k] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Gaurao V. Dhoke
- RWTH Aachen University, Chair of Biotechnology, Worringer Weg 3, D-52074 Aachen, Germany
| | - Mehdi D. Davari
- RWTH Aachen University, Chair of Biotechnology, Worringer Weg 3, D-52074 Aachen, Germany
| | - Ulrich Schwaneberg
- RWTH Aachen University, Chair of Biotechnology, Worringer Weg 3, D-52074 Aachen, Germany
| | - Marco Bocola
- RWTH Aachen University, Chair of Biotechnology, Worringer Weg 3, D-52074 Aachen, Germany
| |
Collapse
|
21
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 301] [Impact Index Per Article: 33.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
22
|
Sokkar P, Boulanger E, Thiel W, Sanchez-Garcia E. Hybrid Quantum Mechanics/Molecular Mechanics/Coarse Grained Modeling: A Triple-Resolution Approach for Biomolecular Systems. J Chem Theory Comput 2015; 11:1809-18. [DOI: 10.1021/ct500956u] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Pandian Sokkar
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Eliot Boulanger
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Elsa Sanchez-Garcia
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
23
|
Werner HJ, Knizia G, Krause C, Schwilk M, Dornbach M. Scalable Electron Correlation Methods I.: PNO-LMP2 with Linear Scaling in the Molecular Size and Near-Inverse-Linear Scaling in the Number of Processors. J Chem Theory Comput 2015; 11:484-507. [DOI: 10.1021/ct500725e] [Citation(s) in RCA: 103] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Gerald Knizia
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Christine Krause
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Max Schwilk
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Mark Dornbach
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
24
|
Nam K. Acceleration of Ab Initio QM/MM Calculations under Periodic Boundary Conditions by Multiscale and Multiple Time Step Approaches. J Chem Theory Comput 2014; 10:4175-83. [PMID: 26588116 DOI: 10.1021/ct5005643] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Development of multiscale ab initio quantum mechanical and molecular mechanical (AI-QM/MM) method for periodic boundary molecular dynamics (MD) simulations and their acceleration by multiple time step approach are described. The developed method achieves accuracy and efficiency by integrating the AI-QM/MM level of theory and the previously developed semiempirical (SE) QM/MM-Ewald sum method [J. Chem. Theory Comput. 2005, 1, 2] extended to the smooth particle-mesh Ewald (PME) summation method. In the developed methods, the total energy of the simulated system is evaluated at the SE-QM/MM-PME level of theory to include long-range QM/MM electrostatic interactions, which is then corrected on the fly using the AI-QM/MM level of theory within the real space cutoff. The resulting energy expression enables decomposition of total forces applied to each atom into forces determined at the low-level SE-QM/MM method and correction forces at the AI-QM/MM level, to integrate the system using the reversible reference system propagator algorithm. The resulting method achieves a substantial speed-up of the entire calculation by minimizing the number of time-consuming energy and gradient evaluations at the AI-QM/MM level. Test calculations show that the developed multiple time step AI-QM/MM method yields MD trajectories and potential of mean force profiles comparable to single time step QM/MM results. The developed method, together with message passing interface (MPI) parallelization, accelerates the present AI-QM/MM MD simulations about 30-fold relative to the speed of single-core AI-QM/MM simulations for the molecular systems tested in the present work, making the method less than one order slower than the SE-QM/MM methods under periodic boundary conditions.
Collapse
Affiliation(s)
- Kwangho Nam
- Department of Chemistry and Computational Life Science Cluster (CLiC), Umeå University , 901 87, Umeå, Sweden
| |
Collapse
|
25
|
Thiel W. Computational Catalysis-Past, Present, and Future. Angew Chem Int Ed Engl 2014; 53:8605-13. [DOI: 10.1002/anie.201402118] [Citation(s) in RCA: 75] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2014] [Indexed: 01/03/2023]
|
26
|
|
27
|
Ferrario V, Ebert C, Svendsen A, Besenmatter W, Gardossi L. An integrated platform for automatic design and screening of virtual mutants based on 3D-QSAR analysis. ACTA ACUST UNITED AC 2014. [DOI: 10.1016/j.molcatb.2013.12.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
28
|
Quantum mechanical modeling: a tool for the understanding of enzyme reactions. Biomolecules 2013; 3:662-702. [PMID: 24970187 PMCID: PMC4030948 DOI: 10.3390/biom3030662] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2013] [Revised: 09/17/2013] [Accepted: 09/19/2013] [Indexed: 01/16/2023] Open
Abstract
Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.
Collapse
|
29
|
Hediger MR, De Vico L, Rannes JB, Jäckel C, Besenmatter W, Svendsen A, Jensen JH. In silico screening of 393 mutants facilitates enzyme engineering of amidase activity in CalB. PeerJ 2013; 1:e145. [PMID: 24010022 PMCID: PMC3757469 DOI: 10.7717/peerj.145] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2013] [Accepted: 08/06/2013] [Indexed: 11/25/2022] Open
Abstract
Our previously presented method for high throughput computational screening of mutant activity (Hediger et al., 2012) is benchmarked against experimentally measured amidase activity for 22 mutants of Candida antarctica lipase B (CalB). Using an appropriate cutoff criterion for the computed barriers, the qualitative activity of 15 out of 22 mutants is correctly predicted. The method identifies four of the six most active mutants with ≥3-fold wild type activity and seven out of the eight least active mutants with ≤0.5-fold wild type activity. The method is further used to screen all sterically possible (386) double-, triple- and quadruple-mutants constructed from the most active single mutants. Based on the benchmark test at least 20 new promising mutants are identified.
Collapse
Affiliation(s)
- Martin R Hediger
- Department of Chemistry, University of Copenhagen , Copenhagen , Denmark
| | | | | | | | | | | | | |
Collapse
|
30
|
Zhang J, Dolg M. Third-Order Incremental Dual-Basis Set Zero-Buffer Approach: An Accurate and Efficient Way To Obtain CCSD and CCSD(T) Energies. J Chem Theory Comput 2013; 9:2992-3003. [DOI: 10.1021/ct400284d] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jun Zhang
- Department of Chemistry, University of Cologne, Cologne, Germany
| | - Michael Dolg
- Department of Chemistry, University of Cologne, Cologne, Germany
| |
Collapse
|
31
|
Steinmann C, Fedorov DG, Jensen JH. Mapping enzymatic catalysis using the effective fragment molecular orbital method: towards all ab initio biochemistry. PLoS One 2013; 8:e60602. [PMID: 23593259 PMCID: PMC3625203 DOI: 10.1371/journal.pone.0060602] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Accepted: 02/28/2013] [Indexed: 12/02/2022] Open
Abstract
We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol−1 for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively.
Collapse
Affiliation(s)
- Casper Steinmann
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmitri G. Fedorov
- Nanosystems Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
- * E-mail:
| |
Collapse
|
32
|
Paasche A, Schirmeister T, Engels B. Benchmark Study for the Cysteine-Histidine Proton Transfer Reaction in a Protein Environment: Gas Phase, COSMO, QM/MM Approaches. J Chem Theory Comput 2013; 9:1765-77. [PMID: 26587634 DOI: 10.1021/ct301082y] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Proton transfer reactions are of crucial interest for the investigation of proteins. We have investigated the accuracy of commonly used quantum chemical methods for the description of proton transfer reactions in different environments (gas phase, COSMO, QM/MM) using the proton transfer between the catalytic dyad residues cysteine 145 and histidine 41 of SARS coronavirus main protease as a case study. The test includes thermodynamic, kinetic, and structural properties. The study comprises computationally demanding ab initio approaches (HF, CC2, MP2, SCS-CC2, SCS-MP2, CCSD(T)), popular density functional theories (BLYP, B3LYP, M06-2X), and semiempirical methods (MNDO/d, AM1, RM1, PM3, PM6). The approximated coupled cluster approach LCCSD(T) is taken as a reference method. We find that the robustness of the tested methods with respect to the environment correlates well with the level of theory. As an example HF, CC2, MP2, and their SCS variants show similar errors for gas phase, COSMO, or QM/MM computations. In contrast for semiempirical methods, the errors strongly diversify if one goes from gas phase to COSMO or QM/MM. Particular problems are observed for the recent semiempirical methods PM6 and RM1, which show the best performance for gas phase calculations but possess larger errors in conjunction with COSMO. Finally, a combination of SCS-MP2 and B3LYP or M06-2X allows reliable estimates about remaining errors.
Collapse
Affiliation(s)
- Alexander Paasche
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| | - Tanja Schirmeister
- Institut für Pharmazie und Biochemie, Johannes Gutenberg-Universität Mainz, Staudinger Weg 5, 55128 Mainz, Germany
| | - Bernd Engels
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
33
|
Hediger MR, De Vico L, Svendsen A, Besenmatter W, Jensen JH. A computational methodology to screen activities of enzyme variants. PLoS One 2012; 7:e49849. [PMID: 23284627 PMCID: PMC3524253 DOI: 10.1371/journal.pone.0049849] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2012] [Accepted: 10/14/2012] [Indexed: 11/22/2022] Open
Abstract
We present a fast computational method to efficiently screen enzyme activity. In the presented method, the effect of mutations on the barrier height of an enzyme-catalysed reaction can be computed within 24 hours on roughly 10 processors. The methodology is based on the PM6 and MOZYME methods as implemented in MOPAC2009, and is tested on the first step of the amide hydrolysis reaction catalyzed by the Candida Antarctica lipase B (CalB) enzyme. The barrier heights are estimated using adiabatic mapping and shown to give barrier heights to within 3 kcal/mol of B3LYP/6-31G(d)//RHF/3-21G results for a small model system. Relatively strict convergence criteria (0.5 kcal/(molÅ)), long NDDO cutoff distances within the MOZYME method (15 Å) and single point evaluations using conventional PM6 are needed for reliable results. The generation of mutant structures and subsequent setup of the semiempirical calculations are automated so that the effect on barrier heights can be estimated for hundreds of mutants in a matter of weeks using high performance computing.
Collapse
Affiliation(s)
- Martin R. Hediger
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | | | | | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
- * E-mail:
| |
Collapse
|
34
|
Kurashige Y, Yang J, Chan GKL, Manby FR. Optimization of orbital-specific virtuals in local Møller-Plesset perturbation theory. J Chem Phys 2012; 136:124106. [PMID: 22462834 DOI: 10.1063/1.3696962] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
We present an orbital-optimized version of our orbital-specific-virtuals second-order Møller-Plesset perturbation theory (OSV-MP2). The OSV model is a local correlation ansatz with a small basis of virtual functions for each occupied orbital. It is related to the Pulay-Saebø approach, in which domains of virtual orbitals are drawn from a single set of projected atomic orbitals; but here the virtual functions associated with a particular occupied orbital are specifically tailored to the correlation effects in which that orbital participates. In this study, the shapes of the OSVs are optimized simultaneously with the OSV-MP2 amplitudes by minimizing the Hylleraas functional or approximations to it. It is found that optimized OSVs are considerably more accurate than the OSVs obtained through singular value decomposition of diagonal blocks of MP2 amplitudes, as used in our earlier work. Orbital-optimized OSV-MP2 recovers smooth potential energy surfaces regardless of the number of virtuals. Full optimization is still computationally demanding, but orbital optimization in a diagonal or Kapuy-type MP2 approximation provides an attractive scheme for determining accurate OSVs.
Collapse
Affiliation(s)
- Yuki Kurashige
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi, 444-8585, Japan.
| | | | | | | |
Collapse
|
35
|
Cho KB, Kim EJ, Seo MS, Shaik S, Nam W. Correlating DFT-calculated energy barriers to experiments in nonheme octahedral Fe(IV)O species. Chemistry 2012; 18:10444-53. [PMID: 22714711 DOI: 10.1002/chem.201200096] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2012] [Revised: 04/20/2012] [Indexed: 11/08/2022]
Abstract
The experimentally measured bimolecular reaction rate constant, k(2), should in principle correlate with the theoretically calculated rate-limiting free energy barrier, ΔG(≠), through the Eyring equation, but it fails quite often to do so due to the inability of current computational methods to account in a precise manner for all the factors contributing to ΔG(≠). This is further aggravated by the exponential sensitivity of the Eyring equation to these factors. We have taken herein a pragmatic approach for C-H activation reactions of 1,4-cyclohexadiene with a variety of octahedral nonheme Fe(IV)O complexes. The approach consists of empirically determining two constants that would aid in predicting experimental k(2) values uniformly from theoretically calculated electronic energy (ΔE(≠)) values. Shown in this study is the predictive power as well as insights into energy relationships in Fe(IV)O C-H activation reactions. We also find that the difference between ΔG(≠) and ΔE(≠) converges at slow reactions, in a manner suggestive of changes in the importance of the triplet spin state weight in the overall reaction.
Collapse
Affiliation(s)
- Kyung-Bin Cho
- Department of Bioinspired Science, Ewha Womans University, 120-750 Seoul, Korea
| | | | | | | | | |
Collapse
|
36
|
Yang J, Chan GKL, Manby FR, Schütz M, Werner HJ. The orbital-specific-virtual local coupled cluster singles and doubles method. J Chem Phys 2012; 136:144105. [DOI: 10.1063/1.3696963] [Citation(s) in RCA: 150] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
|
37
|
Wong KY, Gu H, Zhang S, Piccirilli JA, Harris ME, York DM. Characterization of the reaction path and transition states for RNA transphosphorylation models from theory and experiment. Angew Chem Int Ed Engl 2011; 51:647-51. [PMID: 22076983 DOI: 10.1002/anie.201104147] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2011] [Revised: 10/06/2011] [Indexed: 12/31/2022]
Affiliation(s)
- Kin-Yiu Wong
- BioMaPS Institute for Quantitative Biology, Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 610 Taylor Road, Piscataway, NJ 08854, USA.
| | | | | | | | | | | |
Collapse
|
38
|
Wong KY, Gu H, Zhang S, Piccirilli JA, Harris ME, York DM. Characterization of the Reaction Path and Transition States for RNA Transphosphorylation Models from Theory and Experiment. Angew Chem Int Ed Engl 2011. [DOI: 10.1002/ange.201104147] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
39
|
Werner HJ, Schütz M. An efficient local coupled cluster method for accurate thermochemistry of large systems. J Chem Phys 2011; 135:144116. [DOI: 10.1063/1.3641642] [Citation(s) in RCA: 232] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Adler TB, Werner HJ. An explicitly correlated local coupled cluster method for calculations of large molecules close to the basis set limit. J Chem Phys 2011; 135:144117. [DOI: 10.1063/1.3647565] [Citation(s) in RCA: 106] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
|
41
|
Werner HJ, Knowles PJ, Knizia G, Manby FR, Schütz M. Molpro: a general-purpose quantum chemistry program package. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.82] [Citation(s) in RCA: 2174] [Impact Index Per Article: 167.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
42
|
Local Approximations for an Efficient and Accurate Treatment of Electron Correlation and Electron Excitations in Molecules. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2011. [DOI: 10.1007/978-90-481-2853-2_14] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
43
|
Efficient Explicitly Correlated Coupled-Cluster Approximations. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2010. [DOI: 10.1007/978-90-481-2885-3_21] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
44
|
Mata RA. Application of high level wavefunction methods in quantum mechanics/molecular mechanics hybrid schemes. Phys Chem Chem Phys 2010; 12:5041-52. [DOI: 10.1039/b918608e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
45
|
Sumowski CV, Ochsenfeld C. A Convergence Study of QM/MM Isomerization Energies with the Selected Size of the QM Region for Peptidic Systems. J Phys Chem A 2009; 113:11734-41. [DOI: 10.1021/jp902876n] [Citation(s) in RCA: 89] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Chris Vanessa Sumowski
- Theoretische Chemie, Universität Tübingen, Auf der Morgenstelle 8, D-72076 Tübingen, Germany
| | - Christian Ochsenfeld
- Theoretische Chemie, Universität Tübingen, Auf der Morgenstelle 8, D-72076 Tübingen, Germany
| |
Collapse
|
46
|
Adler TB, Werner HJ. Local explicitly correlated coupled-cluster methods: Efficient removal of the basis set incompleteness and domain errors. J Chem Phys 2009; 130:241101. [DOI: 10.1063/1.3160675] [Citation(s) in RCA: 140] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
47
|
Abstract
Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.
Collapse
Affiliation(s)
- Hans Martin Senn
- Department of Chemistry, WestCHEM and University of Glasgow, Glasgow G12 8QQ, UK.
| | | |
Collapse
|
48
|
Adler TB, Werner HJ, Manby FR. Local explicitly correlated second-order perturbation theory for the accurate treatment of large molecules. J Chem Phys 2009; 130:054106. [DOI: 10.1063/1.3040174] [Citation(s) in RCA: 151] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
|
50
|
Marsman M, Grüneis A, Paier J, Kresse G. Second-order Mo̸ller–Plesset perturbation theory applied to extended systems. I. Within the projector-augmented-wave formalism using a plane wave basis set. J Chem Phys 2009; 130:184103. [DOI: 10.1063/1.3126249] [Citation(s) in RCA: 175] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|