Predicting the NMR spectra of nucleotides by DFT calculations: cyclic uridine monophosphate

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ 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 ¹ H and ¹³ 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.

Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment

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 ¹H and 327 ¹³C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating ¹H and ¹³C 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 ¹H and 1.09 ppm for ¹³C). 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.

Computational 1 H NMR: Part 3. Biochemical studies

This is the third and the last part of three closely interrelated reviews dealing with computation of ¹ H nuclear magnetic resonance chemical shifts and ¹ H-¹ H spin-spin coupling constants. Present review deals with the computation of these parameters in biologically active natural products, carbohydrates, and other molecules of biological origin focusing on stereochemical applications of computational ¹ H nuclear magnetic resonance to these objects.

An automated framework for NMR chemical shift calculations of small organic molecules

When using nuclear magnetic resonance (NMR) to assist in chemical identification in complex samples, researchers commonly rely on databases for chemical shift spectra. However, authentic standards are typically depended upon to build libraries experimentally. Considering complex biological samples, such as blood and soil, the entirety of NMR spectra required for all possible compounds would be infeasible to ascertain due to limitations of available standards and experimental processing time. As an alternative, we introduce the in silico Chemical Library Engine (ISiCLE) NMR chemical shift module to accurately and automatically calculate NMR chemical shifts of small organic molecules through use of quantum chemical calculations. ISiCLE performs density functional theory (DFT)-based calculations for predicting chemical properties—specifically NMR chemical shifts in this manuscript—via the open source, high-performance computational chemistry software, NWChem. ISiCLE calculates the NMR chemical shifts of sets of molecules using any available combination of DFT method, solvent, and NMR-active nuclei, using both user-selected reference compounds and/or linear regression methods. Calculated NMR chemical shifts are provided to the user for each molecule, along with comparisons with respect to a number of metrics commonly used in the literature. Here, we demonstrate ISiCLE using a set of 312 molecules, ranging in size up to 90 carbon atoms. For each, calculation of NMR chemical shifts have been performed with 8 different levels of DFT theory, and with solvation effects using the implicit solvent Conductor-like Screening Model. The DFT method dependence of the calculated chemical shifts have been systematically investigated through benchmarking and subsequently compared to experimental data available in the literature. Furthermore, ISiCLE has been applied to a set of 80 methylcyclohexane conformers, combined via Boltzmann weighting and compared to experimental values. We demonstrate that our protocol shows promise in the automation of chemical shift calculations and, ultimately, the expansion of chemical shift libraries.

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.

Chemical Constituents Identified from Fruit Body of Cordyceps bassiana and Their Anti-Inflammatory Activity

Cordyceps bassiana is one of Cordyceps species with anti-oxidative, anti-cancer, anti-inflammatory, anti-diabetic, anti-obesity, anti-angiogenic, and anti-nociceptive activities. This mushroom has recently demonstrated to have an ability to reduce 2,4-dinitrofluorobenzene-induced atopic dermatitis symptoms in NC/Nga mice. In this study, we further examined phytochemical properties of this mushroom by column chromatography and HPLC analysis. By chromatographic separation and spectroscopic analysis, 8 compounds, such as 1,9-dimethylguanine (1), adenosine (2), uridine (3), nicotinamide (4), 3-methyluracil (5), 1,7-dimethylxanthine (6), nudifloric acid (7), and mannitol (8) were identified from 6 different fractions and 4 more subfractions. Through evaluation of their anti-inflammatory activities using reporter gene assay and mRNA analysis, compound 1 was found to block luciferase activity induced by NF-κB and AP-1, suppress the mRNA levels of cyclooxygenase (COX)-2 and tumor necrosis factor (TNF)-α. Therefore, our data strongly suggests that compound 1 acts as one of major principles in Cordyceps bassiana with anti-inflammatory and anti-atopic dermatitis activities.

1H NMR spectra of alcohols and diols in chloroform: DFT/GIAO calculation of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in chloroform have been computed on the basis of density functional theory, the solvent being included by the integral-equation-formalism polarisable continuum model of Gaussian 09. Relative energies of all conformers are calculated at the Perdew, Burke and Ernzerhof (PBE)0/6-311+G(d,p) level, and NMR shifts by the gauge-including atomic orbital method with the PBE0/6-311+G(d,p) geometry and the cc-pVTZ basis set. The 208 computed CH proton NMR shifts for 34 alcohols correlate very well with the experimental values, with a gradient of 1.00 ± 0.01 and intercept close to zero; the overall root mean square difference (RMSD) is 0.08 ppm. Shifts for CH protons of diols in chloroform are well correlated with the theoretical values for (isotropic) benzene, with similar gradient and intercept (1.02 ± 0.01, -0.13 ppm), but the overall RMSD is slightly higher, 0.12 ppm. This approach generally gives slightly better results than the CHARGE model of Abraham et al. The shifts of unsaturated alcohols in benzene have been re-examined with Gaussian 09, but the overall fit for CH protons is not improved, and OH proton shifts are worse. Shifts of vinyl protons in alkenols are systematically overestimated, and the correlation of computed shifts against the experimental data for unsaturated alcohols follows a quadratic equation. Splitting the 20 compounds studied into two sets, and applying empirical scaling based on the quadratic for the first set to the second set, gives an RMSD of 0.10 ppm. A multi-standard approach gives a similar result.

1H NMR spectra of butane-1,4-diol and other 1,4-diols: DFT calculation of shifts and coupling constants

The proton nuclear magnetic resonance (NMR) spectra of butane-1,4-diol, pentane-1,4-diol, (S,S)-hexane-2,5-diol, 2,5-dimethylhexane-2,5-diol and cyclohexane-1,4-diols (cis and trans) in benzene and some other solvents have been analysed. The conformer distribution and the NMR shifts of these diols in benzene have been computed on the basis of the density functional theory, the solvent being included by means of the integral-equation-formalism polarizable continuum model implemented in Gaussian 09. Relative Gibbs energies of all conformers are calculated at the Perdew, Burke and Ernzerhof (PBE)0/6-311+G(d,p) level and NMR shifts by the gauge-including atomic orbital method with the PBE0/6-311+G(d,p) geometry and the cc-pVTZ basis set. Vicinal three-bond coupling constants for the acyclic diols are calculated from the relative conformer populations, the geometries and generalized Karplus equations developed by Altona's group; these correlate well with the experimental values. The solvent dependence of coupling constants for butane-1,4-diol is attributed to conformational change. Coupling constants for the rigid cyclohexane-1,4-diols do not change with solvent and are readily explained in terms of their geometries. The NMR shifts of hydrogen-bonded protons in individual conformers of alkane-1,n-diols show a very rough correlation with the OH · · · OH distances. The computed overall NMR shifts for CH protons in 1,2-diols, 1,3-diols and 1,4-diols are systematically high but correlate very well with the experimental values, with a gradient of 1.07 ± 0.01; those for OH protons correlate less well.

1H NMR spectra of alkane-1,3-diols in benzene: GIAO/DFT shift calculations

The proton nuclear magnetic resonance (NMR) spectra of propane-1,3-diol, 2-methylpropane-1,3-diol, 2,2-dimethylpropane-1,3-diol, butane-1,3-diol, 3-methylbutane-1,3-diol, pentane-2,4-diols (dl and meso), 2-methylpentane-2,4-diol and cyclohexane-1,3-diols (cis and trans) in benzene have been analysed. The conformer distribution and the NMR shifts of these diols have been computed on the basis of density functional theory, the solvent being included by means of the integral equation formalism phase continuum model (IEFPCM) implemented in Gaussian 09. Relative Gibbs energies of all conformers are calculated at the Perdew, Burke and Ernzerhof (PBE)0/6-311 + G(d,p) level, and NMR shifts by the gauge-including atomic orbital method with the PBE0/6-311 + G(d,p) geometry and the cc-pVTZ basis set. Vicinal coupling constants for 1,2- and 1,3-diols are rationalised in terms of relative conformer populations and geometries. The NMR shifts of hydrogen-bonded protons in individual conformers of alkane-1,n-diols show a very rough correlation with the OH⋯OH distances. The computed overall NMR shifts for CH protons in 1,2- and 1,3-diols are systematically high but correlate very well with the experimental values, with a gradient of 1.07 ± 0.01. Some values for nonequivalent methylene protons in 1,3-diols are reversed, calculation giving enhanced values for the proton anti to the COH bonds. Errors in the NMR shifts computed for the OH protons of nonsymmetrical diols appear to be related to relative populations of conformers where one or other of the OH groups is the donor. Some results based on the second-order Møller-Plesset approach, the Becke three-parameter Lee-Yang-Parr method and on the IEFPCM solvation model implemented in Gaussian 03 are included.

Empirical and DFT GIAO quantum-mechanical methods of (13)C chemical shifts prediction: competitors or collaborators?

The accuracy of (13)C chemical shift prediction by both DFT GIAO quantum-mechanical (QM) and empirical methods was compared using 205 structures for which experimental and QM-calculated chemical shifts were published in the literature. For these structures, (13)C chemical shifts were calculated using HOSE code and neural network (NN) algorithms developed within our laboratory. In total, 2531 chemical shifts were analyzed and statistically processed. It has been shown that, in general, QM methods are capable of providing similar but inferior accuracy to the empirical approaches, but quite frequently they give larger mean average error values. For the structural set examined in this work, the following mean absolute errors (MAEs) were found: MAE(HOSE) = 1.58 ppm, MAE(NN) = 1.91 ppm and MAE(QM) = 3.29 ppm. A strategy of combined application of both the empirical and DFT GIAO approaches is suggested. The strategy could provide a synergistic effect if the advantages intrinsic to each method are exploited.

Theoretical investigation on multinuclear NMR chemical shifts of some tris(trifluoromethyl)boron complexes

Tris(trifluoromethyl)boron complexes have unusual properties and may find applications in many fields of chemistry, biology, and physics. To gain insight into their NMR properties, the isotropic 11B, 13C, and 19F NMR chemical shifts of a series of tris(trifluoromethyl)boron complexes were systematically studied using the gauge-included atomic orbitals (GIAO) method at the levels of B3LYP/6-31 + G(d,p)//B3LYP/6-31G* and B3LYP/6-311 + G(d,p)//B3LYP/6-311 + G(d,p). Solvent effects were taken into account by polarizable continuum models (PCM). The calculated results were compared with the experimental values. The reason that the structurally inequivalent fluorine atoms in a specific species give a same chemical shift in experimental measurements is attributed to the fast rotation of CF3 group around the B-C(F3) bond because of the low energy barrier. The calculated 11B, 13C(F3), and 19F chemical shifts are in good agreement with the experimental measurements, while the deviations of calculated 13C(X, X = O, N) chemical shifts are slightly large. For the latter, the average absolute deviations of the results from B3LYP/6-311 + G(d,p)//B3LYP/6-311 + G(d,p) are smaller than those from B3LYP/6-31 + G(d,p)//B3LYP/6-31G*, and the inclusion of PCM reduces the deviation values. The calculated 19F and 11B chemical shieldings of (CF3)3BCO are greatly dependent on the optimized structures, while the influence of structural parameters on the calculated 13C chemical shieldings is minor.