Proximity and the heteroatom effects on the carbon-13 chemical shifts of methyl-substituted phenols,anilines and thiophenols

Upfield substituent-induced ¹³C chemical shifts for aryl carbons of polymethyl substituted benzenes, phenols, anilines and thiophenols were investigated as a function of the proximity between substituents X and CH₃ (X = CH₃, NH₂, OH and SH). The results indicate that the induced shifts of the substituted aryl carbons are, in general, independent of the polar substituent but depend on the number of adjacent substituted aryl carbons. A ?2.0 ppm upfield shift was found for a substituted aryl carbon adjacent to one substituted aryl carbon and a ?3.8 ppm upfield shift for a substituted aryl carbon bound by two substituted aryl carbons. It is suggested that the near additivity of the upfield shifts is the result of changes in the bond order between the aromatic ring carbons in the region of the substituted aryl carbons due to distortion of the ring. The ¹³C chemical shifts of the methyl substituents for methyl substituted phenols, anilines and thiophenols were determined, and it was found that the values could be predicted from the additivity parameters reported for the analogous methylbenzenes plus an additional pair-interaction term associated with the through-space electronic influence of the heteroatom.     

Effects of amino acid phi,psi propensities and secondary structure interactions in modulating H alpha chemical shifts in peptide and protein beta-sheet.

H alpha chemical shifts are often used as indicators of secondary structure formation in protein structural analysis and peptide folding studies. On the basis of NMR analysis of model beta-sheet and alpha-helical peptides, together with a statistical analysis of protein structures for which NMR data are available, we show that although the gross pattern of H alpha chemical shifts reflects backbone torsion angles, longer range effects from distant amino acids are the dominant factor determining experimental chemical shifts in beta-sheets of peptides and proteins. These show context-dependent variations that aid structural assignment and highlight anomalous shifts that may be of structural significance and provide insights into beta-sheet stability.     

Carbon-13 NMR shielding in the twenty common amino acids: comparisons with experimental results in proteins

We have used ab initio quantum chemical techniques to compute the (13)C(alpha) and (13)C(beta) shielding surfaces for the 14 amino acids not previously investigated (R. H. Havlin et al., J. Am. Chem. Soc. 1997, 119, 11951-11958) in their most popular conformations. The spans (Omega = sigma(33) - sigma(11)) of all the tensors reported here are large ( approximately 34 ppm) and there are only very minor differences between helical and sheet residues. This is in contrast to the previous report in which Val, Ile and Thr were reported to have large ( approximately 12 ppm) differences in Omega between helical and sheet geometries. Apparently, only the beta-branched (beta-disubstituted) amino acids have such large CSA span (Omega) differences; however, there are uniformly large differences in the solution-NMR-determined CSA (Deltasigma = sigma(orth) - sigma(par)) between helices and sheets in all amino acids considered. This effect is overwhelmingly due to a change in shielding tensor orientation. With the aid of such shielding tensor orientation information, we computed Deltasigma values for all of the amino acids in calmodulin/M13 and ubiquitin. For ubiquitin, we find only a 2.7 ppm rmsd between theory and experiment for Deltasigma over an approximately 45 ppm range, a 0.96 slope, and an R(2) = 0.94 value when using an average solution NMR structure. We also report C(beta) shielding tensor results for these same amino acids, which reflect the small isotropic chemical shift differences seen experimentally, together with similar C(beta) shielding tensor magnitudes and orientations. In addition, we describe the results of calculations of C(alpha), C(beta), C(gamma)1, C(gamma)2, and C(delta) shifts in the two isoleucine residues in bovine pancreatic trypsin inhibitor and the four isoleucines in a cytochrome c and demonstrate that the side chain chemical shifts are strongly influenced by chi(2) torsion angle effects. There is very good agreement between theory and experiment using either X-ray or average solution NMR structures. Overall, these results show that both C(alpha) backbone chemical shift anisotropy results as well as backbone and side chain (13)C isotropic shifts can now be predicted with good accuracy by using quantum chemical methods, which should facilitate solution structure determination/refinement using such shielding tensor surface information.     

Solid-state NMR and quantum chemical investigations of 13Calpha shielding tensor magnitudes and orientations in peptides: determining phi and psi torsion angles

We report the experimental determination of the (13)C(alpha) chemical shift tensors of Ala, Leu, Val, Phe, and Met in a number of polycrystalline peptides with known X-ray or de novo solid-state NMR structures. The 700 Hz dipolar coupling between (13)C(alpha) and its directly bonded (14)N permits extraction of both the magnitude and the orientation of the shielding tensor with respect to the C(alpha)-N bond vector. The chemical shift anisotropy (CSA) is recoupled under magic-angle spinning using the SUPER technique (Liu et al., J. Magn. Reson. 2002, 155, 15-28) to yield quasi-static chemical shift powder patterns. The tensor orientation is extracted from the (13)C-(14)N dipolar modulation of the powder line shapes. The magnitudes and orientations of the experimental (13)C(alpha) chemical shift tensors are found to be in good accord with those predicted from quantum chemical calculations. Using these principal values and orientations, supplemented with previously measured tensor orientations from (13)C-(15)N and (13)C-(1)H dipolar experiments, we are able to predict the (phi, psi, chi(1)) angles of Ala and Val within 5.8 degrees of the crystallographic values. This opens up a route to accurate determination of torsion angles in proteins based on shielding tensor magnitude and orientation information using labeled compounds, as well as the structure elucidation of noncrystalline organic compounds using natural abundance (13)C NMR techniques.     

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.     

Toward direct determination of conformations of protein building units from multidimensional NMR experiments VI: chemical shift analysis of his to gain 3D structure and protonation state information

NMR--chemical shift structure correlations were investigated by using GIAO-RB3LYP/6-311++G(2d,2p) formalism. Geometries and chemical shifts (CSI values) of 103 different conformers of N'-formyl-L-histidinamide were determined including both neutral and charged protonation forms. Correlations between amino acid torsional angle values and chemical shifts were investigated for the first time for an aromatic and polar amino acid residue whose side chain may carry different charges. Linear correlation coefficients of a significant level were determined between chemical shifts and dihedral angles for CSI[(1)H(alpha)]/phi, CSI[(13)C(alpha)]/phi, and CSI[(13)C(alpha)]/psi. Protonation of the imidazole ring induces the upfield shift of CSI[(13)C(alpha)] for positively charged histidines and an opposite effect for the negative residue. We investigated the correspondence of theoretical and experimental (13)C(alpha), (13)C(beta), and (1)H(alpha) chemical shifts and the nine basic conformational building units characteristic for proteins. These three chemical shift values allow the identification of conformational building units at 80% accuracy. These results enable the prediction of additional regular secondary structural elements (e.g., polyProlineII, inverse gamma-turns) and loops beyond the assignment of chemical shifts to alpha-helices and beta-pleated sheets. Moreover, the location of the His residue can be further specified in a beta-sheet. It is possible to determine whether the appropriate residue is located at the middle or in a first/last beta-strand within a beta-sheet based on calculated CSI values. Thus, the attractive idea of establishing local residue specific backbone folding parameters in peptides and proteins by employing chemical shift information (e.g., (1)H(alpha) and (13)C(alpha)) obtained from selected heteronuclear correlation NMR experiments (e.g., 2D-HSQC) is reinforced.     

Relationships between 31P chemical shift tensors and conformation of nucleic acid backbone: a DFT study

Density functional theory (DFT) has been applied to study the conformational dependence of 31P chemical shift tensors in B-DNA. The gg and gt conformations of backbone phosphate groups representing BI- and BII-DNA have been examined. Calculations have been carried out on static models of dimethyl phosphate (dmp) and dinucleoside-3',5'-monophosphate with bases replaced by hydrogen atoms in vacuo as well as in an explicit solvent. Trends in 31P chemical shift anisotropy (CSA) tensors with respect to the backbone torsion angles alpha, zeta, beta, and epsilon are presented. Although these trends do not change qualitatively upon solvation, quantitative changes result in the reduction of the chemical shift anisotropy. For alpha and zeta in the range from 270 degrees to 330 degrees and from 240 degrees to 300 degrees , respectively, the delta22 and delta33 principal components vary within as much as 30 ppm, showing a marked dependence on backbone conformation. The calculated 31P chemical shift tensor principal axes deviate from the axes of O-P-O bond angles by at most 5 degrees . For solvent models, our results are in a good agreement with experimental estimates of relative gg and gt isotropic chemical shifts. Solvation also brings the theoretical deltaiso of the gg conformation closer to the experimental gg data of barium diethyl phosphate.     

Investigation of substituent effects of chalcones by 13C n.m.r. spectroscopy

¹³C chemical shifts for 23 para- and meta-substituted chalcones of the types 1 and 2 have been determined. The aromatic shieldings are compared with previous results for other aromatic derivatives. Correlations of the ¹³C chemical shifts of vinyl carbons and carbonyl carbons as well as ring carbons with Hammett σ parameters, π electron densities and the reactivity parameters of Swain and Lupton provide a consistent picture of electronic effects transmitted through the carbon framework of the compounds studied.     

β-Sheet 13C structuring shifts appear only at the H-bonded sites of hairpins

The (13)C chemical shifts measured for designed β-hairpins indicate that the structuring shifts (upfield for Cα and C', downfield for Cβ) previously reported as diagnostic for β-structuring in proteins appear only at the H-bonded strand residues. The resulting periodicity of structuring shift magnitudes is not, however, a consequence of H-bonding status; rather, it reflects a previously unrecognized alternation in the backbone torsion angles of β-strands. This feature of hairpins is also likely to be present in proteins. The study provides reference values for the expectation shifts for (13)C sites in β-structures that should prove useful in the characterization of the folding equilibria of β-sheet models.     

Conformation and configuration of tertiary amines via GIAO-derived (13)C NMR chemical shifts and a multiple independent variable regression analysis.

The (13)C chemical shifts of six tertiary amines of unambiguous conformational structure are compared to predicted (13)C NMR chemical shifts obtained via empirically scaled GIAO shieldings for geometries from MM3 molecular mechanics calculations. An average deviation, absolute value of Deltadelta(av), of 0.8 ppm and a maximum deviation, absolute value of Deltadelta(max), of 2.8 ppm between predicted and experimental (13)C shifts of the six tertiary amines of unambiguous structure are found. In several cases of tertiary amines subject to rapid exchange, where experimental (13)C shifts at room temperature are weighted averages of multiple conformers, a comparison of calculated (13)C shifts of all reasonable MM3 predicted conformers with experimental (13)C shifts via a multiple independent variable regression analysis provides an efficient method of determining the major and minor conformers. The examples presented are 2-methyl-2-azabicyclo[2.2.1]heptane and 1,6-diazabicyclo[4.3.1]decane, which each have two expected contributing structures, and 2-(diethylamino)propane and 1,8-diazabicyclo[6.3.1]dodecane, where ten and seven low-energy conformers, respectively, are predicted by MM3 calculations.     

Experimental and theoretical investigation of the 13C and 15N chemical shift tensors in melanostatin-exploring the chemical shift tensor as a structural probe

The determination of backbone conformations in powdered peptides using 13C and 15N shift tensor information is explored. The 13C and 15N principal shift values in natural abundance 13C and 15N melanostatin (L-Pro-L-Leu-Gly amide) are measured using the FIREMAT technique. Furthermore, the orientation of the C-N bond in the 13C shift principal axis system for the backbone carbons is obtained from the presence of the 13C-14N dipolar coupling. The Ramachandran angles for the title compound are obtained from solid-state NMR data by comparing the experimentally determined shift tensor information to systematic theoretical shielding calculations on N-formyl-L-amino acid-amide models. The effects of geometry optimization and neglect of intermolecular interactions on the theoretical shielding values in the model compounds are investigated. The sets of NMR derived Ramachandran angles are assembled in a set of test structures that are compared to the available single-crystal X-ray structure. Shift tensor calculations on the test structures and the X-ray structure are used to further assess the importance of intermolecular interactions when the shift tensor is used as a structural probe in powdered peptides.     

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.     

A solid state 13C NMR, crystallographic, and quantum chemical investigation of chemical shifts and hydrogen bonding in histidine dipeptides

We report the first solid-state NMR, crystallographic, and quantum chemical investigation of the origins of the 13C NMR chemical shifts of the imidazole group in histidine-containing dipeptides. The chemical shift ranges for Cgamma and Cdelta2 seen in eight crystalline dipeptides were very large (12.7-13.8 ppm); the shifts were highly correlated (R2= 0.90) and were dominated by ring tautomer effects and intermolecular interactions. A similar correlation was found in proteins, but only for buried residues. The imidazole 13C NMR chemical shifts were predicted with an overall rms error of 1.6-1.9 ppm over a 26 ppm range, by using quantum chemical methods. Incorporation of hydrogen bond partner molecules was found to be essential in order to reproduce the chemical shifts seen experimentally. Using AIM (atoms in molecules) theory we found that essentially all interactions were of a closed shell nature and the hydrogen bond critical point properties were highly correlated with the N...H...O (average R2= 0.93) and Nepsilon2...H...N (average R2= 0.98) hydrogen bond lengths. For Cepsilon1, the 13C chemical shifts were also highly correlated with each of these properties (at the Nepsilon2 site), indicating the dominance of intermolecular interactions for Cepsilon1. These results open up the way to analyzing 13C NMR chemical shifts, tautomer states (from Cdelta2, Cepsilon1 shifts), and hydrogen bond properties (from Cepsilon1 shifts) of histidine residue in proteins and should be applicable to imidazole-containing drug molecules bound to proteins, as well.     

Chi(1) rotamer populations and angles of mobile surface side chains are accurately predicted by a torsion angle database potential of mean force

The equilibrium angles and distributions of chi(1) rotamers for mobile surface side chains of the small, 63-residue, B1 domain of protein L have been calculated from the static crystal structure by rigid body/torsion angle simulated annealing using a torsion angle database potential of mean force and compared to those deduced by Monte Carlo analysis of side chain residual dipolar couplings measured in solution. Good agreement between theory and experiment is observed, indicating that for side chains undergoing rotamer averaging that is fast on the chemical shift time scale, the equilibrium angles and distribution of chi(1) rotamers are largely determined by the backbone phi/psi torsion angles.     

Performance of revised STO(1M)-3G basis set for prediction of 5-fluorocytosine chemical shifts

Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied. The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted ¹³C and ¹H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for ¹³C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.     

Determination of peptide backbone torsion angles using double-quantum dipolar recoupling solid-state NMR spectroscopy

Several approaches for utilizing dipolar recoupling solid-state NMR (ssNMR) techniques to determine local structure at high resolution in peptides and proteins have been developed. However, many of these techniques measure only one torsion angle or are accurate for only certain classes of secondary structure. Additionally, the efficiency with which these dipolar recoupling experiments suppress the deleterious effects of chemical shift anisotropy (CSA) at high magnetic field strengths varies. Dipolar recoupling with a windowless sequence (DRAWS) has proven to be an effective pulse sequence for exciting double-quantum (DQ) coherences between adjacent carbonyl carbons along the peptide backbone. By allowing this DQ coherence to evolve, it is possible to measure the relative orientations of the CSA tensors and subsequently use this information to determine the Ramachandran torsion angles phi and psi. Here, we explore the accuracies of the assumptions made in interpreting DQ-DRAWS data and demonstrate their fidelity in measuring torsion angles corresponding to a variety of secondary structures irrespective of hydrogen-bonding patterns. It is shown how a simple choice of isotopic labels and experimental conditions allows accurate measurement of backbone secondary structures without any prior knowledge. This approach is considerably more sensitive for determining structure in helices and has comparable accuracy for beta-sheet and extended conformations relative to other methods. We also illustrate the ability of DQ-DRAWS to distinguish between structures in heterogeneous samples.     

(67)Zn NMR chemical shifts and electric field gradients in zinc complexes: a quantum chemical investigation

We have used quantum chemical methods to predict 67Zn NMR chemical shifts as well as quadrupole coupling constants (CQ) in a series of biomimetic and inorganic zinc complexes. The 67Zn chemical shifts are predicted with an R2 = 0.975, corresponding to a 24.3 ppm or 6.7% error over the entire 365 ppm 67Zn chemical shift range. The 67Zn CQ values are predicted with an R2 = 0.991, corresponding to a 1.17 MHz or 3.0% error over the entire 38.75 MHz range. The 67Zn NMR shifts in a series of complexes containing N,O ligands are, in general, highly correlated with the number of oxygen ligands. The ability to compute 67Zn NMR shifts as well as CQ values opens up the possibility of using both of these properties in structure determination or refinement in proteins.     

Determination of multiple torsion-angle constraints in U-(13)C,(15)N-labeled peptides: 3D (1)H-(15)N-(13)C-(1)H dipolar chemical shift NMR spectroscopy in rotating solids

We demonstrate constraint of peptide backbone and side-chain conformation with 3D (1)H-(15)N-(13)C-(1)H dipolar chemical shift, magic-angle spinning NMR experiments. In these experiments, polarization is transferred from (15)N[i] by ramped SPECIFIC cross polarization to the (13)C(alpha)[i], (13)C(beta)[i], and (13)C(alpha)[i - 1] resonances and evolves coherently under the correlated (1)H-(15)N and (1)H-(13)C dipolar couplings. The resulting set of frequency-labeled (15)N(1)H-(13)C(1)H dipolar spectra depend strongly upon the molecular torsion angles phi[i], chi1[i], and psi[i - 1]. To interpret the data with high precision, we considered the effects of weakly coupled protons and differential relaxation of proton coherences via an average Liouvillian theory formalism for multispin clusters and employed average Hamiltonian theory to describe the transfer of (15)N polarization to three coupled (13)C spins ((13)C(alpha)[i], (13)C(beta)[i], and (13)C(alpha)[i - 1]). Degeneracies in the conformational solution space were minimized by combining data from multiple (15)N(1)H-(13)C(1)H line shapes and analogous data from other 3D (1)H-(13)C(alpha)-(13)C(beta)-(1)H (chi1), (15)N-(13)C(alpha)-(13)C'-(15)N (psi), and (1)H-(15)N[i]-(15)N[i + 1]-(1)H (phi, psi) experiments. The method is demonstrated here with studies of the uniformly (13)C,(15)N-labeled solid tripeptide N-formyl-Met-Leu-Phe-OH, where the combined data constrains a total of eight torsion angles (three phi, three chi1, and two psi): phi(Met) = -146 degrees, psi(Met) = 159 degrees, chi1(Met) = -85 degrees, phi(Leu) = -90 degrees, psi(Leu) = -40 degrees, chi1(Leu) = -59 degrees, phi(Phe) = -166 degrees, and chi1(Phe) = 56 degrees. The high sensitivity and dynamic range of the 3D experiments and the data analysis methods provided here will permit immediate application to larger peptides and proteins when sufficient resolution is available in the (15)N-(13)C chemical shift correlation spectra.