1
|
Wei W, Siegbahn PEM, Liao R. Mechanism of the Dinuclear Iron Enzymep‐Aminobenzoate N‐oxygenase from Density Functional Calculations. ChemCatChem 2018. [DOI: 10.1002/cctc.201801072] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Wen‐Jie Wei
- 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 EngineeringHuazhong University of Science and Technology Wuhan 430074 P. R. China
| | - Per E. M. Siegbahn
- Department of Organic Chemistry, Arrhenius LaboratoryStockholm University Stockholm SE-10691 Sweden
| | - 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 EngineeringHuazhong University of Science and Technology Wuhan 430074 P. R. China
| |
Collapse
|
2
|
Abbat S, Jaladanki CK, Bharatam PV. Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis. J Mol Graph Model 2018; 87:76-88. [PMID: 30508692 DOI: 10.1016/j.jmgm.2018.11.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2018] [Revised: 11/16/2018] [Accepted: 11/19/2018] [Indexed: 11/18/2022]
Abstract
The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP+ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP+ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.
Collapse
Affiliation(s)
- Sheenu Abbat
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Chaitanya K Jaladanki
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Prasad V Bharatam
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India; Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India.
| |
Collapse
|
3
|
Wang WJ, Wei WJ, Liao RZ. Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase. Phys Chem Chem Phys 2018; 20:15784-15794. [DOI: 10.1039/c8cp02683a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.
Collapse
Affiliation(s)
- Wen-Juan Wang
- 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
| | - Wen-Jie Wei
- 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
| | - 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
| |
Collapse
|
4
|
Krámos B, Menyhárd DK, Oláh J. Direct hydride shift mechanism and stereoselectivity of P450nor confirmed by QM/MM calculations. J Phys Chem B 2012; 116:872-85. [PMID: 22148861 DOI: 10.1021/jp2080918] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Nitric oxide reductase (P450(nor)) found in Fusarium oxysporum catalyzes the reduction of nitric oxide to N(2)O in a multistep process. The reducing agent, NADH, is bound in the distal pocket of the enzyme, and direct hydride transfer occurs from NADH to the nitric oxide bound heme enzyme, forming intermediate I. Here we studied the possibility of hydride transfer from NADH to both the nitrogen and oxygen of the heme-bound nitric oxide, using quantum chemical and combined quantum mechanics/molecular mechanics (QM/MM) calculations, on two different protein models, representing both possible stereochemistries, a syn- and an anti-NADH arrangement. All calculations clearly favor hydride transfer to the nitrogen of nitric oxide, and the QM-only barrier and kinetic isotope effects are good agreement with the experimental values of intermediate I formation. We obtained higher barriers in the QM/MM calculations for both pathways, but hydride transfer to the nitrogen of nitric oxide is still clearly favored. The barriers obtained for the syn, Pro-R conformation of NADH are lower and show significantly less variation than the barriers obtained in the case of anti conformation. The effect of basis set and wide range of functionals on the obtained results are also discussed.
Collapse
Affiliation(s)
- Balázs Krámos
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Budapest, Hungary
| | | | | |
Collapse
|
5
|
Ribeiro RF, Marenich AV, Cramer CJ, Truhlar DG. Use of Solution-Phase Vibrational Frequencies in Continuum Models for the Free Energy of Solvation. J Phys Chem B 2011; 115:14556-62. [DOI: 10.1021/jp205508z] [Citation(s) in RCA: 675] [Impact Index Per Article: 51.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Raphael F. Ribeiro
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| | - Aleksandr V. Marenich
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
6
|
David C, Enescu M. Free Energy Calculations on Disulfide Bridges Reduction in Proteins by Combining ab Initio and Molecular Mechanics Methods. J Phys Chem B 2010; 114:3020-7. [DOI: 10.1021/jp910340t] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Catalina David
- Laboratoire de Chimie Physique et Rayonnement, UMR CEA E4, University of Franche-Comte, 16 route de Gray, 25030 Besancon, France
| | - Mironel Enescu
- Laboratoire de Chimie Physique et Rayonnement, UMR CEA E4, University of Franche-Comte, 16 route de Gray, 25030 Besancon, France
| |
Collapse
|
7
|
Zhang R, Lev B, Cuervo JE, Noskov SY, Salahub DR. A Guide to QM/MM Methodology and Applications. ADVANCES IN QUANTUM CHEMISTRY 2010. [DOI: 10.1016/s0065-3276(10)59010-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
|
8
|
Alonso H, Cummins PL, Gready JE. Methyltetrahydrofolate:corrinoid/iron−sulfur Protein Methyltransferase (MeTr): Protonation State of the Ligand and Active-Site Residues. J Phys Chem B 2009; 113:14787-96. [DOI: 10.1021/jp900181g] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hernán Alonso
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Peter L. Cummins
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
9
|
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
|
10
|
|
11
|
Condic-Jurkic K, Zipse H, Smith DM. A compound QM/MM procedure: Comparative performance on a pyruvate formate-lyase model system. J Comput Chem 2009; 31:1024-35. [DOI: 10.1002/jcc.21389] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
12
|
Kannappan B, Gready JE. Redefinition of rubisco carboxylase reaction reveals origin of water for hydration and new roles for active-site residues. J Am Chem Soc 2008; 130:15063-80. [PMID: 18855361 DOI: 10.1021/ja803464a] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Crystallographic, mutagenesis, kinetic, and computational studies on Rubisco over three decades have revealed much about its catalytic mechanism and the role played by several active-site residues. However, key questions remain unanswered. Specific details of the carboxylase and oxygenase mechanisms, required to underpin the rational re-engineering of Rubisco, are still speculative. Here we address critical gaps in knowledge with a definitive comprehensive computational investigation of the mechanism of carboxylase activity at the Rubisco active site. Density functional theory calculations (B3LYP/6-31G(d,p)) were performed on active-site fragment models of a size up to 77 atoms, not previously possible computationally. All amino acid residues suspected to play roles in the acid-base chemistry in the multistep reaction, and interacting directly with the central Mg (2+) atom and the reactive moiety of substrate and intermediates, were included. The results provide a firm basis for us to propose a novel mechanism for the entire sequence of reactions in the carboxylase catalysis and to define precise roles for the active-site residues, singly and in concert. In this mechanism, the carbamylated LYS201 plays a more limited role than previously proposed but is crucial for initiating the reaction by acting as a base in the enolization. We suggest a wider role for HIS294, with involvement in the carboxylation, hydration, and C2-C3 bond-scission steps, consistent with the suggestion of Harpel et al. (1998) but contrary to the consensus view of Cleland et al. (1998). In contrast to the common assumption that the water molecule for the hydration step comes from within the active site, we propose that the Mg-coordinated water is not dissociated at the start of the gas-addition reaction but rather remains coordinated and is used for the hydration of the C3 carbon atom. New roles are also proposed for LYS175, GLU204, and HIS294. The mechanism suggests roles in the gas-addition step for residues in three spatially distinct regions of the active site, HIS294 and LYS334 in the C-terminal domain of the large subunit (LSU), but also hitherto unsuspected roles for a cluster of three residues (ASN123, GLU60, and TYR20) in the N-terminal domain of the partner LSU of the dimer containing the active site. Our new mechanism is supported by existing experimental data, provides new convincing interpretations of previously puzzling data, and allows new insights into mutational strategies for improving Rubisco activity.
Collapse
Affiliation(s)
- Babu Kannappan
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, Canberra, ACT 0200, Australia
| | | |
Collapse
|