1
|
Świderek K, Martí S, Arafet K, Moliner V. Computational study of the mechanism of a polyurethane esterase A (PueA) from Pseudomonas chlororaphis. Faraday Discuss 2024; 252:323-340. [PMID: 38836643 DOI: 10.1039/d4fd00022f] [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: 06/06/2024]
Abstract
The effective management of plastic waste has become a global imperative, given our reliance on a linear model in which plastics are manufactured, used once, and then discarded. This has led to the pervasive accumulation of plastic debris in landfills and environmental contamination. Recognizing this issue, numerous initiatives are underway to address the environmental repercussions associated with plastic disposal. In this study, we investigate the possible molecular mechanism of polyurethane esterase A (PueA), which has been previously identified as responsible for the degradation of a polyester polyurethane (PU) sample in Pseudomonas chlororaphis, as an effort to develop enzymatic biodegradation solutions. After generating the unsolved 3D structure of the protein by AlphaFold2 from its known genome, the enzymatic hydrolysis of the same model PU compound previously used in experiments has been explored employing QM/MM molecular dynamics simulations. This required a preliminary analysis of the 3D structure of the apo-enzyme, identifying the putative active site, and the search for the optimal protein-substrate binding site. Finally, the resulting free energy landscape indicates that wild-type PueA can degrade PU chains, although with low-level activity. The reaction takes place by a characteristic four-step path of the serine hydrolases, involving an acylation followed by a diacylation step. Energetics and structural analysis of the evolution of the active site along the reaction suggests that PueA can be considered a promising protein scaffold for further development to achieve efficient biodegradation of PU.
Collapse
Affiliation(s)
- Katarzyna Świderek
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain.
| | - Sergio Martí
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain.
| | - Kemel Arafet
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain.
| | - Vicent Moliner
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain.
| |
Collapse
|
2
|
Akintola O, Farren-Dai M, Ren W, Bhosale S, Britton R, Świderek K, Moliner V, Bennet AJ. Glycoside Hydrolase Catalysis: Do Substrates and Mechanism-Based Covalent Inhibitors React via Matching Transition States? ACS Catal 2022. [DOI: 10.1021/acscatal.2c04027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Affiliation(s)
- Oluwafemi Akintola
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| | - Marco Farren-Dai
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| | - Weiwu Ren
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| | - Sandeep Bhosale
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| | - Robert Britton
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| | - Katarzyna Świderek
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071Castellón, Spain
| | - Vicent Moliner
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071Castellón, Spain
| | - Andrew J. Bennet
- Department of Chemistry, Simon Fraser University, Burnaby, British ColumbiaV5A 1S6, Canada
| |
Collapse
|
3
|
Serrano-Aparicio N, Ferrer S, Świderek K. Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies. Pharmaceuticals (Basel) 2022; 15:ph15050531. [PMID: 35631358 PMCID: PMC9143130 DOI: 10.3390/ph15050531] [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/31/2022] [Revised: 04/21/2022] [Accepted: 04/22/2022] [Indexed: 12/02/2022] Open
Abstract
20S proteasome is a main player in the protein degradation pathway in the cytosol, thus intervening in multiple pivotal cellular processes. Over the years the proteasome has emerged as a crucial target for the treatment of many diseases such as neurodegenerative diseases, cancer, autoimmune diseases, developmental disorders, cystic fibrosis, diabetes, cardiac diseases, atherosclerosis, and aging. In this work, the mechanism of proteasome covalent inhibition with bisbenzyl-protected homobelactosin C (hBelC) was explored using quantum mechanics/molecular mechanics (QM/MM) methods. Molecular dynamic simulations were used to describe key interactions established between the hBelC and its unique binding mode in the primed site of the β5 subunit. The free energy surfaces were computed to characterize the kinetics and thermodynamics of the inhibition process. This study revealed that although the final inhibition product for hBelC is formed according to the same molecular mechanism as one described for hSalA, the free energy profile of the reaction pathway differs significantly from the one previously reported for γ-lactam-β-lactone containing inhibitors in terms of the height of the activation barrier as well as the stabilization of the final product. Moreover, it was proved that high stabilization of the covalent adduct formed between β5-subunit and hBelC, together with the presence of aminocarbonyl side chain in the structure of the inhibitor which prevents the hydrolysis of the ester bond from taking place, determines its irreversible character.
Collapse
|
4
|
Serrano-Aparicio N, Moliner V, Świderek K. On the Origin of the Different Reversible Characters of Salinosporamide A and Homosalinosporamide A in the Covalent Inhibition of the Human 20S Proteasome. ACS Catal 2021. [DOI: 10.1021/acscatal.1c02614] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Natalia Serrano-Aparicio
- Biocomp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Biocomp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| | - Katarzyna Świderek
- Biocomp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
5
|
Kim B, Snyder R, Nagaraju M, Zhou Y, Ojeda-May P, Keeton S, Hege M, Shao Y, Pu J. Reaction Path-Force Matching in Collective Variables: Determining Ab Initio QM/MM Free Energy Profiles by Fitting Mean Force. J Chem Theory Comput 2021; 17:4961-4980. [PMID: 34283604 PMCID: PMC9064116 DOI: 10.1021/acs.jctc.1c00245] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
First-principles determination of free energy profiles for condensed-phase chemical reactions is hampered by the daunting costs associated with configurational sampling on ab initio quantum mechanical/molecular mechanical (AI/MM) potential energy surfaces. Here, we report a new method that enables efficient AI/MM free energy simulations through mean force fitting. In this method, a free energy path in collective variables (CVs) is first determined on an efficient reactive aiding potential. Based on the configurations sampled along the free energy path, correcting forces to reproduce the AI/MM forces on the CVs are determined through force matching. The AI/MM free energy profile is then predicted from simulations on the aiding potential in conjunction with the correcting forces. Such cycles of correction-prediction are repeated until convergence is established. As the instantaneous forces on the CVs sampled in equilibrium ensembles along the free energy path are fitted, this procedure faithfully restores the target free energy profile by reproducing the free energy mean forces. Due to its close connection with the reaction path-force matching (RP-FM) framework recently introduced by us, we designate the new method as RP-FM in collective variables (RP-FM-CV). We demonstrate the effectiveness of this method on a type-II solution-phase SN2 reaction, NH3 + CH3Cl (the Menshutkin reaction), simulated with an explicit water solvent. To obtain the AI/MM free energy profiles, we employed the semiempirical AM1/MM Hamiltonian as the base level for determining the string minimum free energy pathway, along which the free energy mean forces are fitted to various target AI/MM levels using the Hartree-Fock (HF) theory, density functional theory (DFT), and the second-order Møller-Plesset perturbation (MP2) theory as the AI method. The forces on the bond-breaking and bond-forming CVs at both the base and target levels are obtained by force transformation from Cartesian to redundant internal coordinates under the Wilson B-matrix formalism, where the linearized FM is facilitated by the use of spline functions. For the Menshutkin reaction tested, our FM treatment greatly reduces the deviations on the CV forces, originally in the range of 12-33 to ∼2 kcal/mol/Å. Comparisons with the experimental and benchmark AI/MM results, tests of the new method under a variety of simulation protocols, and analyses of the solute-solvent radial distribution functions suggest that RP-FM-CV can be used as an efficient, accurate, and robust method for simulating solution-phase chemical reactions.
Collapse
Affiliation(s)
- Bryant Kim
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Ryan Snyder
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Mulpuri Nagaraju
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Yan Zhou
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Pedro Ojeda-May
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Seth Keeton
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Mellisa Hege
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of
Oklahoma, 101 Stephenson Pkwy, Norman, OK 73019
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana
University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN
46202
| |
Collapse
|
6
|
Alvarado O, García-Meseguer R, Ruiz-Pernía JJ, Tuñon I, Delgado EJ. Mechanistic study of the biosynthesis of R-phenylacetylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations. Arch Biochem Biophys 2021; 707:108849. [PMID: 33832752 DOI: 10.1016/j.abb.2021.108849] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cβ atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4 kcal mol-1 at 27 °C. From this point, the reaction continues with the abstraction of Hβ atom of the HEThDP intermediate by the Oβ atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH+ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9 kcal mol-1 at 27 °C.
Collapse
Affiliation(s)
- Omar Alvarado
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile; Departamento de Química, Facultad de Ciencias, Universidad del Bío-Bío, Avenida Collao 1202, Concepción, Chile
| | - Rafael García-Meseguer
- School of Mathematics, University of Bristol, Bristol, UK; Department of Physical Chemistry, Universitat de Valencia, 46100, Burjassot, Spain
| | | | - Iñaki Tuñon
- Department of Physical Chemistry, Universitat de Valencia, 46100, Burjassot, Spain
| | - Eduardo J Delgado
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| |
Collapse
|
7
|
Serrano-Aparicio N, Moliner V, Świderek K. Nature of Irreversible Inhibition of Human 20S Proteasome by Salinosporamide A. The Critical Role of Lys–Asp Dyad Revealed from Electrostatic Effects Analysis. ACS Catal 2021. [DOI: 10.1021/acscatal.0c05313] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Katarzyna Świderek
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
8
|
Alvarado O, García-Meseguer R, Ruiz-Pernía JJ, Tuñon I, Delgado EJ. Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations. Arch Biochem Biophys 2021; 701:108807. [PMID: 33587902 DOI: 10.1016/j.abb.2021.108807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 02/06/2021] [Accepted: 02/08/2021] [Indexed: 11/26/2022]
Abstract
The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cβ atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4kcal mol-1 at 27 °C. From this point, the reaction continues with the abstraction of Hβ atom of the HEThDP intermediate by the Oβ atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH+ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9kcal mol-1 at 27 °C.
Collapse
Affiliation(s)
- Omar Alvarado
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile; Departamento de Química, Facultad de Ciencias, Universidad del Bío-Bío, Avenida Collao 1202, Concepción, Chile
| | - Rafael García-Meseguer
- School of Mathematics, University of Bristol, Bristol, United Kingdom; Department of Physical Chemistry, Universitat de Valencia, 46100, Burjassot, Spain
| | | | - Iñaki Tuñon
- Department of Physical Chemistry, Universitat de Valencia, 46100, Burjassot, Spain
| | - Eduardo J Delgado
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| |
Collapse
|
9
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
10
|
De Raffele D, Martí S, Moliner V. Understanding the Directed Evolution of De Novo Retro-Aldolases from QM/MM Studies. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01165] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Daria De Raffele
- Departament de Química Física i Analítica; Universitat Jaume I., 12071 Castellón, Spain
| | - Sergio Martí
- Departament de Química Física i Analítica; Universitat Jaume I., 12071 Castellón, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica; Universitat Jaume I., 12071 Castellón, Spain
| |
Collapse
|
11
|
Galmés MÀ, García-Junceda E, Świderek K, Moliner V. Exploring the Origin of Amidase Substrate Promiscuity in CALB by a Computational Approach. ACS Catal 2019. [DOI: 10.1021/acscatal.9b04002] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Miquel À Galmés
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Eduardo García-Junceda
- Departamento de Química Orgánica Biológica, Instituto de Química Orgánica General, CSIC Juan de la Cierva 3, 28006 Madrid, Spain
| | - Katarzyna Świderek
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
12
|
Ahsan M, Senapati S. Water Plays a Cocatalytic Role in Epoxide Ring Opening Reaction in Aspartate Proteases: A QM/MM Study. J Phys Chem B 2019; 123:7955-7964. [PMID: 31468966 DOI: 10.1021/acs.jpcb.9b04575] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Aspartate proteases are potential targets for various diseases, and many of their inhibitors are FDA-approved drugs. However, these peptidomimetic and reversibly bound drugs become ineffective upon prolonged use. Attempts have been made to design and synthesize various nonpeptidic epoxide-based irreversible inhibitors to combat the drug-resistance enigma. Here, we study the mechanism of epoxide ring opening in two widely studied aspartate proteases, HIV-1 protease and pepsin. Our results from QM/MM molecular dynamics show that the epoxide ring opening in aspartate proteases follow a two-step mechanism with the formation of an oxyanion intermediate, stabilized by a set of water molecules in the protein active site. These water molecules by virtue of "low-barrier hydrogen bonds" with the epoxide ring reduce the intrinsic reaction barrier while remaining structurally unperturbed, thus playing a cocatalytic role. We validated our results by reproducing the experimentally observed protease/pepsin-epoxide covalent complexes as end products. The observed stability of our oxyanion intermediate in a four-water-coordinated state is also consistent with the reported stable state of the hydroxide ion in water as OH-(H2O)4. Our study could pave the way for the design of new class "HIV protease irreversible inhibitors" from the acquired knowledge of the structures of intermediate and transition states traced during the explored reaction mechanism.
Collapse
Affiliation(s)
- Mohd Ahsan
- Department of Biotechnology and BJM School of Biosciences , Indian Institute of Technology Madras , Chennai 600036 , India
| | - Sanjib Senapati
- Department of Biotechnology and BJM School of Biosciences , Indian Institute of Technology Madras , Chennai 600036 , India
| |
Collapse
|
13
|
De Raffele D, Martí S, Moliner V. QM/MM Theoretical Studies of a de Novo Retro-Aldolase Design. ACS Catal 2019. [DOI: 10.1021/acscatal.8b04457] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Daria De Raffele
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Sergio Martí
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
14
|
Theoretical study of the inhibition mechanism of human 20S proteasome by dihydroeponemycin. Eur J Med Chem 2019; 164:399-407. [DOI: 10.1016/j.ejmech.2018.12.062] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Revised: 12/12/2018] [Accepted: 12/24/2018] [Indexed: 01/10/2023]
|
15
|
Wu J, Shen L, Yang W. Internal force corrections with machine learning for quantum mechanics/molecular mechanics simulations. J Chem Phys 2018; 147:161732. [PMID: 29096448 DOI: 10.1063/1.5006882] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ab initio quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulation is a useful tool to calculate thermodynamic properties such as potential of mean force for chemical reactions but intensely time consuming. In this paper, we developed a new method using the internal force correction for low-level semiempirical QM/MM molecular dynamics samplings with a predefined reaction coordinate. As a correction term, the internal force was predicted with a machine learning scheme, which provides a sophisticated force field, and added to the atomic forces on the reaction coordinate related atoms at each integration step. We applied this method to two reactions in aqueous solution and reproduced potentials of mean force at the ab initio QM/MM level. The saving in computational cost is about 2 orders of magnitude. The present work reveals great potentials for machine learning in QM/MM simulations to study complex chemical processes.
Collapse
Affiliation(s)
- Jingheng Wu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Lin Shen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
16
|
Świderek K, Tuñón I, Williams IH, Moliner V. Insights on the Origin of Catalysis on Glycine N-Methyltransferase from Computational Modeling. J Am Chem Soc 2018; 140:4327-4334. [PMID: 29460630 PMCID: PMC6613375 DOI: 10.1021/jacs.7b13655] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor-acceptor distance (DAD): Is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The computed activation free energies and KIEs are in good agreement with experimental data, but the mutations do not meaningfully affect the DAD: Compression cannot explain the experimental variations on KIEs. On the contrary, electrostatic properties in the active site correlate with the catalytic activity of wild type and mutants. The plasticity of the enzyme moderates the effects of the mutations, explaining the rather small degree of variation in KIEs and reactivities.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Departament de Química Física i Analítica; Universitat Jaume I, 12071 Castellón (Spain)
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjasot (Spain)
| | - Ian H. Williams
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| | - Vicent Moliner
- Departament de Química Física i Analítica; Universitat Jaume I, 12071 Castellón (Spain)
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| |
Collapse
|
17
|
Abstract
Approaches to determine chlorine kinetic isotope effects (Cl-KIEs) on enzymatic dehalogenations are discussed and illustrated by representative examples. Three aspects are considered. First methodology for experimental measurement of Cl-KIEs, with stress being on FAB-IRMS technique developed in our laboratory, is described. Subsequently, we concentrate our discussion on the consequences of reaction complexity in the interpretation of experimental values, a problem especially important in cases of polychlorinated reactants. The most fruitful studies of enzymatic dehalogenations by Cl-KIEs require their theoretical evaluation, hence the computational focus of the second part of this chapter.
Collapse
|
18
|
Bordes I, Castillo R, Moliner V. Theoretical Study of the Phosphoryl Transfer Reaction from ATP to Dha Catalyzed by DhaK from Escherichia coli. J Phys Chem B 2017; 121:8878-8892. [PMID: 28850238 DOI: 10.1021/acs.jpcb.7b04862] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Protein kinases, representing one of the largest protein families involved in almost all aspects of cell life, have become one of the most important targets for the development of new drugs to be used in, for instance, cancer treatments. In this article an exhaustive theoretical study of the phosphoryl transfer reaction from adenosine triphosphate (ATP) to dihydroxyacetone (Dha) catalyzed by DhaK from Escherichia coli (E. coli) is reported. Two different mechanisms, previously proposed for the phosphoryl transfer from ATP to the hydroxyl side chain of specific serine, threonine, or tyrosine residues, have been explored based on the generation of free energy surfaces (FES) computed with hybrid QM/MM potentials. The results suggest that the substrate-assisted phosphoryl and proton-transfer mechanism is kinetically more favorable than the mechanism where an aspartate would be activating the Dha. Although the details of the mechanisms appear to be dramatically dependent on the level of theory employed in the calculations (PM3/MM, B3LYP:PM3/MM, or B3LYP/MM), the transition states (TSs) for the phosphoryl transfer step appear to be described as a concerted step with different degrees of synchronicity in the breaking and forming bonds process in both explored mechanisms. Residues of the active site belonging to different subunits of the protein, such as Gly78B, Thr79A, Ser80A, Arg178B, and one Mg2+ cation, would be stabilizing the transferred phosphate in the TS. Asp109A would have a structural role by posing the Dha and other residues of the active site in the proper orientation. The information derived from our calculations not only reveals the role of the enzyme and the particular residues of its active site, but it can assist in the rational design of new more specific inhibitors.
Collapse
Affiliation(s)
- I Bordes
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - R Castillo
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - V Moliner
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain.,Department of Chemistry, University of Bath , Bath BA2 7AY, United Kingdom
| |
Collapse
|
19
|
Świderek K, Arafet K, Kohen A, Moliner V. Benchmarking Quantum Mechanics/Molecular Mechanics (QM/MM) Methods on the Thymidylate Synthase-Catalyzed Hydride Transfer. J Chem Theory Comput 2017; 13:1375-1388. [PMID: 28192669 PMCID: PMC5371049 DOI: 10.1021/acs.jctc.6b01032] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Given the ubiquity of hydride-transfer reactions in enzyme-catalyzed processes, identifying the appropriate computational method for evaluating such biological reactions is crucial to perform theoretical studies of these processes. In this paper, the hydride-transfer step catalyzed by thymidylate synthase (TSase) is studied by examining hybrid quantum mechanics/molecular mechanics (QM/MM) potentials via multiple semiempirical methods and the M06-2X hybrid density functional. Calculations of protium and tritium transfer in these reactions across a range of temperatures allowed calculation of the temperature dependence of kinetic isotope effects (KIE). Dynamics and quantum-tunneling effects are revealed to have little effect on the reaction rate, but are significant in determining the KIEs and their temperature dependence. A good agreement with experiments is found, especially when computed for RM1/MM simulations. The small temperature dependence of quantum tunneling corrections and the quasiclassical contribution term cancel each other, while the recrossing transmission coefficient seems to be temperature-independent over the interval of 5-40 °C.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
- Institute of Applied Radiation Chemistry, Lodz University of Technology , 90-924 Lodz, Poland
| | - Kemel Arafet
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| | - Amnon Kohen
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| |
Collapse
|
20
|
Smyrnova D, Zinovjev K, Tuñón I, Ceulemans A. Thermal Isomerization Mechanism in Dronpa and Its Mutants. J Phys Chem B 2016; 120:12820-12825. [PMID: 28002952 DOI: 10.1021/acs.jpcb.6b10859] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The photoswitching speed of the reversibly switchable fluorescent proteins (RSFPs) from the family of green fluorescent proteins (GFPs) changes upon mutation which is of direct importance for various high-resolution techniques. Dronpa is one of the most used RSFPs. Its point mutants rsFastLime (Dronpa V157G) and rsKame (Dronpa V157L) exhibit a striking difference in their photoswitching speed. Here the QM/MM on-the-fly string method is used in order to explore the details of the thermal isomerization mechanism. The four principal ways in which isomerization may occur have been scrutinized for each of the three proteins. It has been shown that thermal isomerization occurs via a one-bond-flip mechanism in all three proteins, although, in rsKame, where the chromophore is constrained more, the activation free energy difference between hula-twist and one-bond-flip is significantly smaller. Functional mode analysis has been applied to examine the motions of the amino acids during the isomerization. It clearly identifies the importance of Val/Leu 157 as well as the amino acids in the α-helix during the isomerization.
Collapse
Affiliation(s)
- Daryna Smyrnova
- Quantum Chemistry and Physical Chemistry Division, Department of Chemistry, KU Leuven , Celestijnenlaan 200F, 3001 Heverlee, Belgium
| | - Kirill Zinovjev
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain
| | - Arnout Ceulemans
- Quantum Chemistry and Physical Chemistry Division, Department of Chemistry, KU Leuven , Celestijnenlaan 200F, 3001 Heverlee, Belgium
| |
Collapse
|
21
|
Krzemińska A, Moliner V, Świderek K. Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease. J Am Chem Soc 2016; 138:16283-16298. [PMID: 27935692 DOI: 10.1021/jacs.6b06856] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.
Collapse
Affiliation(s)
- Agnieszka Krzemińska
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology , Zeromskiego 116, 90-924 Lodz, Poland
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| | - Katarzyna Świderek
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology , Zeromskiego 116, 90-924 Lodz, Poland.,Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castelló, Spain
| |
Collapse
|
22
|
Swiderek K, Kohen A, Moliner V. The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism. Phys Chem Chem Phys 2016; 17:30793-804. [PMID: 25868526 DOI: 10.1039/c5cp01239b] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene-dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate-enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results.
Collapse
Affiliation(s)
- Katarzyna Swiderek
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain. and Institute of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA 52242, USA
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain.
| |
Collapse
|
23
|
Aranda J, Zinovjev K, Świderek K, Roca M, Tuñón I. Unraveling the Reaction Mechanism of Enzymatic C5-Cytosine Methylation of DNA. A Combined Molecular Dynamics and QM/MM Study of Wild Type and Gln119 Variant. ACS Catal 2016. [DOI: 10.1021/acscatal.6b00394] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Juan Aranda
- Departament
de Quı́mica Fı́sica, Universitat de València, 46100 Burjassot, Spain
| | - Kirill Zinovjev
- Departament
de Quı́mica Fı́sica, Universitat de València, 46100 Burjassot, Spain
| | - Katarzyna Świderek
- Departament
de Quı́mica Fı́sica, Universitat de València, 46100 Burjassot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Maite Roca
- Departament
de Quı́mica Fı́sica, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament
de Quı́mica Fı́sica, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
24
|
Gillet N, Ruiz-Pernía JJ, de la Lande A, Lévy B, Lederer F, Demachy I, Moliner V. QM/MM study of l-lactate oxidation by flavocytochrome b2. Phys Chem Chem Phys 2016; 18:15609-18. [DOI: 10.1039/c6cp00395h] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Free energy surfaces calculated from a state-of-the-art computational methodology highlight the role of active site residues in l-lactate oxidation by flavocytochrome b2.
Collapse
Affiliation(s)
- N. Gillet
- Laboratoire de Chimie-Physique
- CNRS UMR 8000
- Université Paris Sud
- 91405 Cedex Orsay
- France
| | - J. J. Ruiz-Pernía
- Departament de Química Física i Analítica
- Universitat Jaume I
- 12071 Castellón
- Spain
| | - A. de la Lande
- Laboratoire de Chimie-Physique
- CNRS UMR 8000
- Université Paris Sud
- 91405 Cedex Orsay
- France
| | - B. Lévy
- Laboratoire de Chimie-Physique
- CNRS UMR 8000
- Université Paris Sud
- 91405 Cedex Orsay
- France
| | - F. Lederer
- Laboratoire de Chimie-Physique
- CNRS UMR 8000
- Université Paris Sud
- 91405 Cedex Orsay
- France
| | - I. Demachy
- Laboratoire de Chimie-Physique
- CNRS UMR 8000
- Université Paris Sud
- 91405 Cedex Orsay
- France
| | - V. Moliner
- Departament de Química Física i Analítica
- Universitat Jaume I
- 12071 Castellón
- Spain
| |
Collapse
|
25
|
Świderek K, Moliner V. Computational Studies of Candida Antarctica Lipase B to Test Its Capability as a Starting Point To Redesign New Diels-Alderases. J Phys Chem B 2015; 120:2053-70. [PMID: 26624234 DOI: 10.1021/acs.jpcb.5b10527] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The design of new biocatalysts is a target that is receiving increasing attention. One of the most popular reactions in this regard is the Diels-Alder cycloaddition because of its applications in organic synthesis and the absence of efficient natural enzymes that catalyze it. In this paper, the possibilities of using the highly promiscuous Candida Antarctica lipase B as a protein scaffold to redesign a Diels-Alderase has been explored by means of theoretical quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. Free energy surfaces have been computed for two reactions in the wild-type and in several mutants with hybrid AM1/MM potentials with corrections at M06-2X/MM level. The study of the counterpart reactions in solution has allowed performing comparative analysis that render interesting conclusions. Since the dienophile anchors very well in the oxyanion hole of all tested protein variants, the slight electronic changes from reactant complex to the transition state suggest that mutations should be focused in favoring the formation of reactive conformations of a reactant complex that, in turn, would reduce the energy barrier.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Institute of Applied Radiation Chemistry, Lodz University of Technology , 90-924 Lodz, Poland.,Departamento de Química Física y Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - Vicent Moliner
- Departamento de Química Física y Analítica, Universitat Jaume I , 12071 Castellón, Spain
| |
Collapse
|
26
|
Manna RN, Zinovjev K, Tuñón I, Dybala-Defratyka A. Dehydrochlorination of Hexachlorocyclohexanes Catalyzed by the LinA Dehydrohalogenase. A QM/MM Study. J Phys Chem B 2015; 119:15100-9. [DOI: 10.1021/acs.jpcb.5b07538] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rabindra Nath Manna
- Institute
of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, Lodz 90-924, Poland
| | - Kirill Zinovjev
- Departament
de Química Física, Universitat de Valéncia, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat de Valéncia, 46100 Burjassot, Spain
| | - Agnieszka Dybala-Defratyka
- Institute
of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, Lodz 90-924, Poland
| |
Collapse
|
27
|
Bordes I, Ruiz-Pernía JJ, Castillo R, Moliner V. A computational study of the phosphoryl transfer reaction between ATP and Dha in aqueous solution. Org Biomol Chem 2015; 13:10179-90. [PMID: 26303076 DOI: 10.1039/c5ob01079a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Phosphoryl transfer reactions are ubiquitous in biology, being involved in processes ranging from energy and signal transduction to the replication genetic material. Dihydroxyacetone phosphate (Dha-P), an intermediate of the synthesis of pyruvate and a very important building block in nature, can be generated by converting free dihydroxyacetone (Dha) through the action of the dihydroxyacetone kinase enzyme. In this paper the reference uncatalyzed reaction in solution has been studied in order to define the foundations of the chemical reaction and to determine the most adequate computational method to describe this electronically complex reaction. In particular, the phosphorylation reaction mechanism between adenosine triphosphate (ATP) and Dha in aqueous solution has been studied by means of quantum mechanics/molecular mechanics (QM/MM) Molecular Dynamics (MD) simulations with the QM subset of atoms described with semi-empirical and DFT methods. The results appear to be strongly dependent on the level of calculation, which will have to be taken into account for future studies of the reaction catalyzed by enzymes. In particular, PM3/MM renders lower free energy barriers and a less endergonic process than AM1d/MM and PM6/MM methods. Nevertheless, the concerted pathway was not located with the former combination of potentials.
Collapse
Affiliation(s)
- I Bordes
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | | | | | | |
Collapse
|
28
|
Alves CN, Silva JRA, Roitberg AE. Insights into the mechanism of oxidation of dihydroorotate to orotate catalysed by human class 2 dihydroorotate dehydrogenase: a QM/MM free energy study. Phys Chem Chem Phys 2015; 17:17790-6. [PMID: 26087682 DOI: 10.1039/c5cp02016f] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dihydroorotate dehydrogenase (DHOD) enzyme catalyzes the unique redox reaction in the de novo pyrimidine biosynthesis pathway. In this reaction, the oxidation of dihydroorotate (DHO) to orotate (OA) and reduction of the flavin mononucleotide (FMN) cofactor is catalysed by DHOD. The class 2 DHOD, to which the human enzyme belongs, was experimentally shown to follow a stepwise mechanism but the data did not allow the determination of the order of bond-breaking in a stepwise oxidation of DHO. The goal of this study is to understand the reaction mechanism at the molecular level of class 2 DHOD, which may aid in the design of inhibitors that selectively impact the activity of only certain members of the enzyme family. In this paper, the catalytic mechanism of oxidation of DHO to OA in human DHOD was studied using a hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) approach and Molecular Dynamics (MD) simulations. The free energy barriers calculated reveal that the mechanism in human DHOD occurs via a stepwise reaction pathway. In the first step, a proton is abstracted from the C5 of DHO to the deprotonated Ser215 side chain. Whereas, in the second step, the transfer of the hydride or hydride equivalent from the C6 of DHO to the N5 of FMN, where free energy barrier calculated by the DFT/MM level is 10.84 kcal mol(-1). Finally, a residual decomposition analysis was carried out in order to elucidate the influence of the catalytic region residues during DHO oxidation.
Collapse
Affiliation(s)
- Cláudio Nahum Alves
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, Belém, Brazil. mail:
| | | | | |
Collapse
|
29
|
Świderek K, Tuñón I, Moliner V, Bertran J. Protein Flexibility and Preorganization in the Design of Enzymes. The Kemp Elimination Catalyzed by HG3.17. ACS Catal 2015. [DOI: 10.1021/cs501904w] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Katarzyna Świderek
- Departament
de Química Física, Universitat de València, 46100 Burjasot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjasot, Spain
| | - Vicent Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Joan Bertran
- Departament
de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| |
Collapse
|
30
|
Theoretical studies of the hydrolysis of antibiotics catalyzed by a metallo-β-lactamase. Arch Biochem Biophys 2015; 582:116-26. [PMID: 25622886 DOI: 10.1016/j.abb.2015.01.013] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2014] [Revised: 01/13/2015] [Accepted: 01/14/2015] [Indexed: 12/27/2022]
Abstract
In this paper, hybrid QM/MM molecular dynamics (MD) simulations have been performed to explore the mechanisms of hydrolysis of two antibiotics, Imipenen (IMI), an antibiotic belonging to the subgroup of carbapenems, and the Cefotaxime (CEF), a third-generation cephalosporin antibiotic, in the active site of a mono-nuclear β-lactamase, CphA from Aeromonas hydrophila. Significant different transition state structures are obtained for the hydrolysis of both antibiotics: while the TS of the CEF is an ionic species with negative charge on nitrogen, the IMI TS presents a tetrahedral-like character with negative charge on oxygen atom of the carbonyl group of the lactam ring. Thus, dramatic conformational changes can take place in the cavity of CphA to accommodate different substrates, which would be the origin of its substrate promiscuity. Since CphA shows only activity against carbapenem antibiotic, this study sheds some light into the origin of the selectivity of the different MbL and, as a consequence, into the discovery of specific and potent MβL inhibitors against a broad spectrum of bacterial pathogens. We have finally probed that a re-parametrization of semiempirical methods should be done to properly describe the behavior the metal cation in active site, Zn(2+), when used in QM/MM calculations.
Collapse
|
31
|
Silva JRA, Roitberg AE, Alves CN. A QM/MM free energy study of the oxidation mechanism of dihydroorotate dehydrogenase (class 1A) from Lactococcus lactis. J Phys Chem B 2015; 119:1468-73. [PMID: 25564307 DOI: 10.1021/jp512860r] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The dihydroorotate dehydrogenase (DHOD) class 1A enzyme catalyzes is the only redox enzyme in the biosynthetic pathway (de novo) of pyrimidines where dihydroorotate (DHO) is oxidized to orotate (OA) coupled to reduction of a flavin mononucleotide (FMN) cofactor. The rupture of two DHO C-H bonds can proceed in a concerted or stepwise way. Herein, the catalytic mechanism of DHOD from Lactococcus lactis involving DHO oxidation (first half-reaction) was described using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach and molecular dynamics simulations. The free energy profile obtained from self-consistent charge-density functional tight binding/molecular mechanics calculations (corrected by DFT/MM) reveals that this occurs with the proton abstraction from DHO C5 to Cys130 deprotonated and DHO H6 is transferred to FMN N5 in a concerted mechanism with a very low barrier of 5.64 kcal/mol. Finally, through a residual decomposition analysis, the residues that have the main influence on the redox reaction were identified.
Collapse
Affiliation(s)
- José Rogério A Silva
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará , Belém, PA 66075-110, Brazil
| | | | | |
Collapse
|
32
|
Świderek K, Tuñón I, Martí S, Moliner V. Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase. ACS Catal 2015; 5:1172-1185. [PMID: 25705562 DOI: 10.1021/cs501704f] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In the last decade L-Lactate Dehydrogenase (LDH) has become an extremely useful marker in both clinical diagnosis and in monitoring the course of many human diseases. It has been assumed from the 80s that the full catalytic process of LDH starts with the binding of the cofactor and the substrate followed by the enclosure of the active site by a mobile loop of the protein before the reaction to take place. In this paper we show that the chemical step of the LDH catalyzed reaction can proceed within the open loop conformation, and the different reactivity of the different protein conformations would be in agreement with the broad range of rate constants measured in single molecule spectrometry studies. Starting from a recently solved X-ray diffraction structure that presented an open loop conformation in two of the four chains of the tetramer, QM/MM free energy surfaces have been obtained at different levels of theory. Depending on the level of theory used to describe the electronic structure, the free energy barrier for the transformation of pyruvate into lactate with the open conformation of the protein varies between 12.9 and 16.3 kcal/mol, after quantizing the vibrations and adding the contributions of recrossing and tunneling effects. These values are very close to the experimentally deduced one (14.2 kcal·mol-1) and ~2 kcal·mol-1 smaller than the ones obtained with the closed loop conformer. Calculation of primary KIEs and IR spectra in both protein conformations are also consistent with our hypothesis and in agreement with experimental data. Our calculations suggest that the closure of the active site is mainly required for the inverse process; the oxidation of lactate to pyruvate. According to this hypothesis H4 type LDH enzyme molecules, where it has been propose that lactate is transformed into pyruvate, should have a better ability to close the mobile loop than the M4 type LDH molecules.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Sergio Martí
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
| | - Vicent Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
| |
Collapse
|
33
|
Moonsamy S, Bhakat S, Soliman MES. Dynamic features of apo and bound HIV-Nef protein reveal the anti-HIV dimerization inhibition mechanism. J Recept Signal Transduct Res 2014; 35:346-56. [DOI: 10.3109/10799893.2014.984310] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
34
|
Aranda J, Zinovjev K, Roca M, Tuñón I. Dynamics and Reactivity in Thermus aquaticus N6-Adenine Methyltransferase. J Am Chem Soc 2014; 136:16227-39. [DOI: 10.1021/ja5077124] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Juan Aranda
- Departament de Química
Física, Universitat de València, 46100 Burjassot, Spain
| | - Kirill Zinovjev
- Departament de Química
Física, Universitat de València, 46100 Burjassot, Spain
| | - Maite Roca
- Departament de Química
Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química
Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
35
|
Polyak I, Benighaus T, Boulanger E, Thiel W. Quantum mechanics/molecular mechanics dual Hamiltonian free energy perturbation. J Chem Phys 2014; 139:064105. [PMID: 23947841 DOI: 10.1063/1.4817402] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dual Hamiltonian free energy perturbation (DH-FEP) method is designed for accurate and efficient evaluation of the free energy profile of chemical reactions in quantum mechanical/molecular mechanical (QM/MM) calculations. In contrast to existing QM/MM FEP variants, the QM region is not kept frozen during sampling, but all degrees of freedom except for the reaction coordinate are sampled. In the DH-FEP scheme, the sampling is done by semiempirical QM/MM molecular dynamics (MD), while the perturbation energy differences are evaluated from high-level QM/MM single-point calculations at regular intervals, skipping a pre-defined number of MD sampling steps. After validating our method using an analytic model potential with an exactly known solution, we report a QM/MM DH-FEP study of the enzymatic reaction catalyzed by chorismate mutase. We suggest guidelines for QM/MM DH-FEP calculations and default values for the required computational parameters. In the case of chorismate mutase, we apply the DH-FEP approach in combination with a single one-dimensional reaction coordinate and with a two-dimensional collective coordinate (two individual distances), with superior results for the latter choice.
Collapse
Affiliation(s)
- Iakov Polyak
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | | | | |
Collapse
|
36
|
García-Meseguer R, Zinovjev K, Roca M, Ruiz-Pernía JJ, Tuñón I. Linking Electrostatic Effects and Protein Motions in Enzymatic Catalysis. A Theoretical Analysis of Catechol O-Methyltransferase. J Phys Chem B 2014; 119:873-82. [DOI: 10.1021/jp505746x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
| | - Kirill Zinovjev
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Maite Roca
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Javier J. Ruiz-Pernía
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
37
|
Silva JRA, Roitberg AE, Alves CN. Catalytic mechanism of L,D-transpeptidase 2 from Mycobacterium tuberculosis described by a computational approach: insights for the design of new antibiotics drugs. J Chem Inf Model 2014; 54:2402-10. [PMID: 25149147 DOI: 10.1021/ci5003069] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Tuberculosis is perhaps the most persistent human disease caused by an infections bacterium, Mycobacterium tuberculosis. The L,D-transpeptidase enzyme catalyzes the formation of 3 → 3 peptidoglycan cross-links of the Mtb cell wall and facilitates resistance against classical β-lactams. Herein, the experimentally proposed mechanism for LdtMt2 was studied by performing QM/MM MD simulations. The whole mechanistic process includes two stages: acylation and deacylation. During the acylation step, two steps were observed: the first step is a thiolate/imidazole ion-pair in the zwitterionic form, and the second step is the nucleophilic attack on the carboxyl carbon of the natural substrate accompanied by the breaking of the peptide bond on substrate. In the deacylation step the acyl-enzyme suffers a nucleophilic attack on the carboxyl carbon by the amine group of the second substrate. Our free energy results obtained by PMF analysis reveal that the first step (acylation) is the rate-limiting step in the whole catalytic mechanism in accordance with the experimental proposal. Also, the residues responsible for binding of the substrate and transition state stabilization were identified by energy decomposition methods.
Collapse
Affiliation(s)
- José Rogério A Silva
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará , Belém, PA 66075-110, Brazil
| | | | | |
Collapse
|
38
|
Świderek K, Martí S, Moliner V. Theoretical Study of Primary Reaction of Pseudozyma antarctica Lipase B as the Starting Point To Understand Its Promiscuity. ACS Catal 2014. [DOI: 10.1021/cs401047k] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Katarzyna Świderek
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924 Lodz, Poland
| | - Sergio Martí
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
39
|
Reis M, Alves CN, Lameira J, Tuñón I, Martí S, Moliner V. The catalytic mechanism of glyceraldehyde 3-phosphate dehydrogenase from Trypanosoma cruzi elucidated via the QM/MM approach. Phys Chem Chem Phys 2013; 15:3772-85. [PMID: 23389436 DOI: 10.1039/c3cp43968b] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been identified as a key enzyme involved in glycolysis processes for energy production in the Trypanosoma cruzi parasite. This enzyme catalyses the oxidative phosphorylation of glyceraldehyde 3-phosphate (G3P) in the presence of inorganic phosphate (Pi) and nicotinamide adenosine dinucleotide (NAD+). The catalytic mechanism used by GAPDH has been intensively investigated. However, the individual roles of Pi and the C3 phosphate of G3P (Ps) sites, as well as some residues such as His194 in the catalytic mechanism, remain unclear. In this study, we have employed Molecular Dynamics (MD) simulations within hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to obtain the Potential of Mean Force of the catalytic oxidative phosphorylation mechanism of the G3P substrate used by GAPDH. According to our results, the first stage of the reaction (oxidoreduction) takes place in the Pi site (energetically more favourable), with the formation of oxyanion thiohemiacetal and thioacylenzyme intermediates without acid-base assistance of His194. Analysis of the interaction energy by residues shows that Arg249 has an important role in the ability of the enzyme to bind the G3P substrate, which interacts with NAD+ and other important residues, such as Cys166, Glu109, Thr167, Ser247 and Thr226, in the GAPDH active site. Finally, the inhibition mechanism of the GAPDH enzyme by the 3-(p-nitrophenoxycarboxyl)-3-ethylene propyl dihydroxyphosphonate inhibitor was investigated in order to contribute to the design of new inhibitors of GAPDH from Trypanosoma cruzi.
Collapse
Affiliation(s)
- Mauro Reis
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110, Belém, PA, Brazil
| | | | | | | | | | | |
Collapse
|
40
|
Zinovjev K, Ruiz-Pernía JJ, Tuñón I. Toward an Automatic Determination of Enzymatic Reaction Mechanisms and Their Activation Free Energies. J Chem Theory Comput 2013; 9:3740-9. [PMID: 26584125 DOI: 10.1021/ct400153r] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present a combination of the string method and a path collective variable for the exploration of the free energy surface associated to a chemical reaction in condensed environments. The on-the-fly string method is employed to find the minimum free energy paths on a multidimensional free energy surface defined in terms of interatomic distances, which is a convenient selection to study bond forming/breaking processes. Once the paths have been determined, a reaction coordinate is defined as a measure of the advance of the system along these paths. This reaction coordinate can be then used to trace the reaction Potential of Mean Force from which the activation free energy can be obtained. This combination of methodologies has been here applied to the study, by means of Quantum Mechanics/Molecular Mechanics simulations, of the reaction catalyzed by guanidinoacetate methyltransferase. This enzyme catalyzes the methylation of guanidinoacetate by S-adenosyl-l-methionine, a reaction that involves a methyl transfer and a proton transfer and for which different reaction mechanisms have been proposed.
Collapse
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain.,Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - J Javier Ruiz-Pernía
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain
| |
Collapse
|
41
|
de Farias Silva N, Lameira J, Alves CN, Martí S. Computational study of the mechanism of half-reactions in class 1A dihydroorotate dehydrogenase from Trypanosoma cruzi. Phys Chem Chem Phys 2013; 15:18863-71. [DOI: 10.1039/c3cp52692e] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
42
|
Roca M, Pascual-Ahuir JL, Tuñón I. Reversibility and Diffusion in Mandelythiamin Decarboxylation. Searching Dynamical Effects in Decarboxylation Reactions. J Am Chem Soc 2012; 134:10509-14. [DOI: 10.1021/ja3026455] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maite Roca
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| | | | - Iñaki Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
43
|
Ferrer S, Ruiz-Pernía J, Martí S, Moliner V, Tuñón I, Bertrán J, Andrés J. Hybrid schemes based on quantum mechanics/molecular mechanics simulations goals to success, problems, and perspectives. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2012; 85:81-142. [PMID: 21920322 DOI: 10.1016/b978-0-12-386485-7.00003-x] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The development of characterization techniques, advanced synthesis methods, as well as molecular modeling has transformed the study of systems in a well-established research field. The current research challenges in biocatalysis and biotransformation evolve around enzyme discovery, design, and optimization. How can we find or create enzymes that catalyze important synthetic reactions, even reactions that may not exist in nature? What is the source of enzyme catalytic power? To answer these and other related questions, the standard strategies have evolved from trial-and-error methodologies based on chemical knowledge, accumulated experience, and common sense into a clearly multidisciplinary science that allows one to reach the molecular design of tailor-made enzyme catalysts. This is even more so when one refers to enzyme catalysts, for which the detailed structure and composition are known and can be manipulated to introduce well-defined residues which can be implicated in the chemical rearrangements taking place in the active site. The methods and techniques of theoretical and computational chemistry are becoming more and more important in both understanding the fundamental biological roles of enzymes and facilitating their utilization in biotechnology. Improvement of the catalytic function of enzymes is important from scientific and industrial viewpoints, and to put this fact in the actual perspective as well as the potentialities, we recommend the very recent report of Sanderson [Sanderson, K. (2011). Chemistry: enzyme expertise. Nature 471, 397.]. Great fundamental advances have been made toward the ab initio design of enzyme catalysts based on molecular modeling. This has been based on the molecular mechanistic knowledge of the reactions to be catalyzed, together with the development of advanced synthesis and characterization techniques. The corresponding molecular mechanism can be studied by means of powerful quantum chemical calculations. The catalytic active site can be optimized to improve the transition state analogues (TSA) and to enhance the catalytic activity, even improve the active site to favor a desired direction of some promiscuous enzymes. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. Current computational tools to assist experimentalists for the design and engineering of proteins with desired catalytic properties are described. The interplay between enzyme design, molecular simulations, and experiments will be presented to emphasize the interdisciplinary nature of this research field. This text highlights the recent advances and examples selected from our laboratory are shown, of how the applications of these tools are a first attempt to de novo design of protein active sites. Identification of neutral/advantageous/deleterious mutation platforms can be exploited to penetrate some of Nature's closely guarded secrets of chemical reactivity. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. The first part describes briefly how the molecular modeling is carried out. Then, we discuss the requirements of hybrid quantum mechanical/molecular mechanics molecular dynamics (QM/MM MD) simulations, analyzing what are the basis of these theoretical methodologies, how we can use them with a view to its application in the study of enzyme catalysis, and what are the best methodologies for assessing its catalytic potential. In the second part, we focus on some selected examples, taking as a common guide the chorismate to prephenate rearrangement, studying the corresponding molecular mechanism in vacuo, in solution and in an enzyme environment. In addition, examples involving catalytic antibodies (CAs) and promiscuous enzymes will be presented. Finally, a special emphasis is made to provide some hints about the logical evolution that can be anticipated in this research field. Moreover, it helps in understanding the open directions in this area of knowledge and highlights the importance of computational approaches in discovering specific drugs and the impact on the rational design of tailor-made enzymes.
Collapse
Affiliation(s)
- Silvia Ferrer
- Departamento de Química Física y Analítica, Universitat Jaume I, Castellón, Spain
| | | | | | | | | | | | | |
Collapse
|
44
|
Lameira J, Alves CN, Tuñón I, Martí S, Moliner V. Enzyme Molecular Mechanism as a Starting Point to Design New Inhibitors: A Theoretical Study of O-GlcNAcase. J Phys Chem B 2011; 115:6764-75. [DOI: 10.1021/jp202079e] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jerônimo Lameira
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110 Belém, PA, Brazil
| | - Cláudio Nahum Alves
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110 Belém, PA, Brazil
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, València, Spain
| | - Sergio Martí
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
45
|
Ruiz Pernía JJ, Tuñón I, Williams IH. Computational Simulation of the Lifetime of the Methoxymethyl Cation in Water. A Simple Model for a Glycosyl Cation: When Is an Intermediate an Intermediate? J Phys Chem B 2010; 114:5769-74. [DOI: 10.1021/jp910539j] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- J. Javier Ruiz Pernía
- Department of Chemistry, University of Bath, Bath, BA2 7AY, U.K., Departamento de Química, Universidad de Valencia, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Department of Chemistry, University of Bath, Bath, BA2 7AY, U.K., Departamento de Química, Universidad de Valencia, 46100 Burjassot, Spain
| | - Ian H. Williams
- Department of Chemistry, University of Bath, Bath, BA2 7AY, U.K., Departamento de Química, Universidad de Valencia, 46100 Burjassot, Spain
| |
Collapse
|
46
|
Theoretical QM/MM studies of enzymatic pericyclic reactions. Interdiscip Sci 2010; 2:115-31. [DOI: 10.1007/s12539-010-0095-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2009] [Revised: 12/07/2009] [Accepted: 12/09/2009] [Indexed: 11/25/2022]
|
47
|
Ranaghan KE, Mulholland AJ. Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods. INT REV PHYS CHEM 2010. [DOI: 10.1080/01442350903495417] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
48
|
Tishchenko O, Truhlar DG. Non-Hermitian Multiconfiguration Molecular Mechanics. J Chem Theory Comput 2009; 5:1454-61. [DOI: 10.1021/ct900077g] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Oksana Tishchenko
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
49
|
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
|
50
|
A QM/MM study of the reaction mechanism for the 3′-processing step catalyzed by HIV-1 integrase. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2008.08.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|