1
|
Corbella M, Bravo J, Demkiv AO, Calixto AR, Sompiyachoke K, Bergonzi C, Elias MH, Kamerlin SCL. Catalytic Redundancies and Conformational Plasticity Drives Selectivity and Promiscuity in Quorum Quenching Lactonases. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.01.592096. [PMID: 38746346 PMCID: PMC11092605 DOI: 10.1101/2024.05.01.592096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
Several enzymes from the metallo-β-lactamase-like family of lactonases (MLLs) degrade N- acyl-L-homoserine lactones (AHLs). In doing so, they play a role in a microbial communication system, quorum sensing, which contributes to pathogenicity and biofilm formation. There is currently great interest in designing quorum quenching ( QQ ) enzymes that can interfere with this communication and be used in a range of industrial and biomedical applications. However, tailoring these enzymes for specific targets requires a thorough understanding of their mechanisms and the physicochemical properties that determine their substrate specificities. We present here a detailed biochemical, computational, and structural study of the MLL GcL, which is highly proficient, thermostable, and has broad substrate specificity. Strikingly, we show that GcL does not only accept a broad range of substrates but is also capable of utilizing different reaction mechanisms that are differentially used in function of the substrate structure or the remodeling of the active site via mutations. Comparison of GcL to other lactonases such as AiiA and AaL demonstrates similar mechanistic promiscuity, suggesting this is a shared feature across lactonases in this enzyme family. Mechanistic promiscuity has previously been observed in the lactonase/paraoxonase PON1, as well as with protein tyrosine phosphatases that operate via a dual general-acid mechanism. The apparent prevalence of this phenomenon is significant from both a biochemical and an engineering perspective: in addition to optimizing for specific substrates, it is possible to optimize for specific mechanisms, opening new doors not just for the design of novel quorum quenching enzymes, but also of other mechanistically promiscuous enzymes.
Collapse
|
2
|
Allen AEA, Csányi G. Toward transferable empirical valence bonds: Making classical force fields reactive. J Chem Phys 2024; 160:124108. [PMID: 38526105 DOI: 10.1063/5.0196952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 03/07/2024] [Indexed: 03/26/2024] Open
Abstract
The empirical valence bond technique allows classical force fields to model reactive processes. However, parametrization from experimental data or quantum mechanical calculations is required for each reaction present in the simulation. We show that the parameters present in the empirical valence bond method can be predicted using a neural network model and the SMILES strings describing a reaction. This removes the need for quantum calculations in the parametrization of the empirical valence bond technique. In doing so, we have taken the first steps toward defining a new procedure for enabling reactive atomistic simulations. This procedure would allow researchers to use existing classical force fields for reactive simulations, without performing additional quantum mechanical calculations.
Collapse
Affiliation(s)
- Alice E A Allen
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
3
|
Oanca G, van der Ent F, Åqvist J. Efficient Empirical Valence Bond Simulations with GROMACS. J Chem Theory Comput 2023; 19:6037-6045. [PMID: 37623818 PMCID: PMC10500987 DOI: 10.1021/acs.jctc.3c00714] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Indexed: 08/26/2023]
Abstract
We describe a protocol to perform empirical valence bond (EVB) simulations using GROMACS software. EVB is a fast and reliable method that allows one to determine the reaction free-energy profiles in complex systems, such as enzymes, by employing classical force fields to represent a chemical reaction. Therefore, running EVB simulations is basically as fast as any classical molecular dynamics simulation, and the method uses standard free-energy calculations to map the free-energy change along a given reaction path. To exemplify and validate our EVB implementation, we replicated two cases of our earlier enzyme simulations. One of these addresses the decomposition of the activation free energy into its enthalpic and entropic components, and the other is focused on calculating the overall catalytic effect of the enzyme compared to the same reaction in water. These two examples give virtually identical results to those obtained with programs that were specifically designed for EVB simulations and show that the GROMACS implementation is robust and can be used for very large systems.
Collapse
Affiliation(s)
- Gabriel Oanca
- Department of Cell and Molecular Biology,
Biomedical Center, Uppsala University, Uppsala SE-751 24, Sweden
| | - Florian van der Ent
- Department of Cell and Molecular Biology,
Biomedical Center, Uppsala University, Uppsala SE-751 24, Sweden
| | - Johan Åqvist
- Department of Cell and Molecular Biology,
Biomedical Center, Uppsala University, Uppsala SE-751 24, Sweden
| |
Collapse
|
4
|
Zhang J, Tang Z, Zhang X, Zhu H, Zhao R, Lu Y, Gao J. Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States. J Chem Theory Comput 2023; 19:1777-1789. [PMID: 36917687 DOI: 10.1021/acs.jctc.2c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Zhen Tang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Xiaoyong Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Hong Zhu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China
| | - Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,Institute of Theoretical Chemistry, Jilin University, Changchun, 130023 Jilin, P. R. China
| | - Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
5
|
Hong G, Pachter R. Effects of inter-radical interactions and scavenging radicals on magnetosensitivity: spin dynamics simulations of proposed radical pairs. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2023; 52:27-37. [PMID: 36792823 DOI: 10.1007/s00249-023-01630-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/07/2022] [Revised: 01/22/2023] [Accepted: 01/23/2023] [Indexed: 02/17/2023]
Abstract
Although the magnetosensitivity to weak magnetic fields, such as the geomagnetic field, which was exhibited by radical pairs that are potentially responsible for avian navigation, has been previously investigated by spin dynamics simulations, understanding this behavior for proposed radical pairs in other species is limited. These include, for example, radical pairs formed in the single-cell green alga Chlamydomonas reinhardtii (CraCRY) and in Columba livia (ClCRY4). In addition, the radical pair of FADH• with the one-electron reduced cyclobutane thymine dimer that was shown to be sensitive to weak magnetic fields has been of interest. In this work, we investigated the directional magnetosensitivity of these radical pairs to a weak magnetic field by spin dynamics simulations. We find significant reduction in the magnetosensitivity by inclusion of dipolar and exchange interactions, which can be mitigated by a scavenging radical, as demonstrated for the [FAD•- TyrD•] radical pair in CraCRY, but not for the [FADH• T□T•-] radical pair because of the large exchange coupling. The directional magnetosensitivity of the ClCRY4 [FAD•- TyrE•] radical pair can survive this adverse effect even without the scavenging reaction, possibly motivating further experimental exploration.
Collapse
Affiliation(s)
- Gongyi Hong
- Air Force Research Laboratory, Materials and Manufacturing Directorate, Wright-Patterson Air Force Base, Ohio, 45433, USA
| | - Ruth Pachter
- Air Force Research Laboratory, Materials and Manufacturing Directorate, Wright-Patterson Air Force Base, Ohio, 45433, USA.
| |
Collapse
|
6
|
Menon A, Pascazio L, Nurkowski D, Farazi F, Mosbach S, Akroyd J, Kraft M. OntoPESScan: An Ontology for Potential Energy Surface Scans. ACS OMEGA 2023; 8:2462-2475. [PMID: 36687109 PMCID: PMC9850739 DOI: 10.1021/acsomega.2c06948] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Accepted: 11/30/2022] [Indexed: 06/17/2023]
Abstract
In this work, a new OntoPESScan ontology is developed for the semantic representation of one-dimensional potential energy surface (PES) scans, a central concept in computational chemistry. This ontology is developed in line with knowledge graph principles and The World Avatar (TWA) project. OntoPESScan is linked to other ontologies for chemistry in TWA, including OntoSpecies, which helps uniquely identify species along the PES and access their properties, and OntoCompChem, which allows the association of potential energy surfaces with quantum chemical calculations and the concepts used to derive them. A force-field fitting agent is also developed that makes use of the information in the OntoPESScan ontology to fit force fields to reactive surfaces of interest on the fly by making use of the empirical valence bond methodology. This agent is demonstrated to successfully parametrize two cases, namely, a PES scan on ethanol and a PES scan on a localized π-radical PAH hypothesized to play a role in soot formation during combustion. OntoPESScan is an extension to the capabilities of TWA and, in conjunction with potential further ontological support for molecular dynamics and reactions, will further progress toward an open, continuous, and self-growing knowledge graph for chemistry.
Collapse
Affiliation(s)
- Angiras Menon
- Department
of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, U.K.
| | - Laura Pascazio
- CARES, Cambridge Centre for Advanced Research and Education
in Singapore, 1 Create
Way, CREATE Tower, #05-05, Singapore 138602
| | - Daniel Nurkowski
- CMCL
Innovations, Sheraton House, Castle Park, Cambridge CB3 0AX, U.K.
| | - Feroz Farazi
- Department
of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, U.K.
| | - Sebastian Mosbach
- Department
of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, U.K.
- CARES, Cambridge Centre for Advanced Research and Education
in Singapore, 1 Create
Way, CREATE Tower, #05-05, Singapore 138602
| | - Jethro Akroyd
- Department
of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, U.K.
- CARES, Cambridge Centre for Advanced Research and Education
in Singapore, 1 Create
Way, CREATE Tower, #05-05, Singapore 138602
| | - Markus Kraft
- Department
of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, U.K.
- CARES, Cambridge Centre for Advanced Research and Education
in Singapore, 1 Create
Way, CREATE Tower, #05-05, Singapore 138602
- School
of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, Singapore 637459
- The
Alan Turing Institute, London NW1 2BD, United
Kingdom
| |
Collapse
|
7
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
8
|
Zotti LA, Dednam W, Lombardi EB, Palacios JJ. Constrained DFT for Molecular Junctions. NANOMATERIALS 2022; 12:nano12071234. [PMID: 35407352 PMCID: PMC9002544 DOI: 10.3390/nano12071234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Revised: 03/31/2022] [Accepted: 04/04/2022] [Indexed: 02/01/2023]
Abstract
We have explored the use of constrained density functional theory (cDFT) for molecular junctions based on benzenediamine. By elongating the junction, we observe that the energy gap between the ionization potential and the electronic affinity increases with the stretching distance. This is consistent with the trend expected from the electrostatic screening. A more detailed analysis shows how this influences the charge distribution of both the individual metal layers and the molecular atoms. Overall, our work shows that constrained DFT is a powerful tool for studying screening effects in molecular junctions.
Collapse
Affiliation(s)
- Linda Angela Zotti
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain;
- Correspondence:
| | - Wynand Dednam
- Department of Physics, Science Campus, University of South Africa, Private Bag X6, Florida Park 1710, South Africa; (W.D.); (E.B.L.)
| | - Enrico B. Lombardi
- Department of Physics, Science Campus, University of South Africa, Private Bag X6, Florida Park 1710, South Africa; (W.D.); (E.B.L.)
| | - Juan Jose Palacios
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain;
- Departamento de Física de la Materia Condensada and Instituto Nicolás Cabrera (INC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| |
Collapse
|
9
|
Why Monoamine Oxidase B Preferably Metabolizes N-Methylhistamine over Histamine: Evidence from the Multiscale Simulation of the Rate-Limiting Step. Int J Mol Sci 2022; 23:ijms23031910. [PMID: 35163835 PMCID: PMC8836602 DOI: 10.3390/ijms23031910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 02/04/2022] [Accepted: 02/05/2022] [Indexed: 11/18/2022] Open
Abstract
Histamine levels in the human brain are controlled by rather peculiar metabolic pathways. In the first step, histamine is enzymatically methylated at its imidazole Nτ atom, and the produced N-methylhistamine undergoes an oxidative deamination catalyzed by monoamine oxidase B (MAO-B), as is common with other monoaminergic neurotransmitters and neuromodulators of the central nervous system. The fact that histamine requires such a conversion prior to oxidative deamination is intriguing since MAO-B is known to be relatively promiscuous towards monoaminergic substrates; its in-vitro oxidation of N-methylhistamine is about 10 times faster than that for histamine, yet this rather subtle difference appears to be governing the decomposition pathway. This work clarifies the MAO-B selectivity toward histamine and N-methylhistamine by multiscale simulations of the rate-limiting hydride abstraction step for both compounds in the gas phase, in aqueous solution, and in the enzyme, using the established empirical valence bond methodology, assisted by gas-phase density functional theory (DFT) calculations. The computed barriers are in very good agreement with experimental kinetic data, especially for relative trends among systems, thereby reproducing the observed MAO-B selectivity. Simulations clearly demonstrate that solvation effects govern the reactivity, both in aqueous solution as well as in the enzyme although with an opposing effect on the free energy barrier. In the aqueous solution, the transition-state structure involving histamine is better solvated than its methylated analog, leading to a lower barrier for histamine oxidation. In the enzyme, the higher hydrophobicity of N-methylhistamine results in a decreased number of water molecules at the active side, leading to decreased dielectric shielding of the preorganized catalytic electrostatic environment provided by the enzyme. This renders the catalytic environment more efficient for N-methylhistamine, giving rise to a lower barrier relative to histamine. In addition, the transition state involving N-methylhistamine appears to be stabilized by the surrounding nonpolar residues to a larger extent than with unsubstituted histamine, contributing to a lower barrier with the former.
Collapse
|
10
|
Bensberg M, Türtscher PL, Unsleber JP, Reiher M, Neugebauer J. Solvation Free Energies in Subsystem Density Functional Theory. J Chem Theory Comput 2022; 18:723-740. [PMID: 34985890 DOI: 10.1021/acs.jctc.1c00864] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise the scalability of the model and transferability of the predicted solvent effect, we aim to retain both, for different solutes as well as for different solvents. The key for the transferability is the consistent subsystem decomposition of solute and solvent. The key for the scalability is the performance of subsystem DFT for increasing numbers of subsystems. We investigate molecular dynamics and stationary point sampling of solvent configurations and compare the resulting (Gibbs) free energies to experiment and theoretical methods. We can show that with our hybrid model reaction barriers and reaction energies are accurately reproduced compared to experimental data.
Collapse
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Paul L Türtscher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jan P Unsleber
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
11
|
Stoppelman JP, McDaniel JG. Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM +][OAc -]. J Chem Phys 2021; 155:104112. [PMID: 34525833 DOI: 10.1063/5.0063187] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Reactive molecular dynamics simulations enable a detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the empirical valence bond framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NNs). The physics-based terms ensure the correct asymptotic behavior of electrostatic, polarization, and dispersion interactions and are compatible with existing solvent force fields. NNs are utilized for a versatile description of short-range orbital interactions within the transition state region and accurate rendering of vibrational motion of the reacting complex. We demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium cation with acetate to form 1-ethyl-3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ∼1 kJ mol-1 mean absolute error accuracy within the transition state region for the gas-phase complex. To characterize the solvent modulation of the reaction profile, we compute potentials of mean force for the gas-phase reaction as well as the reaction within a four-ion cluster and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes the formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.
Collapse
Affiliation(s)
- John P Stoppelman
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Jesse G McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
12
|
Joutsuka T, Ando K. Efficient Free-energy Calculation of Proton Transfer by Constrained Density Functional Theory and Geometrically Restrained Molecular Dynamics Simulation. CHEM LETT 2021. [DOI: 10.1246/cl.210132] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Tatsuya Joutsuka
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, Hitachi, Ibaraki 316-8511, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Woman’s Christian University, 2-6-1 Zenpukuji, Suginami-ku, Tokyo 167-8585, Japan
| |
Collapse
|
13
|
Artiukhin DG, Eschenbach P, Matysik J, Neugebauer J. Theoretical Assessment of Hinge-Type Models for Electron Donors in Reaction Centers of Photosystems I and II as well as of Purple Bacteria. J Phys Chem B 2021; 125:3066-3079. [PMID: 33749260 DOI: 10.1021/acs.jpcb.0c10656] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hinge-type molecular models for electron donors in reaction centers of Photosystems I and II and purple bacteria were investigated using a two-state computational approach based on frozen-density embedding (FDE). This methodology, dubbed FDE-diab, is known to avoid consequences of the self-interaction error as far as intermolecular phenomena are concerned, which allows a prediction of qualitatively correct spin densities for large biomolecular systems. The calculated spin density distributions are in a good agreement with available experimental results and demonstrated a very high sensitivity to changes in the relative orientation of cofactors and amino acid protonation states. This allows a validation of the previously proposed hinge-type models providing hints on possible protonation states of axial histidine molecules.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Jörg Matysik
- Institut für Analytische Chemie, Universität Leipzig, Linnéstr. 3, 04103 Leipzig, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
14
|
Biedermann M, Diddens D, Heuer A. rs@md: Introducing Reactive Steps at the Molecular Dynamics Simulation Level. J Chem Theory Comput 2021; 17:1074-1085. [DOI: 10.1021/acs.jctc.0c01189] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Myra Biedermann
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
| | - Diddo Diddens
- Helmholtz-Institute Münster Ionics in Energy Storage (IEK-12), Forschungszentrum Jülich GmbH, Corrensstrasse 46, 48149 Münster, Germany
| | - Andreas Heuer
- Institute of Physical Chemistry, University of Münster, Corrensstraße 28/30, 48149 Münster, Germany
- Helmholtz-Institute Münster Ionics in Energy Storage (IEK-12), Forschungszentrum Jülich GmbH, Corrensstrasse 46, 48149 Münster, Germany
| |
Collapse
|
15
|
Hong G, Pachter R, Essen LO, Ritz T. Electron transfer and spin dynamics of the radical-pair in the cryptochrome from Chlamydomonas reinhardtii by computational analysis. J Chem Phys 2020; 152:065101. [PMID: 32061221 DOI: 10.1063/1.5133019] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
In an effort to elucidate the origin of avian magnetoreception, it was postulated that a radical-pair formed in a cryptochrome upon light activation provided the basis for the mechanism that enables an inclination compass sensitive to the geomagnetic field. Photoreduction in this case involves formation of a flavin adenine dinucleotide (FAD)-tryptophan (TRP) radical-pair, following electron transfer within a conserved TRP triad in the cryptochrome. Recently, an animal-like cryptochrome from Chlamydomonas reinhardtii (CraCRY) was analyzed, demonstrating the role of a fourth aromatic residue, which serves as a terminal electron donor in the photoreduction pathway, resulting in the creation of a more distal radical-pair and exhibiting fast electron transfer. In this work, we investigated the electron transfer in CraCRY with a combination of free energy molecular dynamics (MD) simulations, frozen density functional theory, and QM/MM MD simulations, supporting the suggestion of a proton coupled electron transfer mechanism. Spin dynamics simulations discerned details on the dependence of the singlet yield on the direction of the external magnetic field for the [FAD•- TYRH•+] and [FAD•- TYR•] radical-pairs in CraCRY, in comparison with the previously modeled [FAD•- TRPH•+] radical-pair.
Collapse
Affiliation(s)
- Gongyi Hong
- Air Force Research Laboratory, Materials and Manufacturing Directorate, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Ruth Pachter
- Air Force Research Laboratory, Materials and Manufacturing Directorate, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Lars-Oliver Essen
- Department of Chemistry, Center for Synthetic Microbiology, Philipps University, Marburg 35032, Germany
| | - Thorsten Ritz
- Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| |
Collapse
|
16
|
Joutsuka T, Ando K. Constrained Density Functional Theory Molecular Dynamics Simulation of Deprotonation in Aqueous Silicic Acid. J Phys Chem B 2020; 124:8323-8330. [PMID: 32856917 DOI: 10.1021/acs.jpcb.0c05096] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The solubility of silicic acid in water is so low that the molecular mechanism behind the physical properties such as pKa remains poorly understood, despite the importance in fields such as chemistry and geology. Theoretical calculations provide the molecular-level information on such a rare chemical species yet face difficulties in the selection of reaction coordinates and the rare-event sampling of proton transfers (PTs) in condensed phases. We propose a straightforward calculation scheme of pKa with ab initio molecular dynamics (MD) simulation and the constrained density functional theory (CDFT), which provides structural and dynamical properties such as radial distribution functions, vibrational spectra, and reaction paths. The calculated reaction free energies of deprotonations agreed with experiments within a few kcal/mol. Analysis of the solvation structure shows that, after deprotonation, the hydronium ion D3O+ repels away from the deprotonated silicic acid SiO(OD)3- without forming a stable contact-ion pair. The calculated vibrational spectra are consistent with the spectroscopic measurements, and the dynamical analysis of the reaction path quantifies the couplings of the OD stretch vibrations of silicic acid and water to the PT reaction in terms of the vertical energy-gap coordinate that includes both the solute and the solvent components.
Collapse
Affiliation(s)
- Tatsuya Joutsuka
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, Hitachi 316-8511 Ibaraki, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Womans Christian University, 2-6-1 Zenpukuji, Suginami-ku, Tokyo 167-8585, Japan
| |
Collapse
|
17
|
Scivetti I, Sen K, Elena AM, Todorov I. Reactive Molecular Dynamics at Constant Pressure via Nonreactive Force Fields: Extending the Empirical Valence Bond Method to the Isothermal-Isobaric Ensemble. J Phys Chem A 2020; 124:7585-7597. [PMID: 32820921 DOI: 10.1021/acs.jpca.0c05461] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Empirical Valence Bond (EVB) method offers a suitable framework to obtain reactive potentials through the coupling of nonreactive force fields. In this formalism, most of the implemented coupling terms are built using functional forms that depend on spatial coordinates, while parameters are fitted against reference data to model the change of chemistry between the participating nonreactive states. In this work, we demonstrate that the use of such coupling terms precludes the computation of the stress tensor for condensed phase systems and prevents the possibility to carry out EVB molecular dynamics in the isothermal-isobaric (NPT) ensemble. Alternatively, we make use of coupling terms that depend on the energy gaps, defined as the energy differences between the participating nonreactive force fields, and derive a general expression for the EVB stress tensor suitable for computation. Implementation of this new methodology is tested for a model of a single reactive malonaldehyde solvated in nonreactive water. Mass densities and probability distributions for the values of the energy gaps computed in the NPT ensemble reveal a negligible role of the reactive potential in the limit of low concentrated solutions, thus corroborating for the first time the validity of approximations based on the canonical NVT ensemble, customarily adopted for EVB simulations. The presented formalism also aims to contribute to future implementations and extensions of the EVB method to research the limit of highly concentrated solutions.
Collapse
Affiliation(s)
- Ivan Scivetti
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K.,Department of Chemistry, University of Liverpool, Liverpool L69 3BX, U.K
| | - Kakali Sen
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Alin M Elena
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Ilian Todorov
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| |
Collapse
|
18
|
Hydride Abstraction as the Rate-Limiting Step of the Irreversible Inhibition of Monoamine Oxidase B by Rasagiline and Selegiline: A Computational Empirical Valence Bond Study. Int J Mol Sci 2020; 21:ijms21176151. [PMID: 32858935 PMCID: PMC7503497 DOI: 10.3390/ijms21176151] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Revised: 08/11/2020] [Accepted: 08/21/2020] [Indexed: 12/23/2022] Open
Abstract
Monoamine oxidases (MAOs) catalyze the degradation of a very broad range of biogenic and dietary amines including many neurotransmitters in the brain, whose imbalance is extensively linked with the biochemical pathology of various neurological disorders, and are, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. Still, despite this practical significance, the precise molecular mechanism underlying the irreversible MAO inhibition with clinically used propargylamine inhibitors rasagiline and selegiline is still not unambiguously determined, which hinders the rational design of improved inhibitors devoid of side effects current drugs are experiencing. To address this challenge, we present empirical valence bond QM/MM simulations of the rate-limiting step of the MAO inhibition involving the hydride anion transfer from the inhibitor α-carbon onto the N5 atom of the flavin adenin dinucleotide (FAD) cofactor. The proposed mechanism is strongly supported by the obtained free energy profiles, which confirm a higher reactivity of selegiline over rasagiline, while the calculated difference in the activation Gibbs energies of ΔΔG‡ = 3.1 kcal mol-1 is found to be in very good agreement with that from the measured literature kinact values that predict a 1.7 kcal mol-1 higher selegiline reactivity. Given the similarity with the hydride transfer mechanism during the MAO catalytic activity, these results verify that both rasagiline and selegiline are mechanism-based irreversible inhibitors and offer guidelines in designing new and improved inhibitors, which are all clinically employed in treating a variety of neuropsychiatric and neurodegenerative conditions.
Collapse
|
19
|
Koner D, Salehi SM, Mondal P, Meuwly M. Non-conventional force fields for applications in spectroscopy and chemical
reaction dynamics. J Chem Phys 2020; 153:010901. [DOI: 10.1063/5.0009628] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Seyedeh Maryam Salehi
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Padmabati Mondal
- Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra
Pradesh, India
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel,
Switzerland and Department of Chemistry, Brown University, Providence, Rhode Island, USA
| |
Collapse
|
20
|
Ortega DE, Matute RA. Description of the Reaction Intermediate Stabilization for the Zimmerman Di-π-methane Rearrangement on the Basis of a Parametric Diabatic Analysis. J Phys Chem A 2020; 124:3573-3580. [DOI: 10.1021/acs.jpca.0c01145] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Daniela E. Ortega
- Centro Integrativo de Biología y Química Aplicada (CIBQA), Universidad Bernardo OHiggins, Santiago 8370854, Chile
| | - Ricardo A. Matute
- Centro Integrativo de Biología y Química Aplicada (CIBQA), Universidad Bernardo OHiggins, Santiago 8370854, Chile
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
21
|
Ortega DE, Ormazábal-Toledo R, Contreras R, Matute RA. Theoretical insights into the E1cB/E2 mechanistic dichotomy of elimination reactions. Org Biomol Chem 2019; 17:9874-9882. [DOI: 10.1039/c9ob02004g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Modeling the E1cB, E2, and E1cB/E2 borderline mechanisms in terms of carbanion stabilization.
Collapse
Affiliation(s)
- Daniela E. Ortega
- Centro Integrativo de Biología y Química Aplicada (CIBQA)
- Universidad Bernardo O Higgins
- Santiago 8370854
- Chile
| | - Rodrigo Ormazábal-Toledo
- Centro Integrativo de Biología y Química Aplicada (CIBQA)
- Universidad Bernardo O Higgins
- Santiago 8370854
- Chile
| | - Renato Contreras
- Departamento de Química
- Facultad de Ciencias
- Universidad de Chile
- Santiago
- Chile
| | - Ricardo A. Matute
- Centro Integrativo de Biología y Química Aplicada (CIBQA)
- Universidad Bernardo O Higgins
- Santiago 8370854
- Chile
- Division of Chemistry and Chemical Engineering
| |
Collapse
|
22
|
Kulkarni Y, Kamerlin SCL. Computational physical organic chemistry using the empirical valence bond approach. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.07.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
23
|
Meuwly M. Reactive molecular dynamics: From small molecules to proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1386] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Markus Meuwly
- Department of Chemistry University of Basel Basel Switzerland
- Department of Chemistry Brown University Providence Rhode Island
| |
Collapse
|
24
|
Mosquera MA, Ratner MA, Schatz GC. Locally coupled open subsystems: A formalism for affordable electronic structure calculations featuring fractional charges and size consistency. J Chem Phys 2018; 149:034105. [DOI: 10.1063/1.5038557] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Martín A. Mosquera
- Department of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, USA
| | - Mark A. Ratner
- Department of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, USA
| | - George C. Schatz
- Department of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, USA
| |
Collapse
|
25
|
Hong G, Pachter R, Ritz T. Coupling Drosophila melanogaster Cryptochrome Light Activation and Oxidation of the Kvβ Subunit Hyperkinetic NADPH Cofactor. J Phys Chem B 2018; 122:6503-6510. [PMID: 29847128 DOI: 10.1021/acs.jpcb.8b03493] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Motivated by the observations on the involvement of light-induced processes in the Drosophila melanogaster cryptochrome (DmCry) in regulation of the neuronal firing rate, which is achieved by a redox-state change of its voltage-dependent K+ channel Kvβ subunit hyperkinetic (Hk) reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor, we propose in this work two hypothetical pathways that may potentially enable such coupling. In the first pathway, triggered by blue-light-induced formation of a radical pair [FAD•-TRP•+] in DmCry, the hole (TRP•+) may hop to Hk, for example, through a tryptophan chain and oxidize NADPH, possibly leading to inhibition of the N-terminus inactivation in the K+ channel. In a second possible pathway, DmCry's FAD•- is reoxidized by molecular oxygen, producing H2O2, which then diffuses to Hk and oxidizes NADPH. In this work, by applying a combination of quantum and empirical-based methods for free-energy calculations, we find that the oxidation of NADPH by TRP•+ or H2O2 and the reoxidation of FAD•- by O2 are thermodynamically feasible. Our results may have an implication in identifying a magnetic sensing signal transduction pathway, specifically upon Drosophila's Hk NADPH cofactor oxidation, with a subsequent inhibition of the K+ channel N-terminus inactivation gate, permitting K+ flux.
Collapse
Affiliation(s)
- Gongyi Hong
- Air Force Research Laboratory , Wright-Patterson Air Force Base , Dayton , Ohio 45433-7702 , United States
| | - Ruth Pachter
- Air Force Research Laboratory , Wright-Patterson Air Force Base , Dayton , Ohio 45433-7702 , United States
| | - Thorsten Ritz
- Department of Physics and Astronomy , University of California , Irvine , California 92697-4575 , United States
| |
Collapse
|
26
|
Artiukhin DG, Neugebauer J. Frozen-density embedding as a quasi-diabatization tool: Charge-localized states for spin-density calculations. J Chem Phys 2018; 148:214104. [DOI: 10.1063/1.5023290] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Denis G. Artiukhin
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
27
|
Maršavelski A, Petrović D, Bauer P, Vianello R, Kamerlin SCL. Empirical Valence Bond Simulations Suggest a Direct Hydride Transfer Mechanism for Human Diamine Oxidase. ACS OMEGA 2018; 3:3665-3674. [PMID: 30023875 PMCID: PMC6044848 DOI: 10.1021/acsomega.8b00346] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/26/2018] [Accepted: 03/19/2018] [Indexed: 06/08/2023]
Abstract
Diamine oxidase (DAO) is an enzyme involved in the regulation of cell proliferation and the immune response. This enzyme performs oxidative deamination in the catabolism of biogenic amines, including, among others, histamine, putrescine, spermidine, and spermine. The mechanistic details underlying the reductive half-reaction of the DAO-catalyzed oxidative deamination which leads to the reduced enzyme cofactor and the aldehyde product are, however, still under debate. The catalytic mechanism was proposed to involve a prototropic shift from the substrate-Schiff base to the product-Schiff base, which includes the rate-limiting cleavage of the Cα-H bond by the conserved catalytic aspartate. Our detailed mechanistic study, performed using a combined quantum chemical cluster approach with empirical valence bond simulations, suggests that the rate-limiting cleavage of the Cα-H bond involves direct hydride transfer to the topaquinone cofactor-a mechanism that does not involve the formation of a Schiff base. Additional investigation of the D373E and D373N variants supported the hypothesis that the conserved catalytic aspartate is indeed essential for the reaction; however, it does not appear to serve as the catalytic base, as previously suggested. Rather, the electrostatic contributions of the most significant residues (including D373), together with the proximity of the Cu2+ cation to the reaction site, lower the activation barrier to drive the chemical reaction.
Collapse
Affiliation(s)
- Aleksandra Maršavelski
- Computational
Organic Chemistry and Biochemistry Group, Division of Organic Chemistry
and Biochemistry, Ruđer Bošković
Institute, Bijenička
cesta 54, 10000 Zagreb, Croatia
- Department
of Chemistry, Faculty of Science, University
of Zagreb, Horvatovac
102a, 10000 Zagreb, Croatia
- Department
of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Dušan Petrović
- Department
of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Paul Bauer
- Department
of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
- Department
of Biophysics, SciLifeLab, KTH Royal Institute
of Technology, S-10691 Stockholm, Sweden
| | - Robert Vianello
- Computational
Organic Chemistry and Biochemistry Group, Division of Organic Chemistry
and Biochemistry, Ruđer Bošković
Institute, Bijenička
cesta 54, 10000 Zagreb, Croatia
| | | |
Collapse
|
28
|
Barrozo A, Liao Q, Esguerra M, Marloie G, Florián J, Williams NH, Kamerlin SCL. Computer simulations of the catalytic mechanism of wild-type and mutant β-phosphoglucomutase. Org Biomol Chem 2018; 16:2060-2073. [DOI: 10.1039/c8ob00312b] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
β-Phosphoglucomutase (β-PGM) has served as an important model system for understanding biological phosphoryl transfer.
Collapse
Affiliation(s)
- Alexandre Barrozo
- Science for Life Laboratory
- Department of Cell and Molecular Biology
- Uppsala University
- S-75124 Uppsala
- Sweden
| | - Qinghua Liao
- Science for Life Laboratory
- Department of Cell and Molecular Biology
- Uppsala University
- S-75124 Uppsala
- Sweden
| | - Mauricio Esguerra
- Science for Life Laboratory
- Department of Cell and Molecular Biology
- Uppsala University
- S-75124 Uppsala
- Sweden
| | - Gaël Marloie
- Science for Life Laboratory
- Department of Cell and Molecular Biology
- Uppsala University
- S-75124 Uppsala
- Sweden
| | - Jan Florián
- Department of Chemistry and Biochemistry
- Loyola University Chicago
- Chicago
- USA
| | | | | |
Collapse
|
29
|
Purg M, Kamerlin SCL. Empirical Valence Bond Simulations of Organophosphate Hydrolysis: Theory and Practice. Methods Enzymol 2018; 607:3-51. [DOI: 10.1016/bs.mie.2018.06.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
30
|
Antipov SV, Bhattacharyya S, El Hage K, Xu ZH, Meuwly M, Rothlisberger U, Vaníček J. Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061509. [PMID: 29376107 PMCID: PMC5758379 DOI: 10.1063/1.4996559] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 10/11/2017] [Indexed: 06/07/2023]
Abstract
Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.
Collapse
Affiliation(s)
- Sergey V Antipov
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Swarnendu Bhattacharyya
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Krystel El Hage
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Oanca G, Stare J, Mavri J. How fast monoamine oxidases decompose adrenaline? Kinetics of isoenzymes A and B evaluated by empirical valence bond simulation. Proteins 2017; 85:2170-2178. [PMID: 28836294 DOI: 10.1002/prot.25374] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Revised: 08/18/2017] [Accepted: 08/21/2017] [Indexed: 12/22/2022]
Abstract
This work scrutinizes kinetics of decomposition of adrenaline catalyzed by monoamine oxidase (MAO) A and B enzymes, a process controlling the levels of adrenaline in the central nervous system and other tissues. Experimental kinetic data for MAO A and B catalyzed decomposition of adrenaline are reported only in the form of the maximum reaction rate. Therefore, we estimated the experimental free energy barriers form the kinetic data of closely related systems using regression method, as was done in our previous study. By using multiscale simulation on the Empirical Valence Bond (EVB) level, we studied the chemical reactivity of the MAO A catalyzed decomposition of adrenaline and we obtained a value of activation free energy of 17.3 ± 0.4 kcal/mol. The corresponding value for MAO B is 15.7 ± 0.7 kcal/mol. Both values are in good agreement with the estimated experimental barriers of 16.6 and 16.0 kcal/mol for MAO A and MAO B, respectively. The fact that we reproduced the kinetic data and preferential catalytic effect of MAO B over MAO A gives additional support to the validity of the proposed hydride transfer mechanism. Furthermore, we demonstrate that adrenaline is preferably involved in the reaction in a neutral rather than in a protonated form due to considerably higher barriers computed for the protonated adrenaline substrate. The results are discussed in the context of chemical mechanism of MAO enzymes and possible applications of multiscale simulation to rationalize the effects of MAO activity on adrenaline level.
Collapse
Affiliation(s)
- Gabriel Oanca
- Department of Computational Biochemistry and Drug Design, National Institute of Chemistry, Ljubljana, Slovenia.,Faculty of Physics, Alexandru Ioan Cuza University of Iasi, Iasi, Romania
| | - Jernej Stare
- Department of Computational Biochemistry and Drug Design, National Institute of Chemistry, Ljubljana, Slovenia
| | - Janez Mavri
- Department of Computational Biochemistry and Drug Design, National Institute of Chemistry, Ljubljana, Slovenia
| |
Collapse
|
32
|
Synthesis, characterization and photophysical behavior of heteroleptic ruthenium-complexed ladder-like structured polysilsesquioxanes. Macromol Res 2017. [DOI: 10.1007/s13233-017-5110-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
33
|
Multiscale simulation of monoamine oxidase catalyzed decomposition of phenylethylamine analogs. Eur J Pharmacol 2017; 817:46-50. [PMID: 28583428 DOI: 10.1016/j.ejphar.2017.05.061] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Revised: 03/31/2017] [Accepted: 05/31/2017] [Indexed: 01/01/2023]
Abstract
Phenylethylamine (PEA) is an endogenous amphetamine and its levels are increased by physical activity. As other biogenic monoamines, it is decomposed by monoamine oxidase (MAO) enzymes. The chemical mechanism of MAO, and flavoenzymes in general, is a subject of heated debate. We have previously shown that the rate-limiting step of MAO catalysis involves a hydride transfer from the substrate methylene group vicinal to the amino group to the N5 atom of the lumiflavin co-factor moiety. By using multiscale simulation on the Empirical Valence Bond (EVB) level, we studied the chemical reactivity of the monoamine oxidase B catalyzed decomposition of PEA and its two derivatives: p-chloro-β-methylphenylamine (p-CMP) and p-methoxy-β-methylphenethylamine (p-MMP). We calculated activation free energies of 17.1kcal/mol (PEA), 18.4kcal/mol (p-MMP) and 20.0kcal/mol (p-CMP), which are in excellent agreement with the experimental values of 16.7kcal/mol for PEA and 18.3kcal/mol for p-MMP, while the experimental value for p-CMP is not available. This gives strong support to the validity of our hydride transfer mechanism for both MAO A and B isoforms. The results are discussed in the context of the interplay between MAO point mutations and neuropsychiatric disorders.
Collapse
|
34
|
Hori Y, Ida T, Mizuno M. Potential energy construction in the diabatic picture for quantum mechanical rate constants of intermolecular proton transfer. Phys Chem Chem Phys 2017. [DOI: 10.1039/c7cp03024j] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We propose a simple method for potential construction in the diabatic picture and the estimation of thermal rate constants for intermolecular proton transfer reactions using quantum dynamics simulations carried out on the constructed potentials.
Collapse
Affiliation(s)
- Yuta Hori
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| | - Tomonori Ida
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| | - Motohiro Mizuno
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| |
Collapse
|
35
|
Poberžnik M, Purg M, Repič M, Mavri J, Vianello R. Empirical Valence Bond Simulations of the Hydride-Transfer Step in the Monoamine Oxidase A Catalyzed Metabolism of Noradrenaline. J Phys Chem B 2016; 120:11419-11427. [PMID: 27734680 DOI: 10.1021/acs.jpcb.6b09011] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines, such as dopamine, serotonin, and noradrenaline (NA), which is why they have been extensively implicated in the etiology and course of various neurodegenerative disorders and, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. The precise chemical mechanism through which MAOs regulate the amine concentration, which is vital for the development of novel inhibitors, is still not unambiguously determined in the literature. In this work, we present atomistic empirical valence bond simulations of the rate-limiting step of the MAO-A-catalyzed NA (norepinephrine) degradation, involving hydride transfer from the substrate α-methylene group to the flavin moiety of the flavin adenine dinucleotide prosthetic group, employing the full dimensionality and thermal fluctuations of the hydrated enzyme, with extensive configurational sampling. We show that MAO-A lowers the free energy of activation by 14.3 kcal mol-1 relative to that of the same reaction in aqueous solution, whereas the calculated activation free energy of ΔG‡ = 20.3 ± 1.6 kcal mol-1 is found to be in reasonable agreement with the correlated experimental value of 16.5 kcal mol-1. The results presented here strongly support the fact that both MAO-A and MAO-B isoforms function by the same hydride-transfer mechanism. We also considered a few point mutations of the "aromatic cage" tyrosine residue (Tyr444Phe, Tyr444Leu, Tyr444Trp, Tyr444His, and Tyr444Glu), and the calculated changes in the reaction barriers are in agreement with the experimental values, thus providing further support to the proposed mechanism.
Collapse
Affiliation(s)
- Matic Poberžnik
- Department of Physical and Organic Chemistry, Jožef Stefan Institute , Jamova cesta 39, SI-1000 Ljubljana, Slovenia
| | - Miha Purg
- Department of Cell and Molecular Biology, Uppsala Biomedical Centre , Husargatan 3, S-75124 Uppsala, Sweden
| | - Matej Repič
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry , Hajdrihova ulica 19, SI-1000 Ljubljana, Slovenia
| | - Janez Mavri
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry , Hajdrihova ulica 19, SI-1000 Ljubljana, Slovenia
| | - Robert Vianello
- Computational Organic Chemistry and Biochemistry Group, Ruđer Bošković Institute , Bijenička cesta 54, HR-10000 Zagreb, Croatia
| |
Collapse
|
36
|
Empirical Force Fields for Mechanistic Studies of Chemical Reactions in Proteins. Methods Enzymol 2016. [PMID: 27498633 DOI: 10.1016/bs.mie.2016.05.045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
Following chemical reactions in atomistic detail is one of the most challenging aspects of current computational approaches to chemistry. In this chapter the application of adiabatic reactive MD (ARMD) and its multistate version (MS-ARMD) are discussed. Both methods allow to study bond-breaking and bond-forming processes in chemical and biological processes. Particular emphasis is put on practical aspects for applying the methods to investigate the dynamics of chemical reactions. The chapter closes with an outlook of possible generalizations of the methods discussed.
Collapse
|
37
|
Vianello R, Domene C, Mavri J. The Use of Multiscale Molecular Simulations in Understanding a Relationship between the Structure and Function of Biological Systems of the Brain: The Application to Monoamine Oxidase Enzymes. Front Neurosci 2016; 10:327. [PMID: 27471444 PMCID: PMC4945635 DOI: 10.3389/fnins.2016.00327] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2016] [Accepted: 06/28/2016] [Indexed: 01/17/2023] Open
Abstract
HIGHLIGHTS Computational techniques provide accurate descriptions of the structure and dynamics of biological systems, contributing to their understanding at an atomic level.Classical MD simulations are a precious computational tool for the processes where no chemical reactions take place.QM calculations provide valuable information about the enzyme activity, being able to distinguish among several mechanistic pathways, provided a carefully selected cluster model of the enzyme is considered.Multiscale QM/MM simulation is the method of choice for the computational treatment of enzyme reactions offering quantitative agreement with experimentally determined reaction parameters.Molecular simulation provide insight into the mechanism of both the catalytic activity and inhibition of monoamine oxidases, thus aiding in the rational design of their inhibitors that are all employed and antidepressants and antiparkinsonian drugs. Aging society and therewith associated neurodegenerative and neuropsychiatric diseases, including depression, Alzheimer's disease, obsessive disorders, and Parkinson's disease, urgently require novel drug candidates. Targets include monoamine oxidases A and B (MAOs), acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and various receptors and transporters. For rational drug design it is particularly important to combine experimental synthetic, kinetic, toxicological, and pharmacological information with structural and computational work. This paper describes the application of various modern computational biochemistry methods in order to improve the understanding of a relationship between the structure and function of large biological systems including ion channels, transporters, receptors, and metabolic enzymes. The methods covered stem from classical molecular dynamics simulations to understand the physical basis and the time evolution of the structures, to combined QM, and QM/MM approaches to probe the chemical mechanisms of enzymatic activities and their inhibition. As an illustrative example, the later will focus on the monoamine oxidase family of enzymes, which catalyze the degradation of amine neurotransmitters in various parts of the brain, the imbalance of which is associated with the development and progression of a range of neurodegenerative disorders. Inhibitors that act mainly on MAO A are used in the treatment of depression, due to their ability to raise serotonin concentrations, while MAO B inhibitors decrease dopamine degradation and improve motor control in patients with Parkinson disease. Our results give strong support that both MAO isoforms, A and B, operate through the hydride transfer mechanism. Relevance of MAO catalyzed reactions and MAO inhibition in the context of neurodegeneration will be discussed.
Collapse
Affiliation(s)
- Robert Vianello
- Computational Organic Chemistry and Biochemistry Group, Ruđer Bošković InstituteZagreb, Croatia
| | - Carmen Domene
- Department of Chemistry, King's College LondonLondon, UK
- Chemistry Research Laboratory, University of OxfordOxford, UK
| | - Janez Mavri
- Department of Computational Biochemistry and Drug Design, National Institute of ChemistryLjubljana, Slovenia
| |
Collapse
|
38
|
Bauer P, Carlsson ÅJ, Amrein BA, Dobritzsch D, Widersten M, Kamerlin SCL. Conformational diversity and enantioconvergence in potato epoxide hydrolase 1. Org Biomol Chem 2016; 14:5639-51. [PMID: 27049844 PMCID: PMC5315018 DOI: 10.1039/c6ob00060f] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2016] [Accepted: 03/31/2016] [Indexed: 01/11/2023]
Abstract
Potato epoxide hydrolase 1 (StEH1) is a biocatalytically important enzyme that exhibits rich enantio- and regioselectivity in the hydrolysis of chiral epoxide substrates. In particular, StEH1 has been demonstrated to enantioconvergently hydrolyze racemic mixes of styrene oxide (SO) to yield (R)-1-phenylethanediol. This work combines computational, crystallographic and biochemical analyses to understand both the origins of the enantioconvergent behavior of the wild-type enzyme, as well as shifts in activities and substrate binding preferences in an engineered StEH1 variant, R-C1B1, which contains four active site substitutions (W106L, L109Y, V141K and I155V). Our calculations are able to reproduce both the enantio- and regioselectivities of StEH1, and demonstrate a clear link between different substrate binding modes and the corresponding selectivity, with the preferred binding modes being shifted between the wild-type enzyme and the R-C1B1 variant. Additionally, we demonstrate that the observed changes in selectivity and the corresponding enantioconvergent behavior are due to a combination of steric and electrostatic effects that modulate both the accessibility of the different carbon atoms to the nucleophilic side chain of D105, as well as the interactions between the substrate and protein amino acid side chains and active site water molecules. Being able to computationally predict such subtle effects for different substrate enantiomers, as well as to understand their origin and how they are affected by mutations, is an important advance towards the computational design of improved biocatalysts for enantioselective synthesis.
Collapse
Affiliation(s)
- P. Bauer
- Science for Life Laboratory , Department of Cell and Molecular Biology , Uppsala University , BMC Box 596 , S-751 24 Uppsala , Sweden .
| | - Å. Janfalk Carlsson
- Department of Chemistry-BMC , Uppsala University , BMC Box 576 , S-751 23 Uppsala , Sweden . ;
| | - B. A. Amrein
- Science for Life Laboratory , Department of Cell and Molecular Biology , Uppsala University , BMC Box 596 , S-751 24 Uppsala , Sweden .
| | - D. Dobritzsch
- Department of Chemistry-BMC , Uppsala University , BMC Box 576 , S-751 23 Uppsala , Sweden . ;
| | - M. Widersten
- Department of Chemistry-BMC , Uppsala University , BMC Box 576 , S-751 23 Uppsala , Sweden . ;
| | - S. C. L. Kamerlin
- Science for Life Laboratory , Department of Cell and Molecular Biology , Uppsala University , BMC Box 596 , S-751 24 Uppsala , Sweden .
| |
Collapse
|
39
|
Hernández-Fernández F, Pavanello M, Visscher L. Effect of metallation, substituents and inter/intra-molecular polarization on electronic couplings for hole transport in stacked porphyrin dyads. Phys Chem Chem Phys 2016; 18:21122-32. [DOI: 10.1039/c6cp00516k] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Hole transport properties for stacked porphyrins are analyzed with density functional theory.
Collapse
Affiliation(s)
- F. Hernández-Fernández
- Theoretical Chemistry
- Faculty of Sciences
- VU University Amsterdam
- 1081 HV Amsterdam
- The Netherlands
| | - M. Pavanello
- Department of Chemistry
- Rutgers University
- Newark
- USA
| | - L. Visscher
- Theoretical Chemistry
- Faculty of Sciences
- VU University Amsterdam
- 1081 HV Amsterdam
- The Netherlands
| |
Collapse
|
40
|
Danielsson J, Meuwly M. Atomistic Simulation of Adiabatic Reactive Processes Based on Multi-State Potential Energy Surfaces. J Chem Theory Comput 2015; 4:1083-93. [PMID: 26636362 DOI: 10.1021/ct800066q] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The adiabatic reactive molecular dynamics (ARMD) method provides a framework to study chemical reactions using molecular dynamics simulations with minimal computational overhead. Here, ARMD is generalized to an arbitrary reactive process between two states in which reactants and products can be treated by an atomistic force field. The implementation is described, and the method is applied to two systems: the kinetics of NO rebinding to myoglobin (Mb) as a validation system and the conformational transition in neuroglobin (Ngb) which explores the full functionality of ARMD. For MbNO, the nonexponential kinetics observed both in experiment and earlier ARMD studies is reproduced. Furthermore, the sensitivity of the results with respect to the asymptotic separation between the two potential energy surfaces (NO bound and unbound) is studied.
Collapse
Affiliation(s)
- Jonas Danielsson
- Chemistry Department, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Chemistry Department, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
41
|
|
42
|
Wesolowski TA, Shedge S, Zhou X. Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. Chem Rev 2015; 115:5891-928. [DOI: 10.1021/cr500502v] [Citation(s) in RCA: 219] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Sapana Shedge
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Xiuwen Zhou
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| |
Collapse
|
43
|
Shurki A, Derat E, Barrozo A, Kamerlin SCL. How valence bond theory can help you understand your (bio)chemical reaction. Chem Soc Rev 2015; 44:1037-52. [DOI: 10.1039/c4cs00241e] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This Tutorial Review describes modern VB theory, using a didactic example of an SN2 reaction to illustrate its immediate applications.
Collapse
Affiliation(s)
- Avital Shurki
- Institute for Drug Research
- School of Pharmacy
- The Hebrew University of Jerusalem
- Jerusalem
- Israel
| | - Etienne Derat
- Institut Parisien de Chimie Moléculaire
- UMR 8232
- Paris
- France
| | - Alexandre Barrozo
- Science for Life Laboratory
- Department of Cell and Molecular Biology
- Uppsala University
- Uppsala
- Sweden
| | | |
Collapse
|
44
|
Repič M, Vianello R, Purg M, Duarte F, Bauer P, Kamerlin SCL, Mavri J. Empirical valence bond simulations of the hydride transfer step in the monoamine oxidase B catalyzed metabolism of dopamine. Proteins 2014; 82:3347-55. [PMID: 25220264 DOI: 10.1002/prot.24690] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2014] [Accepted: 09/10/2014] [Indexed: 01/05/2023]
Abstract
Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines such as dopamine, serotonin and noradrenaline. In this work, we present a comprehensive study of the rate-limiting step of dopamine degradation by MAO B, which consists in the hydride transfer from the methylene group of the substrate to the flavin moiety of the FAD prosthetic group. This article builds on our previous quantum chemical study of the same reaction using a cluster model (Vianello et al., Eur J Org Chem 2012; 7057), but now considering the full dimensionality of the hydrated enzyme with extensive configurational sampling. We show that MAO B is specifically tuned to catalyze the hydride transfer step from the substrate to the flavin moiety of the FAD prosthetic group and that it lowers the activation barrier by 12.3 kcal mol⁻¹ compared to the same reaction in aqueous solution, a rate enhancement of more than nine orders of magnitude. Taking into account the deprotonation of the substrate prior to the hydride transfer reaction, the activation barrier in the enzyme is calculated to be 16.1 kcal mol⁻¹, in excellent agreement with the experimental value of 16.5 kcal mol⁻¹. Additionally, we demonstrate that the protonation state of the active site residue Lys296 does not have an influence on the hydride transfer reaction.
Collapse
Affiliation(s)
- Matej Repič
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry, Hajdrihova 19, SI-1000, Ljubljana, Slovenia
| | | | | | | | | | | | | |
Collapse
|
45
|
McKenzie RH. A diabatic state model for double proton transfer in hydrogen bonded complexes. J Chem Phys 2014; 141:104314. [DOI: 10.1063/1.4895007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Jacob CR, Neugebauer J. Subsystem density-functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1175] [Citation(s) in RCA: 229] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Christoph R. Jacob
- Center for Functional Nanostructures and Institute of Physical Chemistry; Karlsruhe Institute of Technology (KIT); Karlsruhe Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut; Westfälische Wilhelms-Universität Münster; Münster Germany
| |
Collapse
|
47
|
Shurki A, Sharir-Ivry A. Valence Bond-Based Hybrid Quantum Mechanics Molecular Mechanics Approaches and Proper Inclusion of the Effect of the Surroundings. Isr J Chem 2014. [DOI: 10.1002/ijch.201400038] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
48
|
Ramos P, Pavanello M. Quantifying Environmental Effects on the Decay of Hole Transfer Couplings in Biosystems. J Chem Theory Comput 2014; 10:2546-56. [DOI: 10.1021/ct400921r] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
49
|
Solovyeva A, Pavanello M, Neugebauer J. Describing long-range charge-separation processes with subsystem density-functional theory. J Chem Phys 2014; 140:164103. [DOI: 10.1063/1.4871301] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|
50
|
Migliore A, Polizzi NF, Therien M, Beratan DN. Biochemistry and theory of proton-coupled electron transfer. Chem Rev 2014; 114:3381-465. [PMID: 24684625 PMCID: PMC4317057 DOI: 10.1021/cr4006654] [Citation(s) in RCA: 350] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2013] [Indexed: 02/01/2023]
Affiliation(s)
- Agostino Migliore
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Nicholas F. Polizzi
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Michael
J. Therien
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - David N. Beratan
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|