Natural-abundance solid-state 2H NMR spectroscopy at high magnetic field

High-resolution solid-state (2)H NMR spectroscopy provides a method for measuring (1)H NMR chemical shifts in solids and is advantageous over the direct measurement of high-resolution solid-state (1)H NMR spectra, as it requires only the application of routine magic angle sample spinning (MAS) and routine (1)H decoupling methods, in contrast to the requirement for complex pulse sequences for homonuclear (1)H decoupling and ultrafast MAS in the case of high-resolution solid-state (1)H NMR. However, a significant obstacle to the routine application of high-resolution solid-state (2)H NMR is the very low natural abundance of (2)H, with the consequent problem of inherently low sensitivity. Here, we explore the feasibility of measuring (2)H MAS NMR spectra of various solids with natural isotopic abundances at high magnetic field (850 MHz), focusing on samples of amino acids, peptides, collagen, and various organic solids. The results show that high-resolution solid-state (2)H NMR can be used successfully to measure isotropic (1)H chemical shifts in favorable cases, particularly for mobile functional groups, such as methyl and -N(+)H(3) groups, and in some cases phenyl groups. Furthermore, we demonstrate that routine (2)H MAS NMR measurements can be exploited for assessing the relative dynamics of different functional groups in a molecule and for assessing whole-molecule motions in the solid state. The magnitude and field-dependence of second-order shifts due to the (2)H quadrupole interaction are also investigated, on the basis of analysis of simulated and experimental (1)H and (2)H MAS NMR spectra of fully deuterated and selectively deuterated samples of the α polymorph of glycine at two different magnetic field strengths.     

NMR nomenclature: nuclear spin properties and conventions for chemical shifts. IUPAC Recommendations 2001. International Union of Pure and Applied Chemistry. Physical Chemistry Division. Commission on Molecular Structure and Spectroscopy

An Addendum has been published for this article in Magnetic Resonance in Chemistry 40(9) 2002, 622. A unified scale is recommended for reporting the NMR chemical shifts of all nuclei relative to the ¹H resonance of tetramethylsilane. The unified scale is designed to provide a precise ratio, Ξ, of the resonance frequency of a given nuclide to that of the primary reference, the ¹H resonance of tetramethylsilane (TMS) in dilute solution (volume fraction, ? < 1%) in chloroform. Referencing procedures are discussed, including matters of practical application of the unified scale. Special attention is paid to recommended reference samples, and values of Ξ for secondary references on the unified scale are listed, many of which are the results of new measurements. Some earlier recommendations relating to the reporting of chemical shifts are endorsed. The chemical shift, δ, is redefined to avoid previous ambiguities but to leave practical usage unchanged. Relations between the unified scale and recently published recommendations for referencing in aqueous solutions (for specific use in biochemical work) are discussed, as well as the special effects of working in the solid state with magic‐angle spinning. In all, nine new recommendations relating to chemical shifts are made. Standardized nuclear spin data are also presented in tabular form for the stable (and some unstable) isotopes of all elements with non‐zero quantum numbers. The information given includes quantum numbers, isotopic abundances, magnetic moments, magnetogyric ratios and receptivities, together with quadrupole moments and linewidth factors (where appropriate). Copyright 2001 IUPAC. Reprinted with permission from Pure Appl. Chem. 2001; 73: 1795.     

Calculating nuclear magnetic resonance chemical shifts in solvated systems

The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.     

Standardization of chemical shifts of TMS and solvent signals in NMR solvents

The standard for chemical shift is dilute tetramethylsilane (TMS) in CDCl3, but many measurements are made relative to TMS in other solvents, the proton resonance of the solvent peak or relative to the lock frequency. Here, the chemical shifts of TMS and the proton and deuterium chemical shifts of the solvent signals of several solvents are measured over a wide temperature range. This allows for the use of TMS or the solvent and lock signal as a secondary reference for other NMR signals, as compared with dilute TMS in CDCl3 at a chosen temperature; 25 degrees C is chosen here. An accuracy of 0.02 ppm is achievable for dilute solutions, provided that the interaction with the solvent is not very strong. The proton chemical shift of residual water is also reported where appropriate.     

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.     

Solvent effects and dynamic averaging of 195Pt NMR shielding in cisplatin derivatives

The influences of solvent effects and dynamic averaging on the (195)Pt NMR shielding and chemical shifts of cisplatin and three cisplatin derivatives in aqueous solution were computed using explicit and implicit solvation models. Within the density functional theory framework, these simulations were carried out by combining ab initio molecular dynamics (aiMD) simulations for the phase space sampling with all-electron relativistic NMR shielding tensor calculations using the zeroth-order regular approximation. Structural analyses support the presence of a solvent-assisted "inverse" or "anionic" hydration previously observed in similar square-planar transition-metal complexes. Comparisons with computationally less demanding implicit solvent models show that error cancellation is ubiquitous when dealing with liquid-state NMR simulations. After aiMD averaging, the calculated chemical shifts for the four complexes are in good agreement with experiment, with relative deviations between theory and experiment of about 5% on average (1% of the Pt(II) chemical shift range).     

Characterization of a new hexasodium diphosphopentamolybdate hydrate, Na6[P2Mo5O23]x7H2O, by 23Na MQMAS NMR spectroscopy and X-ray powder diffraction

A novel hexasodium disphosphopentamolybdate hydrate, Na6[P2Mo5O23]x7H2O, has been identified using X-ray powder diffraction, 1H, 23Na, and 31P magic-angle spinning (MAS) NMR, and 23Na multiple-quantum (MQ) MAS NMR. Powder XRD reveals that the hydrate belongs to the triclinic spacegroup P1 with cell dimensions a = 10.090(3) A, b = 15.448(5) A, c = 8.460(4) A, alpha = 101.45(6) degrees, beta = 104.09(2) degrees, gamma = 90.71(5) degrees, and Z = 2. The number of water molecules of crystallization has been determined on the basis of a quantitative evaluation of the 1H MAS NMR spectrum, the crystallographic unit cell volume, and a hydrogen content analysis. The 23Na MQMAS NMR spectra of Na6[P2Mo5O23]x7H2O, obtained at three different magnetic fields, clearly resolve resonances from six different sodium sites and allow a determination of the second-order quadrupolar effect parameters and isotropic chemical shifts for the individual resonances. These data are used to determine the quadrupole coupling parameters (CQ and eta Q) from simulations of the complex line shapes of the central transitions, observed in 23Na MAS NMR spectra at the three magnetic fields. This analysis illustrates the advantages of combining MQMAS and MAS NMR at moderate and high magnetic fields for a precise determination of quadrupole coupling parameters and isotropic chemical shifts for multiple sodium sites in inorganic systems. 31P MAS NMR demonstrates the presence of two distinct P sites in the asymmetric unit of Na6[P2Mo5O23].7H2O while the 31P chemical shielding anisotropy parameters, determined for this hydrate and for Na6[P2Mo5O23]x13H2O, show that these two hydrates can easily be distinguished using 31P MAS NMR.     

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.     

2H Natural-abundance MAS NMR spectroscopy: an alternative approach to obtain 1H chemical shifts in solids

(2)H NMR was examined as an approach to determine (1)H chemical shifts in solids. For high-resolution observation, the line width due to (2)H quadrupole interaction and chemical-shift anisotropy was removed by magic-angle spinning and that due to (1)H-(2)H dipolar interactions by (1)H decoupling. Further, we showed that the sensitivity can be enhanced by applying (1)H to (2)H cross polarization and by adding spinning-sideband spectra. These make it possible to obtain (2)H natural-abundance MAS spectra revealing highly resolved (2)H signals. The second-order quadrupole effects of (2)H are also examined.     

MP2 calculation of 77Se NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the Second‐order Möller‐Plesset perturbation theory (MP2) calculation of ⁷⁷Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390–466 ppm) dominates over diamagnetic term (192–198 ppm) resulting in a total shielding relativistic correction of about 230–260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70–80 ppm. In the ⁷⁷Se NMR chemical shifts scale, relativistic corrections are about 20–30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of ⁷⁷Se NMR chemical shifts of the medium‐sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of ⁷⁷Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%. Copyright © 2015 John Wiley & Sons, Ltd.     

Unexpected effects of third-order cross-terms in heteronuclear spin systems under simultaneous radio-frequency irradiation and magic-angle spinning NMR

We recently noted [R. K. Harris, P. Hodgkinson, V. Zorin, J.-N. Dumez, B. Elena, L. Emsley, E. Salager, and R. Stein, Magn. Reson. Chem. 48, S103 (2010)] anomalous shifts in apparent (1)H chemical shifts in experiments using (1)H homonuclear decoupling sequences to acquire high-resolution (1)H NMR spectra for organic solids under magic-angle spinning (MAS). Analogous effects were also observed in numerical simulations of model (13)C,(1)H spin systems under homonuclear decoupling and involving large (13)C,(1)H dipolar couplings. While the heteronuclear coupling is generally assumed to be efficiently suppressed by sample spinning at the magic angle, we show that under conditions typically used in solid-state NMR, there is a significant third-order cross-term from this coupling under the conditions of simultaneous MAS and homonuclear decoupling for spins directly bonded to (1)H. This term, which is of the order of 100 Hz under typical conditions, explains the anomalous behaviour observed on both (1)H and (13)C spins, including the fast dephasing observed in (13)C{(1)H} heteronuclear spin-echo experiments under (1)H homonuclear decoupling. Strategies for minimising the impact of this effect are also discussed.     

Analyzing Discrepancies in Chemical-Shift Predictions of Solid Pyridinium Fumarates

Highly accurate chemical-shift predictions in molecular solids are behind the success and rapid development of NMR crystallography. However, unusually large errors of predicted hydrogen and carbon chemical shifts are sometimes reported. An understanding of these deviations is crucial for the reliability of NMR crystallography. Here, recently reported large deviations of predicted hydrogen and carbon chemical shifts of a series of solid pyridinium fumarates are thoroughly analyzed. The influence of the geometry optimization protocol and of the computational level of NMR calculations on the accuracy of predicted chemical shifts is investigated. Periodic calculations with GGA, meta-GGA and hybrid functionals are employed. Furthermore, molecular corrections at the coupled-cluster singles-and-doubles (CCSD) level are calculated. The effect of nuclear delocalization on the structure and NMR shielding is also investigated. The geometry optimization with a computationally demanding hybrid functional leads to a substantial improvement in proton chemical-shift predictions.     

Molecular insight into the electrostatic membrane surface potential by 14n/31p MAS NMR spectroscopy: nociceptin-lipid association

Exploiting naturally abundant (14)N and (31)P nuclei by high-resolution MAS NMR (magic angle spinning nuclear magnetic resonance) provides a molecular view of the electrostatic potential present at the surface of biological model membranes, the electrostatic charge distribution across the membrane interface, and changes that occur upon peptide association. The spectral resolution in (31)P and (14)N MAS NMR spectra is sufficient to probe directly the negatively charged phosphate and positively charged choline segment of the electrostatic P(-)-O-CH(2)-CH(2)-N(+)(CH(3))(3) headgroup dipole of zwitterionic DMPC (dimyristoylphosphatidylcholine) in mixed-lipid systems. The isotropic shifts report on the size of the potential existing at the phosphate and ammonium group within the lipid headgroup while the chemical shielding anisotropy ((31)P) and anisotropic quadrupolar interaction ((14)N) characterize changes in headgroup orientation in response to surface potential. The (31)P/(14)N isotropic chemical shifts for DMPC show opposing systematic changes in response to changing membrane potential, reflecting the size of the electrostatic potential at opposing ends of the P(-)-N(+) dipole. The orientational response of the DMPC lipid headgroup to electrostatic surface variations is visible in the anisotropic features of (14)N and (31)P NMR spectra. These features are analyzed in terms of a modified "molecular voltmeter" model, with changes in dynamic averaging reflecting the tilt of the C(beta)-N(+)(CH)(3) choline and PO(4)(-) segment. These properties have been exploited to characterize the changes in surface potential upon the binding of nociceptin to negatively charged membranes, a process assumed to proceed its agonistic binding to its opoid G-protein coupled receptor.     

NMR methods for studying the structure and dynamics of oncogenic and antihistaminic peptides in biomembranes

We present several applications of both wide-line and magic angle spinning (MAS) solid-state NMR of bicelles in which are embedded fragments of a tyrosine kinase receptor or enkephalins. The magnetically orientable bicelle membranes are shown to be of particular interest for studying the functional properties of lipids and proteins in a state that is very close to their natural environment. Quadrupolar, dipolar and chemical shielding interactions can be used to determine minute alterations of internal membrane dynamics and the orientation of peptides with respect to the membrane plane. MAS of bicelles can in turn lead to high-resolution proton spectra of hydrated membranes. Using deuterium-proton contrast methods one can then obtain pseudo-high-resolution proton spectra of peptides or proteins embedded in deuterated membranes and determine their atomic 3D structure using quasi-conventional liquid-state NMR methods.     

17O and 29Si NMR parameters of MgSiO3 phases from high-resolution solid-state NMR spectroscopy and first-principles calculations

The 29Si and 17O NMR parameters of six polymorphs of MgSiO3 were determined through a combination of high-resolution solid-state NMR and first-principles gauge including projector augmented wave (GIPAW) formalism calculations using periodic boundary conditions. MgSiO3 is an important component of the Earth's mantle that undergoes structural changes as a function of pressure and temperature. For the lower pressure polymorphs (ortho-, clino-, and protoenstatite), all oxygen species in the 17O high-resolution triple-quantum magic angle spinning (MAS) NMR spectra were resolved and assigned. These assignments differ from those tentatively suggested in previous work on the basis of empirical experimental correlations. The higher pressure polymorphs of MgSiO3 (majorite, akimotoite, and perovskite) are stabilized at pressures corresponding to the Earth's transition zone and lower mantle, with perovskite being the major constituent at depths >660 km. We present the first 17O NMR data for these materials and confirm previous 29Si work in the literature. The use of high-resolution multiple-quantum MAS (MQMAS) and satellite-transition MAS (STMAS) experiments allows us to resolve distinct oxygen species, and full assignments are suggested. The six polymorphs exhibit a wide variety of structure types, providing an ideal opportunity to consider the variation of NMR parameters (both shielding and quadrupolar) with local structure, including changes in coordination number, local geometry (bond distances and angles), and bonding. For example, we find that, although there is a general correlation of increasing 17O chemical shift with increasing Si-O bond length, the shift observed also depends upon the exact coordination environment.     

Effects of different GIAO and CSGT models and basis sets on 2-aryl-1,3,4-oxadiazole derivatives

The direct molecular structure implementations of the gage-including atomic orbital (GIAO), individual gages for atoms in molecules (IGAIM) and continuous set of gage transformations (CSGT) methods for calculating nuclear magnetic shielding tensors at both the Hartree-Fock (HF) and density functional (B3LYP) levels of theory with 6-31G(d), 6-311G(d), 6-31++G(d,p), 6-311++G(d,p), and 6-311++G(df,pd) basis sets are presented. Dependence on the ¹H and ¹³C NMR chemical shifts on the choice of method and basis set have been investigated. Also, these chemical shifts of 2-aryl-1,3,4-oxadiazoles 5a–g have been performed related to dihedral angles (C4–C3–C2–O) of two conformers. The optimized molecular geometries and ¹H and ¹³C chemical shift values of 2-aryl-1,3,4-oxadiazoles 5a–g in the ground state have been obtained. The linear correlation coefficients of ¹³C NMR chemical shifts for these molecules were given. The new nuclear magnetic shielding tensors of tetramethylsilane (TMS) were calculated. The data of 2-aryl-1,3,4-oxadiazole derivatives display significant molecular structure and NMR analysis. Also, these provide the basis for future design of efficient materials having the 1,3,4-oxadiazole core.     

Prediction of the 1H and 13C NMR spectra of alpha-D-glucose in water by DFT methods and MD simulations

We have applied computational protocols based on DFT and molecular dynamics simulations to the prediction of the alkyl 1H and 13C chemical shifts of alpha-d-glucose in water. Computed data have been compared with accurate experimental chemical shifts obtained in our laboratory. 13C chemical shifts do not show a marked solvent effect. In contrast, the results for 1H chemical shifts provided by structures optimized in the gas phase are only fair and point out that it is necessary to take into account both the flexibility of the glucose structure and the strong effect exerted by solvent water thereupon. Thus, molecular dynamics simulations were carried out to model both the internal geometry as well as the influence of solvent molecules on the conformational distribution of the solute. Snapshots from the simulation were used as input to DFT NMR calculations with varying degrees of sophistication. The most important factor that affects the accuracy of computed 1H chemical shifts is the solute geometry; the effect of the solvent on the shielding constants can be reasonably accounted for by self-consistent reaction field models without the need of explicitly including solvent molecules in the NMR property calculation.     

Separate-local-field NMR spectroscopy on half-integer quadrupolar nuclei

New approaches to the characterization of resonances in the solid-state NMR spectroscopy of half-integer quadrupolar nuclei are explored, on the basis of the acquisition of heteronuclear separate-local-field spectra on rotating solids. In their two-dimensional version, these experiments correlate for each chemical site a second-order quadrupolar MAS powder pattern with the dipolar MAS sideband pattern to nearby heteronuclei. As 3D NMR sequences, such 2D anisotropic correlation spectra become separated for inequivalent chemical sites along a third, isotropic dimension. Extending in such manner separate-local-field NMR approaches to quadrupoles facilitates the assignment of inequivalent resonances to specific structural environments, and provides new tools for the investigation of dynamics in solids. Details about these 2D and 3D NMR experiments are given, and their application is illustrated with 1H-23Na recoupling experiments on mononucleotides possessing multiple bound cations.     

A multinuclear solid-state NMR study of alkali metal ions in tetraphenylborate salts, M[BPh4] (M = Na, K, Rb and Cs): what is the NMR signature of cation-pi interactions?

We report a multinuclear solid-state ( (23)Na, (39)K, (87)Rb, (133)Cs) NMR study of tetraphenylborate salts, M[BPh 4] (M = Na, K, Rb, Cs). These compounds are isostructural in the solid state with the alkali metal ion surrounded by four phenyl groups resulting in strong cation-pi interactions. From analyses of solid-state NMR spectra obtained under stationary and magic-angle spinning (MAS) conditions at 11.75 and 21.15 T, we have obtained the quadrupole coupling constants, C Q, and the chemical shift tensor parameters for the alkali metal ions in these compounds. We found that the observed quadrupole coupling constant for M (+) in M[BPh 4] is determined by a combination of nuclear quadrupole moment, Sternheimer antishielding factor, and unit cell dimensions. On the basis of a comparison between computed paramagnetic and diamagnetic contributions to the total chemical shielding values for commonly found cation-ligand interactions, we conclude that cation-pi interactions give rise to significantly lower paramagnetic shielding contributions than other cation-ligand interactions. As a result, highly negative chemical shifts are expected to be the NMR signature for cations interacting exclusively with pi systems.     

High-resolution NMR correlation spectra of disordered solids

We show how high-resolution NMR spectra can be obtained for solids for which the spectra are normally broadened due to structural disorder. The method relies on correlations in the chemical shifts between pairs of coupled spins. It is found experimentally that there are strong correlations in the chemical shifts between neighboring spins in both phosphorus-31 and carbon-13 spectra. These correlations can be exploited not only to provide resolution in two-dimensional spectra, but also to yield "chains" of correlated chemical shifts, constituting a valuable new source of structural information for disordered materials.