Computational studies of 13C NMR chemical shifts of saccharides

The (13)C NMR chemical shifts for alpha-D-lyxofuranose, alpha-D-lyxopyranose (1)C(4), alpha-D-lyxopyranose (4)C(1), alpha-D-glucopyranose (4)C(1), and alpha-D-glucofuranose have been studied at ab initio and density-functional theory levels using TZVP quality basis set. The methods were tested by calculating the nuclear magnetic shieldings for tetramethylsilane (TMS) at different levels of theory using large basis sets. Test calculations on the monosaccharides showed B3LYP(TZVP) and BP86(TZVP) to be cost-efficient levels of theory for calculation of NMR chemical shifts of carbohydrates. The accuracy of the molecular structures and chemical shifts calculated at the B3LYP(TZVP) level is comparable to those obtained at the MP2(TZVP) level. Solvent effects were considered by surrounding the saccharides by water molecules and also by employing a continuum solvent model. None of the applied methods to consider solvent effects was successful. The B3LYP(TZVP) and MP2(TZVP)(13)C NMR chemical shift calculations yielded without solvent and rovibrational corrections an average deviation of 5.4 ppm and 5.0 ppm between calculated and measured shifts. A closer agreement between calculated and measured chemical shifts can be obtained by using a reference compound that is structurally reminiscent of saccharides such as neat methanol. An accurate shielding reference for carbohydrates can be constructed by adding an empirical constant shift to the calculated chemical shifts, deduced from comparisons of B3LYP(TZVP) or BP86(TZVP) and measured chemical shifts of monosaccharides. The systematic deviation of about 3 ppm for O(1)H chemical shifts can be designed to hydrogen bonding, whereas solvent effects on the (1)H NMR chemical shifts of C(1)H were found to be small. At the B3LYP(TZVP) level, the barrier for the torsional motion of the hydroxyl group at C(6) in alpha-D-glucofuranose was calculated to 7.5 kcal mol(-1). The torsional displacement was found to introduce large changes of up to 10 ppm to the (13)C NMR chemical shifts yielding uncertainties of about +/-2 ppm in the chemical shifts.     

NMR spectra of salicylohydroxamic acid in DMSO-d6 solution: a DFT study

The ¹H, ¹³C and ¹H, ¹³C COSY NMR spectra of salicylohydroxamic acid (sha) were measured in DMSO-d₆ solution. The B3LYP GIAO method with the 6-311++G(d,p) basis set was chosen to reproduce the experimental spectra. All possible zusammen and entgegen conformers of monomeric sha were computed. After geometry optimisation (B3LYP/6-311++G(d,p)) only nine independent models of the molecule were shown to be stable. Additionally, the NMR chemical shifts of the Onsager model of the most stable monomer were calculated. The computed chemical shifts for the labile protons for all aforementioned geometries meaningfully underestimated experimental results suggesting the existence of the H-bonded structure of sha in DMSO solution. The most probable two dimeric structures along with two solvent-bounded aggregates were subsequently calculated at the same level of theory. The best agreement was obtained for sha H-bonded with two DMSO molecules (confirmed by the absence of concentration effect). The relative error not exceeding 10 and 4% for chemical shifts in ¹H and ¹³C NMR spectra of sha–(DMSO)₂, respectively, showed that the applied method with the B3LYP/6-311++G(d,p) basis set was efficient to predict the NMR shifts of a compound with strong H-bonds. Thus, this allows to assign properly NMR resonances to specific structure formed in DMSO solution.     

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1) H and (13) C?NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13) C and (1) H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.     

Orientational ordering of a nematic liquid crystal and its mixture with its chain-perfluorinated analogue

The orientation of different segments of 4'-cyanophenyl 4-heptylbenzoate (7CPB) has been investigated using ¹³C NMR. The method of proton-encoded local field (PELF) spectroscopy in combination with off-magic-angle spinning (OMAS) of the sample was used. High resolution 2D spectra were obtained, from which the order parameters were calculated. Linear relations between the obtained order parameters and anisotropic chemical shifts determined by 1D ¹³C NMR were established and semi-empirical parameters were achieved. A 1:2 mixture of 7CPB and its chain-perfluorinated analogue (7PFCPB) showed interesting phase behaviour with the change of temperature. It was studied by the use of ¹³C NMR and polarizing optical microscopy. The order parameters of 7CPB in the smectic A phase of the mixture were calculated using the semi-empirical parameters obtained from the 2D NMR method.     

1H chemical shifts in NMR: Part 22-Prediction of the 1H chemical shifts of alcohols, diols and inositols in solution, a conformational and solvation investigation

The (1)H NMR spectra of a number of alcohols, diols and inositols are reported and assigned in CDCl(3), D(2)O and DMSO-d(6) (henceforth DMSO) solutions. These data were used to investigate the effects of the OH group on the (1)H chemical shifts in these molecules and also the effect of changing the solvent. Inspection of the (1)H chemical shifts of those alcohols which were soluble in both CDCl(3) and D(2)O shows that there is no difference in the chemical shifts in the two solvents, provided that the molecules exist in the same conformation in the two solvents. In contrast, DMSO gives rise to significant and specific solvation shifts. The (1)H chemical shifts of these compounds in the three solvents were analysed using the CHARGE model. This model incorporates the electric field, magnetic anisotropy and steric effects of the functional group for long-range protons together with functions for the calculation of the two- and three-bond effects. The long-range effect of the OH group was quantitatively explained without the inclusion of either the C--O bond anisotropy or the C--OH electric field. Differential beta and gamma effects for the 1,2-diol group needed to be included to obtain accurate chemical shift predictions. For DMSO solution the differential solvent shifts were calculated in CHARGE on the basis of a similar model, incorporating two-bond, three-bond and long-range effects. The analyses of the (1)H spectra of the inositols and their derivatives in D(2)O and DMSO solution also gave the ring (1)H,(1)H coupling constants and for DMSO solution the CH--OH couplings and OH chemical shifts. The (1)H,(1)H coupling constants were calculated in the CHARGE program by an extension of the cos(2)phi equation to include the orientation effects of electronegative atoms and the CH--OH couplings by a simple cos(2)phi equation. Comparison of the observed and calculated couplings confirmed the proposed conformations of myo-inositol, chiro-inositol, quebrachitol and allo-inositol. The OH chemical shifts were also calculated in the CHARGE program. Comparison of the observed and calculated OH chemical shifts and CH.OH couplings suggested the existence of intramolecular hydrogen bonding in a myo-inositol derivative.     

双氮桥联1,10-菲咯啉大环化合物分子构型的从头算研究

应用规范不变原子轨道法(GIAO)在RHF/6-31G**和B3LYP/6-31G**水平上计算了质子化双氮桥联1,10-菲咯啉大环化合物(H4HAPP2+)C2h和C2h构型的1HNMR,并用TDDFT法计算了H4HAPP2+电子光谱.结果表明,B3LYP/6-31G*优化的C2h构型为较优构型,经谐振频率验证无虚频,C2h构型是H4HAPP2+合理的对称性构型.     

Accurate prediction of proton chemical shifts. II. Peptide analogues

Proton chemical shifts of eight cyclic amide molecules were measured in DMSO and D2O solutions. The magnetic shieldings of the corresponding aliphatic, aromatic, and amide protons were calculated by Hartree-Fock and DFT, using the 6-311G**, 6-311++G**, and TZVP basis sets. For aliphatic protons, all of these methods reproduce the experimental values in DMSO solutions excellently after linear regression. The Hartree-Fock method tends to give slightly better agreement than DFT. The best performance is given by the HF/6-311G** method, with an rms deviation of 0.068 ppm. The deviations from experimental chemical shifts in D2O solutions are only slightly larger than those in DMSO solutions. This suggests that we can use the calculated gas phase proton chemical shifts directly to predict experimental data in various solvents, including water. For amide protons, which exchange with water and form hydrogen bonds with DMSO, only modest agreement is obtained, as expected. The present studies confirm that the GIAO approach can reach high accuracy for the relative chemical shifts of aliphatic and aromatic protons at a low cost. Such calculations may provide constraints for the conformational analysis of proteins and other macromolecules.     

Experimental and theoretical NMR and IR studies of the side-chain orientation effects on the backbone conformation of dehydrophenylalanine residue

Conformation of N-acetyl-(E)-dehydrophenylalanine N', N'-dimethylamide (Ac-(E)-ΔPhe-NMe(2)) in solution, a member of (E)-α, β-dehydroamino acids, was studied by NMR and infrared spectroscopy and the results were compared with those obtained for (Z) isomer. To support the spectroscopic interpretation, the Φ, Ψ potential energy surfaces were calculated at the MP2/6-31 + G(d,p) level of theory in chloroform solution modeled by the self-consistent reaction field-polarizable continuum model method. All minima were fully optimized by the MP2 method and their relative stabilities were analyzed in terms of π-conjugation, internal H-bonds and dipole interactions between carbonyl groups. The obtained NMR spectral features were compared with theoretical nuclear magnetic shieldings, calculated using Gauge Independent Atomic Orbitals (GIAO) approach and rescaled to theoretical chemical shifts using benzene as reference. The calculated indirect nuclear spin-spin coupling constants were compared with available experimental parameters.     

Quantum chemical and nuclear magnetic resonance spectral studies on molecular properties and electronic structure of berberine and berberrubine

The structural and electronic properties of berberine and berberrubine have been studied extensively using density functional theory (DFT) employing B3LYP exchange correlation. The geometries of these molecules have been fully optimized at the B3LYP/6-311G** level. The chemical shift of 1H and 13C resonances in NMR spectra of these molecules have been calculated using the gauge invariant atomic model (GIAO) method as implemented in Gaussian 98. One- and two-dimensional HSQC (1H-13C), HMBC (1H-13C) and ROESY (1H-1H) spectra were recorded at 500 MHz for the berberine molecule in D(2)O solution. All proton and carbon resonances were unambiguously assigned, and inter-proton distances obtained from ten observed NOE contacts. A restrained molecular dynamics (RMD) approach was used to get the optimized solution structure of berberine. The structure of berberine and berberrubine molecules was also obtained using the ROESY data available in literature. Comparison of the calculated NMR chemical shifts with the experimental values revealed that DFT methods produce very good results for both proton and carbon chemical shifts. The importance of the basis sets to the calculated NMR parameters is discussed. It has been found that calculated structure and chemical shifts in the gas phase predicted with B3LYP/6-311G** are in very good agreement with the present experimental data and the measured values reported earlier.     

Effects of Quantum Nuclear Delocalisation on NMR Parameters from Path Integral Molecular Dynamics

The influence of nuclear delocalisation on NMR chemical shifts in molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory calculations of shielding tensors. Nuclear quantum effects are shown to explain previously observed systematic deviations in correlations between calculated and experimental chemical shifts, with particularly large PIMD‐induced changes (up to 23 ppm) observed for carbon atoms in methyl groups. The PIMD approach also enables isotope substitution effects on chemical shifts and J couplings to be predicted in excellent agreement with experiment for both isolated molecules and molecular crystals. An approach based on convoluting calculated shielding or coupling surfaces with probability distributions of selected bond distances and valence angles obtained from PIMD simulations is used to calculate isotope effects.     

Structure-NMR chemical shift relationships for novel functionalized derivatives of quinoxalines

(1)H, (13)C and (15)N NMR chemical shifts for a variety of novel quinoxalines were determined by different 2D methods and were calculated using the GIAO DFT approach. Comparison with experimental data shows good correlations in the case of (1)H, (13)C and (15)N chemical shifts. Different combinations of basis sets were tested. In non-polar solvents quinoxalines exist as dimers owing to strong hydrogen bonding. Calculations for dimers improve the correlation between experiment and theory. Additive empirical methods for estimating chemical shifts have drawbacks and have to be used with a great care for this type of compound.     

29/28, 30/28Silicon, 34/32sulfur and 80/77selenium isotope-induced fluorine chemical shifts through two bonds

Isotope effects on fluorine chemical shifts induced by heteroatoms bonded covalently to a carbon atom bearing fluorine atoms were studied. For each compound, the isotope-induced chemical shifts 2delta19F(X) through two bonds were measured for the heteroatom (X = 29/28Si, 30/28Si, 34/32S and 80/77Se). The 1delta19F(13/12C) values for the carbon bonded to the fluorine atoms were also recorded. Examination of the 19F NMR data showed homogeneity of the isotope-induced chemical shifts along the rows of the periodic table and regularity down the columns (from 10 to 15 ppb per mass unit for the second row to 0.4 ppb for the fourth row). It became negligible for atoms of the fifth row.     

Quantum chemical computational studies on 5-(4-bromophenylamino)-2-methylsulfanylmethyl-2H-1,2,3-triazol-4-carboxylic acid ethyl ester

The optimized molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) (1)H and (13)C NMR shift values of 5-(4-bromophenylamino)-2-methylsulfanylmethyl-2H-1,2,3-triazol-4-carboxylic acid ethyl ester have been calculated by using Hartree-Fock (HF) and density functional method (DFT/B3LYP) with 6-31G(d), 6-31G(d,p) and LANL2DZ basis sets. The optimized molecular geometric parameters were presented and compared with the data obtained from X-ray diffraction. In order to fit the calculated harmonic wavenumbers to the experimentally observed ones, scaled quantum mechanics force field (SQM FF) methodology was proceeded. Correlation factors between the experimental and calculated (1)H chemical shift values of the title compound in vacuum and in CHCl(3) solution by using the conductor-like screening continuum solvation model (COSMO) were reported. The calculated results showed that the optimized geometry well reproduces the crystal structure. The theoretical vibrational frequencies and chemical shifts are in very good agreement with the experimental data. In solvent media the energetic behavior of the title compound was also examined by using the B3LYP method with the 6-31G(d) basis set, applying the COSMO model. The obtained results indicated that the total energy of the title compound decreases with increasing polarity of the solvent. Furthermore, molecular electrostatic potential (MEP), natural bond orbital (NBO) and frontier molecular orbitals (FMOs) of the title compound were performed by the B3LYP/LANL2DZ method, and also thermodynamic parameters for the title compound were calculated at all the HF and B3LYP levels.     

Ab initio calculations of the O1s XPS spectra of ZnO and Zn oxo compounds

O1s core level binding energies of oxygen atoms in bulk ZnO, at different ZnO surfaces, and in some Zn oxo compounds were calculated by means of wave function based quantum chemical ab initio methods. Initial and final state effects were obtained by Koopmans' theorem and at the DeltaSCF level, respectively. After correction for scalar relativistic effects and electron correlation, the calculated XPS peak positions are in excellent agreement with the available experimental data for all systems included in the present study. The O1s core level shifts between an isolated H2O molecule and the Zn oxo compounds or ZnO, as well as between oxygen atoms in bulk ZnO and at various ZnO surfaces, can be understood by means of Madelung potentials and electronic relaxation or screening. XPS spectra were calculated for various cluster models which are designed to describe different possibilities of stabilizing the polar O-terminated ZnO(0001) surface by the adsorption of H atoms. The experimental spectra are only compatible with the theoretical results for the fully hydroxylated H-ZnO(0001) surface exhibiting a (1x1) surface structure.     

Spectroscopic, electronic structure and natural bond orbital analysis of o-fluoronitrobenzene and p-fluoronitrobenzene: a comparative study

Experimental FTIR, FT-Raman and FT-NMR spectroscopic studies of o-fluoronitrobenzene and p-fluoronitrobenzene have been carried out. A detailed quantum chemical calculations have been performed using DFT/B3LYP method with 6-311++G** and 6-31G** basis sets. Complete vibrational analyses of the compounds were performed. The temperature dependence of thermodynamic properties has been analysed. The atomic charges, electronic exchange interaction and charge delocalisation of the molecule have been performed by natural bond orbital (NBO) analysis. Molecular electrostatic surface potential (MESP), total electron density distribution and frontier molecular orbitals (FMOs) are constructed at B3LYP/6-311++G** level to understand the electronic properties. The charge density distribution and site of chemical reactivity of the molecules have been obtained by mapping electron density isosurface with electrostatic potential surfaces (ESP). The electronic properties, HOMO and LUMO energies were measured by time-dependent TD-DFT approach. (1)H and (13)C NMR spectra were recorded and (1)H and (13)C nuclear magnetic resonance chemical shifts of the molecule were calculated. The (1)H and (13)C nuclear magnetic resonance (NMR) chemical shifts of the molecules in chloroform solvent and in gas phase were calculated by using the Gauge-Independent Atomic Orbital (GIAO) method and are found to be in good agreement with experimental values. The theoretical parameters obtained at B3LYP levels have been compared with the experimental values.     

Visualization of homoaromaticity in cations, neutral molecules and anions by spatial magnetic properties (through space NMR shieldings)—an H/C NMR chemical shift study

Prototypes for homoaromaticity in cations, neutral molecules, and anions are theoretically studied at the MP2 level of theory. For the global minimum structures on the potential energy surface both ¹H/¹³C chemical shifts and spatial magnetic properties as through space NMR shieldings (TSNMRS) were calculated by the GIAO perturbation method. The TSNMRS are visualized as iso-chemical-shielding surfaces (ICSS) of different sign and size. Coincident experimental and computed ¹H/¹³C chemical shifts afforded the possibility to decide from the TSNMRSs at hand on both the existence and the size of homoaromaticity in the molecules studied.     

Ab initio/GIAO-CCSD(T) study of structures, energies, and 13C NMR chemical shifts of C4H7(+) and C5H9(+) ions: relative stability and dynamic aspects of the cyclopropylcarbinyl vs bicyclobutonium ions

The structures and energies of the carbocations C 4H 7 (+) and C 5H 9 (+) were calculated using the ab initio method. The (13)C NMR chemical shifts of the carbocations were calculated using the GIAO-CCSD(T) method. The pisigma-delocalized bisected cyclopropylcarbinyl cation, 1 and nonclassical bicyclobutonium ion, 2 were found to be the minima for C 4H 7 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level the structure 2 is 0.4 kcal/mol more stable than the structure 1. The (13)C NMR chemical shifts of 1 and 2 were calculated by the GIAO-CCSD(T) method. Based on relative energies and (13)C NMR chemical shift calculations, an equilibrium involving the 1 and 2 in superacid solutions is most likely responsible for the experimentally observed (13)C NMR chemical shifts, with the latter as the predominant equilibrating species. The alpha-methylcyclopropylcarbinyl cation, 4, and nonclassical bicyclobutonium ion, 5, were found to be the minima for C 5H 9 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level ion 5 is 5.9 kcal/mol more stable than the structure 4. The calculated (13)C NMR chemical shifts of 5 agree rather well with the experimental values of C 5H 9 (+).     

Configuration and E/Z interconversion mechanism of O(S)-allyl-S(O)-methyl-N-(acridin-9-yl) iminothiocarbonate

The configuration and dynamic behavior of O-allyl-S-methyl-N-(acridin-9-yl)iminothiocarbonate (1) and its S-allyl-O-methyl regioisomer (2) were studied using quantum chemical calculations and by applying a novel graphical method to scatter maps obtained from MD simulations for evaluation of an NOE-weighted internuclear distance (r(NOE)). Energy calculations indicated that the Z configuration was predominant for each compound and, further, this was supported both by the calculated chemical shifts and the r(NOE). Both N-inversion- and rotation-type transition-state structures were also calculated for the E/Z isomerization process, the results indicating that the preferred interconversion mechanism for 1 is N-inversion, but contrastingly, interconversion via rotation is equally as probable as N-inversion for 2. This supports the notion that one or the other or both pathways can be active and each system needs to be assessed on a case-by-case basis.     

Towards Accurate Predictions of Proton NMR Spectroscopic Parameters in Molecular Solids

The factors contributing to the accuracy of quantum-chemical calculations for the prediction of proton NMR chemical shifts in molecular solids are systematically investigated. Proton chemical shifts of six solid amino acids with hydrogen atoms in various bonding environments (CH, CH₂, CH₃, OH, SH and NH₃) were determined experimentally using ultra-fast magic-angle spinning and proton-detected 2D NMR experiments. The standard DFT method commonly used for the calculations of NMR parameters of solids is shown to provide chemical shifts that deviate from experiment by up to 1.5 ppm. The effects of the computational level (hybrid DFT functional, coupled-cluster calculation, inclusion of relativistic spin-orbit coupling) are thoroughly discussed. The effect of molecular dynamics and nuclear quantum effects are investigated using path-integral molecular dynamics (PIMD) simulations. It is demonstrated that the accuracy of the calculated proton chemical shifts is significantly better when these effects are included in the calculations.     

Oxoperoxo vanadium(V) complexes of L-lactic acid: density functional theory study of structure and NMR chemical shifts

Various combinations of density functionals and pseudopotentials with associated valence basis-sets are compared for reproducing the known solid-state structure of [V 2O 2(OO) 2 l-lact 2] (2-) cis . Gas-phase optimizations at the B3LYP/SBKJC level have been found to provide a structure that is close to that seen in the solid state by X-ray diffraction. Although this may result in part from error compensation, this optimized structure allowed satisfactory reproduction of solution multinuclear NMR chemical shifts of the complex in all-electron DFT-IGLO calculations (UDFT-IGLO-PW91 level), suggesting that it is probably close to that found in solution. This combination of approaches has subsequently been used to optimize the structures of the vanadium oxoperoxo complexes [V 2O 3(OO) l-lact 2] (2-) cis , [V 2O 3(OO) l-lact 2] (2-) trans , and [VO(OO)( l-lact)(H 2O)] (-) cis . The (1)H, (13)C, (51)V, and (17)O NMR chemical shifts for these complexes have been calculated and compared with the experimental solution chemical shifts. Excellent agreement is seen with the (13)C chemical shifts, while somewhat inferior agreement is found for (1)H shifts. The (51)V and (17)O chemical shifts of the dioxo vanadium centers are well reproduced, with differences between theoretical and experimental shifts ranging from 22.9 to 35.6 ppm and from 25.1 to 43.7 ppm, respectively. Inferior agreement is found for oxoperoxo vanadium centers, with differences varying from 137.3 to 175.0 ppm for (51)V shifts and from 148.7 to 167.0 ppm for (17)O(oxo) shifts. The larger errors are likely to be due to overestimated peroxo O-O distances. The chosen methodology is able to predict and analyze a number of interesting structural features for vanadium(V) oxoperoxocomplexes of alpha-hydroxycarboxylic acids.