1
|
Nottoli M, Vanich E, Cupellini L, Scalmani G, Pelosi C, Lipparini F. Importance of Polarizable Embedding for Computing Optical Rotation: The Case of Camphor in Ethanol. J Phys Chem Lett 2024:7992-7999. [PMID: 39078659 DOI: 10.1021/acs.jpclett.4c01550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
Solvation effects on optical rotation are notoriously challenging to model for computational chemistry, as the specific rotatory power of a molecule can vary wildly going from apolar to polar or even protic solvents. To address such a problem, we present a polarizable embedding implementation of an electric and magnetic response property based on density functional theory and the AMOEBA polarizable force field, and apply such an implementation to the study of the optical rotation of camphor in ethanol. By comparing a continuum model, and electrostatic and polarizable embedding QM/MM models, we observe that accounting for the environment's polarization gives rise to not only a different quantitative prediction, in very good agreement with experiments for the QM/AMOEBA model, but also to a very different qualitative picture, with the values of the optical rotation computed along a classical molecular dynamics trajectory with electrostatic embedding being statistically uncorrelated to the ones obtained with the polarizable description.
Collapse
Affiliation(s)
- Michele Nottoli
- Institute of Applied Analysis and Numerical Simulation, Universität Stuttgart, Pfaffenwaldring 57, D-70569, Stuttgart, Germany
| | - Edoardo Vanich
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street Building 40, Wallingford, Connecticut 06492, United States
| | - Chiara Pelosi
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| |
Collapse
|
2
|
Calcagno F, Maryasin B, Garavelli M, Avagliano D, Rivalta I. Modeling solvent effects and convergence of 31P-NMR shielding calculations with COBRAMM. J Comput Chem 2024; 45:1562-1575. [PMID: 38514234 DOI: 10.1002/jcc.27338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 02/15/2024] [Accepted: 02/22/2024] [Indexed: 03/23/2024]
Abstract
Solvent effects on 31P-NMR parameters for triphenylphosphine oxide and triphenylphosphine in chloroform have been extensively investigated by testing different solvation models. The solvent is described implicitly, mixed implicitly/explicitly, and using full explicit models. Polarizable continuum model (PCM), molecular dynamic (MD) simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations are used to disclose the effects of solute/solvent interactions and, more generally, the role of the embedding in NMR simulations. The results show the beneficial effect of carrying out QM/MM optimizations on top of geometries directly extracted from classical MD simulations, used to ensure representative conformational sampling. The nuclear shielding convergence has been tested against a different number of snapshots and with the inclusion of solvent shells into the QM region. An automated MD//QM/MM//GIAO protocol, implemented in the COBRAMM package, is here proposed and tested on trimethyl phosphite showing that our approach boosts the convergence of nuclear shielding satisfactorily. The present work aims to be a stepping-stone to assess proper QM/MM computational strategies in simulating chemical shifts in non-homogeneous systems like supramolecular and biological systems.
Collapse
Affiliation(s)
- Francesco Calcagno
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
- Center for Chemical Catalysis - C3, University of Bologna, Bologna, Italy
| | - Boris Maryasin
- Institute of Organic Chemistry, University of Vienna, Vienna, Austria
- Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Marco Garavelli
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
| | - Davide Avagliano
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
| | - Ivan Rivalta
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
- Center for Chemical Catalysis - C3, University of Bologna, Bologna, Italy
- ENSL, CNRS, Laboratoire de Chimie UMR 5182, Lyon, France
| |
Collapse
|
3
|
Carrasco-Busturia D, Olsen JMH. Polarizable Embedding Potentials through Molecular Fractionation with Conjugate Caps Including Hydrogen Bonds. J Chem Theory Comput 2023; 19:6510-6520. [PMID: 37665268 DOI: 10.1021/acs.jctc.3c00613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
Polarizable embedding (PE) refers to classical embedding approaches, such as those used in quantum mechanics/molecular mechanics (QM/MM), that allow mutual polarization between the quantum and classical regions. The quality of the embedding potential is critical to provide accurate results, e.g., for spectroscopic properties and dynamical processes. High-quality embedding-potential parameters can be obtained by dividing the classical region into smaller fragments and deriving the parameters from ab initio calculations on the fragments. For solvents and other systems composed of small molecules, the fragments can be individual molecules, while a more complicated fragmentation procedure is needed for larger molecules, such as proteins and nucleic acids. One such fragmentation strategy is the molecular fractionation with conjugate caps (MFCC) approach. As is widely known, hydrogen bonds play a key role in many biomolecular systems, e.g., in proteins, where they are responsible for the secondary structure. In this work, we assess the effects of including hydrogen-bond fragmentation in the MFCC procedure [MFCC(HB)] for deriving the embedding-potential parameters. The MFCC(HB) extension is evaluated on several molecular systems, ranging from small model systems to proteins, directly in terms of molecular electrostatic potentials and embedding potentials and indirectly in terms of selected properties of chromophores embedded in water and complex protein environments.
Collapse
|
4
|
Kvedaravičiūtė S, Carrasco-Busturia D, Møller KB, Olsen JMH. Polarizable Embedding without Artificial Boundary Polarization. J Chem Theory Comput 2023; 19:5122-5141. [PMID: 37458793 DOI: 10.1021/acs.jctc.3c00434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/09/2023]
Abstract
We present a fully self-consistent polarizable embedding (PE) model that does not suffer from unphysical boundary polarization. This is achieved through the use of the minimum-image convention (MIC) in the induced electrostatics. It is a simple yet effective approach that includes a more physically accurate description of the polarization throughout the molecular system. Using PE with MIC (PE-MIC), we shed new light on the limitations of commonly employed cutoff models, such as the droplet model, when used in PE calculations. Specifically, we investigate the effects of the unphysical polarization at the outer boundary by comparing induced dipoles and the associated electrostatic potentials, as well as some optical properties of solute-solvent and biomolecular systems. We show that the magnitude of the inaccuracies caused by the unphysical polarization depends on multiple parameters: the nature of the quantum subsystem and of the environment, the cutoff model and distance, and the calculated property.
Collapse
Affiliation(s)
| | | | - Klaus B Møller
- DTU Chemistry, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | | |
Collapse
|
5
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
6
|
Jørgensen FK, Reinholdt P, Hedegård ED, Kongsted J. Nuclear Magnetic Shielding Constants with the Polarizable Density Embedding Model. J Chem Theory Comput 2022; 18:7384-7393. [PMID: 36332108 DOI: 10.1021/acs.jctc.2c00829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We extend the polarizable density embedding (PDE) model to support the calculation of nuclear magnetic resonance (NMR) shielding constants using gauge-including atomic orbitals (GIAOs) within a density functional theory (DFT) framework. The PDE model divides the total system into fragments, describing some by quantum mechanics (QM) and the others through an embedding model. The PDE model uses anisotropic polarizabilities, inter-fragment two-electron Coulomb integrals, and a non-local repulsion operator to emulate the QM effects. The terms involving Coulomb integrals are straightforwardly extended with GIAOs. In contrast, we consider two approaches to handle the gauge dependency of the non-local operator, employing either simple symmetrization or a gauge transformation. We find the latter approach to be most stable with respect to increasing the basis set size of the QM region. We examine the accuracy of the PDE model for calculating NMR shielding constants on several solutes in a water solution. The performance is compared with the classical polarizable embedding (PE) model in addition to supermolecular reference calculations. Based on these systems, we address the basis set convergence characteristics and the QM region size requirements. Furthermore, we investigate the performance of the PDE model for a system with significant electron spill-out. In many cases, we find that the PDE model outperforms the PE model, especially regarding the accuracy of nuclear shielding constants when using small QM region sizes and in systems with significant electron spill-out.
Collapse
Affiliation(s)
- Frederik Kamper Jørgensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Erik Donovan Hedegård
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| |
Collapse
|
7
|
Palivec V, Pohl R, Kaminský J, Martinez-Seara H. Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment. J Chem Theory Comput 2022; 18:4373-4386. [PMID: 35687789 DOI: 10.1021/acs.jctc.2c00127] [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/28/2022]
Abstract
Determining the structure of saccharides in their native environment is crucial to understanding their function and more accurately targeting their utilization. Nuclear magnetic resonance observables such as the nuclear Overhauser effect or spin-spin coupling constants are routinely utilized to study saccharides in their native water environment. However, while highly sensitive to the local environment, chemical shifts are mostly overlooked, despite being commonly measured for compounds identification. Although chemical shifts carry considerable structural information, their direct association with structure is notoriously difficult. This is mostly due to the similarity in the chemical nature of most saccharides causing similar physicochemical environments close to sugar C and H atoms, resulting in comparable chemical shifts. The rise of computational power allows one to compute reliable chemical shifts and use them to determine atomistic details of these sugars in solution. However, any prediction is severely limited by the computational protocol used and its accuracy. In this work, we studied a set of 31 saccharides on which we evaluated various computational protocols to calculate the total number of 375 1H and 327 13C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating 1H and 13C chemical shifts that we recommend for further use. These protocols can help with the interpretation of experimental spectra, but we also show that they are also capable of structure prediction independently. This is possible because of the low mean absolute deviations of calculated shifts from the experiment (0.06 ppm for 1H and 1.09 ppm for 13C). We explore different solvation methods, basis sets, and optimization schemes to reach such accuracy. A correct sampling of the conformation phase space of flexible sugar molecules is also key to obtaining accurately converged theoretical chemical shifts. The linear regression method was applied to convert the calculated isotropic nuclear magnetic shielding constants to simulated chemical shifts comparable with the experiment. The achieved level of accuracy can help in utilizing chemical shifts for elucidating the 3D atomistic structure of saccharides in aqueous solutions. All linear regression parameters obtained on our extensive set of sugars for all the tested protocols can be reutilized in future works.
Collapse
Affiliation(s)
- Vladimír Palivec
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Radek Pohl
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| |
Collapse
|
8
|
Unzueta PA, Beran GJO. Polarizable continuum models provide an effective electrostatic embedding model for fragment-based chemical shift prediction in challenging systems. J Comput Chem 2020; 41:2251-2265. [PMID: 32748418 DOI: 10.1002/jcc.26388] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 06/04/2020] [Accepted: 07/04/2020] [Indexed: 12/25/2022]
Abstract
Ab initio nuclear magnetic resonance chemical shift prediction provides an important tool for interpreting and assigning experimental spectra, but it becomes computationally prohibitive in large systems. The computational costs can be reduced considerably by fragmentation of the large system into a series of contributions from many smaller subsystems. However, the presence of charged functional groups and the need to partition the system across covalent bonds create complications in biomolecules that typically require the use of large fragments and careful descriptions of the electrostatic environment. The present work shows how a model that combines chemical shielding contributions from non-overlapping monomer and dimer fragments embedded in a polarizable continuum model provides a simple, easy-to-implement, and computationally inexpensive approach for predicting chemical shifts in complex systems. The model's performance proves rather insensitive to the continuum dielectric constant, making the selection of the optimal embedding dielectric less critical. The PCM-embedded fragment model is demonstrated to perform well across systems ranging from molecular crystals to proteins.
Collapse
Affiliation(s)
- Pablo A Unzueta
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| | - Gregory J O Beran
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| |
Collapse
|
9
|
Dietschreit JCB, Wagner A, Le TA, Klein P, Schindelin H, Opatz T, Engels B, Hellmich UA, Ochsenfeld C. Predicting
19
F NMR Chemical Shifts: A Combined Computational and Experimental Study of a Trypanosomal Oxidoreductase–Inhibitor Complex. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202000539] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Johannes C. B. Dietschreit
- Theoretical Chemistry Department of Chemistry University of Munich (LMU) Butenandtstr. 7 81377 Munich Germany
| | - Annika Wagner
- Dept. Chemistry Section Biochemistry Johannes Gutenberg-Universität Mainz 55128 Mainz Germany
- Centre for Biomolecular Magnetic Resonance (BMRZ) Goethe-University Frankfurt Max-von-Laue Str. 9 60438 Frankfurt Germany
| | - T. Anh Le
- Institute for Physical and Theoretical Chemistry University of Würzburg Emil-Fischer-Straße 42 97074 Würzburg Germany
| | - Philipp Klein
- Dept. Chemistry Section Organic Chemistry Johannes Gutenberg-Universität Mainz 55128 Mainz Germany
| | - Hermann Schindelin
- Institute of Structural Biology Rudolf Virchow Center for Experimental Biomedicine University of Würzburg 97080 Würzburg Germany
| | - Till Opatz
- Dept. Chemistry Section Organic Chemistry Johannes Gutenberg-Universität Mainz 55128 Mainz Germany
| | - Bernd Engels
- Institute for Physical and Theoretical Chemistry University of Würzburg Emil-Fischer-Straße 42 97074 Würzburg Germany
| | - Ute A. Hellmich
- Dept. Chemistry Section Biochemistry Johannes Gutenberg-Universität Mainz 55128 Mainz Germany
- Centre for Biomolecular Magnetic Resonance (BMRZ) Goethe-University Frankfurt Max-von-Laue Str. 9 60438 Frankfurt Germany
| | - Christian Ochsenfeld
- Theoretical Chemistry Department of Chemistry University of Munich (LMU) Butenandtstr. 7 81377 Munich Germany
- Max Planck Institute for Solid State Research 70569 Stuttgart Germany
| |
Collapse
|
10
|
Dietschreit JCB, Wagner A, Le TA, Klein P, Schindelin H, Opatz T, Engels B, Hellmich UA, Ochsenfeld C. Predicting 19 F NMR Chemical Shifts: A Combined Computational and Experimental Study of a Trypanosomal Oxidoreductase-Inhibitor Complex. Angew Chem Int Ed Engl 2020; 59:12669-12673. [PMID: 32239740 PMCID: PMC7496126 DOI: 10.1002/anie.202000539] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2020] [Revised: 03/22/2020] [Indexed: 02/02/2023]
Abstract
The absence of fluorine from most biomolecules renders it an excellent probe for NMR spectroscopy to monitor inhibitor-protein interactions. However, predicting the binding mode of a fluorinated ligand from a chemical shift (or vice versa) has been challenging due to the high electron density of the fluorine atom. Nonetheless, reliable 19 F chemical-shift predictions to deduce ligand-binding modes hold great potential for in silico drug design. Herein, we present a systematic QM/MM study to predict the 19 F NMR chemical shifts of a covalently bound fluorinated inhibitor to the essential oxidoreductase tryparedoxin (Tpx) from African trypanosomes, the causative agent of African sleeping sickness. We include many protein-inhibitor conformations as well as monomeric and dimeric inhibitor-protein complexes, thus rendering it the largest computational study on chemical shifts of 19 F nuclei in a biological context to date. Our predicted shifts agree well with those obtained experimentally and pave the way for future work in this area.
Collapse
Affiliation(s)
| | - Annika Wagner
- Dept. ChemistrySection BiochemistryJohannes Gutenberg-Universität Mainz55128MainzGermany
- Centre for Biomolecular Magnetic Resonance (BMRZ)Goethe-University FrankfurtMax-von-Laue Str. 960438FrankfurtGermany
| | - T. Anh Le
- Institute for Physical and Theoretical ChemistryUniversity of WürzburgEmil-Fischer-Straße 4297074WürzburgGermany
| | - Philipp Klein
- Dept. ChemistrySection Organic ChemistryJohannes Gutenberg-Universität Mainz55128MainzGermany
| | - Hermann Schindelin
- Institute of Structural BiologyRudolf Virchow Center for Experimental BiomedicineUniversity of Würzburg97080WürzburgGermany
| | - Till Opatz
- Dept. ChemistrySection Organic ChemistryJohannes Gutenberg-Universität Mainz55128MainzGermany
| | - Bernd Engels
- Institute for Physical and Theoretical ChemistryUniversity of WürzburgEmil-Fischer-Straße 4297074WürzburgGermany
| | - Ute A. Hellmich
- Dept. ChemistrySection BiochemistryJohannes Gutenberg-Universität Mainz55128MainzGermany
- Centre for Biomolecular Magnetic Resonance (BMRZ)Goethe-University FrankfurtMax-von-Laue Str. 960438FrankfurtGermany
| | - Christian Ochsenfeld
- Theoretical ChemistryDepartment of ChemistryUniversity of Munich (LMU)Butenandtstr. 781377MunichGermany
- Max Planck Institute for Solid State Research70569StuttgartGermany
| |
Collapse
|
11
|
Olsen JMH, Reine S, Vahtras O, Kjellgren E, Reinholdt P, Hjorth Dundas KO, Li X, Cukras J, Ringholm M, Hedegård ED, Di Remigio R, List NH, Faber R, Cabral Tenorio BN, Bast R, Pedersen TB, Rinkevicius Z, Sauer SPA, Mikkelsen KV, Kongsted J, Coriani S, Ruud K, Helgaker T, Jensen HJA, Norman P. Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems. J Chem Phys 2020; 152:214115. [PMID: 32505165 DOI: 10.1063/1.5144298] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.
Collapse
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Simen Reine
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Erik Kjellgren
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Karen Oda Hjorth Dundas
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Xin Li
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Janusz Cukras
- Department of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Magnus Ringholm
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Erik D Hedegård
- Division of Theoretical Chemistry, Lund University, SE-223 62 Lund, Sweden
| | - Roberto Di Remigio
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Rasmus Faber
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | | | - Radovan Bast
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Thomas Bondo Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Kenneth Ruud
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Trygve Helgaker
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| |
Collapse
|
12
|
Imamura K, Yamazaki T, Yokogawa D, Higashi M, Sato H. Nuclear magnetic shielding of molecule in solution based on reference interaction site model self-consistent field with spatial electron density distribution. J Chem Phys 2020; 152:194102. [DOI: 10.1063/5.0008903] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Takeshi Yamazaki
- 1QB Information Technologies (1QBit), 200-1285 West Pender Street, Vancouver, British Columbia V6E 4B1, Canada
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyoto 615-8520, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyoto 615-8520, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
13
|
Cabral BJ. Born-Oppenheimer molecular dynamics, hydrogen bond interactions and magnetic properties of liquid hydrogen cyanide. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.09.092] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
14
|
Caprasecca S, Cupellini L, Jurinovich S, Loco D, Lipparini F, Mennucci B. A polarizable QM/MM description of environment effects on NMR shieldings: from solvated molecules to pigment–protein complexes. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2264-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
15
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
16
|
|
17
|
Galamba N, Cabral BJC. Magnetic properties and core electron binding energies of liquid water. J Chem Phys 2018; 148:044510. [PMID: 29390846 DOI: 10.1063/1.5011037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The magnetic properties and the core and inner valence electron binding energies of liquid water are investigated. The adopted methodology relies on the combination of molecular dynamics and electronic structure calculations. Born-Oppenheimer molecular dynamics with the Becke and Lee-Yang-Parr functionals for exchange and correlation, respectively, and includes an empirical correction (BLYP-D3) functional and classical molecular dynamics with the TIP4P/2005-F model were carried out. The Keal-Tozer functional was applied for predicting magnetic shielding and spin-spin coupling constants. Core and inner valence electron binding energies in liquid water were calculated with symmetry adapted cluster-configuration interaction. The relationship between the magnetic shielding constant σ(17O), the role played by the oxygen atom as a proton acceptor and donor, and the tetrahedral organisation of liquid water are investigated. The results indicate that the deshielding of the oxygen atom in water is very dependent on the order parameter (q) describing the tetrahedral organisation of the hydrogen bond network. The strong sensitivity of magnetic properties on changes of the electronic density in the nuclei environment is illustrated by a correlation between σ(17O) and the energy gap between the 1a1[O1s] (core) and the 2a1 (inner valence) orbitals of water. Although several studies discussed the eventual connection between magnetic properties and core electron binding energies, such a correlation could not be clearly established. Here, we demonstrate that for liquid water this correlation exists although involving the gap between electron binding energies of core and inner valence orbitals.
Collapse
Affiliation(s)
- N Galamba
- Centro de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Edifício C8, Campo Grande, 1749-016 Lisboa, Portugal
| | - Benedito J C Cabral
- BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
| |
Collapse
|
18
|
Tupikina EY, Efimova AA, Denisov GS, Tolstoy PM. NMR Chemical Shift of a Helium Atom as a Probe for Electronic Structure of FH, F–, (FHF)−, and FH2+. J Phys Chem A 2017; 121:9654-9662. [DOI: 10.1021/acs.jpca.7b10189] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- E. Yu. Tupikina
- Department
of Physics and ‡Center for Magnetic Resonance, St. Petersburg State University, St. Petersburg 198504, Russia
| | - A. A. Efimova
- Department
of Physics and ‡Center for Magnetic Resonance, St. Petersburg State University, St. Petersburg 198504, Russia
| | - G. S. Denisov
- Department
of Physics and ‡Center for Magnetic Resonance, St. Petersburg State University, St. Petersburg 198504, Russia
| | | |
Collapse
|
19
|
Zhao D, Song R, Li W, Ma J, Dong H, Li S. Accurate Prediction of NMR Chemical Shifts in Macromolecular and Condensed-Phase Systems with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:5231-5239. [DOI: 10.1021/acs.jctc.7b00380] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dongbo Zhao
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Ruiheng Song
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Wei Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Hao Dong
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|
20
|
Nørby MS, Olsen JMH, Steinmann C, Kongsted J. Modeling Electronic Circular Dichroism within the Polarizable Embedding Approach. J Chem Theory Comput 2017; 13:4442-4451. [DOI: 10.1021/acs.jctc.7b00712] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Morten S. Nørby
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | | | - Casper Steinmann
- Department
of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Jacob Kongsted
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
21
|
Hedegård ED, Bast R, Kongsted J, Olsen JMH, Jensen HJA. Relativistic Polarizable Embedding. J Chem Theory Comput 2017; 13:2870-2880. [DOI: 10.1021/acs.jctc.7b00162] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
| | - Radovan Bast
- High
Performance Computing Group, UiT The Arctic University of Norway, Tromsø 9037, Norway
| | - Jacob Kongsted
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | | | - Hans Jørgen Aagaard Jensen
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
22
|
Sundholm D, Rauhalahti M, Özcan N, Mera-Adasme R, Kussmann J, Luenser A, Ochsenfeld C. Nuclear Magnetic Shieldings of Stacked Aromatic and Antiaromatic Molecules. J Chem Theory Comput 2017; 13:1952-1962. [PMID: 28287722 DOI: 10.1021/acs.jctc.6b01250] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dage Sundholm
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FIN-00014 Helsinki, Finland
| | - Markus Rauhalahti
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FIN-00014 Helsinki, Finland
| | - Nergiz Özcan
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FIN-00014 Helsinki, Finland
| | - Raúl Mera-Adasme
- Departamento
de Ciencias del Ambiente, Universidad de Santiago de Chile (USACH), Av. Libertador Bernardo O’Higgins 3363, 9170022 Estación Central, Chile
| | - Jörg Kussmann
- Department
of Chemistry, University of Munich (LMU), München D-81377, Germany
| | - Arne Luenser
- Department
of Chemistry, University of Munich (LMU), München D-81377, Germany
| | | |
Collapse
|
23
|
Steinmann C, Bratholm LA, Olsen JMH, Kongsted J. Automated Fragmentation Polarizable Embedding Density Functional Theory (PE-DFT) Calculations of Nuclear Magnetic Resonance (NMR) Shielding Constants of Proteins with Application to Chemical Shift Predictions. J Chem Theory Comput 2017; 13:525-536. [PMID: 27992211 DOI: 10.1021/acs.jctc.6b00965] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Full-protein nuclear magnetic resonance (NMR) shielding constants based on ab initio calculations are desirable, because they can assist in elucidating protein structures from NMR experiments. In this work, we present NMR shielding constants computed using a new automated fragmentation (J. Phys. Chem. B 2009, 113, 10380-10388) approach in the framework of polarizable embedding density functional theory. We extend our previous work to give both basis set recommendations and comment on how large the quantum mechanical region should be to successfully compute 13C NMR shielding constants that are comparable with experiment. The introduction of a probabilistic linear regression model allows us to substantially reduce the number of snapshots that are needed to make comparisons with experiment. This approach is further improved by augmenting snapshot selection with chemical shift predictions by which we can obtain a representative subset of snapshots that gives the smallest predicted error, compared to experiment. Finally, we use this subset of snapshots to calculate the NMR shielding constants at the PE-KT3/pcSseg-2 level of theory for all atoms in the protein GB3.
Collapse
Affiliation(s)
- Casper Steinmann
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, United Kingdom.,Department of Physics, Chemistry and Pharmacy, University of Southern Denmark , DK-5230 Odense M, Denmark
| | | | | | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark , DK-5230 Odense M, Denmark
| |
Collapse
|
24
|
Bennie SJ, van der Kamp MW, Pennifold RCR, Stella M, Manby FR, Mulholland AJ. A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase. J Chem Theory Comput 2016; 12:2689-97. [PMID: 27159381 DOI: 10.1021/acs.jctc.6b00285] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.
Collapse
Affiliation(s)
- Simon J Bennie
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
- School of Biochemistry, Biomedical Sciences Building, University Walk, University of Bristol , Bristol BS8 1TD, U.K
| | - Robert C R Pennifold
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Martina Stella
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Frederick R Manby
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| |
Collapse
|
25
|
Chikayama E, Shimbo Y, Komatsu K, Kikuchi J. The Effect of Molecular Conformation on the Accuracy of Theoretical (1)H and (13)C Chemical Shifts Calculated by Ab Initio Methods for Metabolic Mixture Analysis. J Phys Chem B 2016; 120:3479-87. [PMID: 26963288 DOI: 10.1021/acs.jpcb.5b12748] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
NMR spectroscopy is a powerful method for analyzing metabolic mixtures. The information obtained from an NMR spectrum is in the form of physical parameters, such as chemical shifts, and construction of databases for many metabolites will be useful for data interpretation. To increase the accuracy of theoretical chemical shifts for development of a database for a variety of metabolites, the effects of sets of conformations (structural ensembles) and the levels of theory on computations of theoretical chemical shifts were systematically investigated for a set of 29 small molecules in the present study. For each of the 29 compounds, 101 structures were generated by classical molecular dynamics at 298.15 K, and then theoretical chemical shifts for 164 (1)H and 123 (13)C atoms were calculated by ab initio quantum chemical methods. Six levels of theory were used by pairing Hartree-Fock, B3LYP (density functional theory), or second order Møller-Plesset perturbation with 6-31G or aug-cc-pVDZ basis set. The six average fluctuations in the (1)H chemical shift were ±0.63, ± 0.59, ± 0.70, ± 0.62, ± 0.75, and ±0.66 ppm for the structural ensembles, and the six average errors were ±0.34, ± 0.27, ± 0.32, ± 0.25, ± 0.32, and ±0.25 ppm. The results showed that chemical shift fluctuations with changes in the conformation because of molecular motion were larger than the differences between computed and experimental chemical shifts for all six levels of theory. In conclusion, selection of an appropriate structural ensemble should be performed before theoretical chemical shift calculations for development of an accurate database for a variety of metabolites.
Collapse
Affiliation(s)
- Eisuke Chikayama
- RIKEN Center for Sustainable Resource Science , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan.,Department of Information Systems, Niigata University of International and Information Studies , 3-1-1 Mizukino, Nishi-ku, Niigata, Niigata 950-2292, Japan
| | - Yudai Shimbo
- NEC Solution Innovators, Ltd. , 2-2-41 Ekimae, Kashiwazaki, Niigata 945-0055, Japan
| | - Keiko Komatsu
- RIKEN Center for Sustainable Resource Science , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Jun Kikuchi
- RIKEN Center for Sustainable Resource Science , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan.,Graduate School of Medical Life Science, Yokohama City University , 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan.,Graduate School of Bioagricultural Sciences, Nagoya University , 1 Furo-cho, Chikusa-ku, Nagoya, Aichi 464-0810, Japan
| |
Collapse
|
26
|
List NH, Olsen JMH, Kongsted J. Excited states in large molecular systems through polarizable embedding. Phys Chem Chem Phys 2016; 18:20234-50. [DOI: 10.1039/c6cp03834d] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.
Collapse
Affiliation(s)
- Nanna Holmgaard List
- Department of Physics, Chemistry and Pharmacy
- University of Southern Denmark
- 5230 Odense M
- Denmark
| | | | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy
- University of Southern Denmark
- 5230 Odense M
- Denmark
| |
Collapse
|
27
|
Hartman JD, Monaco S, Schatschneider B, Beran GJO. Fragment-based (13)C nuclear magnetic resonance chemical shift predictions in molecular crystals: An alternative to planewave methods. J Chem Phys 2015; 143:102809. [PMID: 26374002 DOI: 10.1063/1.4922649] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We assess the quality of fragment-based ab initio isotropic (13)C chemical shift predictions for a collection of 25 molecular crystals with eight different density functionals. We explore the relative performance of cluster, two-body fragment, combined cluster/fragment, and the planewave gauge-including projector augmented wave (GIPAW) models relative to experiment. When electrostatic embedding is employed to capture many-body polarization effects, the simple and computationally inexpensive two-body fragment model predicts both isotropic (13)C chemical shifts and the chemical shielding tensors as well as both cluster models and the GIPAW approach. Unlike the GIPAW approach, hybrid density functionals can be used readily in a fragment model, and all four hybrid functionals tested here (PBE0, B3LYP, B3PW91, and B97-2) predict chemical shifts in noticeably better agreement with experiment than the four generalized gradient approximation (GGA) functionals considered (PBE, OPBE, BLYP, and BP86). A set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided based on these benchmark calculations. Statistical cross-validation procedures are used to demonstrate the robustness of these fits.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Stephen Monaco
- The Pennsylvania State University, The Eberly Campus, 2201 University Dr, Lemont Furnace, Pennsylvania 15456, USA
| | - Bohdan Schatschneider
- The Pennsylvania State University, The Eberly Campus, 2201 University Dr, Lemont Furnace, Pennsylvania 15456, USA
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
28
|
List NH, Beerepoot MTP, Olsen JMH, Gao B, Ruud K, Jensen HJA, Kongsted J. Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone. J Chem Phys 2015; 142:034119. [PMID: 25612701 DOI: 10.1063/1.4905909] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.
Collapse
Affiliation(s)
- Nanna Holmgaard List
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, Odense M, Odense DK-5230 Denmark
| | - Maarten T P Beerepoot
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Jógvan Magnus Haugaard Olsen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, Odense M, Odense DK-5230 Denmark
| | - Bin Gao
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, Odense M, Odense DK-5230 Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, Odense M, Odense DK-5230 Denmark
| |
Collapse
|
29
|
Steinmann C, Kongsted J. Electronic Energy Transfer in Polarizable Heterogeneous Environments: A Systematic Investigation of Different Quantum Chemical Approaches. J Chem Theory Comput 2015; 11:4283-93. [DOI: 10.1021/acs.jctc.5b00470] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Casper Steinmann
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| |
Collapse
|
30
|
Hartman JD, Neubauer TJ, Caulkins BG, Mueller LJ, Beran GJO. Converging nuclear magnetic shielding calculations with respect to basis and system size in protein systems. JOURNAL OF BIOMOLECULAR NMR 2015; 62:327-40. [PMID: 25993979 PMCID: PMC4512207 DOI: 10.1007/s10858-015-9947-2] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2015] [Accepted: 05/12/2015] [Indexed: 05/25/2023]
Abstract
Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van der Waals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size.
Collapse
Affiliation(s)
- Joshua D. Hartman
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, Tel.: +1-951-827-7869
| | - Thomas J. Neubauer
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, Tel.: +1-951-827-7869
| | - Bethany G. Caulkins
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, Tel.: +1-951-827-7869
| | - Leonard J. Mueller
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, Tel.: +1-951-827-7869
| | - Gregory J. O. Beran
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, Tel.: +1-951-827-7869
| |
Collapse
|
31
|
de la Lande A, Gillet N, Chen S, Salahub DR. Progress and challenges in simulating and understanding electron transfer in proteins. Arch Biochem Biophys 2015; 582:28-41. [PMID: 26116376 DOI: 10.1016/j.abb.2015.06.016] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2015] [Revised: 06/15/2015] [Accepted: 06/22/2015] [Indexed: 11/19/2022]
Abstract
This Review presents an overview of the most common numerical simulation approaches for the investigation of electron transfer (ET) in proteins. We try to highlight the merits of the different approaches but also the current limitations and challenges. The article is organized into three sections. Section 2 deals with direct simulation algorithms of charge migration in proteins. Section 3 summarizes the methods for testing the applicability of the Marcus theory for ET in proteins and for evaluating key thermodynamic quantities entering the reaction rates (reorganization energies and driving force). Recent studies interrogating the validity of the theory due to the presence of non-ergodic effects or of non-linear responses are also described. Section 4 focuses on the tunneling aspects of electron transfer. How can the electronic coupling between charge transfer states be evaluated by quantum chemistry approaches and rationalized? What interesting physics regarding the impact of protein dynamics on tunneling can be addressed? We will illustrate the different sections with examples taken from the literature to show what types of system are currently manageable with current methodologies. We also take care to recall what has been learned on the biophysics of ET within proteins thanks to the advent of atomistic simulations.
Collapse
Affiliation(s)
- Aurélien de la Lande
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France.
| | - Natacha Gillet
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France
| | - Shufeng Chen
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France
| | - Dennis R Salahub
- Department of Chemistry, CMS - Centre for Molecular Simulation and IQST - Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive N.W., Calgary, Alberta T2N 1N4, Canada.
| |
Collapse
|
32
|
Olsen JMH, List NH, Kristensen K, Kongsted J. Accuracy of Protein Embedding Potentials: An Analysis in Terms of Electrostatic Potentials. J Chem Theory Comput 2015; 11:1832-42. [DOI: 10.1021/acs.jctc.5b00078] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Laboratory
of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Nanna Holmgaard List
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Kasper Kristensen
- qLEAP
Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| | - Jacob Kongsted
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
33
|
Olsen JMH, Steinmann C, Ruud K, Kongsted J. Polarizable Density Embedding: A New QM/QM/MM-Based Computational Strategy. J Phys Chem A 2015; 119:5344-55. [DOI: 10.1021/jp510138k] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Laboratory of
Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Casper Steinmann
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Kenneth Ruud
- Centre for Theoretical
and Computational Chemistry, Department of Chemistry, University of Tromsø—The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Jacob Kongsted
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| |
Collapse
|