1
|
Schattenberg C, Kaupp M. Implementation and First Evaluation of Strong-Correlation-Corrected Local Hybrid Functionals for the Calculation of NMR Shieldings and Shifts. J Phys Chem A 2024; 128:2253-2271. [PMID: 38456430 PMCID: PMC10961831 DOI: 10.1021/acs.jpca.3c08507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2023] [Revised: 02/01/2024] [Accepted: 02/09/2024] [Indexed: 03/09/2024]
Abstract
Local hybrid functionals containing strong-correlation factors (scLHs) and range-separated local hybrids (RSLHs) have been integrated into an efficient coupled-perturbed Kohn-Sham implementation for the calculation of nuclear shielding constants. Several scLHs and the ωLH22t RSLH have then been evaluated for the first time for the extended NS372 benchmark set of main-group shieldings and shifts and the TM70 benchmark of 3d transition-metal shifts. The effects of the strong-correlation corrections have been analyzed with respect to the spatial distribution of the sc-factors, which locally diminish exact-exchange admixture at certain regions in a molecule. The scLH22t, scLH23t-mBR, and scLH23t-mBR-P functionals, which contain a "damped" strong-correlation factor to retain the excellent performance of the underlying LH20t functional for weakly correlated situations, tend to make smaller corrections to shieldings and shifts than the "undamped" scLH22ta functional. While the latter functional can also deteriorate agreement with the reference data in certain weakly correlated cases, it provides overall better performance, in particular for systems where static correlation is appreciable. This pertains only to a minority of systems in the NS372 main-group test set but to many more systems in the TM70 transition-metal test set, in particular for high-oxidation-state complexes, e.g., Cr(+VI) complexes and other systems with stretched bonds. Another undamped scLH, the simpler LDA-based scLH21ct-SVWN-m, also tends to provide significant improvements in many cases. The differences between the functionals and species can be rationalized on the basis of one-dimensional plots of the strong-correlation factors, augmented by isosurface plots of the fractional orbital density (FOD). Position-dependent exact-exchange admixture is thus shown to provide substantial flexibility in treating response properties like NMR shifts for both weakly and strongly correlated systems.
Collapse
Affiliation(s)
- Caspar
Jonas Schattenberg
- Research
Unit of Structural Chemistry & Computational Biophysics, Leibniz-Forschungsinstitut für Molekulare Pharmakologie
(FMP), Robert-Roessle-Str.
10, 13125 Berlin, Germany
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
2
|
Kleine Büning JB, Grimme S. Computation of CCSD(T)-Quality NMR Chemical Shifts via Δ-Machine Learning from DFT. J Chem Theory Comput 2023. [PMID: 37262324 DOI: 10.1021/acs.jctc.3c00165] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
NMR spectroscopy undoubtedly plays a central role in determining molecular structures across different chemical disciplines, and the accurate computational prediction of NMR parameters is highly desirable. In this work, a new Δ-machine learning approach is presented to correct DFT-computed NMR chemical shifts using input features from the calculation and in addition highly accurate reference data at the CCSD(T)/pcSseg-2 level of theory with a basis set extrapolation scheme. The model is trained on a data set containing 1000 optimized and geometrically distorted structures of small organic molecules comprising most elements of the first three periods and containing data for 7090 1H and 4230 13C NMR chemical shifts. Applied to the PBE0/pcSseg-2 method, the mean absolute deviation (MAD) on the internal NMR shift test set is reduced by 81% for 1H and 92% for 13C at virtually no additional computational cost. For 12 different DFT functional and basis set combinations, the MAD of the ML-corrected NMR shifts ranges from 0.021 to 0.039 ppm (1H) and from 0.38 to 1.07 ppm (13C). Importantly, the new method consistently outperforms the simple and widely used linear regression correction technique. This behavior is reproduced on three different external benchmark sets, confirming the generality and robustness of the correction scheme, which can easily be applied in DFT-based spectral simulations.
Collapse
Affiliation(s)
- Julius B Kleine Büning
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
3
|
Bartalucci E, Malär AA, Mehnert A, Kleine Büning JB, Günzel L, Icker M, Börner M, Wiebeler C, Meier BH, Grimme S, Kersting B, Wiegand T. Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State. Angew Chem Int Ed Engl 2023; 62:e202217725. [PMID: 36630178 DOI: 10.1002/anie.202217725] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 01/11/2023] [Accepted: 01/11/2023] [Indexed: 01/12/2023]
Abstract
The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.
Collapse
Affiliation(s)
- Ettore Bartalucci
- Max-Planck-Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany.,Institute of Technical and Macromolecular Chemistry, RWTH Aachen University, Worringerweg 2, 52074, Aachen, Germany
| | | | - Anne Mehnert
- Institute of Inorganic Chemistry, Leipzig University, Johannisallee 29, 04103, Leipzig, Germany
| | - Julius B Kleine Büning
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4, 53115, Bonn, Germany
| | - Lennart Günzel
- Institute of Inorganic Chemistry, Leipzig University, Johannisallee 29, 04103, Leipzig, Germany
| | - Maik Icker
- Institute of Organic Chemistry, Leipzig University Linnéstraße 3, 04103, Leipzig, Germany
| | - Martin Börner
- Institute of Inorganic Chemistry, Leipzig University, Johannisallee 29, 04103, Leipzig, Germany
| | - Christian Wiebeler
- Institute of Analytic Chemistry, Leipzig University, Linnéstraße 3, 04103, Leipzig, Germany.,Wilhelm-Ostwald-Institute for Physical and Theoretical Chemistry, Leipzig University, Linnéstraße 2, 04103, Leipzig, Germany
| | - Beat H Meier
- Physical Chemistry, ETH Zurich, 8093, Zurich, Switzerland
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4, 53115, Bonn, Germany
| | - Berthold Kersting
- Institute of Inorganic Chemistry, Leipzig University, Johannisallee 29, 04103, Leipzig, Germany
| | - Thomas Wiegand
- Max-Planck-Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany.,Institute of Technical and Macromolecular Chemistry, RWTH Aachen University, Worringerweg 2, 52074, Aachen, Germany.,previous address: Physical Chemistry, ETH Zurich, 8093, Zurich, Switzerland
| |
Collapse
|
4
|
Jensen F. Basis Set Extrapolation of Vibrational Frequencies. J Phys Chem A 2023; 127:2859-2863. [PMID: 36943253 DOI: 10.1021/acs.jpca.3c00064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2023]
Abstract
We investigate three different approaches for extrapolating harmonic vibrational frequencies to the complete basis set limit, by direct extrapolation of the frequencies, by calculation of frequencies based on extrapolated Hessians, or based on the Hessian from optimization of the extrapolated energy surface. For regular molecules, the three extrapolation procedures yield essentially identical results, but for loosely bound complexes, the frequencies derived from extrapolated Hessians yield unpredictable behavior. None of the basis set extrapolations, however, provide any significant improvement over the results upon which the extrapolation is based.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
5
|
Rzepiela K, Kaminský J, Buczek A, Broda MA, Kupka T. Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements? MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27238230. [PMID: 36500321 PMCID: PMC9737175 DOI: 10.3390/molecules27238230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 11/17/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022]
Abstract
The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H3PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.
Collapse
Affiliation(s)
- Kacper Rzepiela
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 166 10 Prague, Czech Republic
- Correspondence: (J.K.); (T.K.)
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Małgorzata A. Broda
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
- Correspondence: (J.K.); (T.K.)
| |
Collapse
|
6
|
Experimental and theoretical characterization of chelidonic acid structure. Struct Chem 2022. [DOI: 10.1007/s11224-022-02026-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Abstract
Abstract
Chelidonic acid (4-oxo-4H-pyran-2,6-dicarboxylic acid) is present in plants of Papaveraceae family, especially in Chelidonium majus. Due to its anticancer, antibacterial, hepatoprotective, and antioxidant properties, it has been used in medical treatments. In this work, the X-ray structure of methanol solvate of chelidonic acid was determined. Layers of chelidonic acid are held by hydrogen bonds via COOH and C = O fragments and additionally bridged by methanol. The formed H-bond network between two acid units is different from typical –COOH dimers observed, e.g., in crystals of isophtalic acid. The molecular structure of 2,6-dimethyl-γ-pyrone (2Me4PN) and chelidonic acid, a 2,6-dicarboxylic derivate of γ-pyrone (4PN), was verified in silico using density functional theory (DFT-B3LYP) combined with large correlation-consistent basis sets. The impact of –CH3 and –COOH substituents on 4PN ring structure, dipole moments, geometric/magnetic indexes of aromaticity, and NBO charges was assessed following unconstrained geometry optimization in the gas phase, chloroform, methanol, DMSO, and water with solvent effect introduced using the polarized continuous model (PCM). H-bond network formed in chelidonic acid–methanol complex was analyzed and their interaction energy estimated. Theoretical modeling enabled prediction of accurate structural parameters, dipole moments, and geometric/magnetic indexes of aromaticity of the studied 4PN, 2Me4PN, and chelidonic acid molecules.
Collapse
|
7
|
Krivdin LB. Computational 1 H and 13 C NMR in structural and stereochemical studies. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2022; 60:733-828. [PMID: 35182410 DOI: 10.1002/mrc.5260] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 02/14/2022] [Accepted: 02/16/2022] [Indexed: 06/14/2023]
Abstract
Present review outlines the advances and perspectives of computational 1 H and 13 C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational 1 H and 13 C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| |
Collapse
|
8
|
Makulski W, Aucar JJ, Aucar GA. Ammonia: the molecule for establishing 14N and 15N absolute shielding scales and a source of information on nuclear magnetic moments. J Chem Phys 2022; 157:084306. [DOI: 10.1063/5.0096523] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Multinuclear NMR studies of the gaseous mixtures 3He/14NH3 and 3He/15NH3 are reported. Precise analysis of the 3He, 14N, 15N and 1H resonance frequencies show linear dependence on the gas density. Extrapolation of these results to the zero-pressure limit gives ν0(1H), ν0(14N) and ν0(15N) resonance frequencies of the isolated ammonia molecule at 300K. The analogous value for 3He atoms in gaseous mixtures ν0(3He) was measured as well. The application of a new scheme to introduce the most important electronic effects on NMR shieldings, together with highly accurate quantum chemical calculations allow the 14/15N and 1H shielding of the isolated ammonia molecule to be obtained with the greatest accuracy and precision. For the first time, these studies were carried out on ammonia within the so-called four-component relativistic framework. The NMR frequency comparison method provides an approach for determining the 14N and 15N nuclear magnetic moments. The new shielding parameters in ammonia were used for reevaluation of the entire nitrogen absolute shielding scale. Additionally, the absolute shielding values of several gaseous compounds and secondary reference substances in liquids were presented. It was established that 14N and 15N absolute shielding constants in 14NH3 and 15NH3 are very similar, and only differ by less than 0.01 ppm, which is not usually measurable in NMR experiments. Precise calculations of 14N and 15N dipole moments were also made from these accurate shielding values.
Collapse
Affiliation(s)
| | - Juan J. Aucar
- Natural and Exact Science Faculty, Northeastern University of Argentina, Argentina
| | - Gustavo A. Aucar
- Natural and Exact Science Faculty, Northeastern University of Argentina, Argentina
| |
Collapse
|
9
|
Jaworski A, Hedin N. Electron correlation and vibrational effects in predictions of paramagnetic NMR shifts. Phys Chem Chem Phys 2022; 24:15230-15244. [PMID: 35703010 DOI: 10.1039/d2cp01206e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electronic structure calculations are fundamentally important for the interpretation of nuclear magnetic resonance (NMR) spectra from paramagnetic systems that include organometallic and inorganic compounds, catalysts, or metal-binding sites in proteins. Prediction of induced paramagnetic NMR shifts requires knowledge of electron paramagnetic resonance (EPR) parameters: the electronic g tensor, zero-field splitting D tensor, and hyperfine A tensor. The isotropic part of A, called the hyperfine coupling constant (HFCC), is one of the most troublesome properties for quantum chemistry calculations. Yet, even relatively small errors in calculations of HFCC tend to propagate into large errors in the predicted NMR shifts. The poor quality of A tensors that are currently calculated using density functional theory (DFT) constitutes a bottleneck in improving the reliability of interpretation of the NMR spectra from paramagnetic systems. In this work, electron correlation effects in calculations of HFCCs with a hierarchy of ab initio methods were assessed, and the applicability of different levels of DFT approximations and the coupled cluster singles and doubles (CCSD) method was tested. These assessments were performed for the set of selected test systems comprising an organic radical, and complexes with transition metal and rare-earth ions, for which experimental data are available. Severe deficiencies of DFT were revealed but the CCSD method was able to deliver good agreement with experimental data for all systems considered, however, at substantial computational costs. We proposed a more computationally tractable alternative, where the A was computed with the coupled cluster theory exploiting locality of electron correlation. This alternative is based on the domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) method. In this way the robustness and reliability of the coupled cluster theory were incorporated into the modern formalism for the prediction of induced paramagnetic NMR shifts, and became applicable to systems of chemical interest. This approach was verified for the bis(cyclopentadienyl)vanadium(II) complex (Cp2V; vanadocene), and the metal-binding site of the Zn2+ → Co2+ substituted superoxide dismutase (SOD) metalloprotein. Excellent agreement with experimental NMR shifts was achieved, which represented a substantial improvement over previous theoretical attempts. The effects of vibrational corrections to orbital shielding and hyperfine tensor were evaluated and discussed within the second-order vibrational perturbation theory (VPT2) framework.
Collapse
Affiliation(s)
- Aleksander Jaworski
- Department of Materials and Environmental Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden.
| | - Niklas Hedin
- Department of Materials and Environmental Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden.
| |
Collapse
|
10
|
Rusakov YY, Rusakova IL. New pecS- n ( n = 1, 2) basis sets for quantum chemical calculations of the NMR chemical shifts of H, C, N, and O nuclei. J Chem Phys 2022; 156:244112. [DOI: 10.1063/5.0096907] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper demonstrates the performance of our previously suggested property-energy consistent method on the example of the generation of effective basis sets, pecS-1 and pecS-2, suited for the calculation of hydrogen, carbon, nitrogen, and oxygen chemical shifts. The new basis sets were successfully approbated in the GIAO-DFT calculations of the chemical shifts of 35 molecules using six different functionals. The pecS-1 basis set demonstrated very good accuracy, which makes this small basis set an effective means for the large-scale computations. At the same time, the pecS-2 basis set also gave very accurate results, thus putting it on a par with the other commensurate basis sets suited for the chemical shifts calculations.
Collapse
Affiliation(s)
- Yuriy Yu. Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Irina L. Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
11
|
Palivec V, Pohl R, Kaminský J, Martinez-Seara H. Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment. J Chem Theory Comput 2022; 18:4373-4386. [PMID: 35687789 DOI: 10.1021/acs.jctc.2c00127] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Determining the structure of saccharides in their native environment is crucial to understanding their function and more accurately targeting their utilization. Nuclear magnetic resonance observables such as the nuclear Overhauser effect or spin-spin coupling constants are routinely utilized to study saccharides in their native water environment. However, while highly sensitive to the local environment, chemical shifts are mostly overlooked, despite being commonly measured for compounds identification. Although chemical shifts carry considerable structural information, their direct association with structure is notoriously difficult. This is mostly due to the similarity in the chemical nature of most saccharides causing similar physicochemical environments close to sugar C and H atoms, resulting in comparable chemical shifts. The rise of computational power allows one to compute reliable chemical shifts and use them to determine atomistic details of these sugars in solution. However, any prediction is severely limited by the computational protocol used and its accuracy. In this work, we studied a set of 31 saccharides on which we evaluated various computational protocols to calculate the total number of 375 1H and 327 13C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating 1H and 13C chemical shifts that we recommend for further use. These protocols can help with the interpretation of experimental spectra, but we also show that they are also capable of structure prediction independently. This is possible because of the low mean absolute deviations of calculated shifts from the experiment (0.06 ppm for 1H and 1.09 ppm for 13C). We explore different solvation methods, basis sets, and optimization schemes to reach such accuracy. A correct sampling of the conformation phase space of flexible sugar molecules is also key to obtaining accurately converged theoretical chemical shifts. The linear regression method was applied to convert the calculated isotropic nuclear magnetic shielding constants to simulated chemical shifts comparable with the experiment. The achieved level of accuracy can help in utilizing chemical shifts for elucidating the 3D atomistic structure of saccharides in aqueous solutions. All linear regression parameters obtained on our extensive set of sugars for all the tested protocols can be reutilized in future works.
Collapse
Affiliation(s)
- Vladimír Palivec
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Radek Pohl
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo námĕstí 542/2, Prague 6 CZ166 10, Czech Republic
| |
Collapse
|
12
|
Bharadwaj VS, Westawker LP, Crowley MF. Towards Elucidating Structure–Spectra Relationships in Rhamnogalacturonan II: Computational Protocols for Accurate 13C and 1H Shifts for Apiose and Its Borate Esters. Front Mol Biosci 2022; 8:756219. [PMID: 35141275 PMCID: PMC8820409 DOI: 10.3389/fmolb.2021.756219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Accepted: 11/24/2021] [Indexed: 11/13/2022] Open
Abstract
Apiose is a naturally occurring, uncommon branched-chain pentose found in plant cell walls as part of the complex polysaccharide Rhamnogalacturonan II (RG-II). The structural elucidation of the three-dimensional structure of RG-II by nuclear magnetic resonance (NMR) spectroscopy is significantly complicated by the ability of apiose to cross-link via borate ester linkages to form RG-II dimers. Here, we developed a computational approach to gain insight into the structure–spectra relationships of apio–borate complexes in an effort to complement experimental assignments of NMR signals in RG-II. Our protocol involved structure optimizations using density functional theory (DFT) followed by isotropic magnetic shielding constant calculations using the gauge-invariant atomic orbital (GIAO) approach to predict chemical shifts. We evaluated the accuracy of 23 different functional–basis set (FBS) combinations with and without implicit solvation for predicting the experimental 1H and 13C shifts of a methyl apioside and its three borate derivatives. The computed NMR predictions were evaluated on the basis of the overall shift accuracy, relative shift ordering, and the ability to distinguish between dimers and monomers. We demonstrate that the consideration of implicit solvation during geometry optimizations in addition to the magnetic shielding constant calculations greatly increases the accuracy of NMR chemical shift predictions and can correctly reproduce the ordering of the 13C shifts and yield predictions that are, on average, within 1.50 ppm for 13C and 0.12 ppm for 1H shifts for apio–borate compounds.
Collapse
|
13
|
Hägele G, Nazarski RB, Schmitz A, Xing S, Janiak C. 1H NMR spectra, structure, and conformational exchange of S- n-alkyl-tetrahydrothiophenium cations of some ionic liquids. PHOSPHORUS SULFUR 2022. [DOI: 10.1080/10426507.2021.2023149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Gerhard Hägele
- Institut für Anorganische und Strukturchemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Ryszard B. Nazarski
- Department of Physical Chemistry, Faculty of Chemistry, Theoretical and Structural Chemistry Group, University of Lodz, Łódź, Poland
| | - Alexa Schmitz
- Institut für Anorganische und Strukturchemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Shanghua Xing
- Institut für Anorganische und Strukturchemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Christoph Janiak
- Institut für Anorganische und Strukturchemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
14
|
Crittenden DL. A new double-reference correction scheme for accurate and efficient computation of NMR chemical shieldings. Phys Chem Chem Phys 2022; 24:27055-27063. [DOI: 10.1039/d2cp03992c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
Our novel correction procedure yields high-accuracy DFT predictions of absolute NMR shieldings and enables outliers due to relativistic effects or manifestly inadequate modelling of electron correlation to be easily and unambiguously identified.
Collapse
Affiliation(s)
- Deborah L. Crittenden
- School of Physical and Chemical Sciences, University of Canterbury, Christchurch, 8140, New Zealand
| |
Collapse
|
15
|
Schattenberg CJ, Kaupp M. Extended Benchmark Set of Main-Group Nuclear Shielding Constants and NMR Chemical Shifts and Its Use to Evaluate Modern DFT Methods. J Chem Theory Comput 2021; 17:7602-7621. [PMID: 34797677 DOI: 10.1021/acs.jctc.1c00919] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
An extended theoretical benchmark set, NS372, for light main-group nuclear shieldings and NMR shifts has been constructed based on high-level GIAO-CCSD(T)/pcSseg-3//CCSD(T)/cc-pVQZ reference data. After removal of the large static-correlation cases O3, F3-, and BH from the statistical evaluations for the 17O, 19F, and 11B subsets, the benchmark comprises overall 372 shielding values in 117 molecules with a wide range of electronic-structure situations, containing 124 1H, 14 11B, 93 13C, 43 15N, 31 17O, 47 19F, 14 31P, and 6 33S shielding constants. The CCSD(T)/pcSseg-3 data are shown to be close to the basis-set and method limit and thus provide an excellent benchmark to evaluate more approximate methods, such as density functional approaches. This dataset has been used to evaluate Hartree-Fock (HF) and MP2, and a wide range of exchange-correlation functionals from local density approximation (LDA) to generalized gradient approximations (GGAs) and meta-GGAs (focusing on their current-density functional implementations), as well as global hybrid, range-separated hybrid, local hybrid, and double-hybrid functionals. Starting with absolute shielding constants, the DSD-PBEP86 double hybrid is confirmed to provide the highest accuracy, with an aggregate relative mean absolute error (rel. MAE) of only 0.9%, followed by MP2 (1.1%). MP2 and double hybrids only show larger errors for a few systems with the largest static-correlation effects. The double-hybrid B2GP-PLYP, the two local hybrids cLH12ct-SsirPW92 and cLH12ct-SsifPW92, and the current-density functional meta-GGA cB97M-V follow closely behind (all 1.5%), as do some further functionals, cLH20t and cMN15-L (both 1.6%), as well as B2PLYP and KT3 (both 2.0%). Functionals on the lower rungs of the usual ladder offer the advantage of lower computational cost and access to larger molecules. Closer examination also reveals the best-performing methods for individual nuclei in the test set. Different ways of treating τ-dependent functionals are evaluated. When moving from absolute shielding constants to chemical shifts, some of the methods can benefit from systematic error compensation, and the overall error range somewhat narrows. Further methods now achieve the 2% threshold of relative MAEs, including functionals based on TPSS (TPSSh, cmPSTS).
Collapse
Affiliation(s)
- Caspar Jonas Schattenberg
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
16
|
Nguyen TT. 1H/ 13C chemical shift calculations for biaryls: DFT approaches to geometry optimization. ROYAL SOCIETY OPEN SCIENCE 2021; 8:210954. [PMID: 34631126 PMCID: PMC8479412 DOI: 10.1098/rsos.210954] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Accepted: 09/01/2021] [Indexed: 06/13/2023]
Abstract
Twelve common density functional methods and seven basis sets for geometry optimization were evaluated on the accuracy of 1H/13C NMR chemical shift calculations for biaryls. For these functionals, 1H shifts calculations for gas phase optimized geometries were significantly less accurate than those for in-solution optimized structures, while 13C results were not strongly influenced by geometry optimization methods and solvent effects. B3LYP, B3PW91, mPW1PW91 and ωB97XD were the best-performing functionals with lowest errors; among seven basis sets, DGDZVP2 and 6-31G(d,p) outperformed the others. The combination of these functionals and basis sets resulted in high accuracy with CMAEmin = 0.0327 ppm (0.76%) and 0.888 ppm (0.58%) for 1H and 13C, respectively. The selected functionals and basis set were validated when consistently producing optimized structures with high accuracy results for 1H and 13C chemical shift calculations of two other biaryls. This study highly recommends the IEFPCM/B3LYP, B3PW91, mPW1PW91 or ωB97XD/DGDZVP2 or 6-31G(d,p) level of theory for the geometry optimization step, especially the solvent incorporation, which would lead to high accuracy 1H/13C calculation. This work would assist in the fully structural assignments of biaryls and provide insights into in-solution biaryl conformations.
Collapse
Affiliation(s)
- Thien T Nguyen
- Institute of Research and Development, Duy Tan University, Da Nang 550000, Vietnam
- Faculty of Pharmacy, College of Medicine and Pharmacy, Duy Tan University, Da Nang 550000, Vietnam
| |
Collapse
|
17
|
DFT Calculations of 1H NMR Chemical Shifts of Geometric Isomers of Conjugated Linolenic Acids, Hexadecatrienyl Pheromones, and Model Triene-Containing Compounds: Structures in Solution and Revision of NMR Assignments. Molecules 2021; 26:molecules26113477. [PMID: 34200468 PMCID: PMC8201138 DOI: 10.3390/molecules26113477] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Revised: 06/02/2021] [Accepted: 06/04/2021] [Indexed: 11/29/2022] Open
Abstract
A DFT study of the 1H NMR chemical shifts, δ(1H), of geometric isomers of 18:3 conjugated linolenic acids (CLnAs), hexadecatrienyl pheromones, and model triene-containing compounds is presented, using standard functionals (B3LYP and PBE0) as well as corrections for dispersion interactions (B3LYP-D3, APFD, M06–2X and ωB97XD). The results are compared with literature experimental δ(1H) data in solution. The closely spaced “inside” olefinic protons are significantly more deshielded due to short-range through-space H…H steric interactions and appear close to or even beyond δ-values of aromatic systems. Several regularities of the computational δ(1H) of the olefinic protons of the conjugated double bonds are reproduced very accurately for the lowest-energy DFT-optimized single conformer for all functionals used and are in very good agreement with experimental δ(1H) in solution. Examples are provided of literature studies in which experimental resonance assignments deviate significantly from DFT predictions and, thus, should be revised. We conclude that DFT calculations of 1H chemical shifts of trienyl compounds are powerful tools (i) for the accurate prediction of δ(1H) even with less demanding functionals and basis sets; (ii) for the unequivocal identification of geometric isomerism of conjugated trienyl systems that occur in nature; (iii) for tackling complex problems of experimental resonance assignments due to extensive signal overlap; and (iv) for structure elucidation in solution.
Collapse
|
18
|
Colherinhas G. Updating atomic charge parameters of aliphatic amino acids: a quest to improve the performance of molecular modeling via sequential molecular dynamics and DFT-GIAO-NMR calculations. Phys Chem Chem Phys 2021; 23:8413-8425. [PMID: 33876005 DOI: 10.1039/d1cp00183c] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
In this work, we observe the behavior of the dipole moment, atomic charges, solute-solvent interactions and NMR spectroscopy of aliphatic amino acids in a water solution via the computational simulations of classical molecular dynamics and DFT quantum calculations. Our results indicate that the convergence of the atomic charge of the solute, from an iterative process, together with the dipole moment of the amino acid, alters the lifetime of hydrogen bonds present in the first solvation shell, resulting in the modification of its structure and dynamics. Using GIAO-DFT-NMR calculations, we assessed the impact of these structural solute-solvent modifications on the magnetic shielding constants of the solute carbon atoms. In this sense, we evaluate the importance of an update in parameters that describe atomic charges present in the CHARMM36 force field.
Collapse
Affiliation(s)
- Guilherme Colherinhas
- Departamento de Física, CEPAE, Universidade Federal de Goiás, 74690-900, Goiânia, GO, Brazil.
| |
Collapse
|
19
|
Lozynski M, Rusinska-Roszak D. Finding the direct energy-structure correlations in intramolecular aromaticity assisted hydrogen bonding (AAHB). J Mol Graph Model 2021; 105:107884. [PMID: 33725643 DOI: 10.1016/j.jmgm.2021.107884] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2020] [Revised: 02/25/2021] [Accepted: 02/25/2021] [Indexed: 11/25/2022]
Abstract
A predictive model for intramolecular hydrogen bond energy (EHB) calculation of polyaromatic ortho-hydroxyaldehydes based on a set of small, functionalized hydrocarbons is developed. The complete data set of 18 compounds was used for this study. The model is based on one of four optional categories of molecular descriptors: geometric, spectroscopic, bond order and topological indices. The model of Wiberg bond indices (WBIs) as descriptors of the CC involved bond based on stepwise regression has acceptable prediction abilities for 14 structures of ortho-hydroxyformylobenzo[a]pyrene derivatives already at the semiempirical level. The presented correlation enables a significantly more rapid and quantitative description of the hydrogen bonding strength than the much more time-consuming MTA method. Thus, WBIs are shown to provide a reliable means for fast prescreening of the energy of chelate hydrogen bonds potentially for any polyaromatic derivatives.
Collapse
Affiliation(s)
- Marek Lozynski
- Faculty of Chemical Technology, Poznan University of Technology, Berdychowo 4, 60-965, Poznan, Poland
| | - Danuta Rusinska-Roszak
- Faculty of Chemical Technology, Poznan University of Technology, Berdychowo 4, 60-965, Poznan, Poland.
| |
Collapse
|
20
|
Hansen PE, Saeed BA, Rutu RS, Kupka T. One-bond 1 J( 15 N,H) coupling constants at sp 2 -hybridized nitrogen of Schiff bases, enaminones and similar compounds: A theoretical study. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:750-762. [PMID: 32415997 DOI: 10.1002/mrc.5052] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2020] [Revised: 05/07/2020] [Accepted: 05/11/2020] [Indexed: 06/11/2023]
Abstract
1 J(15 N,H) coupling constants for enaminones and NH-forms of intramolecularly hydrogen-bonded Schiff bases as model compounds for sp2 -hybridized nitrogen atoms are evaluated using density functional theory (DFT) to find the optimal functionals and basis sets. Ammonia is used as a test molecule and its one-bond coupling constant is compared with experiment. A methylamine Schiff base of a truncated molecule of gossypol is used for checking the performance of selected B3LYP, O3LYP, PBE, BHandH, and APFD density functionals and standard, modified, and dedicated basis sets for coupling constants. Both in vacuum and in chloroform, modeled by the simple continuum model of solvent, the modified basis sets predict significantly better the 1 J(15 N,H) value in ammonia and in the methylamine Schiff base of a truncated molecule of gossypol than the standard basis sets. This procure is then used on a broad set of intramolecularly hydrogen-bonded molecules, and a good correlation between calculated and experimental one-bond NH coupling constants is obtained. The 1 J(15 N,H) couplings are slightly overestimated. The calculated data show for hydrogen-bonded NH interatomic distances that the calculated values depend on the NH bond lengths. The shorter the bond lengths, the larger the 1 J(15 N,H). A useful correlation between 1 J(15 N,H) and NH bond length is derived that enables realistic predictions of one-bond NH coupling constants. The calculations reproduce experimentally observed trends for the studied molecules.
Collapse
Affiliation(s)
- Poul Erik Hansen
- Department of Science and Environment, Roskilde University, Roskilde, Denmark
| | - Bahjat A Saeed
- Department of Chemistry, College of Education for Pure Sciences, University of Basrah, Iraq
| | - Rita S Rutu
- Department of Pharmaceutical Chemistry, College of Pharmacy, University of Basrah, Iraq
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| |
Collapse
|
21
|
Kupka T, Broda MA, Wieczorek PP. What is the form of muscimol from fly agaric mushroom (Amanita muscaria) in water? An insight from NMR experiment supported by molecular modeling. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:584-593. [PMID: 31912552 DOI: 10.1002/mrc.4990] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Revised: 12/30/2019] [Accepted: 01/02/2020] [Indexed: 06/10/2023]
Abstract
The biologically active alkaloid muscimol is present in fly agaric mushroom (Amanita muscaria), and its structure and action is related to human neurotransmitter γ -aminobutyric acid (GABA). The current study reports on determination of muscimol form present in water solution using multinuclear 1 H and 13 C nuclear magnetic resonance (NMR) experiments supported by density functional theory molecular modeling. The structures of three forms of free muscimol molecule both in the gas phase and in the presence of water solvent, modeled by polarized continuous model, and nuclear magnetic isotropic shieldings, the corresponding chemical shifts, and indirect spin-spin coupling constants were calculated. Several J-couplings observed in proton and carbon NMR spectra, not available before, are reported. The obtained experimental spectra, supported by theoretical calculations, favor the zwitterion form of muscimol in water. This structure differs from NH isomer, previously determined in dimethyl sulfoxide (DMSO) solution. In addition, positions of signals C3 and C5 are reversed in both solvents.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| | | | | |
Collapse
|
22
|
Schattenberg CJ, Reiter K, Weigend F, Kaupp M. An Efficient Coupled-Perturbed Kohn–Sham Implementation of NMR Chemical Shift Computations with Local Hybrid Functionals and Gauge-Including Atomic Orbitals. J Chem Theory Comput 2020; 16:931-943. [DOI: 10.1021/acs.jctc.9b00944] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Caspar Jonas Schattenberg
- Technische Universität Berlin, Institut für Chemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Kevin Reiter
- Karlsruher Institut für Technologie (KIT), Institut für Nanotechnologie, Postfach 3640, D-76021 Karlsruhe, Germany
| | - Florian Weigend
- Karlsruher Institut für Technologie (KIT), Institut für Nanotechnologie, Postfach 3640, D-76021 Karlsruhe, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
23
|
Krivdin LB. Computational 1 H NMR: Part 2. Chemical applications. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:5-14. [PMID: 31125992 DOI: 10.1002/mrc.4896] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/24/2019] [Revised: 05/12/2019] [Accepted: 05/18/2019] [Indexed: 06/09/2023]
Abstract
This is the second one of three closely interrelated reviews dealing with computation of 1 H NMR chemical shifts and 1 H-1 H spin-spin coupling constants prepared for Magnetic Resonance in Chemistry. Presented in this review are some basic notes and illustrative examples of how modern computational 1 H NMR could be used for structural elucidation and stereoelectronic studies of the medium-sized organic molecules involving saturated, unsaturated, aromatic, and heteroaromatic compounds together with their functional derivatives and coordination complexes to get deeper insight into their stereochemical structure and stereodynamic behavior.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
- Chair of Chemistry, Angarsk State Technical University, Tchaikovsky St. 60, 665835, Angarsk, Russia
| |
Collapse
|
24
|
Krivdin LB. Computational 1 H NMR: Part 1. Theoretical background. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2019; 57:897-914. [PMID: 30963636 DOI: 10.1002/mrc.4873] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2019] [Revised: 03/31/2019] [Accepted: 04/01/2019] [Indexed: 06/09/2023]
Abstract
This is the first one of the three closely interrelated reviews to be published in Magnetic Resonance in Chemistry dealing with accordingly theoretical background, chemical applications, and biochemical studies of and by means of computational 1 H NMR. Presented in the first part of the review is a general outline of the modern theoretical methods and accuracy factors of computational 1 H NMR involving locally dense basis set schemes, solvent effects, vibrational corrections, and relativistic effects performed at the density functional theory and/or nonempirical levels. This review is dedicated to Prof. Stephan Sauer in view of his invaluable contribution to the field of computational nuclear magnetic resonance.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
- Angarsk State Technical University, Angarsk, Russia
| |
Collapse
|
25
|
Srishailam K, Reddy BV, Rao GR. Investigation of torsional potentials, hindered rotation, molecular structure and vibrational properties of some biphenyl carboxaldehydes using spectroscopic techniques and density functional formalism. J Mol Struct 2019. [DOI: 10.1016/j.molstruc.2019.06.064] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
26
|
Performance of polarization-consistent vs. correlation-consistent basis sets for CCSD(T) prediction of water dimer interaction energy. J Mol Model 2019; 25:313. [DOI: 10.1007/s00894-019-4200-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2019] [Accepted: 08/30/2019] [Indexed: 10/26/2022]
Abstract
Abstract
Detailed study of Jensen’s polarization-consistent vs. Dunning’s correlation-consistent basis set families performance on the extrapolation of raw and counterpoise-corrected interaction energies of water dimer using coupled cluster with single, double, and perturbative correction for connected triple excitations (CCSD(T)) in the complete basis set (CBS) limit are reported. Both 3-parameter exponential and 2-parameter inverse-power fits vs. the cardinal number of basis set, as well as the number of basis functions were analyzed and compared with one of the most extensive CCSD(T) results reported recently. The obtained results for both Jensen- and Dunning-type basis sets underestimate raw interaction energy by less than 0.136 kcal/mol with respect to the reference value of − 4.98065 kcal/mol. The use of counterpoise correction further improves (closer to the reference value) interaction energy. Asymptotic convergence of 3-parameter fitted interaction energy with respect to both cardinal number of basis set and the number of basis functions are closer to the reference value at the CBS limit than other fitting approaches considered here. Separate fits of Hartree-Fock and correlation interaction energy with 3-parameter formula additionally improved the results, and the smallest CBS deviation from the reference value is about 0.001 kcal/mol (underestimated) for CCSD(T)/aug-cc-pVXZ calculations. However, Jensen’s basis set underestimates such value to 0.012 kcal/mol. No improvement was observed for using the number of basis functions instead of cardinal number for fitting.
Collapse
|
27
|
Kupka T, Mnich A, Broda MA. Performance of revised STO(1M)-3G basis set for prediction of 5-fluorocytosine chemical shifts. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2019; 57:489-498. [PMID: 31013546 DOI: 10.1002/mrc.4879] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2019] [Revised: 04/14/2019] [Accepted: 04/15/2019] [Indexed: 06/09/2023]
Abstract
Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied. The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted 13 C and 1 H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for 13 C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Adrianna Mnich
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| |
Collapse
|
28
|
Impact of the ΔPhe configuration on the Boc-Gly-ΔPhe-NHMe conformation: experiment and theory. Struct Chem 2019. [DOI: 10.1007/s11224-019-01387-w] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
29
|
Kupka T, Gajda Ł, Stobiński L, Kołodziej Ł, Mnich A, Buczek A, Broda MA. Local aromaticity mapping in the vicinity of planar and nonplanar molecules. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2019; 57:359-372. [PMID: 31034627 DOI: 10.1002/mrc.4880] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Accepted: 04/23/2019] [Indexed: 06/09/2023]
Abstract
We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H2 Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, 3 He, and 7 Li+ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and 3 He chemical shifts in the middle of neutral C60 and C606- anion agreed well with the reported experimental NMR values for He@C60 and He@C606- .
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Łukasz Gajda
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Leszek Stobiński
- Faculty of Process and Chemical Engineering, Warsaw University of Technology, 1, Waryńskiego Street, 00-645, Warsaw, Poland
- nanoMAT Ltd., Waryńskiego 1 (FL. 4, RM. 401), 00-645, Warsaw, Poland
| | | | - Adrianna Mnich
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, Opole, Poland
| | | |
Collapse
|
30
|
Krivdin LB. Computational protocols for calculating 13C NMR chemical shifts. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2019; 112-113:103-156. [PMID: 31481156 DOI: 10.1016/j.pnmrs.2019.05.004] [Citation(s) in RCA: 58] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2019] [Revised: 05/19/2019] [Accepted: 05/20/2019] [Indexed: 06/10/2023]
Abstract
The most recent results dealing with the computation of 13C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
31
|
Lacerda EG, Kamounah FS, Coutinho K, Sauer SPA, Hansen PE, Hammerich O. Computational Prediction of 1 H and 13 C NMR Chemical Shifts for Protonated Alkylpyrroles: Electron Correlation and Not Solvation is the Salvation. Chemphyschem 2018; 20:78-91. [PMID: 30452112 DOI: 10.1002/cphc.201801066] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2018] [Indexed: 12/18/2022]
Abstract
Prediction of chemical shifts in organic cations is known to be a challenge. In this article we meet this challenge for α-protonated alkylpyrroles, a class of compounds not yet studied in this context, and present a combined experimental and theoretical study of the 13 C and 1 H chemical shifts in three selected pyrroles. We have investigated the importance of the solvation model, basis set, and quantum chemical method with the goal of developing a simple computational protocol, which allows prediction of 13 C and 1 H chemical shifts with sufficient accuracy for identifying such compounds in mixtures. We find that density functional theory with the B3LYP functional is not sufficient for reproducing all 13 C chemical shifts, whereas already the simplest correlated wave function model, Møller-Plesset perturbation theory (MP2), leads to almost perfect agreement with the experimental data. Treatment of solvent effects generally improves the agreement with experiment to some extent and can in most cases be accomplished by a simple polarizable continuum model. The only exception is the NH proton, which requires inclusion of explicit solvent molecules in the calculation.
Collapse
Affiliation(s)
- Evanildo G Lacerda
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen, Denmark.,Instituto de Física da Universidade de São Paulo, CP 66318, 05314-970, São Paulo, SP, Brazil
| | - Fadhil S Kamounah
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen, Denmark.,Department of Science and Environment, Roskilde University, Universitetsvej 1, DK-4000, Roskilde, Denmark
| | - Kaline Coutinho
- Instituto de Física da Universidade de São Paulo, CP 66318, 05314-970, São Paulo, SP, Brazil
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen, Denmark
| | - Poul Erik Hansen
- Department of Science and Environment, Roskilde University, Universitetsvej 1, DK-4000, Roskilde, Denmark
| | - Ole Hammerich
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen, Denmark
| |
Collapse
|
32
|
Abstract
We investigate by explicit parameter optimization to what extent basis sets of polarized double-ζ quality can introduce compensating errors in five different density functional methods. It is shown that minor changes in the contraction coefficients of the valence functions in the basis sets can have a significant impact and allow different density functional methods to achieve very similar performances. This holds for nuclear magnetic shielding constants and for isomerization energies, barrier heights, and noncovalent interactions. It is furthermore shown that errors due to neglect of vibrational and solvent effects can be absorbed in the combined method and basis set errors. These findings hold for data sets consisting of 50-150 data points. This raises the question of whether the common practice of identifying combinations of density functional methods and basis sets that have a good performance against a selected set of reference data should be considered as data fitting in the combined parameter space spanned by the method and basis set.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry , Aarhus University , DK-8000 Aarhus , Denmark
| |
Collapse
|
33
|
Cuny J, Jolibois F, Gerber IC. Evaluation of Gas-to-Liquid 17O Chemical Shift of Water: A Test Case for Molecular and Periodic Approaches. J Chem Theory Comput 2018; 14:4041-4051. [DOI: 10.1021/acs.jctc.8b00243] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques (LCPQ/IRSAMC), Université de Toulouse and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Franck Jolibois
- LPCNO, Université Fédérale de Toulouse Midi-Pyrénées, INSA-CNRS-UPS, 135 avenue de Rangueil, 31077 Cedex 4 Toulouse, France
| | - Iann C. Gerber
- LPCNO, Université Fédérale de Toulouse Midi-Pyrénées, INSA-CNRS-UPS, 135 avenue de Rangueil, 31077 Cedex 4 Toulouse, France
| |
Collapse
|
34
|
Nishitani S, Maekawa Y, Sakata T. Understanding the Molecular Structure of the Sialic Acid-Phenylboronic Acid Complex by using a Combined NMR Spectroscopy and DFT Study: Toward Sialic Acid Detection at Cell Membranes. ChemistryOpen 2018; 7:513-519. [PMID: 30003005 PMCID: PMC6031860 DOI: 10.1002/open.201800071] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2018] [Indexed: 12/03/2022] Open
Abstract
The origin of the unusually high stability of the sialic acid (SA) and phenylboronic acid (PBA) complex was investigated by a combined nuclear magnetic resonance (NMR) spectroscopy and density functional theory (DFT) study. SA is a glycan-terminating monosaccharide, and its importance as a clinical target has long been recognized. Inspired by the fact that the binding properties of SA-PBA complexation are anomalously high relative to those of typical monosaccharides, great effort has been made to build a clinical platform with the use of PBA as a SA-selective receptor. Although a number of applications have been reported in recent years, the ability of PBA to recognize SA-terminating surface glycans selectively is still unclear, because high-affinity SA-PBA complexation might not occur in a physiological environment. In particular, different forms of SA (α- and β-pyranose) were not considered in detail. To answer this question, the combined NMR spectroscopy/DFT study revealed that the advantageous binding properties of the SA-PBA complex arise from ester bonding involving the α-carboxylate moieties (C1 and C2) of β-SA but not α-SA. Moreover, the facts that the C2 atom is blocked by a glycoside bond in a physiological environment and that α-SA basically exists on membrane-bound glycans in a physiological environment lead to the conclusion that PBA cannot selectively recognize the SA unit to discriminate specific types of cells. Our results have a significant impact on the field of SA-based cell recognition.
Collapse
Affiliation(s)
- Shoichi Nishitani
- Department of Materials EngineeringGraduate School of EngineeringThe University of Tokyo7-3-1 Hongo, Bunkyo-kuTokyo113–8656Japan
| | - Yuki Maekawa
- Department of Materials EngineeringGraduate School of EngineeringThe University of Tokyo7-3-1 Hongo, Bunkyo-kuTokyo113–8656Japan
| | - Toshiya Sakata
- Department of Materials EngineeringGraduate School of EngineeringThe University of Tokyo7-3-1 Hongo, Bunkyo-kuTokyo113–8656Japan
| |
Collapse
|
35
|
Kupka T, Buczek A, Broda MA, Gajda Ł, Ignasiak M. Convergence of nuclear magnetic shieldings and one-bond 1 J( 11 B 1 H) indirect spin-spin coupling constants in small boron molecules. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2018; 56:338-351. [PMID: 29361201 DOI: 10.1002/mrc.4714] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2017] [Revised: 01/10/2018] [Accepted: 01/13/2018] [Indexed: 06/07/2023]
Abstract
Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH3 molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual 11 B and 1 H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced 11 B shielding in BH molecule.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Łukasz Gajda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Monika Ignasiak
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| |
Collapse
|
36
|
Xin D, Jones PJ, Gonnella NC. DiCE: Diastereomeric in Silico Chiral Elucidation, Expanded DP4 Probability Theory Method for Diastereomer and Structural Assignment. J Org Chem 2018; 83:5035-5043. [DOI: 10.1021/acs.joc.8b00338] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dongyue Xin
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut 06877, United States
| | - Paul-James Jones
- Information Technology, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut 06877, United States
| | - Nina C. Gonnella
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut 06877, United States
| |
Collapse
|
37
|
Roggatz CC, Lorch M, Benoit DM. Influence of Solvent Representation on Nuclear Shielding Calculations of Protonation States of Small Biological Molecules. J Chem Theory Comput 2018; 14:2684-2695. [PMID: 29566332 DOI: 10.1021/acs.jctc.7b01020] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In this study, we assess the influence of solvation on the accuracy and reliability of isotropic nuclear magnetic shielding calculations for amino acids in comparison to experimental data. We focus particularly on the performance of solvation methods for different protonation states, as biological molecules occur almost exclusively in aqueous solution and are subject to protonation with pH. We identify significant shortcomings of current implicit solvent models and present a hybrid solvation approach that improves agreement with experimental data by taking into account the presence of direct interactions between amino acid protonation state and water molecules.
Collapse
|
38
|
Gajda Ł, Kupka T, Broda MA, Leszczyńska M, Ejsmont K. Method and basis set dependence of the NICS indexes of aromaticity for benzene. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2018; 56:265-275. [PMID: 29211311 DOI: 10.1002/mrc.4690] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Revised: 11/21/2017] [Accepted: 11/24/2017] [Indexed: 05/24/2023]
Abstract
The role of theory level in prediction of benzene magnetic indexes of aromaticity is analysed and compared with calculated nuclear magnetic shieldings of 3 He used as NMR probe. Three closely related nucleus-independent chemical shift (NICS) based indexes were calculated for benzene at SCF-HF, MP2, and DFT levels of theory and the impact of basis set on these quantities was studied. The changes of benzene NICS(0), NICS(1), and NICS(1)zz parameters calculated using SCF-HF, MP2 and several density functionals were within 1 to 3 ppm. Similar deviations between magnetic indexes of aromaticity were observed for values calculated with selected basis sets. Only very small effect of polar solvent on benzene aromaticity was predicted. The 3 He nuclear magnetic isotropic shielding (σ) and its zz-components (σzz ) of helium atom approaching the centre of benzene ring from above produced similar curves versus benzene-He distance to NICS parameters calculated for similarly moving Bq ghost atom. We also propose an experimental verification of NICS calculations by designing the 3 He NMR measurement for benzene saturated with helium gas or in low temperature matrices.
Collapse
Affiliation(s)
- Łukasz Gajda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | | | - Krzysztof Ejsmont
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| |
Collapse
|
39
|
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
|
40
|
Marć M, Kupka T, Wieczorek PP, Namieśnik J. Computational modeling of molecularly imprinted polymers as a green approach to the development of novel analytical sorbents. Trends Analyt Chem 2018. [DOI: 10.1016/j.trac.2017.10.020] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
41
|
How reliable are Minnesota density functionals for modeling phosphorus–hydrogen NMR spin–spin coupling constants? Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2182-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
42
|
Krivdin LB. Calculation of 15N NMR chemical shifts: Recent advances and perspectives. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2017; 102-103:98-119. [PMID: 29157495 DOI: 10.1016/j.pnmrs.2017.08.001] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Revised: 08/21/2017] [Accepted: 08/21/2017] [Indexed: 06/07/2023]
Abstract
Recent advances in computation of 15N NMR chemical shifts are reviewed, concentrating mainly on practical aspects of computational protocols and accuracy factors. The review includes the discussion of the level of theory, the choice of density functionals and basis sets together with taking into account solvent effects, rovibrational corrections and relativistic effects. Computational aspects of 15N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
43
|
Xin D, Sader CA, Fischer U, Wagner K, Jones PJ, Xing M, Fandrick KR, Gonnella NC. Systematic investigation of DFT-GIAO 15N NMR chemical shift prediction using B3LYP/cc-pVDZ: application to studies of regioisomers, tautomers, protonation states and N-oxides. Org Biomol Chem 2017; 15:928-936. [PMID: 28050610 DOI: 10.1039/c6ob02450e] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The calculation of 15N NMR chemical shifts has been systematically investigated using density functional theory-gauge including/invariant atomic orbitals (DFT-GIAO) approximation at the B3LYP/cc-pVDZ level of theory. General linear regression terms for 15N chemical shift predictions were calculated for nitromethane and liquid ammonia references in DMSO. Both aliphatic and aromatic nitrogens were studied using a diverse set of molecular scaffolds. Statistical error analysis between experiment and prediction revealed that, with the exception of primary amines, 95% of linear scaled N-15 chemical shifts are within a ±9.56 ppm range. Comparison of the 15N calculated isotropic chemical shifts with the experimentally determined chemical shifts provided accurate assignment of the correct structure in cases where experimental data was ambiguous or inconclusive. Application of 15N prediction proved to be highly effective in identifying the correct regio-isomer, oxidation state, protonation state and preferred tautomer in solution.
Collapse
Affiliation(s)
- Dongyue Xin
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA.
| | - Charles Avery Sader
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA.
| | - Udo Fischer
- Boehringer Ingelheim Pharma GmbH & Co.KG, Germany
| | - Klaus Wagner
- Boehringer Ingelheim Pharma GmbH & Co.KG, Germany
| | - Paul-James Jones
- Information Technology, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA
| | - Minli Xing
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA.
| | - Keith R Fandrick
- Chemical Development, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA
| | - Nina C Gonnella
- Material and Analytical Sciences, Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, Connecticut, USA.
| |
Collapse
|
44
|
Faber R, Buczek A, Kupka T, Sauer SPA. On the convergence of zero-point vibrational corrections to nuclear shieldings and shielding anisotropies towards the complete basis set limit in water. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1210831] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Rasmus Faber
- Department of Chemistry, University of Copenhagen, Copenhagen Ø, Denmark
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| | | |
Collapse
|
45
|
Kupka T, Stachów M, Stobiński L, Kaminský J. DFT study of zigzag (n, 0) single-walled carbon nanotubes: 13C NMR chemical shifts. J Mol Graph Model 2016; 67:14-9. [DOI: 10.1016/j.jmgm.2016.04.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Revised: 04/14/2016] [Accepted: 04/19/2016] [Indexed: 10/21/2022]
|
46
|
Buczek A, Makowski M, Jewgiński M, Latajka R, Kupka T, Broda MA. Toward engineering efficient peptidomimetics. Screening conformational landscape of two modified dehydroaminoacids. Biopolymers 2016; 101:28-40. [PMID: 23606332 DOI: 10.1002/bip.22264] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2013] [Accepted: 04/12/2013] [Indexed: 11/11/2022]
Abstract
Effective peptidomimetics should posses structural rigidity and appropriate interaction pattern leading to potential spatial and electronic matching to the target receptor site. Rational design of such small bioactive molecules could push chemical synthesis and molecular modeling toward faster progress in medicinal chemistry. Conformational properties of N-t-butoxycarbonyl-glycine-(E/Z)-dehydrophenylalanine N',N'-dimethylamides (Boc-Gly-(E/Z)-ΔPhe-NMe2 ) in chloroform were studied by NMR and IR spectroscopy. The experimental findings were supported by extensive calculations at DFT(B3LYP, M06-2X) and MP2 levels of theory and the β-turn tendency for both isomers of the studied dipeptide were determined in vacuum and in solution. The theoretical data and experimental IR results were used as an additional information for the NMR-based determination of the detailed solution conformations of the peptides. The obtained results reveal that N-methylation of C-terminal amide group changes dramatically the conformational properties of studied dehydropeptides. Theoretical conformational analysis reveals that the tendency to adopt β-turn conformations is much weaker for the N-methylated Z isomer (Boc-Gly-(Z)-ΔPhe-NMe2 ), both in vacuum and in polar environment. On the contrary, N-methylated E isomer (Boc-Gly-(E)-ΔPhe-NMe2 ) can easily adopt β-turn conformation, but the backbone torsion angles (φ1, ψ1, φ2, ψ2) are off the limits for common β-turn types.
Collapse
Affiliation(s)
- Aneta Buczek
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| | | | | | | | | | | |
Collapse
|
47
|
Lummiss JAM, Perras FA, McDonald R, Bryce DL, Fogg DE. Sterically Driven Olefin Metathesis: The Impact of Alkylidene Substitution on Catalyst Activity. Organometallics 2016. [DOI: 10.1021/acs.organomet.5b00984] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Justin A. M. Lummiss
- Centre for Catalysis Research & Innovation and Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, ON, Canada, K1N 6N5
| | - Frédéric A. Perras
- Centre for Catalysis Research & Innovation and Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, ON, Canada, K1N 6N5
| | - Robert McDonald
- X-ray
Crystallographic Laboratory, Department of Chemistry, University of Alberta, Edmonton, AB, Canada, T6G 2G2
| | - David L. Bryce
- Centre for Catalysis Research & Innovation and Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, ON, Canada, K1N 6N5
| | - Deryn E. Fogg
- Centre for Catalysis Research & Innovation and Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, ON, Canada, K1N 6N5
| |
Collapse
|
48
|
Faber R, Kaminsky J, Sauer SPA. Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters. GAS PHASE NMR 2016. [DOI: 10.1039/9781782623816-00218] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.
Collapse
Affiliation(s)
- Rasmus Faber
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| | - Jakub Kaminsky
- Department of Molecular Spectroscopy, Institute of Organic Chemistry and Biochemistry 166 10 Prague Czech Republic
| | - Stephan P. A. Sauer
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| |
Collapse
|
49
|
Buczek A, Kupka T, Broda MA, Żyła A. Predicting the structure and vibrational frequencies of ethylene using harmonic and anharmonic approaches at the Kohn-Sham complete basis set limit. J Mol Model 2016; 22:42. [PMID: 26800989 PMCID: PMC4723629 DOI: 10.1007/s00894-015-2902-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Accepted: 12/22/2015] [Indexed: 11/25/2022]
Abstract
In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested. The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies. In this study, the performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, and 6) was tested. For the first time, we demonstrated regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit ![]()
Collapse
Affiliation(s)
- Aneta Buczek
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Adriana Żyła
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.,Department of Biosystematics, University of Opole, Oleska 22, 45-052, Opole, Poland.,Institute of Physics, Adam Mickiewicz University in Poznan, 65, Umultowska Street, 61-614, Poznań, Poland
| |
Collapse
|
50
|
Kupka T, Wieczorek PP. Theoretical and experimental NMR studies on muscimol from fly agaric mushroom (Amanita muscaria). SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2016; 153:216-225. [PMID: 26312739 DOI: 10.1016/j.saa.2015.08.026] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/22/2015] [Revised: 07/13/2015] [Accepted: 08/12/2015] [Indexed: 06/04/2023]
Abstract
In this article we report results of combined theoretical and experimental NMR studies on muscimol, the bioactive alkaloid from fly agaric mushroom (Amanita muscaria). The assignment of (1)H and (13)C NMR spectra of muscimol in DMSO-d6 was supported by additional two-dimensional heteronuclear correlated spectra (2D NMR) and gauge independent atomic orbital (GIAO) NMR calculations using density functional theory (DFT). The effect of solvent in theoretical calculations was included via polarized continuum model (PCM) and the hybrid three-parameter B3LYP density functional in combination with 6-311++G(3df,2pd) basis set enabled calculation of reliable structures of non-ionized (neutral) molecule and its NH and zwitterionic forms in the gas phase, chloroform, DMSO and water. GIAO NMR calculations, using equilibrium and rovibrationally averaged geometry, at B3LYP/6-31G* and B3LYP/aug-cc-pVTZ-J levels of theory provided muscimol nuclear magnetic shieldings. The theoretical proton and carbon chemical shifts were critically compared with experimental NMR spectra measured in DMSO. Our results provide useful information on its structure in solution. We believe that such data could improve the understanding of basic features of muscimol at atomistic level and provide another tool in studies related to GABA analogs.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, Opole University, Oleska 48, 45-052 Opole, Poland.
| | - Piotr P Wieczorek
- Faculty of Chemistry, Opole University, Oleska 48, 45-052 Opole, Poland.
| |
Collapse
|