1
|
Unzueta PA, Greenwell CS, Beran GJO. Predicting Density Functional Theory-Quality Nuclear Magnetic Resonance Chemical Shifts via Δ-Machine Learning. J Chem Theory Comput 2021; 17:826-840. [DOI: 10.1021/acs.jctc.0c00979] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Pablo A. Unzueta
- Department of Chemistry, University of California, Riverside, Riverside, California 92521, United States
| | - Chandler S. Greenwell
- Department of Chemistry, University of California, Riverside, Riverside, California 92521, United States
| | - Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, Riverside, California 92521, United States
| |
Collapse
|
2
|
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
|
3
|
Hartman J, Day GM, Beran GJO. Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions. CRYSTAL GROWTH & DESIGN 2016; 16:6479-6493. [PMID: 27829821 PMCID: PMC5095663 DOI: 10.1021/acs.cgd.6b01157] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Revised: 09/09/2016] [Indexed: 05/10/2023]
Abstract
Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state 13C and 15N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.
Collapse
Affiliation(s)
- Joshua
D. Hartman
- Department
of Chemistry, University of California, Riverside, California 92521 United States
| | - Graeme M. Day
- School
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gregory J. O. Beran
- Department
of Chemistry, University of California, Riverside, California 92521 United States
- E-mail:
| |
Collapse
|
4
|
Hartman JD, Kudla RA, Day GM, Mueller LJ, Beran GJO. Benchmark fragment-based (1)H, (13)C, (15)N and (17)O chemical shift predictions in molecular crystals. Phys Chem Chem Phys 2016; 18:21686-709. [PMID: 27431490 DOI: 10.1039/c6cp01831a] [Citation(s) in RCA: 78] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The performance of fragment-based ab initio(1)H, (13)C, (15)N and (17)O chemical shift predictions is assessed against experimental NMR chemical shift data in four benchmark sets of molecular crystals. Employing a variety of commonly used density functionals (PBE0, B3LYP, TPSSh, OPBE, PBE, TPSS), we explore the relative performance of cluster, two-body fragment, and combined cluster/fragment models. The hybrid density functionals (PBE0, B3LYP and TPSSh) generally out-perform their generalized gradient approximation (GGA)-based counterparts. (1)H, (13)C, (15)N, and (17)O isotropic chemical shifts can be predicted with root-mean-square errors of 0.3, 1.5, 4.2, and 9.8 ppm, respectively, using a computationally inexpensive electrostatically embedded two-body PBE0 fragment model. Oxygen chemical shieldings prove particularly sensitive to local many-body effects, and using a combined cluster/fragment model instead of the simple two-body fragment model decreases the root-mean-square errors to 7.6 ppm. These fragment-based model errors compare favorably with GIPAW PBE ones of 0.4, 2.2, 5.4, and 7.2 ppm for the same (1)H, (13)C, (15)N, and (17)O test sets. Using these benchmark calculations, a set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided and their robustness assessed using statistical cross-validation. We demonstrate the utility of these approaches and the reported scaling parameters on applications to 9-tert-butyl anthracene, several histidine co-crystals, benzoic acid and the C-nitrosoarene SnCl2(CH3)2(NODMA)2.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside, California 92521, USA.
| | | | | | | | | |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Maryasin B, Olbrich M, Trauner D, Ochsenfeld C. Calculated Nuclear Magnetic Resonance Spectra of Polytwistane and Related Hydrocarbon Nanorods. J Chem Theory Comput 2015; 11:1020-6. [DOI: 10.1021/ct5011505] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Boris Maryasin
- Department of Chemistry, University of Munich (LMU), and Munich Center for Integrated Protein Science (CIPSM), Butenandtstrasse 5-13, D-81377 Munich, Germany
| | - Martin Olbrich
- Department of Chemistry, University of Munich (LMU), and Munich Center for Integrated Protein Science (CIPSM), Butenandtstrasse 5-13, D-81377 Munich, Germany
| | - Dirk Trauner
- Department of Chemistry, University of Munich (LMU), and Munich Center for Integrated Protein Science (CIPSM), Butenandtstrasse 5-13, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Department of Chemistry, University of Munich (LMU), and Munich Center for Integrated Protein Science (CIPSM), Butenandtstrasse 5-13, D-81377 Munich, Germany
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
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
|
13
|
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
|