1
|
Wang J, Ren T, Sun G, Zhang N, Zhao L, Zhong R. Mechanism of AGT-Mediated Repair of dG-dC Cross-Links in the Drug Resistance to Chloroethylnitrosoureas: Molecular Docking, MD Simulation, and ONIOM (QM/MM) Investigation. J Chem Inf Model 2024; 64:3411-3429. [PMID: 38511939 DOI: 10.1021/acs.jcim.3c01958] [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: 03/22/2024]
Abstract
Chloroethylnitrosoureas (CENUs) are important chemotherapies applied in the treatment of cancer. They exert anticancer activity by inducing DNA interstrand cross-links (ICLs) via the formation of two O6-alkylguanine intermediates, O6-chloroethylguanine (O6-ClEtG) and N1,O6-ethanoguanine (N1,O6-EtG). However, O6-alkylguanine-DNA alkyltransferase (AGT), a DNA-repair enzyme, can restore the O6-alkylguanine damages and thereby obstruct the formation of ICLs (dG-dC cross-link). In this study, the inhibitory mechanism of ICL formation was investigated to elucidate the drug resistance of CENUs mediated by AGT in detail. Based on the structures of the substrate-enzyme complexes obtained from docking and MD simulations, two ONIOM (QM/MM) models with different sizes of the QM region were constructed. The model with a larger QM region, which included the substrate (O6-ClEtG or N1,O6-EtG), a water molecule, and five residues (Tyr114, Cys145, His146, Lys165, and Glu172) in the active pocket of AGT, accurately described the repairing reaction and generated the results coinciding with the experimental outcomes. The repair process consists of two sequential steps: hydrogen transfer to form a thiolate anion on Cys145 and alkyl transfer from the O6 site of guanine (the rate-limiting step). The repair of N1,O6-EtG was more favorable than that of O6-ClEtG from both kinetics and thermodynamics aspects. Moreover, the comparison of the repairing process with the formation of dG-dC cross-link and the inhibition of AGT by O6-benzylguanine (O6-BG) showed that the presence of AGT could effectively interrupt the formation of ICLs leading to drug resistance, and the inhibition of AGT by O6-BG that was energetically more favorable than the repair of O6-ClEtG could not prevent the repair of N1,O6-EtG. Therefore, it is necessary to completely eliminate AGT activity before CENUs medication to enhance the chemotherapeutic effectiveness. This work provides reasonable explanations for the supposed mechanism of AGT-mediated drug resistance of CENUs and will assist in the development of novel CENU chemotherapies and their medication strategies.
Collapse
Affiliation(s)
- Jiaojiao Wang
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| | - Ting Ren
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| | - Guohui Sun
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| | - Na Zhang
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| | - Lijiao Zhao
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| | - Rugang Zhong
- Beijing Key Laboratory of Environmental and Viral Oncology, College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
2
|
Li X, Duan R, Wang Y, Qu LB, Li Z, Wei D. Insights into N-Heterocyclic Carbene-Catalyzed Oxidative α-C(sp 3)-H Activation of Aliphatic Aldehydes and Cascade [2 + 2] Cycloaddition with Ketimines. J Org Chem 2019; 84:6117-6125. [PMID: 31012584 DOI: 10.1021/acs.joc.9b00295] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Predicting the chemoselectivity of [2 + 2] cyclizations is an important challenge in organic chemistry. Herein, we provided a valuable case for this issue. Density functional theory calculations were performed to systematically study the possible mechanisms and origin of selectivities for the N-heterocyclic carbene (NHC)-catalyzed oxidative α-C(sp3)-H activation of aliphatic aldehydes and the cascade [2 + 2] cycloaddition with ketimines. The [2 + 2] cycloaddition of azolium enolate intermediates to the C═N bond, rather than the C═O bond of ketimine, is revealed to be determined by chemo- and stereoselectivity. By comparing the energy gap between the frontier molecular orbitals (FMOs) of the two reacting parts involved in the [2 + 2] cycloaddition transition states, we propose a new strategy to determine the origin of the reaction chemoselectivity. Moreover, the local nucleophilic index can efficiently predict the active site of ketimines. Further analyses illustrate that NHC can increase the nucleophilicity of aldehydes and the acidity of the α-C(sp3)-H bond, and 3,3',5,5'-tetra- tert-butyl diphenoquinone (DQ) acts as an oxidant and promotes α-C(sp3)-H bond deprotonation. This work is useful not only for understanding the NHC-catalyzed oxidative [2 + 2] annulation but also for developing new applications of the FMO theory in organocatalytic cyclizations.
Collapse
Affiliation(s)
- Xue Li
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| | - Ruihong Duan
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| | - Yanyan Wang
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| | - Ling-Bo Qu
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| | - Zhongjun Li
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| | - Donghui Wei
- The College of Chemistry and Molecular Engineering , Zhengzhou University , 100 Science Avenue , Zhengzhou , Henan 450001 , P. R. China
| |
Collapse
|
3
|
Mamun AA, Toda MJ, Kozlowski PM. Can photolysis of the Co C bond in coenzyme B12-dependent enzymes be used to mimic the native reaction? JOURNAL OF PHOTOCHEMISTRY AND PHOTOBIOLOGY B-BIOLOGY 2019; 191:175-184. [DOI: 10.1016/j.jphotobiol.2018.12.018] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Revised: 12/19/2018] [Accepted: 12/20/2018] [Indexed: 12/22/2022]
|
4
|
Qiu S, Azofra LM, MacFarlane DR, Sun C. Hydrogen bonding effect between active site and protein environment on catalysis performance in H 2-producing [NiFe] hydrogenases. Phys Chem Chem Phys 2018; 20:6735-6743. [PMID: 29457815 DOI: 10.1039/c7cp07685a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H+/e- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S- in cysteine residues shows a larger influence on the energy profile than that via the CN- ligands on Fe.
Collapse
Affiliation(s)
- Siyao Qiu
- School of Chemistry, Faculty of Science, Monash University, Clayton, VIC 3800, Australia.
| | | | | | | |
Collapse
|
5
|
Wick CR, Smith DM. Modeling the Reactions Catalyzed by Coenzyme B 12 Dependent Enzymes: Accuracy and Cost-Quality Balance. J Phys Chem A 2018; 122:1747-1755. [PMID: 29389127 DOI: 10.1021/acs.jpca.7b11798] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The reactions catalyzed by coenzyme B12 dependent enzymes are formally initiated by the homolytic cleavage of a carbon-cobalt bond and a subsequent or concerted H-atom-transfer reaction. A reasonable model chemistry for describing those reactions should, therefore, account for an accurate description of both reactions. The inherent limitation due to the necessary system size renders the coenzyme B12 system a suitable candidate for DFT or hybrid QM/MM methods; however, the accurate description of both homolytic Co-C cleavage and H-atom-transfer reactions within this framework is challenging and can lead to controversial results with varying accuracy. We present an assessment study of 16 common density functionals applied to prototypical model systems for both reactions. H-abstraction reactions were modeled on the basis of four reference reactions designed to resemble a broad range of coenzyme B12 reactions. The Co-C cleavage reaction is treated by an ONIOM(QM/MM) setup that is in excellent agreement with solution-phase experimental data and is as accurate as full DFT calculations on the complete model system. We find that the meta-GGAs TPSS-D3 and M06L-D3 and the meta-hybrid M06-D3 give the best overall performance with MUEs for both types of reactions below 10 kJ mol-1. Our recommended model chemistry allows for a fast and accurate description of coenzyme B12 chemistry that is readily applicable to study the reactions in an enzymatic framework.
Collapse
Affiliation(s)
- Christian R Wick
- Division of Physical Chemistry, Group for Computational Life Sciences, Ruđer Bošković Institute , Bijenička cesta 54, 10000 Zagreb, Croatia
| | - David M Smith
- Division of Physical Chemistry, Group for Computational Life Sciences, Ruđer Bošković Institute , Bijenička cesta 54, 10000 Zagreb, Croatia.,Center for Computational Chemistry, FAU Erlangen-Nürnberg , Nägelsbachstrasse 25, 91052 Erlangen, Germany
| |
Collapse
|
6
|
Yu MJ, Chen SL. From NAD + to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase. Chemistry 2017; 23:7545-7557. [PMID: 28374531 DOI: 10.1002/chem.201700405] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2017] [Indexed: 02/02/2023]
Abstract
Lactate racemase (LarA), a new nickel enzyme discovered recently, catalyzes the racemization between d- and l-lactates with a novel nickel pincer cofactor (Ni-PTTMN) derived from nicotinic acid. In this study, by using DFT and a 200-atom active-site model, LarA is revealed to employ a modified proton-coupled hydride-transfer mechanism in which a hydride is transferred to a cofactor pyridine carbon from the substrate α-carbon along with proton transfer from the substrate hydroxy group to a histidine, and then moved back from the opposite side. Tyr294 and Lys298 provide significant acceleration effects by orientating substrates and stabilizing the negative charge developing at the substrate hydroxy oxygen. The barrier was determined to be 12.0 kcal mol-1 , which reveals enhanced racemase activity relative to the LarA reaction using NAD+ -like cofactors. Compared with NAD+ , Ni-PTTMN has a stronger hydride-addition reactivity in moderate and high environmental polarity and may fit perfectly the moderately polar active site of LarA.
Collapse
Affiliation(s)
- Ming-Jia Yu
- Key Laboratory of Cluster Science of the Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, 100081, P.R. China
| | - Shi-Lu Chen
- Key Laboratory of Cluster Science of the Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, 100081, P.R. China
| |
Collapse
|
7
|
Liao RZ, Chen SL, Siegbahn PEM. Unraveling the Mechanism and Regioselectivity of the B12-Dependent Reductive Dehalogenase PceA. Chemistry 2016; 22:12391-9. [DOI: 10.1002/chem.201601575] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2016] [Indexed: 01/09/2023]
Affiliation(s)
- Rong-Zhen Liao
- Key Laboratory of Material Chemistry for Energy Conversion and Storage; Ministry of Education; Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica; Hubei Key Laboratory of Materials Chemistry and Service Failure; School of Chemistry and Chemical Engineering; Huazhong University of Science and Technology; Wuhan 430074 P. R. China
| | - Shi-Lu Chen
- School of Chemistry; Beijing Institute of Technology; Beijing 100081 P. R. China
| | - Per E. M. Siegbahn
- Department of Organic Chemistry; Arrhenius Laboratory; Stockholm University; 10691 Stockholm Sweden
| |
Collapse
|
8
|
Spataru T, Fernandez F. The Nature of the Co-C Bond Cleavage Processes in Methylcob(II)Alamin and Adenosylcob(III)Alamin. CHEMISTRY JOURNAL OF MOLDOVA 2016. [DOI: 10.19261/cjm.2016.11(1).01] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
|
9
|
Mao X, Sun Y, Pei S. A theoretical investigation into thiophenic derivative cracking mechanism over acidic and cation-exchanged beta zeolites. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.10.019] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
10
|
Liao RZ, Chen SL, Siegbahn PEM. Which Oxidation State Initiates Dehalogenation in the B12-Dependent Enzyme NpRdhA: CoII, CoI, or Co0? ACS Catal 2015. [DOI: 10.1021/acscatal.5b01502] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Rong-Zhen Liao
- Key
Laboratory of Material Chemistry for Energy Conversion and Storage,
Ministry of Education, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, People’s Republic of China
| | - Shi-Lu Chen
- School
of Chemistry, Beijing Institute of Technology, Beijing 100081, People’s Republic of China
| | - Per E. M. Siegbahn
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
11
|
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: 296] [Impact Index Per Article: 32.9] [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
|
12
|
Chung LW, Sameera WMC, Ramozzi R, Page AJ, Hatanaka M, Petrova GP, Harris TV, Li X, Ke Z, Liu F, Li HB, Ding L, Morokuma K. The ONIOM Method and Its Applications. Chem Rev 2015; 115:5678-796. [PMID: 25853797 DOI: 10.1021/cr5004419] [Citation(s) in RCA: 734] [Impact Index Per Article: 81.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Lung Wa Chung
- †Department of Chemistry, South University of Science and Technology of China, Shenzhen 518055, China
| | - W M C Sameera
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Romain Ramozzi
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Alister J Page
- §Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, Australia
| | - Miho Hatanaka
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Galina P Petrova
- ∥Faculty of Chemistry and Pharmacy, University of Sofia, Bulgaria Boulevard James Bourchier 1, 1164 Sofia, Bulgaria
| | - Travis V Harris
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan.,⊥Department of Chemistry, State University of New York at Oswego, Oswego, New York 13126, United States
| | - Xin Li
- #State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Zhuofeng Ke
- ∇School of Chemistry and Chemical Engineering, Sun Yat-sen University, Guangzhou 510275, China
| | - Fengyi Liu
- ○Key Laboratory of Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, China
| | - Hai-Bei Li
- ■School of Ocean, Shandong University, Weihai 264209, China
| | - Lina Ding
- ▲School of Pharmaceutical Sciences, Zhengzhou University, 100 Kexue Avenue, Zhengzhou, Henan 450001, China
| | - Keiji Morokuma
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| |
Collapse
|
13
|
Conrad KS, Jordan CD, Brown KL, Brunold TC. Spectroscopic and Computational Studies of Cobalamin Species with Variable Lower Axial Ligation: Implications for the Mechanism of Co–C Bond Activation by Class I Cobalamin-Dependent Isomerases. Inorg Chem 2015; 54:3736-47. [DOI: 10.1021/ic502665x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Karen S. Conrad
- Department
of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Christopher D. Jordan
- Department
of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Kenneth L. Brown
- Department of Chemistry and Biochemistry, Ohio University, Athens, Ohio 45701, United States
| | - Thomas C. Brunold
- Department
of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
14
|
The entropic contributions in vitamin B12 enzymes still reflect the electrostatic paradigm. Proc Natl Acad Sci U S A 2015; 112:4328-33. [PMID: 25805820 DOI: 10.1073/pnas.1503828112] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
The catalytic power of enzymes containing coenzyme B12 has been, in some respects, the "last bastion" for the strain hypothesis. Our previous study of this system established by a careful sampling that the major part of the catalytic effect is due to the electrostatic interaction between the ribose of the ado group and the protein and that the strain contribution is very small. This finding has not been sufficiently appreciated due to misunderstandings of the power of the empirical valence bond (EVB) calculations and the need of sufficient sampling. Furthermore, some interesting new experiments point toward entropic effects as the source of the catalytic power, casting doubt on the validity of the electrostatic idea, at least, in the case of B12 enzymes. Here, we focus on the observation of the entropic effects and on analyzing their origin. We clarify that our EVB approach evaluates free energies rather than enthalpies and demonstrate by using the restraint release (RR) approach that the observed entropic contribution to the activation barrier is of electrostatic origin. Our study illustrates the power of the RR approach by evaluating the entropic contributions to catalysis and provides further support to our paradigm for the origin of the catalytic power of B12 enzymes. Overall, our study provides major support to our electrostatic preorganization idea and also highlights the basic requirements from ab initio quantum mechanics/molecular mechanics calculations of activation free energies of enzymatic reactions.
Collapse
|
15
|
Using potential energy surface scans to examine the bond dissociation energies of trans-ReOS2N2 and [ReOS3N]1− model complexes. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.08.017] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
16
|
Kepp KP. Co-C dissociation of adenosylcobalamin (coenzyme B12): role of dispersion, induction effects, solvent polarity, and relativistic and thermal corrections. J Phys Chem A 2014; 118:7104-17. [PMID: 25116644 DOI: 10.1021/jp503607k] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantum-chemical cluster modeling is challenged in the limit of large, soft systems by the effects of dispersion and solvent, and well as other physical interactions. Adenosylcobalamin (AdoCbl, coenzyme B12), as one of the most complex cofactors in life, constitutes such a challenge. The cleavage of its unique organometallic Co-C bond has inspired multiple studies of this cofactor. This paper reports the fully relaxed potential energy surface of Co-C cleavage of AdoCbl, including for the first time all side-chain interactions with the dissociating Ado group. Various methods and corrections for dispersion, relativistic effects, solvent polarity, basis set superposition error, and thermal and vibrational effects were investigated, totaling more than 550 single-point energies for the large model. The results show immense variability depending on method, including solvation, functional type, and dispersion, challenging the conceived accuracy of methods used for such systems. In particular, B3LYP-D3 seems to severely underestimate the Co-C bond strength, consistent with previous results, and BP86 remains accurate for cobalamins when dispersion interactions are accounted for.
Collapse
Affiliation(s)
- Kasper P Kepp
- DTU Chemistry, Technical University of Denmark , Building 206, Kgs. Lyngby, DK-2800, Denmark
| |
Collapse
|
17
|
Pang J, Scrutton NS, Sutcliffe MJ. Quantum Mechanics/Molecular Mechanics Studies on the Mechanism of Action of Cofactor Pyridoxal 5′-Phosphate in Ornithine 4,5-Aminomutase. Chemistry 2014; 20:11390-401. [DOI: 10.1002/chem.201402759] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2014] [Indexed: 02/02/2023]
|
18
|
Blomberg MRA, Borowski T, Himo F, Liao RZ, Siegbahn PEM. Quantum chemical studies of mechanisms for metalloenzymes. Chem Rev 2014; 114:3601-58. [PMID: 24410477 DOI: 10.1021/cr400388t] [Citation(s) in RCA: 431] [Impact Index Per Article: 43.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Margareta R A Blomberg
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University , SE-106 91 Stockholm, Sweden
| | | | | | | | | |
Collapse
|
19
|
Mechanistic investigation of methanol to propene conversion catalyzed by H-beta zeolite: a two-layer ONIOM study. J Mol Model 2013; 19:5407-22. [DOI: 10.1007/s00894-013-2030-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2013] [Accepted: 10/07/2013] [Indexed: 10/26/2022]
|
20
|
Kellie JL, Wetmore SD. Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases. CAN J CHEM 2013. [DOI: 10.1139/cjc-2012-0506] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.
Collapse
Affiliation(s)
- Jennifer L. Kellie
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| |
Collapse
|
21
|
DFT Studies of Trans and Cis Influences in the Homolysis of the Co–C Bond in Models of the Alkylcobalamins. J Phys Chem A 2013; 117:3057-68. [DOI: 10.1021/jp311788t] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
22
|
Li Z, Jiang YY, Yeagley AA, Bour JP, Liu L, Chruma JJ, Fu Y. Mechanism of the Pd-catalyzed decarboxylative allylation of α-imino esters: decarboxylation via free carboxylate ion. Chemistry 2012; 18:14527-38. [PMID: 23002017 DOI: 10.1002/chem.201201425] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2012] [Revised: 07/09/2012] [Indexed: 11/09/2022]
Abstract
The Pd-catalyzed decarboxylative allylation of α-(diphenylmethylene)imino esters (1) or allyl diphenylglycinate imines (2) is an efficient method to construct new C(sp(3))-C(sp(3)) bonds. The detailed mechanism of this reaction was studied by theoretical calculations [ONIOM(B3LYP/LANL2DZ+p:PM6)] combined with experimental observations. The overall catalytic cycle was found to consist of three steps: oxidative addition, decarboxylation, and reductive allylation. The oxidative addition of 1 to [(dba)Pd(PPh(3))(2)] (dba = dibenzylideneacetone) produces an allylpalladium cation and a carboxylate anion with a low activation barrier of +9.1 kcal mol(-1). The following rate-determining decarboxylation proceeds via a solvent-exposed α-imino carboxylate anion rather than an O-ligated allylpalladium carboxylate with an activation barrier of +22.7 kcal mol(-1). The 2-azaallyl anion generated by this decarboxylation attacks the face of the allyl ligand opposite to the Pd center in an outer-sphere process to produce major product 3, with a lower activation barrier than that of the minor product 4. A positive linear Hammett correlation [ρ = 1.10 for the PPh(3) ligand] with the observed regioselectivity (3 versus 4) supports an outer-sphere pathway for the allylation step. When Pd combined with the bis(diphenylphosphino)butane (dppb) ligand is employed as a catalyst, the decarboxylation still proceeds via the free carboxylate anion without direct assistance of the cationic Pd center. Consistent with experimental observations, electron-withdrawing substituents on 2 were calculated to have lower activation barriers for decarboxylation and, thus, accelerate the overall reaction rates.
Collapse
Affiliation(s)
- Zhe Li
- Department of Applied Chemistry, China Agricultural University, Beijing 100193, PR China
| | | | | | | | | | | | | |
Collapse
|
23
|
Yoshida T, Hitaoka S, Mashima A, Sugimoto T, Matoba H, Chuman H. Combined QM/MM (ONIOM) and QSAR approach to the study of complex formation of matrix metalloproteinase‑9 with a series of biphenylsulfonamides−LERE-QSAR analysis (V). J Phys Chem B 2012; 116:10283-9. [PMID: 22845734 DOI: 10.1021/jp305476x] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We previously proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as an ab initio fragment molecular orbital ones. In the present work, we applied LERE-QSAR to complex formation of matrix metalloproteinase-9 (MMP-9) with a series of substituted biphenylsulfonamides. The results shows that the overall free-energy change accompanying complex formation is due to predominantly the contribution from the electrostatic interaction with the zinc atom in the active site of MMP-9. Carbonic anhydrase (CA) belongs to the zinc-containing protease family. In contrast to the current case of MMP-9, the overall free-energy change during complex formation of CA with a series of benzenesulfonamides is due to the contributions from the solvation and dissociation free-energy changes, as previously reported. Comparison of the two sets of results indicates quantitative differences in the relative contributions of free-energy components to the overall free-energy change between the two data sets, corresponding with those in the respective classical QSAR equations. The LERE-QSAR procedure was demonstrated to quantitatively reveal differences in the binding mechanisms between the two cases involving similar but different zinc-containing proteins at the electronic and atomic levels.
Collapse
Affiliation(s)
- Tatsusada Yoshida
- Institute of Health Biosciences, The University of Tokushima Graduate School , 1-78 Shomachi, Tokushima 770-8505, Japan
| | | | | | | | | | | |
Collapse
|
24
|
Kozlowski PM, Kumar M, Piecuch P, Li W, Bauman NP, Hansen JA, Lodowski P, Jaworska M. The Cobalt–Methyl Bond Dissociation in Methylcobalamin: New Benchmark Analysis Based on Density Functional Theory and Completely Renormalized Coupled-Cluster Calculations. J Chem Theory Comput 2012; 8:1870-94. [DOI: 10.1021/ct300170y] [Citation(s) in RCA: 83] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pawel M. Kozlowski
- Department of Chemistry, University
of Louisville,
2320 South Brook St., Louisville, Kentucky 40292, United States
| | - Manoj Kumar
- Department of Chemistry, University
of Louisville,
2320 South Brook St., Louisville, Kentucky 40292, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University,
578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Wei Li
- Department of Chemistry, Michigan State University,
578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Nicholas P. Bauman
- Department of Chemistry, Michigan State University,
578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Jared A. Hansen
- Department of Chemistry, Michigan State University,
578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Piotr Lodowski
- Institute
of Chemistry, University of Silesia, Szkolna
9, PL-40 006 Katowice, Poland
| | - Maria Jaworska
- Institute
of Chemistry, University of Silesia, Szkolna
9, PL-40 006 Katowice, Poland
| |
Collapse
|
25
|
Adenosylcobalamin enzymes: theory and experiment begin to converge. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2012; 1824:1154-64. [PMID: 22516318 DOI: 10.1016/j.bbapap.2012.03.012] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/01/2011] [Revised: 02/04/2012] [Accepted: 03/27/2012] [Indexed: 11/21/2022]
Abstract
Adenosylcobalamin (coenzyme B(12)) serves as the cofactor for a group of enzymes that catalyze unusual rearrangement or elimination reactions. The role of the cofactor as the initiator of reactive free radicals needed for these reactions is well established. Less clear is how these enzymes activate the coenzyme towards homolysis and control the radicals once generated. The availability of high resolution X-ray structures combined with detailed kinetic and spectroscopic analyses have allowed several adenosylcobalamin enzymes to be computationally modeled in some detail. Computer simulations have generally obtained good agreement with experimental data and provided valuable insight into the mechanisms of these unusual reactions. Importantly, atomistic modeling of the enzymes has allowed the role of specific interactions between protein, substrate and coenzyme to be explored, leading to mechanistic predictions that can now be tested experimentally. This article is part of a Special Issue entitled: Radical SAM enzymes and Radical Enzymology.
Collapse
|
26
|
Pang J, Li X, Morokuma K, Scrutton NS, Sutcliffe MJ. Large-Scale Domain Conformational Change Is Coupled to the Activation of the Co–C Bond in the B12-Dependent Enzyme Ornithine 4,5-Aminomutase: A Computational Study. J Am Chem Soc 2012; 134:2367-77. [DOI: 10.1021/ja210417k] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
| | - Xin Li
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | | | | |
Collapse
|
27
|
Bucher D, Sandala GM, Durbeej B, Radom L, Smith DM. The Elusive 5′-Deoxyadenosyl Radical in Coenzyme-B12-Mediated Reactions. J Am Chem Soc 2012; 134:1591-9. [DOI: 10.1021/ja207809b] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Denis Bucher
- School of Chemistry and ARC Centre of Excellence
for Free Radical Chemistry
and Biotechnology, University of Sydney, Sydney, NSW 2006, Australia
| | - Gregory M. Sandala
- School of Chemistry and ARC Centre of Excellence
for Free Radical Chemistry
and Biotechnology, University of Sydney, Sydney, NSW 2006, Australia
- Division of Organic
Chemistry and Biochemistry, Ruđer Bošković Institute, 10002 Zagreb, Croatia
| | - Bo Durbeej
- Division of Computational
Physics, IFM Theory and Modelling, Linköping University, SE-581 83 Linköping, Sweden
| | - Leo Radom
- School of Chemistry and ARC Centre of Excellence
for Free Radical Chemistry
and Biotechnology, University of Sydney, Sydney, NSW 2006, Australia
| | - David M. Smith
- Division of Organic
Chemistry and Biochemistry, Ruđer Bošković Institute, 10002 Zagreb, Croatia
- Computer-Chemie-Centrum, University of Erlangen-Nürnberg, 91052 Erlangen, Germany
| |
Collapse
|
28
|
Manzerova J, Krymov V, Gerfen GJ. Investigating the intermediates in the reaction of ribonucleoside triphosphate reductase from Lactobacillus leichmannii: An application of HF EPR-RFQ technology. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2011; 213:32-45. [PMID: 21944735 DOI: 10.1016/j.jmr.2011.08.030] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/18/2011] [Accepted: 08/23/2011] [Indexed: 05/31/2023]
Abstract
In this investigation high-frequency electron paramagnetic resonance spectroscopy (HFEPR) in conjunction with innovative rapid freeze-quench (RFQ) technology is employed to study the exchange-coupled thiyl radical-cob(II)alamin system in ribonucleotide reductase from a prokaryote Lactobacillus leichmannii. The size of the exchange coupling (Jex) and the values of the thiyl radical g tensor are refined, while confirming the previously determined (Gerfen et al. (1996) [20]) distance between the paramagnets. Conclusions relevant to ribonucleotide reductase catalysis and the architecture of the active site are presented. A key part of this work has been the development of a unique RFQ apparatus for the preparation of millisecond quench time RFQ samples which can be packed into small (0.5 mm ID) sample tubes used for CW and pulsed HFEPR--lack of this ability has heretofore precluded such studies. The technology is compatible with a broad range of spectroscopic techniques and can be readily adopted by other laboratories.
Collapse
Affiliation(s)
- Julia Manzerova
- Department of Physiology and Biophysics, Albert Einstein College of Medicine, 1300 Morris Park Ave., Bronx, NY 10461, United States
| | | | | |
Collapse
|
29
|
Rutledge LR, Wetmore SD. Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines. J Am Chem Soc 2011; 133:16258-69. [PMID: 21877721 DOI: 10.1021/ja207181c] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(‡) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.
Collapse
Affiliation(s)
- Lesley R Rutledge
- Department of Chemistry and Biochemistry, University of Lethbridge, Alberta T1K 3M4, Canada
| | | |
Collapse
|
30
|
Hirao H, Morokuma K. [Recent progress in the theoretical studies of structure, function, and reaction of biological molecules]. YAKUGAKU ZASSHI 2011; 131:1151-61. [PMID: 21804318 DOI: 10.1248/yakushi.131.1151] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Essential biomolecular functions often involve electron-related events such as chemical reactions and photoluminescence phenomena. Theoretical description of such electronic processes requires the use of quantum mechanics (QM), but the number of atoms that can be handled with QM is usually smaller than the number of atoms present in a single protein. A reasonable strategy is therefore to give priority to a few tens or hundreds of atoms in the system and deal with them quantum mechanically. Lower-priority atoms influence the event occurring in the higher-priority area; therefore, their effect should also be taken into account. Under these circumstances, a reasonable approach is to apply two or more different theoretical methods to differently prioritized subsystems. QM can be combined, for example, with less accurate yet much less demanding molecular mechanics (MM). Our own N-layered integrated molecular orbital and molecular mechanics (ONIOM) method allows for such hybrid calculations, and our group has been applying it to a wide range of biology-related problems. In this paper, we briefly explain the theoretical background and the procedure for the theoretical investigation of biological systems. Subsequently, we provide an overview of some of our recent studies of metalloenzymes and photobiology-related problems.
Collapse
Affiliation(s)
- Hajime Hirao
- Fukui Institute for Fundamental Chemistry, Kyoto University, Japan
| | | |
Collapse
|
31
|
Chung LW, Hirao H, Li X, Morokuma K. The ONIOM method: its foundation and applications to metalloenzymes and photobiology. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.85] [Citation(s) in RCA: 132] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
32
|
Rommel JB, Kästner J. The fragmentation-recombination mechanism of the enzyme glutamate mutase studied by QM/MM simulations. J Am Chem Soc 2011; 133:10195-203. [PMID: 21612278 DOI: 10.1021/ja202312d] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The radical mechanism of the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is studied with quantum mechanical/molecular mechanical (QM/MM) simulations based on density functional theory (DFT/MM). The hydrogen transfer between the substrate and the cofactor is found to be rate limiting with a barrier of 101.1 kJ mol(-1). A careful comparison to the uncatalyzed reaction in water is performed. The protein influences the reaction predominantly electrostatically and to a lesser degree sterically. Our calculations shed light on the atomistic details of the reaction mechanism. The well-known arginine claw and Glu 171 ( Clostridium cochlearium notation) are found to have the strongest influence on the reaction. However, a catalytic role of Glu 214, Lys 322, Gln 147, Glu 330, Lys 326, and Met 294 is found as well. The arginine claw keeps the intermediates in place and is probably responsible for the enantioselectivity. Glu 171 temporarily accepts a proton from the glutamyl radical intermediate and donates it back at the end of the reaction. We relate our results to experimental data when available. Our simulations lead to further understanding of how glutamate mutase catalyzes the carbon skeleton rearrangement of glutamate.
Collapse
Affiliation(s)
- Judith B Rommel
- Computational Biochemistry Group, Institute of Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | | |
Collapse
|
33
|
Kamachi T, Kouno T, Doitomi K, Yoshizawa K. Generation of adenosyl radical from S-adenosylmethionine (SAM) in biotin synthase. J Inorg Biochem 2011; 105:850-7. [DOI: 10.1016/j.jinorgbio.2011.03.013] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2010] [Revised: 03/16/2011] [Accepted: 03/16/2011] [Indexed: 10/18/2022]
|
34
|
Garden JA, Kennedy AR, Mulvey RE, Robertson SD. Neutral zinc, lower-order zincate and higher-order zincate derivatives of pyrrole: synthesis and structural characterisation of zinc complexes with one, two, three or four pyrrolyl ligands. Dalton Trans 2011; 40:11945-54. [DOI: 10.1039/c1dt11430a] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
35
|
Kawatsu T, Lundberg M, Morokuma K. Protein Free Energy Corrections in ONIOM QM:MM Modeling: A Case Study for Isopenicillin N Synthase (IPNS). J Chem Theory Comput 2010; 7:390-401. [DOI: 10.1021/ct1005592] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tsutomu Kawatsu
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan, and Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Marcus Lundberg
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan, and Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan, and Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
36
|
Chung LW, Li X, Sugimoto H, Shiro Y, Morokuma K. ONIOM Study on a Missing Piece in Our Understanding of Heme Chemistry: Bacterial Tryptophan 2,3-Dioxygenase with Dual Oxidants. J Am Chem Soc 2010; 132:11993-2005. [DOI: 10.1021/ja103530v] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Lung Wa Chung
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan, and Biometal Science Laboratory, RIKEN SPring-8 Center, Harima Institute, Hyogo 679-5148, Japan
| | - Xin Li
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan, and Biometal Science Laboratory, RIKEN SPring-8 Center, Harima Institute, Hyogo 679-5148, Japan
| | - Hiroshi Sugimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan, and Biometal Science Laboratory, RIKEN SPring-8 Center, Harima Institute, Hyogo 679-5148, Japan
| | - Yoshitsugu Shiro
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan, and Biometal Science Laboratory, RIKEN SPring-8 Center, Harima Institute, Hyogo 679-5148, Japan
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan, and Biometal Science Laboratory, RIKEN SPring-8 Center, Harima Institute, Hyogo 679-5148, Japan
| |
Collapse
|
37
|
Rutledge LR, Wetmore SD. The assessment of density functionals for DNA–protein stacked and T-shaped complexes. CAN J CHEM 2010. [DOI: 10.1139/v10-046] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.
Collapse
Affiliation(s)
- Lesley R. Rutledge
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, AB T1K 3M4, Canada
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, AB T1K 3M4, Canada
| |
Collapse
|
38
|
Cho HG, Andrews L. Infrared Spectra of CX2═CoX2 and CX3−CoX Complexes from Reactions of Laser-Ablated Cobalt Atoms with Halomethanes. J Phys Chem A 2010; 114:8056-68. [PMID: 20684577 DOI: 10.1021/jp1026543] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Han-Gook Cho
- Department of Chemistry, University of Incheon, 177 Dohwa-dong, Nam-ku, Incheon, 402-749, South Korea, and Department of Chemistry, University of Virginia, P.O. Box 400319, Charlottesville, Virginia 22904-4319
| | - Lester Andrews
- Department of Chemistry, University of Incheon, 177 Dohwa-dong, Nam-ku, Incheon, 402-749, South Korea, and Department of Chemistry, University of Virginia, P.O. Box 400319, Charlottesville, Virginia 22904-4319
| |
Collapse
|
39
|
Retegan M, Milet A, Jamet H. Comparative Theoretical Studies of the Phosphomonoester Hydrolysis Mechanism by Purple Acid Phosphatases. J Phys Chem A 2010; 114:7110-6. [DOI: 10.1021/jp100478f] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- M. Retegan
- DCM, Equipe Chimie Théorique, UMR CNRS 5250, ICMG, FR CNRS, Université J. Fourier, BP. 53, 38041 Grenoble Cedex, France
| | - A. Milet
- DCM, Equipe Chimie Théorique, UMR CNRS 5250, ICMG, FR CNRS, Université J. Fourier, BP. 53, 38041 Grenoble Cedex, France
| | - H. Jamet
- DCM, Equipe Chimie Théorique, UMR CNRS 5250, ICMG, FR CNRS, Université J. Fourier, BP. 53, 38041 Grenoble Cedex, France
| |
Collapse
|
40
|
Kozlowski PM, Kamachi T, Kumar M, Nakayama T, Yoshizawa K. Theoretical Analysis of the Diradical Nature of Adenosylcobalamin Cofactor−Tyrosine Complex in B12-Dependent Mutases: Inspiring PCET-Driven Enzymatic Catalysis. J Phys Chem B 2010; 114:5928-39. [DOI: 10.1021/jp100573b] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Pawel M. Kozlowski
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, and Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| | - Takashi Kamachi
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, and Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| | - Manoj Kumar
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, and Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| | - Tomonori Nakayama
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, and Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| | - Kazunari Yoshizawa
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, and Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| |
Collapse
|
41
|
|
42
|
Lu Y, Wang Y, Xu Z, Yan X, Luo X, Jiang H, Zhu W. C-X...H contacts in biomolecular systems: how they contribute to protein-ligand binding affinity. J Phys Chem B 2009; 113:12615-21. [PMID: 19708644 DOI: 10.1021/jp906352e] [Citation(s) in RCA: 79] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The hydrogen bond acceptor capability of halogens has long been underappreciated in the field of biology. In this work, we have surveyed structures of protein complexes with halogenated ligands to characterize geometrical preferences of C-X...H contacts and contributions of such interactions to protein-ligand binding affinity. Notably, F...H interactions in biomolecules exhibit a remarkably different behavior as compared to three other kinds of X...H (X = Cl, Br, I) interactions, which has been rationalized by means of ab initio calculations using simple model systems. The C-X...H contacts in biological systems are characterized as weak hydrogen bonding interactions. Furthermore, the electrophile "head on" and nucleophile "side on" interactions of halogens have been extensively investigated through the examination of interactions in protein structures and a two-layer ONIOM-based QM/MM method. In biomolecular systems, C-X...H contacts are recognized as secondary interaction contributions to C-X...O halogen bonds that play important roles in conferring specificity and affinity for halogenated ligands. The results presented here are within the context of their potential applications in drug design, including relevance to the development of accurate force fields for halogens.
Collapse
Affiliation(s)
- Yunxiang Lu
- Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | | | | | | | | | | | | |
Collapse
|
43
|
Morokuma K. Theoretical studies of structure, function and reactivity of molecules--a personal account. PROCEEDINGS OF THE JAPAN ACADEMY. SERIES B, PHYSICAL AND BIOLOGICAL SCIENCES 2009; 85:167-82. [PMID: 19444009 PMCID: PMC3524299 DOI: 10.2183/pjab.85.167] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2008] [Accepted: 03/12/2009] [Indexed: 05/27/2023]
Abstract
Last few decades theoretical/computational studies of structure, function and reactivity of molecules have been contributing significantly in chemistry by explanation of experimental results, better understanding of underlying principles and prediction of the unknown experimental outcome. Accuracy needed in chemistry has long been established, but due to high power dependency of such accurate methods on the molecular size, it has been a major challenge to apply theoretical methods to large molecular systems. In the present article we will review some examples of such applications. One is theoretical study of growth/formation of carbon nanostructures such as fullerenes and carbon nanotubes, using quantum mechanical molecular dynamics method. For growth of single walled carbon nanotube from transition metal cluster, we have demonstrated continued growth of attached nanotube, cap formation and growth from small carbon fragments. For homogeneous catalysis we presented results of studies on N(2) activation by Zr complexes. For biomolecular reactions we use active site and protein models and show that in some catalyses the protein environment is involved in reactions and changes the preferred pathway, and in some other case the effect is modest. The review is concluded with a perspective.
Collapse
Affiliation(s)
- Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto, Japan.
| |
Collapse
|