1
|
Li Z, Zhao B, Zhang H, Zhang Y. High-throughput calculations and machine learning modeling of 17O NMR in non-magnetic oxides. Faraday Discuss 2024. [PMID: 39258887 DOI: 10.1039/d4fd00128a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
The only NMR-active oxygen isotope, oxygen-17 (17O), serves as a sensitive probe due to its large chemical shift range, the electric field gradient at the oxygen site, and the quadrupolar interaction. Consequently, 17O solid-state NMR offers unique insights into local structures and finds significant applications in the studies of disorder, reactivity, and host-guest chemistry. Despite recent advances in sensitivity enhancement, isotopic labeling, and NMR crystallography, the application of 17O solid-state NMR is still hindered by low natural abundance, costly enrichment, and challenges in handling spectrum signals. Density functional theory calculations and machine learning techniques offer an alternative approach to mapping the local crystal structures to NMR parameters. However, the lack of high-quality data remains a challenge, despite the establishment of some datasets. In this study, we implement and execute a high-throughput workflow combining AiiDA and CASTEP to evaluate the NMR parameters. Focusing on non-magnetic oxides, we have chosen over 7100 binary, ternary, and quaternary compounds from the Materials Project database and performed calculations. Furthermore, using various descriptors for the local crystalline environments, we model the 17O NMR parameters using machine learning techniques, further enhancing our ability to predict and understand 17O NMR parameters in oxide crystals.
Collapse
Affiliation(s)
- Zhiyuan Li
- TU Darmstadt, Otto-Berndt-Straße 3, 64287 Darmstadt, Germany.
| | - Bo Zhao
- TU Darmstadt, Otto-Berndt-Straße 3, 64287 Darmstadt, Germany.
| | - Hongbin Zhang
- TU Darmstadt, Otto-Berndt-Straße 3, 64287 Darmstadt, Germany.
| | - Yixuan Zhang
- TU Darmstadt, Otto-Berndt-Straße 3, 64287 Darmstadt, Germany.
| |
Collapse
|
2
|
Drontschenko V, Bangerter FH, Ochsenfeld C. Analytical Second-Order Properties for the Random Phase Approximation: Nuclear Magnetic Resonance Shieldings. J Chem Theory Comput 2023; 19:7542-7554. [PMID: 37863033 DOI: 10.1021/acs.jctc.3c00542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2023]
Abstract
A method for the analytical computation of nuclear magnetic resonance (NMR) shieldings within the direct random phase approximation (RPA) is presented. As a starting point, we use the RPA ground-state energy expression within the resolution-of-the-identity approximation in the atomic-orbital formalism. As has been shown in a recent benchmark study using numerical second derivatives [Glasbrenner, M. J. Chem. Theory Comput. 2022, 18, 192], RPA based on a Hartree-Fock reference shows accuracies comparable to coupled cluster singles and doubles (CCSD) for NMR chemical shieldings. Together with the much lower computational cost of RPA, it has emerged as an accurate method for the computation of NMR shieldings. Therefore, we aim to extend the applicability of RPA NMR to larger systems by introducing analytical second-order derivatives, making it a viable method for the accurate and efficient computation of NMR chemical shieldings.
Collapse
Affiliation(s)
- Viktoria Drontschenko
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Felix H Bangerter
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
3
|
Glasbrenner M, Vogler S, Ochsenfeld C. Efficient low-scaling computation of NMR shieldings at the second-order Møller-Plesset perturbation theory level with Cholesky-decomposed densities and an attenuated Coulomb metric. J Chem Phys 2021; 155:224107. [PMID: 34911319 DOI: 10.1063/5.0069956] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for the computation of nuclear magnetic resonance (NMR) shieldings with second-order Møller-Plesset perturbation theory (MP2) is presented which allows to efficiently compute the entire set of shieldings for a given molecular structure. The equations are derived using Laplace-transformed atomic orbital second-order Møller-Plesset perturbation theory as a starting point. The Z-vector approach is employed for minimizing the number of coupled-perturbed self-consistent-field equations that need to be solved. In addition, the method uses the resolution-of-the-identity approximation with an attenuated Coulomb metric and Cholesky decomposition of pseudo-density matrices. The sparsity in the three-center integrals is exploited with sparse linear algebra approaches, leading to reduced computational cost and memory demands. Test calculations show that the deviations from NMR shifts obtained with canonical MP2 are small if appropriate thresholds are used. The performance of the method is illustrated in calculations on DNA strands and on glycine chains with up to 283 atoms and 2864 basis functions.
Collapse
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Sigurd Vogler
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| |
Collapse
|
4
|
Glasbrenner M, Graf D, Ochsenfeld C. Benchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR Shieldings. J Chem Theory Comput 2021; 18:192-205. [PMID: 34898213 DOI: 10.1021/acs.jctc.1c00866] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A method for computing NMR shieldings with the direct random phase approximation (RPA) and the closely related σ-functionals [Trushin, E.; Thierbach, A.; Görling, A. Toward chemical accuracy at low computational cost: density functional theory with σ-functionals for the correlation energy. J. Chem. Phys. 2021, 154, 014104] is presented, which is based on a finite-difference approach. The accuracy is evaluated in benchmark calculations using high-quality coupled cluster values as a reference. Our results show that the accuracy of the computed NMR shieldings using direct RPA is strongly dependent on the density functional theory reference orbitals and improves with increasing amounts of exact Hartree-Fock exchange in the functional. NMR shieldings computed with direct RPA using a Hartree-Fock reference are significantly more accurate than MP2 shieldings and comparable to CCSD shieldings. Also, the basis set convergence is analyzed and it is shown that at least triple-zeta basis sets are required for reliable results.
Collapse
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany.,Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
5
|
Stauch T, Ganoe B, Wong J, Lee J, Rettig A, Liang J, Li J, Epifanovsky E, Head-Gordon T, Head-Gordon M. Molecular magnetisabilities computed via finite fields: assessing alternatives to MP2 and revisiting magnetic exaltations in aromatic and antiaromatic species. Mol Phys 2021; 119. [DOI: 10.1080/00268976.2021.1990426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Tim Stauch
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Brad Ganoe
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jonathan Wong
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Joonho Lee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jie Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | | | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
6
|
Jinger R, Fliegl H, Bast R, Dimitrova M, Lehtola S, Sundholm D. Spatial Contributions to Nuclear Magnetic Shieldings. J Phys Chem A 2021; 125:1778-1786. [PMID: 33605721 PMCID: PMC8023705 DOI: 10.1021/acs.jpca.0c10884] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2020] [Revised: 02/05/2021] [Indexed: 01/09/2023]
Abstract
We develop a methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities which are calculated from the current density via the Biot-Savart relation. Atomic contributions to nuclear magnetic shielding constants can be estimated within our framework with a Becke partitioning scheme. The new features have been implemented in the GIMIC program and are applied in this work to the study of the 1H and 13C nuclear magnetic shieldings in benzene (C6H6) and cyclobutadiene (C4H4). The new methodology allows a visual inspection of the spatial origins of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constant of a single nucleus, something which has not been hitherto easily accomplished. Analysis of the shielding densities shows that diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions, as the shielding or deshielding is determined by the direction of the current-density flux with respect to the studied nucleus instead of the tropicity. Becke partitioning of the magnetic shieldings shows that the magnetic shielding contributions mainly originate from the studied atom and its nearest neighbors, confirming the localized character of nuclear magnetic shieldings.
Collapse
Affiliation(s)
- Rahul
Kumar Jinger
- Indian
Institute of Science Education and Research, Dr. Homi Bhabha Road, Pashan, Pune 411008, India
| | - Heike Fliegl
- Karlsruhe
Institute of Technology, Institute of Nanotechnology, Hermann-von-Helmholtz Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany
| | - Radovan Bast
- Department
Information Technology, UiT Arctic University
Norway, Tromsø, Norway
| | - Maria Dimitrova
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Susi Lehtola
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
| | - Dage Sundholm
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| |
Collapse
|
7
|
Shi M, Jin X, Wan Z, He X. Automated fragmentation quantum mechanical calculation of 13C and 1H chemical shifts in molecular crystals. J Chem Phys 2021; 154:064502. [PMID: 33588539 DOI: 10.1063/5.0039115] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach was applied to calculate the 13C and 1H nuclear magnetic resonance (NMR) chemical shifts in molecular crystals. Two benchmark sets of molecular crystals were selected to calculate the NMR chemical shifts. Systematic investigation was conducted to examine the convergence of AF-QM/MM calculations and the impact of various density functionals with different basis sets on the NMR chemical shift prediction. The result demonstrates that the calculated NMR chemical shifts are close to convergence when the distance threshold for the QM region is larger than 3.5 Å. For 13C chemical shift calculations, the mPW1PW91 functional is the best density functional among the functionals chosen in this study (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB98, and OPBE), while the OB98 functional is more suitable for the 1H NMR chemical shift prediction of molecular crystals. Moreover, with the B3LYP functional, at least a triple-ζ basis set should be utilized to accurately reproduce the experimental 13C and 1H chemical shifts. The employment of diffuse basis functions will further improve the accuracy for 13C chemical shift calculations, but not for the 1H chemical shift prediction. We further proposed a fragmentation scheme of dividing the central molecule into smaller fragments. By comparing with the results of the fragmentation scheme using the entire central molecule as the core region, the AF-QM/MM calculations with the fragmented central molecule can not only achieve accurate results but also reduce the computational cost. Therefore, the AF-QM/MM approach is capable of predicting the 13C and 1H NMR chemical shifts for molecular crystals accurately and effectively, and could be utilized for dealing with more complex periodic systems such as macromolecular polymers and biomacromolecules. The AF-QM/MM program for molecular crystals is available at https://github.com/shiman1995/NMR.
Collapse
Affiliation(s)
- Man Shi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Zheng Wan
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|
8
|
Dittmer A, Stoychev GL, Maganas D, Auer AA, Neese F. Computation of NMR Shielding Constants for Solids Using an Embedded Cluster Approach with DFT, Double-Hybrid DFT, and MP2. J Chem Theory Comput 2020; 16:6950-6967. [PMID: 32966067 PMCID: PMC7659039 DOI: 10.1021/acs.jctc.0c00067] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
In
this work, we explore the accuracy of post-Hartree–Fock
(HF) methods and double-hybrid density functional theory (DFT) for
the computation of solid-state NMR chemical shifts. We apply an embedded
cluster approach and investigate the convergence with cluster size
and embedding for a series of inorganic solids with long-range electrostatic
interactions. In a systematic study, we discuss the cluster design,
the embedding procedure, and basis set convergence using gauge-including
atomic orbital (GIAO) NMR calculations at the DFT and MP2 levels of
theory. We demonstrate that the accuracy obtained for the prediction
of NMR chemical shifts, which can be achieved for molecular systems,
can be carried over to solid systems. An appropriate embedded cluster
approach allows one to apply methods beyond standard DFT even for
systems for which long-range electrostatic effects are important. We find that an embedded
cluster should include at least one sphere of explicit neighbors around
the nuclei of interest, given that a sufficiently large point charge
and boundary effective potential embedding is applied. Using the pcSseg-3
basis set and GIAOs for the computation of nuclear shielding constants,
accuracies of 1.6 ppm for 7Li, 1.5 ppm for 23Na, and 5.1 ppm for 39K as well as 9.3 ppm for 19F, 6.5 ppm for 35Cl, 7.4 ppm for 79Br, and
7.5 ppm for 25Mg as well as 3.8 ppm for 67Zn
can be achieved with MP2. Comparing various DFT functionals with HF
and MP2, we report the superior quality of results for methods that
include post-HF correlation like MP2 and double-hybrid DFT.
Collapse
Affiliation(s)
- Anneke Dittmer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Dimitrios Maganas
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
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
|
Repisky M, Komorovsky S, Kadek M, Konecny L, Ekström U, Malkin E, Kaupp M, Ruud K, Malkina OL, Malkin VG. ReSpect: Relativistic spectroscopy DFT program package. J Chem Phys 2020; 152:184101. [DOI: 10.1063/5.0005094] [Citation(s) in RCA: 47] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Stanislav Komorovsky
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dubravska cesta 9, SK-84536 Bratislava, Slovakia
| | - Marius Kadek
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Lukas Konecny
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Ulf Ekström
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, N-0315 Oslo, Norway
| | - Elena Malkin
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Martin Kaupp
- Technische Universität Berlin, Institute of Chemistry, Strasse des 17 Juni 135, D-10623 Berlin, Germany
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Olga L. Malkina
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dubravska cesta 9, SK-84536 Bratislava, Slovakia
| | - Vladimir G. Malkin
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dubravska cesta 9, SK-84536 Bratislava, Slovakia
| |
Collapse
|
12
|
Zhang Y, Suo B, Wang Z, Zhang N, Li Z, Lei Y, Zou W, Gao J, Peng D, Pu Z, Xiao Y, Sun Q, Wang F, Ma Y, Wang X, Guo Y, Liu W. BDF: A relativistic electronic structure program package. J Chem Phys 2020; 152:064113. [DOI: 10.1063/1.5143173] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yong Zhang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Bingbing Suo
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Zikuan Wang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Yibo Lei
- Key Laboratory of Synthetic and Natural Functional Molecule Chemistry of Ministry of Education, Shaanxi Key Laboratory of Physico-Inorganic Chemistry, College of Chemistry and Materials Science, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Wenli Zou
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, Hubei 430070, People’s Republic of China
| | - Daoling Peng
- College of Chemistry, South China Normal University, Guangzhou, Guangdong 510006, People’s Republic of China
| | - Zhichen Pu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Qiming Sun
- Tencent America LLC, Palo Alto, California 94306, USA
| | - Fan Wang
- Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Institute of Atomic and Molecular Physics, Sichuan University, Chengdu, Sichuan 610065, People’s Republic of China
| | - Yongtao Ma
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Xiaopeng Wang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| |
Collapse
|
13
|
Gräfenstein J. Efficient calculation of NMR isotopic shifts: Difference-dedicated vibrational perturbation theory. J Chem Phys 2019; 151:244120. [PMID: 31893883 DOI: 10.1063/1.5134538] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
We present difference-dedicated second-order vibrational perturbation theory (VPT2) as an efficient method for the computation of nuclear magnetic resonance (NMR) isotopic shifts, which reflect the geometry dependence of the NMR property in combination with different vibration patterns of two isotopologues. Conventional calculations of isotopic shifts, e.g., by standard VPT2, require scanning the geometry dependence over the whole molecule, which becomes expensive rapidly as the molecule size increases. In DD-VPT2, this scan can be restricted to a small region around the substitution site. At the heart of DD-VPT2 is a set of localized vibration modes common to the two isotopologues and designed such that the difference between the vibration patterns is caught by a small subset of them (usually fewer than 10). We tested the DD-VPT2 method for a series of molecules with increasing size and found that this method provides results with the same quality as VPT2 and in good agreement with the experiment, with computational savings up to 95% and less numerical instabilities. The method is easy to automatize and straightforward to generalize to other molecular properties.
Collapse
Affiliation(s)
- Jürgen Gräfenstein
- Department of Chemistry and Molecular Biology, University of Gothenburg, SE-412 96 Göteborg, Sweden
| |
Collapse
|
14
|
Yuan M, Zhang Y, Qu Z, Xiao Y, Liu W. Sublinear scaling quantum chemical methods for magnetic shieldings in large molecules. J Chem Phys 2019; 150:154113. [DOI: 10.1063/1.5083193] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Minghong Yuan
- Beijing National Center for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yong Zhang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Zhi Qu
- Beijing National Center for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yunlong Xiao
- Beijing National Center for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| |
Collapse
|
15
|
Glasbrenner M, Vogler S, Ochsenfeld C. Linear and sublinear scaling computation of the electronic g-tensor at the density functional theory level. J Chem Phys 2019; 150:024104. [DOI: 10.1063/1.5066266] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| | - Sigurd Vogler
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| |
Collapse
|
16
|
Stoychev GL, Auer AA, Neese F. Efficient and Accurate Prediction of Nuclear Magnetic Resonance Shielding Tensors with Double-Hybrid Density Functional Theory. J Chem Theory Comput 2018; 14:4756-4771. [DOI: 10.1021/acs.jctc.8b00624] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Georgi L. Stoychev
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr 45470, Germany
| | - Alexander A. Auer
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr 45470, Germany
| |
Collapse
|
17
|
Jensen SR, Flå T, Jonsson D, Monstad RS, Ruud K, Frediani L. Magnetic properties with multiwavelets and DFT: the complete basis set limit achieved. Phys Chem Chem Phys 2018; 18:21145-61. [PMID: 27087397 DOI: 10.1039/c6cp01294a] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Multiwavelets are emerging as an attractive alternative to traditional basis sets such as Gaussian-type orbitals and plane waves. One of their distinctive properties is the ability to reach the basis set limit (often a chimera for traditional approaches) reliably and consistently by fixing the desired precision ε. We present our multiwavelet implementation of the linear response formalism, applied to static magnetic properties, at the self-consistent field level of theory (both for Hartree-Fock and density functional theories). We demonstrate that the multiwavelets consistently improve the accuracy of the results when increasing the desired precision, yielding results that have four to five digits precision, thus providing a very useful benchmark which could otherwise only be estimated by extrapolation methods. Our results show that magnetizabilities obtained with the augmented quadruple-ζ basis (aug-cc-pCVQZ) are practically at the basis set limit, whereas absolute nuclear magnetic resonance shielding tensors are more challenging: even by making use of a standard extrapolation method, the accuracy is not substantially improved. In contrast, our results provide a benchmark that: (1) confirms the validity of the extrapolation ansatz; (2) can be used as a reference to achieve a property-specific extrapolation scheme, thus providing a means to obtain much better extrapolated results; (3) allows us to separate functional-specific errors from basis-set ones and thus to assess the level of cancellation between basis set and functional errors often exploited in density functional theory.
Collapse
Affiliation(s)
- Stig Rune Jensen
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway.
| | - Tor Flå
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway. and Department of Mathematics and Statistics, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Dan Jonsson
- High-Performance Computing Group, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Rune Sørland Monstad
- 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.
| | - Luca Frediani
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway.
| |
Collapse
|
18
|
Stoychev GL, Auer AA, Izsák R, Neese F. Self-Consistent Field Calculation of Nuclear Magnetic Resonance Chemical Shielding Constants Using Gauge-Including Atomic Orbitals and Approximate Two-Electron Integrals. J Chem Theory Comput 2018; 14:619-637. [DOI: 10.1021/acs.jctc.7b01006] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Georgi L. Stoychev
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A. Auer
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
19
|
Reiter K, Mack F, Weigend F. Calculation of Magnetic Shielding Constants with meta-GGA Functionals Employing the Multipole-Accelerated Resolution of the Identity: Implementation and Assessment of Accuracy and Efficiency. J Chem Theory Comput 2017; 14:191-197. [DOI: 10.1021/acs.jctc.7b01115] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kevin Reiter
- Institut
für Nanotechnologie, Karlsruher Institut für Technologie, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Fabian Mack
- Institut
für Physikalische Chemie, Abteilung für Theoretische
Chemie, Karlsruher Institut für Technologie, Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Florian Weigend
- Institut
für Nanotechnologie, Karlsruher Institut für Technologie, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institut
für Physikalische Chemie, Abteilung für Theoretische
Chemie, Karlsruher Institut für Technologie, Kaiserstraße 12, 76131 Karlsruhe, Germany
| |
Collapse
|
20
|
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
|
21
|
Kumar C, Kjærgaard T, Helgaker T, Fliegl H. Nuclei-selected atomic-orbital response-theory formulation for the calculation of NMR shielding tensors using density-fitting. J Chem Phys 2016; 145:234108. [DOI: 10.1063/1.4972212] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Chandan Kumar
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Thomas Kjærgaard
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Trygve Helgaker
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Heike Fliegl
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| |
Collapse
|
22
|
Luenser A, Kussmann J, Ochsenfeld C. Computation of indirect nuclear spin–spin couplings with reduced complexity in pure and hybrid density functional approximations. J Chem Phys 2016; 145:124103. [DOI: 10.1063/1.4962260] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Arne Luenser
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU Munich) and Center for Integrated Protein Science Munich (CIPSM, LMU Munich), Butenandtstr. 5–13, D-81377 Munich, Germany
| | - Jörg Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU Munich) and Center for Integrated Protein Science Munich (CIPSM, LMU Munich), Butenandtstr. 5–13, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU Munich) and Center for Integrated Protein Science Munich (CIPSM, LMU Munich), Butenandtstr. 5–13, D-81377 Munich, Germany
| |
Collapse
|
23
|
Sumowski CV, Hanni M, Schweizer S, Ochsenfeld C. Sensitivity of ab Initio vs Empirical Methods in Computing Structural Effects on NMR Chemical Shifts for the Example of Peptides. J Chem Theory Comput 2015; 10:122-33. [PMID: 26579896 DOI: 10.1021/ct400713t] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The structural sensitivity of NMR chemical shifts as computed by quantum chemical methods is compared to a variety of empirical approaches for the example of a prototypical peptide, the 38-residue kaliotoxin KTX comprising 573 atoms. Despite the simplicity of empirical chemical shift prediction programs, the agreement with experimental results is rather good, underlining their usefulness. However, we show in our present work that they are highly insensitive to structural changes, which renders their use for validating predicted structures questionable. In contrast, quantum chemical methods show the expected high sensitivity to structural and electronic changes. This appears to be independent of the quantum chemical approach or the inclusion of solvent effects. For the latter, explicit solvent simulations with increasing number of snapshots were performed for two conformers of an eight amino acid sequence. In conclusion, the empirical approaches neither provide the expected magnitude nor the patterns of NMR chemical shifts determined by the clearly more costly ab initio methods upon structural changes. This restricts the use of empirical prediction programs in studies where peptide and protein structures are utilized for the NMR chemical shift evaluation such as in NMR refinement processes, structural model verifications, or calculations of NMR nuclear spin relaxation rates.
Collapse
Affiliation(s)
- Chris Vanessa Sumowski
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstr. 7, D-81377 Munich, Germany and Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Matti Hanni
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstr. 7, D-81377 Munich, Germany and Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Sabine Schweizer
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstr. 7, D-81377 Munich, Germany and Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstr. 7, D-81377 Munich, Germany and Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstr. 5-13, D-81377 Munich, Germany
| |
Collapse
|
24
|
Hartman JD, Beran GJO. Fragment-Based Electronic Structure Approach for Computing Nuclear Magnetic Resonance Chemical Shifts in Molecular Crystals. J Chem Theory Comput 2015; 10:4862-72. [PMID: 26584373 DOI: 10.1021/ct500749h] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
First-principles chemical shielding tensor predictions play a critical role in studying molecular crystal structures using nuclear magnetic resonance. Fragment-based electronic structure methods have dramatically improved the ability to model molecular crystal structures and energetics using high-level electronic structure methods. Here, a many-body expansion fragment approach is applied to the calculation of chemical shielding tensors in molecular crystals. First, the impact of truncating the many-body expansion at different orders and the role of electrostatic embedding are examined on a series of molecular clusters extracted from molecular crystals. Second, the ability of these techniques to assign three polymorphic forms of the drug sulfanilamide to the corresponding experimental (13)C spectra is assessed. This challenging example requires discriminating among spectra whose (13)C chemical shifts differ by only a few parts per million (ppm) across the different polymorphs. Fragment-based PBE0/6-311+G(2d,p) level chemical shielding predictions correctly assign these three polymorphs and reproduce the sulfanilamide experimental (13)C chemical shifts with 1 ppm accuracy. The results demonstrate that fragment approaches are competitive with the widely used gauge-invariant projector augmented wave (GIPAW) periodic density functional theory calculations.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California , Riverside, California 92521, United States
| | - Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Zhu T, Zhang JZH, He X. Quantum calculation of protein NMR chemical shifts based on the automated fragmentation method. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2015; 827:49-70. [PMID: 25387959 DOI: 10.1007/978-94-017-9245-5_5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
The performance of quantum mechanical methods on the calculation of protein NMR chemical shifts is reviewed based on the recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. By using the Poisson-Boltzmann (PB) model and first solvation water molecules, the influence of solvent effect is also discussed. Benefiting from the fragmentation algorithm, the AF-QM/MM approach is computationally efficient, linear-scaling with a low pre-factor, and thus can be applied to routinely calculate the ab initio NMR chemical shifts for proteins of any size. The results calculated using Density Functional Theory (DFT) show that when the solvent effect is included, this method can accurately reproduce the experimental ¹H NMR chemical shifts, while the ¹³C NMR chemical shifts are less affected by the solvent. However, although the inclusion of solvent effect shows significant improvement for ¹⁵N chemical shifts, the calculated values still have large deviations from the experimental observations. Our study further demonstrates that AF-QM/MM calculated results accurately reflect the dependence of ¹³C(α) NMR chemical shifts on the secondary structure of proteins, and the calculated ¹H chemical shift can be utilized to discriminate the native structure of proteins from decoys.
Collapse
Affiliation(s)
- Tong Zhu
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China
| | | | | |
Collapse
|
28
|
Kussmann J, Luenser A, Beer M, Ochsenfeld C. A reduced-scaling density matrix-based method for the computation of the vibrational Hessian matrix at the self-consistent field level. J Chem Phys 2015; 142:094101. [DOI: 10.1063/1.4908131] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jörg Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Arne Luenser
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Matthias Beer
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| |
Collapse
|
29
|
Maurer M, Ochsenfeld C. Spin component-scaled second-order Møller-Plesset perturbation theory for calculating NMR shieldings. J Chem Theory Comput 2015; 11:37-44. [PMID: 26574201 DOI: 10.1021/ct5007295] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Spin component-scaled and scaled opposite-spin second-order Møller-Plesset perturbation approaches (SCS-MP2 and SOS-MP2) are introduced for calculating NMR chemical shifts in analogy to the well-established scaled approaches for MP2 energies. Gauge-including atomic orbitals (GIAO) are employed throughout this work. The GIAO-SCS-MP2 and GIAO-SOS-MP2 methods typically show superior performance to nonscaled MP2 and are closer to the coupled-cluster singles doubles perturbative triples (CCSD(T))/cc-pVQZ reference values. In addition, the pragmatic use of mixed basis sets for the Hartree-Fock and the correlated part of NMR chemical shift calculations is shown to be beneficial.
Collapse
Affiliation(s)
- Marina Maurer
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstrasse 7, D-81377 Munich, Germany.,Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstrasse 5-13, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU) , Butenandtstrasse 7, D-81377 Munich, Germany.,Center for Integrated Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU) , Butenandtstrasse 5-13, D-81377 Munich, Germany
| |
Collapse
|
30
|
Reid DM, Collins MA. Calculating nuclear magnetic resonance shieldings using systematic molecular fragmentation by annihilation. Phys Chem Chem Phys 2015; 17:5314-20. [DOI: 10.1039/c4cp05116e] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Systematic fragmentation accurately predicts theoretical chemical shieldings.
Collapse
Affiliation(s)
- David M. Reid
- Research School of Chemistry
- The Australian National University
- ACT
- Australia
| | - Michael A. Collins
- Research School of Chemistry
- The Australian National University
- ACT
- Australia
| |
Collapse
|
31
|
Flaig D, Maurer M, Hanni M, Braunger K, Kick L, Thubauville M, Ochsenfeld C. Benchmarking Hydrogen and Carbon NMR Chemical Shifts at HF, DFT, and MP2 Levels. J Chem Theory Comput 2014; 10:572-8. [DOI: 10.1021/ct400780f] [Citation(s) in RCA: 126] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Denis Flaig
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Marina Maurer
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Matti Hanni
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Katharina Braunger
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Leonhard Kick
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Matthias Thubauville
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr.
7, D-81377 Munich, Germany
- Center for Integrated Protein
Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| |
Collapse
|
32
|
Schmidt TC, Paasche A, Grebner C, Ansorg K, Becker J, Lee W, Engels B. QM/MM investigations of organic chemistry oriented questions. Top Curr Chem (Cham) 2014; 351:25-101. [PMID: 22392477 DOI: 10.1007/128_2011_309] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.
Collapse
Affiliation(s)
- Thomas C Schmidt
- Institut für Phys. und Theor. Chemie, Emil-Fischer-Strasse 42, Campus Hubland Nord, 97074, Würzburg, Germany
| | | | | | | | | | | | | |
Collapse
|
33
|
Maurer M, Ochsenfeld C. A linear- and sublinear-scaling method for calculating NMR shieldings in atomic orbital-based second-order Møller-Plesset perturbation theory. J Chem Phys 2013; 138:174104. [PMID: 23656111 DOI: 10.1063/1.4801084] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An atomic-orbital (AO) based formulation for calculating nuclear magnetic resonance chemical shieldings at the second-order Møller-Plesset perturbation theory level is introduced, which provides a basis for reducing the scaling of the computational effort with the molecular size from the fifth power to linear and for a specific nucleus to sublinear. The latter sublinear scaling in the rate-determining steps becomes possible by avoiding global perturbations with respect to the magnetic field and by solving for quantities that involve the local nuclear magnetic spin perturbation instead. For avoiding the calculation of the second-order perturbed density matrix, we extend our AO-based reformulation of the Z-vector method within a density matrix-based scheme. Our pilot implementation illustrates the fast convergence with respect to the required number of Laplace points and the asymptotic scaling behavior in the rate-determining steps.
Collapse
Affiliation(s)
- Marina Maurer
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich, Butenandtstr. 7, D-81377 Munich, Germany
| | | |
Collapse
|
34
|
Flaig D, Ochsenfeld C. An extrapolation method for the efficient calculation of molecular response properties within Born-Oppenheimer molecular dynamics. Phys Chem Chem Phys 2013; 15:9392-6. [PMID: 23666498 DOI: 10.1039/c3cp50204j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The calculation of molecular response properties in dynamic molecular systems is a major challenge that requires sampling over many steps of, e.g., Born-Oppenheimer molecular dynamics (BO-MD) simulations. We present an extrapolation scheme to accelerate such calculations for multiple steps within BO-MD trajectories or equivalently within other sampling methods of conformational space. The extrapolation scheme is related to the one introduced by Pulay and Fogarasi [Chem. Phys. Lett., 2004, 386, 272] for self-consistent field (SCF) energy calculations. We extend the extrapolation to the quantities within our density matrix-based Laplace-transformed coupled perturbed SCF (DL-CPSCF) method that allows for linear-scaling calculations of response properties for large molecular systems. Here, we focus on the example of calculating NMR chemical shifts for which the number of required DL-CPSCF iterations reduces by roughly 40-70%.
Collapse
Affiliation(s)
- Denis Flaig
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | | |
Collapse
|
35
|
Dutt S, Wilch C, Gersthagen T, Talbiersky P, Bravo-Rodriguez K, Hanni M, Sánchez-García E, Ochsenfeld C, Klärner FG, Schrader T. Molecular Tweezers with Varying Anions: A Comparative Study. J Org Chem 2013; 78:6721-34. [DOI: 10.1021/jo4009673] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Som Dutt
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| | - Constanze Wilch
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| | - Thomas Gersthagen
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| | - Peter Talbiersky
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| | - Kenny Bravo-Rodriguez
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1,
45470 Mülheim an der Ruhr, Germany
| | - Matti Hanni
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany and Center for Integrated
Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, 81377
Munich, Germany
| | - Elsa Sánchez-García
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1,
45470 Mülheim an der Ruhr, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry,
Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany and Center for Integrated
Protein Science (CIPSM) at the Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, 81377
Munich, Germany
| | - Frank-Gerrit Klärner
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| | - Thomas Schrader
- Department of Chemistry, University of Duisburg-Essen, Universitätsstr.
7, 45117 Essen, Germany
| |
Collapse
|
36
|
Zhu T, Zhang JZH, He X. Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model. J Chem Theory Comput 2013; 9:2104-14. [PMID: 26583557 DOI: 10.1021/ct300999w] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We have performed a density functional theory (DFT) calculation of the amide proton NMR chemical shift in proteins using a recently developed automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. Systematic investigation was carried out to examine the influence of explicit solvent molecules, cooperative hydrogen bonding effects, density functionals, size of the basis sets, and the local geometry of proteins on calculated chemical shifts. Our result demonstrates that the predicted amide proton ((1)HN) NMR chemical shift in explicit solvent shows remarkable improvement over that calculated with the implicit solvation model. The cooperative hydrogen bonding effect is also shown to improve the accuracy of (1)HN chemical shifts. Furthermore, we found that the OPBE exchange-correlation functional is the best density functional for the prediction of protein (1)HN chemical shifts among a selective set of DFT methods (namely, B3LYP, B3PW91, M062X, M06L, mPW1PW91, OB98, OPBE), and the locally dense basis set of 6-311++G**/4-31G* is shown to be sufficient for (1)HN chemical shift calculation. By taking ensemble averaging into account, (1)HN chemical shifts calculated by the AF-QM/MM approach can be used to validate the performance of various force fields. Our study underscores that the electronic polarization of protein is of critical importance to stabilizing hydrogen bonding, and the AF-QM/MM method is able to describe the local chemical environment in proteins more accurately than most widely used empirical models.
Collapse
Affiliation(s)
- Tong Zhu
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062
| |
Collapse
|
37
|
Hollman DS, Wilke JJ, Schaefer HF. Explicitly correlated atomic orbital basis second order Møller–Plesset theory. J Chem Phys 2013; 138:064107. [DOI: 10.1063/1.4790582] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Kussmann J, Beer M, Ochsenfeld C. Linear-scaling self-consistent field methods for large molecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1138] [Citation(s) in RCA: 69] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
39
|
Maurer SA, Lambrecht DS, Kussmann J, Ochsenfeld C. Efficient distance-including integral screening in linear-scaling Møller-Plesset perturbation theory. J Chem Phys 2013; 138:014101. [DOI: 10.1063/1.4770502] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
40
|
Loibl S, Schütz M. NMR shielding tensors for density fitted local second-order Møller-Plesset perturbation theory using gauge including atomic orbitals. J Chem Phys 2013; 137:084107. [PMID: 22938218 DOI: 10.1063/1.4744102] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An efficient method for the calculation of nuclear magnetic resonance (NMR) shielding tensors is presented, which treats electron correlation at the level of second-order Mo̸ller-Plesset perturbation theory. It uses spatially localized functions to span occupied and virtual molecular orbital spaces, respectively, which are expanded in a basis of gauge including atomic orbitals (GIAOs or London atomic orbitals). Doubly excited determinants are restricted to local subsets of the virtual space and pair energies with an interorbital distance beyond a certain threshold are omitted. Furthermore, density fitting is employed to factorize the electron repulsion integrals. Ordinary Gaussians are employed as fitting functions. It is shown that the errors in the resulting NMR shielding constant, introduced (i) by the local approximation and (ii) by density fitting, are very small or even negligible. The capabilities of the new program are demonstrated by calculations on some extended molecular systems, such as the cyclobutane pyrimidine dimer photolesion with adjacent nucleobases in the native intrahelical DNA double strand (ATTA sequence). Systems of that size were not accessible to correlated ab initio calculations of NMR spectra before. The presented method thus opens the door to new and interesting applications in this area.
Collapse
Affiliation(s)
- Stefan Loibl
- Institute of Physical and Theoretical Chemistry, University of Regensburg, Universitätsstraße 31, D-93040 Regensburg, Germany
| | | |
Collapse
|
41
|
Exner TE, Frank A, Onila I, Möller HM. Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 3. Conformational Sampling and Explicit Solvents Model. J Chem Theory Comput 2012; 8:4818-27. [PMID: 26605634 DOI: 10.1021/ct300701m] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fragment-based quantum chemical calculations are able to accurately calculate NMR chemical shifts even for very large molecules like proteins. But even with systematic optimization of the level of theory and basis sets as well as the use of implicit solvents models, some nuclei like polar protons and nitrogens suffer from poor predictions. Two properties of the real system, strongly influencing the experimental chemical shifts but almost always neglected in the calculations, will be discussed here in great detail: (1) conformational averaging and (2) interactions with first-shell solvent molecules. Classical molecular dynamics simulations in explicit water were carried out for obtaining a representative ensemble including the arrangement of neighboring solvent molecules, which was then subjected to quantum chemical calculations. We could demonstrate with the small test system N-methyl acetamide (NMA) that the calculated chemical shifts show immense variations of up to 6 ppm and 50 ppm for protons and nitrogens, respectively, depending on the snapshot taken from a classical molecular dynamics simulation. Applying the same approach to the HA2 domain of the influenza virus glycoprotein hemagglutinin, a 32-amino-acid-long polypeptide, and comparing averaged values to the experiment, chemical shifts of nonpolar protons and carbon atoms in proteins were calculated with unprecedented accuracy. Additionally, the mean absolute error could be reduced by a factor of 2.43 for polar protons, and reasonable correlations were obtained for nitrogen and carbonyl carbon in contrast to all other studies published so far.
Collapse
Affiliation(s)
- Thomas E Exner
- Department of Chemistry and Zukunftskolleg, University of Konstanz, 78457 Konstanz, Germany.,Theoretical Medicinal Chemistry and Biophysics, Institute of Pharmacy, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Andrea Frank
- Department of Chemistry and Zukunftskolleg, University of Konstanz, 78457 Konstanz, Germany
| | - Ionut Onila
- Theoretical Medicinal Chemistry and Biophysics, Institute of Pharmacy, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Heiko M Möller
- Department of Chemistry and Zukunftskolleg, University of Konstanz, 78457 Konstanz, Germany
| |
Collapse
|
42
|
Flaig D, Beer M, Ochsenfeld C. Convergence of Electronic Structure with the Size of the QM Region: Example of QM/MM NMR Shieldings. J Chem Theory Comput 2012; 8:2260-71. [DOI: 10.1021/ct300036s] [Citation(s) in RCA: 98] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Denis Flaig
- Theoretical
Chemistry, University of Munich (LMU),
Butenandtstrasse 7, D-81377 München, Germany
| | - Matthias Beer
- Theoretical
Chemistry, University of Munich (LMU),
Butenandtstrasse 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Theoretical
Chemistry, University of Munich (LMU),
Butenandtstrasse 7, D-81377 München, Germany
| |
Collapse
|
43
|
Maurer SA, Lambrecht DS, Flaig D, Ochsenfeld C. Distance-dependent Schwarz-based integral estimates for two-electron integrals: Reliable tightness vs. rigorous upper bounds. J Chem Phys 2012; 136:144107. [DOI: 10.1063/1.3693908] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Helgaker T, Coriani S, Jørgensen P, Kristensen K, Olsen J, Ruud K. Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations. Chem Rev 2012; 112:543-631. [DOI: 10.1021/cr2002239] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trygve Helgaker
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy
| | - Poul Jørgensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|