Chiral discrimination via nuclear magnetic shielding polarisabilities from NMR spectroscopy: theoretical study of (R(a))-1,3-dimethylallene, (2R)-2-methyloxirane, and (2R)-N-methyloxaziridine

Three medium-size optically active molecules have been studied to make a guess at candidates suitable for chiral discrimination in an isotropic medium via nuclear magnetic resonance spectroscopy. The criterion for experimental detection is given by the magnitude of the isotropic part of nuclear magnetic shielding polarisability tensors, related to a pseudoscalar of opposite sign for the two enantiomers. The pseudoscalar shielding polarisability at the (17)O nucleus in N-methyloxaziridine, calculated at the Hartree-Fock level, is approximately 7.8 x10(-)(17) mV(-)(1). To obtain an experimentally observable magnetic field induced at the (17)O nucleus in N-methyloxaziridine, electric fields as large as approximately 10(7) - 10(8) Vm(-)(1) should be applied to the probe. The molecular electric dipole moment induced by precession of the magnetic dipole of the (17)O nucleus in a magnetic field of 10 T is, in absolute value, approximately 8.8 x 10(-)(42) Cm. The estimated rf-voltage at a resonance circuit is approximately 10 nV. Smaller values have been estimated for N, C, and H nuclei in 1,3-dimethylallene and 2-methyloxirane.     

Nuclear magnetic shielding constants of liquid water: insights from hybrid quantum mechanics/molecular mechanics models

We present a gauge-origin independent method for the calculation of nuclear magnetic shielding tensors of molecules in a structured and polarizable environment. The method is based on a combination of density functional theory (DFT) or Hartree-Fock wave functions with molecular mechanics. The method is unique in the sense that it includes three important properties that need to be fulfilled in accurate calculations of nuclear magnetic shielding constants: (i) the model includes electron correlation effects, (ii) the model uses gauge-including atomic orbitals to give gauge-origin independent results, and (iii) the effect of the environment is treated self-consistently using a discrete reaction-field methodology. The authors present sample calculations of the isotropic nuclear magnetic shielding constants of liquid water based on a large number of solute-solvent configurations derived from molecular dynamics simulations employing potentials which treat solvent polarization either explicitly or implicitly. For both the (17)O and (1)H isotropic shielding constants the best predicted results compare fairly well with the experimental data, i.e., they reproduce the experimental solvent shifts to within 4 ppm for the (17)O shielding and 1 ppm for the (1)H shielding.     

Density functional theory investigation of hydrogen bonding effects on the oxygen, nitrogen and hydrogen electric field gradient and chemical shielding tensors of anhydrous chitosan crystalline structure

A systematic computational investigation was carried out to characterize the 17O, 14N and 2H electric field gradient, EFG, as well as 17O, 15N, 13C and 1H chemical shielding tensors in the anhydrous chitosan crystalline structure. To include the hydrogen-bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a hexameric cluster. The computations were performed with the B3LYP method and 6-311++G(d,p) and 6-31++G(d,p) standard basis sets using the Gaussian 98 suite of programs. Calculated EFG and chemical shielding tensors were used to evaluate the 17O, 14N and 2H nuclear quadrupole resonance, NQR, and 17O, 15N, 13C and 1H nuclear magnetic resonance, NMR, parameters in the hexameric cluster, which are in good agreement with the available experimental data. The difference between the calculated NQR and NMR parameters of the monomer and hexamer cluster shows how much hydrogen bonding interactions affect the EFG and chemical shielding tensors of each nucleus. These results indicate that both O(3)-H(33)...O(5-3) and N-H(22)...O(6-4) hydrogen bonding have a major influence on NQR and NMR parameters. Also, the quantum chemical calculations indicate that the intra- and intermolecular hydrogen bonding interactions play an essential role in determining the relative orientation of EFG and chemical shielding principal components in the molecular frame axes.     

Leading-order relativistic effects on nuclear magnetic resonance shielding tensors

We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.     

Parity nonconservation contribution to the nuclear magnetic resonance shielding constants of chiral molecules: a four-component relativistic study

A systematic four-component relativistic study of the parity nonconservation (PNC) contribution to the (isotropic) NMR shielding constants of chiral molecules is presented for the P enantiomers of the series H(2)X(2) (X=(17)O,(33)S,(77)Se,(125)Te,(209)Po). The PNC contributions are obtained within a linear response approach at the Hartree-Fock level. A careful design of the basis sets is necessary. The four-component relativistic results based on the Dirac-Coulomb Hamiltonian are compared with the nonrelativistic Levy-Leblond results and those obtained by the spin-free modified Dirac Hamiltonian. The calculations confirm the nonrelativistic scaling law Z(2.4) of the PNC contribution with respect to nuclear charge Z. However, the calculations also show that the overall scaling is significantly modified by relativistic effects. The scalar relativistic effect scales as Z(4.7) for the selected set of molecules, whereas the spin-orbit effect, of opposite sign, scales better than Z(6) and completely dominates the PNC contribution for the heaviest elements. This opens up the intriguing possibility of the experimental observation of PNC effects on NMR parameters of molecules containing heavy atoms. The presented formalism is expected to be valuable in assisting the search for suitable candidate molecules.     

Direct chiral discrimination in NMR spectroscopy

Conventional nuclear magnetic resonance spectroscopy is unable to distinguish between the two mirror-image forms (enantiomers) of a chiral molecule. This is because the NMR spectrum is determined by the chemical shifts and spin–spin coupling constants which – in the absence of a chiral solvent – are identical for the two enantiomers. We discuss how chirality may nevertheless be directly detected in liquid-state NMR spectroscopy: In a chiral molecule, the rotating nuclear magnetic moment induces an electric dipole moment in the direction perpendicular to itself and to the permanent magnetic field of the spectrometer. We present computations of the precessing electric polarization following a π/2 pulse. Our estimates indicate that the electric polarization should be detectable in favourable cases. We also predict that application of an electrostatic field induces a chirally sensitive magnetization oscillating in the direction of the permanent magnetic field. We show that the electric-field-perturbed chemical shift tensor, the nuclear magnetic shielding polarizability, underlies these chiral NMR effects.     

Quadrupole central transition 17O NMR spectroscopy of biological macromolecules in aqueous solution

We demonstrate a general nuclear magnetic resonance (NMR) spectroscopic approach in obtaining high-resolution (17)O (spin-5/2) NMR spectra for biological macromolecules in aqueous solution. This approach, termed quadrupole central transition (QCT) NMR, is based on the multiexponential relaxation properties of half-integer quadrupolar nuclei in molecules undergoing slow isotropic tumbling motion. Under such a circumstance, Redfield's relaxation theory predicts that the central transition, m(I) = +1/2 ? -1/2, can exhibit relatively long transverse relaxation time constants, thus giving rise to relatively narrow spectral lines. Using three robust protein-ligand complexes of size ranging from 65 to 240 kDa, we have obtained (17)O QCT NMR spectra with unprecedented resolution, allowing the chemical environment around the targeted oxygen atoms to be directly probed for the first time. The new QCT approach increases the size limit of molecular systems previously attainable by solution (17)O NMR by nearly 3 orders of magnitude (1000-fold). We have also shown that, when both quadrupole and shielding anisotropy interactions are operative, (17)O QCT NMR spectra display an analogous transverse relaxation optimized spectroscopy type behavior in that the condition for optimal resolution depends on the applied magnetic field. We conclude that, with the currently available moderate and ultrahigh magnetic fields (14 T and higher), this (17)O QCT NMR approach is applicable to a wide variety of biological macromolecules. The new (17)O NMR parameters so obtained for biological molecules are complementary to those obtained from (1)H, (13)C, and (15)N NMR studies.     

Calculation of dipole-shielding polarizabilities (sigma(alphabetagamma)I): the influence of uniform electric field effects on the shielding of backbone nuclei in proteins

A significant contribution to the chemical shielding of a nucleus can arise from uniform electric fields that act to distort the electronic charge distribution surrounding a nucleus and, hence, affect the nuclear shielding. It has been shown by Buckingham (Buckingham, A. D. Can. J. Chem. 1960, 38, 300) that the nuclear magnetic shielding tensor sigmaalphabetaI of a nucleus in the presence of an external weak static uniform electric field E may be expanded using sigmaalphabetaI(E) = sigmaalphabetaI + sigmaalphabetagammaIEgamma + 1/2sigmaalphabetagammadeltaIEgammaEdelta + sigmaalphabeta,gammadeltaIEgammadelta ... The third rank tensor sigmaalphabetagammaI is referred to as the dipole-shielding polarizability and describes the nonlinear response of the electron cloud to first order in E, muI, and B0. We report calculations of sigmaalphabetagammaI for the N, HN, and C' nuclei in N-methyl acetamide (NMA) and show that these tensors can be used to provide considerable insight into the behavior of uniform electric fields upon the shielding of backbone nuclei in proteins. The sigmaalphabetagammaI values for the N, HN, and C' of NMA were calculated using the continuous transformation of the origin of the current density (CTOCD) scheme with the diamagnetic contribution set to zero (CTOCD-DZ). Values are given for the individual tensor components of sigmaalphabetagammaI for each nucleus. To test that the calculations have provided a reasonable estimate for the sigmaalphabetagammaI of N, HN, and C' nuclei in proteins, a pH titration was performed using Hen Lysozyme (HEWL). The pH-induced isotropic shielding changes for the C', N, and HN nuclei in some peptide bonds close to E35 ( approximately <8 A) were extracted from sets of fitted titration curves. Assuming the experimental shielding changes arise solely from uniform electric field effects caused by the deprotonation of E35, without any other pH-induced structural alterations which might lead to a shielding change, the experimental shielding differences were compared to those calculated via the product Agamma(I).Egamma where Agamma(I) = (1/3)sigmaalphaalphagammaIota. The agreement with the experimental data is in many cases reasonable and suggests that, within the Buckingham formalism, the complete sigmaalphabetagammaI tensors reported here will be helpful to resolve the importance of uniform electric fields upon isotropic and anisotropic shielding in proteins and their complexes.     

Methodological aspects in the calculation of parity-violating effects in nuclear magnetic resonance parameters

We examine the quantum chemical calculation of parity-violating (PV) electroweak contributions to the spectral parameters of nuclear magnetic resonance (NMR) from a methodological point of view. Nuclear magnetic shielding and indirect spin-spin coupling constants are considered and evaluated for three chiral molecules, H2O2, H2S2, and H2Se2. The effects of the choice of a one-particle basis set and the treatment of electron correlation, as well as the effects of special relativity, are studied. All of them are found to be relevant. The basis-set dependence is very pronounced, especially at the electron correlated ab initio levels of theory. Coupled-cluster and density-functional theory (DFT) results for PV contributions differ significantly from the Hartree-Fock data. DFT overestimates the PV effects, particularly with nonhybrid exchange-correlation functionals. Beginning from third-row elements, special relativity is of importance for the PV NMR properties, shown here by comparing perturbational one-component and various four-component calculations. In contrast to what is found for nuclear magnetic shielding, the choice of the model for nuclear charge distribution--point charge or extended (Gaussian)--has a significant impact on the PV contribution to the spin-spin coupling constants.     

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH₂O, H₂O, HF, F₂, HCN, SiH₄ and H₂S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS‐2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization‐consistent pcS‐n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS‐2 and aug‐cc‐pVTZ‐J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS‐2, MP2/pcS‐2 and DFT/CBS calculations with pcS‐n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

Nuclear magnetic shielding and quadrupole coupling tensors in liquid water: a combined molecular dynamics simulation and quantum chemical study

Nuclear magnetic resonance (NMR) shielding tensors for the oxygen and hydrogen nuclei, as well as nuclear quadrupole coupling tensors for the oxygen and deuterium nuclei of water in the liquid and gaseous state, are calculated using Hartree-Fock and density functional theory methods, for snapshots sampled from Car-Parrinello molecular dynamics trajectories. Clusters representing local liquid structures and instantaneous configurations of a single molecule representing low-density gas are fed into a quantum chemical program for the calculation of the NMR tensors. The average isotropic and anisotropic tensorial properties of 400 samples in both states, averaged using a common Eckart coordinate frame, are calculated from the data. We report results for the gas-to-liquid chemical shifts of (17)O and (1)H nuclei, as well as the corresponding change in the nuclear quadrupole couplings of (17)O and (2)H. Full thermally averaged shielding and quadrupole coupling tensors are reported for the gaseous and liquid-state water, for the first time in the case of liquid. Electron correlation effects, the difference of classical vs quantum mechanical rovibrational averaging, and different methods of averaging anisotropic properties are discussed.     

Mechanistic study on the opposite migration order of clenbuterol enantiomers in capillary electrophoresis with beta-cyclodextrin and single-isomer heptakis(2,3-diacetyl-6-sulfo)-beta-cyclodextrin

Opposite migration order was observed for the enantiomers of the chiral beta2-adrenergic drug clenbuterol (CL) in capillary electrophoresis (CE) when resolved with native beta-cyclodextrin (beta-CD) and heptakis (2,3-diacetyl-6-sulfo)-beta-CD (HDAS-beta-CD). The possible mechanisms of the affinity reversal of the CL enantiomers depending on the structure of the CD were studied using 1H-nuclear magnetic resonance (1H-NMR) spectrometry and one-dimensional rotating frame nuclear Overhauser and exchange spectrometry (1-D ROESY). Significant differences were observed between the structure of the (+/-)-CL complexes with beta-CD and HDAS-beta-CD.     

First-principles calculations of solid-state (17)O and (29)Si NMR spectra of Mg(2)SiO(4) polymorphs

The nuclear magnetic resonance (NMR) shielding and electric field gradient (EFG) tensors of three polymorphs of Mg(2)SiO(4), forsterite (alpha-Mg(2)SiO(4)), wadsleyite (beta-Mg(2)SiO(4)) and ringwoodite (gamma-Mg(2)SiO(4)), have been calculated using a density functional theory (DFT) approach with a planewave basis set and pseudopotential approximation. These Mg(2)SiO(4) polymorphs are the principal components of the Earth down to depths of 660 km and have been proposed as the hosts of water in the Earth's upper mantle and transition zone. A comparison of our calculations with single-crystal spectroscopic data in the literature for the alpha-polymorph, forsterite, shows that both the magnitude and orientation of the shielding and EFG tensors for O and Si can be obtained with sufficient accuracy to distinguish subtle differences in atomic positions between published structures. We compare calculated (17)O MAS NMR quadrupolar powder lineshapes directly with experimental lineshapes and show that we are able to reproduce them within the precision with which the NMR parameters may be determined from multi-parameter fitting. The relatively small amounts of sample available for the beta- and gamma-polymorphs, arising from the high pressures required for synthesis, has hindered the extraction of NMR parameters in previous work. The application of DFT calculations to these high-pressure polymorphs confirms previous spectral assignments, and provides deeper insight into the empirical correlations and observations reported in the literature. These first-principles methods are highly promising for the determination of local bonding in more complex materials, such as the hydrated forms of Mg(2)SiO(4), by aiding analysis of their multinuclear NMR spectra.     

Phosphoramides: Synthesis,Spectroscopy, and X‐ray Crystallography

A series of new phosphoramides with general formula RP(O)X₂, where R = amino/p‐methylphenoxy and X = amine, were synthesized and characterized by ¹H, ¹³C, ³¹P nuclear magnetic resonance (NMR), and infrared (IR) spectroscopy and elemental analysis. The ³¹P{¹H} NMR spectra show that among compounds 7–9 containing 2‐, 3‐, and 4‐aminopyridinyl moieties, respectively, the shielding order of the P atom decreases as 7 > 9 > 8 . Also, the structure of compound 7 was determined by X‐ray crystallography. In this structure, repeated noncentrosymmetric dimers are formed by two strong intermolecular N(1)‐H(1N)…N(2) and N(3)‐H(3N)…O(1) hydrogen bonds. Taking into account weak intermolecular C(17)‐H(17C)…N(4), C(17)‐H(17E)…N(4), C(2)‐H(2A)…O(2), and also weak aromatic C—H…C interactions, a three‐dimensional polymeric chain is created in the crystalline network. The density functional theory calculations at B3LYP, B3PW91, and M06 levels using the 6–31+G** basis set were in good agreement with the X‐ray crystallography data.     

Study of hydrogen bonds in N-methylacetamide by DFT calculations of oxygen, nitrogen, and hydrogen solid-state NMR parameters

Solid-state nuclear magnetic resonance (NMR) parameters of ¹⁷O, ¹⁴N/¹⁵N, and ²H/¹H nuclei were evaluated in two available neutron crystalline structures of N-methylacetamide (NMA) at 250 and 276 K, NMA-I and NMA-II, respectively. Density functional theory calculations were performed by B3LYP method and 6-311++G** and IGLO-II type basis sets to calculate the electric field gradient (EFG) and chemical shielding (CS) tensors at the sites of mentioned nuclei. In order to investigate hydrogen bonds (HBs) effects on NMR tensors, calculations were performed on four-model systems of NMA: an optimized isolated gas-phase, crystalline monomers, crystalline dimers, and crystalline trimers. Comparing the calculated results reveal the influence of N–H···O=C and C–H···O=C HB types on the NMR tensors which are observable by the evaluated parameters including quadrupole coupling constant, C _Q, and isotropic CS, σ _iso. Furthermore, the results demonstrate more influence of HB on the EFG and CS tensors of NMA at 276 K rather than that of 250 K.     

Molecular magnetic properties of two-copper(II) containing complexes [Cu(II) (1-phenylamidino-O-n-butylurea) en (H2O)]2(2+) and [Cu(II) sulphato-mono (1-phenylamidino-O-methylurea)]2. An EPR study

Electron paramagnetic resonance (EPR) investigations were conducted on [Cu(II) (1-phenylamidino-O-n-butylurea) en (H2O)]2(2+) (1) and [Cu(II) sulphato-mono (1-phenylamidino-O-methylurea)]2 (2) respectively, in the temperature range 300-77K. Fine structure characteristics of S = 1 system, was observed in both complexes with zero field splitting of 0.0525 and 0.0225 cm(-1), respectively, suggesting the formation of dimeric complexes. The presence of the half-field signal (DeltaMs= +/-2), in the complex 1, further confirmed the formation of dimer. The temperature dependence of EPR signal intensity has given evidence for the ferromagnetic (FM) coupling between the two Cu2+ ions. The isotropic exchange interaction constants J, were evaluated from this and were found out to be approximately 57 and approximately 27 cm(-1), respectively, for the complexes 1 and 2. The photoacoustic spectra of these complexes had shown a band around 26,400 cm(-1) characteristic of metal-metal bonding giving an independent support for the existence of dimeric Cu2+ species. The high magnetic moment values at room temperature for complex 1 (2.68 microB) and complex 2 (2.00 microB), obtained from the magnetic susceptibility measurements, support the formation of ferromagnetically coupled Cu2+ dimers.     

Quantum chemical investigation of hyperfine coupling constants on first coordination sphere water molecule of gadolinium(III) aqua complexes

Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H(2)O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H(2)O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive (17)O isotropic hyperfine coupling constant (A(iso)((17)O) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (r(eff)(Gd-O) = 2.72 +/- 0.06 A) of dipolar interaction compared to the mean internuclear distance (r(Gd-O) = 2.56 +/- 0.06 A) obtained from the MD trajectory of [Gd(DOTA)(H(2)O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the (17)O nucleus by approximately 45%. The (1)H isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)((1)H) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.     

Electric quadrupole polarizabilities of nuclear magnetic shielding in some small molecules

Computational procedures, based on (i) the Ramsey common origin approach and (ii) the continuous transformation of the origin of the quantum mechanical current density-diamagnetic zero (CTOCD-DZ), were applied at the Hartree-Fock level to determine electric quadrupole polarizabilities of nuclear magnetic shielding for molecules in the presence of a nonuniform electric field with a uniform gradient. The quadrupole polarizabilities depend on the origin of the coordinate system, but values of the magnetic field induced at a reference nucleus, determined via the CTOCD-DZ approach, are origin independent for any calculations relying on the algebraic approximation, irrespective of size and quality of the (gaugeless) basis set employed. On the other hand, theoretical estimates of the induced magnetic field obtained by single-origin methods are translationally invariant only in the limit of complete basis sets. Calculations of electric quadrupole polarizabilities of nuclear magnetic shielding are reported for H(2), HF, H(2)O, NH(3), and CH(4) molecules.     

Mechanistic study of the enantiomeric recognition of a basic compound with negatively charged single-isomer gamma-cyclodextrin derivatives using capillary electrophoresis,nuclear magnetic resonance spectroscopy,and infrared spectroscopy

The possible mechanisms for the chiral recognition of 2-(R)-N-[1-(6-aminopyridin-2-ylmethyl)piperidin-4-yl]-2-[(1R)-3,3-difluorocyclopentyl]-2-hydroxy-2-phenylacetamide (RR-M3), and its enantiomer (SS-M3) with octakis(2,3-di-O-acetyl-6-sulfo)-gamma-cyclodextrin (ODAS-gamma-CD) and octakis(6-sulfo)-gamma-cyc?pdextrom enantiomer; (OS-gamma-CD), were investigated using capillary electrophoresis (CE), proton ((1)H), fluorine ((19)F) and carbon ((13)C) nuclear magnetic resonance spectroscopy (NMR), and infrared (IR) spectroscopy. Clear evidence for the formation of diastereomeric complexes between the enantiomers and the two CDs was observed. NMR spectra suggest that the phenyl and difluorocyclopentyl rings are involved in the complexation. The phenyl ring on the guest molecule is deeply penetrated into the cavity of OS-gamma-CD, but it is not included into the cavity of ODAS-gamma-CD. The continuous variation plots built based on the (1)H NMR and IR spectra indicate a 1:1 complex stoichiometric ratio of the M3 enantiomers for both CDs. The affinity of the enantiomers for the two CDs is opposite.