Toward direct determination of conformations of protein building units from multidimensional NMR experiments part II: a theoretical case study of formyl-L-valine amide

Chemical shielding anisotropy tensors have been determined for all twenty-seven characteristic conformers of For-L-Val-NH2 using the GIAO-RHF formalism with the 6-31 + G* and TZ2P basis sets. The individual chemical shifts and their conformational averages have been compared to their experimental counterparts taken from the BioMagnetic Resonance Bank (BMRB). At the highest level of theory applied, for all nuclei but the amide proton, deviations between statistically averaged theoretical and experimental chemical shifts are as low as 1-3%. Correlated chemical shift plots of selected nuclei, as function of the respective phi, psi, chi1, and chi2 torsional angles, have been generated. On two-dimensional chemical shift-chemical shift plots, for example, 1H(NH)-15N(NH) and 15N(NH)-13Calpha, regions corresponding to major conformational clusters have been identified, providing a basis for the quantitative identification of conformers from NMR shift data. Experimental NMR resonances of nuclei of valine residues have been deduced from 18 selected proteins, resulting in 93 1Halpha-13Calpha chemical shift pairs. These experimental results have been compared to relevant ab initio values revealing remarkable correlation between the two sets of data. Correlations of 1Halpha and 13Calpha values with backbone conformational parameters (phi and psi) have also been found for all pairs (e.g. 1Halpha/phi and 13Calpha/phi) but 1Halpha/psi. Overall, the appealing idea of establishing backbone folding of proteins by employing chemical shift information alone, obtained from selected multiple-pulse NMR experiments (e.g. 2D-HSQC, 2D-HMQC, and 3D-HNCA), has received further support.     

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.     

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.     

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.     

Tryptophan chemical shift in peptides and proteins: a solid state carbon-13 nuclear magnetic resonance spectroscopic and quantum chemical investigation

We have obtained the carbon-13 nuclear magnetic resonance spectra of a series of tryptophan-containing peptides and model systems, together with their X-ray crystallographic structures, and used quantum chemical methods to predict the (13)C NMR shifts (or shieldings) of all nonprotonated aromatic carbons (C(gamma), C(delta 2) and C(epsilon 2). Overall, there is generally good accord between theory and experiment. The chemical shifts of Trp C(gamma) in several proteins, hen egg white lysozyme, horse myoglobin, horse heart cytochrome c, and four carbonmonoxyhemoglobins, are also well predicted. The overall Trp C(gamma) shift range seen in the peptides and proteins is 11.4 ppm, and individual shifts (or shieldings) are predicted with an rms error of approximately 1.4 ppm (R value = 0.86). Unlike C(alpha) and N(H) chemical shifts, which are primarily a function of the backbone phi,psi torsion angles, the Trp C(gamma) shifts are shown to be correlated with the side-chain torsion angles chi(1) and chi(2) and appear to arise, at least in part, from gamma-gauche interactions with the backbone C' and N(H) atoms. This work helps solve the problem of the chemical shift nonequivalences of nonprotonated aromatic carbons in proteins first identified over 30 years ago and opens up the possibility of using aromatic carbon chemical shift information in structure determination.     

A theoretical case study of type I and type II beta-turns

NMR chemical shielding anisotropy tensors have been computed by employing a medium size basis set and the GIAO-DFT(B3LYP) formalism of electronic structure theory for all of the atoms of type I and type II beta-turn models. The models contain all possible combinations of the amino acid residues Gly, Ala, Val, and Ser, with all possible side-chain orientations where applicable in a dipeptide. The several hundred structures investigated contain either constrained or optimized phi, psi, and chi dihedral angles. A statistical analysis of the resulting large database was performed and multidimensional (2D and 3D) chemical-shift/chemical-shift plots were generated. The (1)H(alpha-13)C(alpha), (13)C(alpha-1)H(alpha-13)C(beta), and (13)C(alpha-1)H(alpha-13)C' 2D and 3D plots have the notable feature that the conformers clearly cluster in distinct regions. This allows straightforward identification of the backbone and side-chain conformations of the residues forming beta-turns. Chemical shift calculations on larger For-(L-Ala)(n)-NH(2) (n=4, 6, 8) models, containing a single type I or type II beta-turn, prove that the simple models employed are adequate. A limited number of chemical shift calculations performed at the highly correlated CCSD(T) level prove the adequacy of the computational method chosen. For all nuclei, statistically averaged theoretical and experimental shifts taken from the BioMagnetic Resonance Bank (BMRB) exhibit good correlation. These results confirm and extend our previous findings that chemical shift information from selected multiple-pulse NMR experiments could be employed directly to extract folding information for polypeptides and proteins.     

Structure distribution in an elastin-mimetic peptide (VPGVG)3 investigated by solid-state NMR

Elastin is an extracellular-matrix protein that imparts elasticity to tissues. We have used solid-state NMR to determine a number of distances and torsion angles in an elastin-mimetic peptide, (VPGVG)3, to understand the structural basis of elasticity. C-H and C-N distances between the V6 carbonyl and the V9 amide segment were measured using 13C-15N and 13C-1H rotational-echo double-resonance experiments. The results indicate the coexistence of two types of intramolecular distances: a third of the molecules have short C-H and C-N distances of 3.3 +/- 0.2 and 4.3 +/- 0.2 A, respectively, while the rest have longer distances of about 7 A. Complementing the distance constraints, we measured the (phi, psi ) torsion angles of the central pentameric unit using dipolar correlation NMR. The -angles of P7 and G8 are predominantly ~150, thus restricting the majority of the peptide to be extended. Combining all torsion angles measured for the five residues, the G8 C chemical shift, and the V6-V9 distances, we obtained a bimodal structure distribution for the PG residues in VPGVG. The minor form is a compact structure with a V6-V9 C=O-HN hydrogen bond and can be either a type II -turn or a previously unidentified turn with Pro (phi = -70, psi= 20 +/- 20) and Gly ( phi= -100 +/- 20, psi = -20 +/- 20). The major form is an extended and distorted beta-strand without a V6-V9 hydrogen bond and differs from the ideal parallel and antiparallel beta-strands. The other three residues in the VPGVG unit mainly adopt antiparallel beta-sheet torsion angles. Since (VPGVG)3 has the same 13C and 15N isotropic and anisotropic chemical shifts as the elastin-mimetic protein (VPGXG)n (X = V and K, n = 195), the observed conformational distribution around Pro and Gly sheds light on the molecular mechanism of elastin elasticity.     

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.     

Variation in protein C(alpha)-related one-bond J couplings

Four types of polypeptide (1)J(C alpha X) couplings are examined, involving the main-chain carbon C(alpha) and either of four possible substituents. A total 3105 values of (1)J(C alpha H alpha), (1)J(C alpha C beta), (1)J(C alpha C'), and (1)J(C alpha N') were collected from six proteins, averaging 143.4 +/- 3.3, 34.9 +/- 2.5, 52.6 +/- 0.9, and 10.7 +/- 1.2 Hz, respectively. Analysis of variances (ANOVA) reveals a variety of factors impacting on (1)J and ranks their relative statistical significance and importance to biomolecular NMR structure refinement. Accordingly, the spread in the (1)J values is attributed, in equal proportions, to amino-acid specific substituent patterns and to polypeptide-chain geometry, specifically torsions phi, psi, and chi(1) circumjacent to C(alpha). The (1)J coupling constants correlate with protein secondary structure. For alpha-helical phi, psi combinations, (1)J(C alpha H alpha) is elevated by more than one standard deviation (147.8 Hz), while both (1)J(C alpha N') and (1)J(C alpha C beta) fall short of their grand means (9.5 and 33.7 Hz). Rare positive phi torsion angles in proteins exhibit concomitant small (1)J(C alpha H alpha) and (1)J(C alpha N') (138.4 and 9.6 Hz) and large (1)J(C alpha C beta) (39.9 Hz) values. The (1)J(C alpha N') coupling varies monotonously over the phi torsion range typical of beta-sheet secondary structure and is largest (13.3 Hz) for phi around -160 degrees. All four coupling types depend on psi and thus help determine a torsion that is notoriously difficult to assess by traditional approaches using (3)J. Influences on (1)J stemming from protein secondary structure and other factors, such as amino-acid composition, are largely independent.     

Packing interactions in hydrated and anhydrous forms of the antibiotic Ciprofloxacin: a solid-state NMR, X-ray diffraction, and computer simulation study

We present an experimental NMR, X-ray diffraction (XRD), and computational study of the supramolecular assemblies of two crystalline forms of Ciprofloxacin: one anhydrate and one hydrate forming water wormholes. The resonance assignment of up to 51 and 54 distinct (13)C and (1)H resonances for the hydrate is reported. The effect of crystal packing, identified by XRD, on the (1)H and (13)C chemical shifts including weak interionic H-bonds, is quantified; (1)H chemical shift changes up to ～-3.5 ppm for CH···π contacts and ～+2 ppm (CH···O((-))); ～+4.7 ppm (((+))NH···O((-))) for H-bonds. Water intake induces chemical shift changes up to 2 and 5 ppm for (1)H and (13)C nuclei, respectively. Such chemical shifts are found to be sensitive detectors of hydration/dehydration in highly insoluble hydrates.     

Solution NMR studies of the A beta(1-40) and A beta(1-42) peptides establish that the Met35 oxidation state affects the mechanism of amyloid formation

The pathogenesis of Alzheimer's disease is characterized by the aggregation and fibrillation of the 40-residue A beta(1-40) and 42-residue A beta(1-42) peptides into amyloid plaques. The structural changes associated with the conversion of monomeric A beta peptide building blocks into multimeric fibrillar beta-strand aggregates remain unknown. Recently, we established that oxidation of the methionine-35 side chain to the sulfoxide (Met35(red) --> Met35(ox)) significantly impedes the rate of aggregation and fibrillation of the A beta peptide. To explore this effect at greater resolution, we carefully compared the (1)H, (15)N, and (13)C NMR chemical shifts of four A beta peptides that had the Met35 reduced or oxidized (A beta(1-40)Met35(red), A beta(1-40)Met35(ox), A beta(1-42)Met35(red), and A beta(1-42)Met35(ox)). With the use of a special disaggregation protocol, the highly aggregation prone A beta peptides could be studied at higher, millimolar concentrations (as required by NMR) in aqueous solution at neutral pH, remaining largely monomeric at 5 degrees C as determined by sedimentation equilibrium studies. The NOE, amide-NH temperature coefficients, and chemical shift indices of the (1)H alpha, (13)C alpha, and (13)C beta established that the four peptides are largely random, extended chain structures, with the Met35(ox) reducing the propensity for beta-strand structure at two hydrophobic regions (Leu17-Ala21 and Ile31-Val36), and turn- or bendlike structures at Asp7-Glu11 and Phe20-Ser26. Additional NMR studies monitoring changes that occur during aging at 37 degrees C established that, along with a gradual loss of signal/noise, the Met35(ox) significantly hindered upfield chemical shift movements of the 2H NMR signals for the His6, His13, and His14 side chains. Taken together, the present NMR studies demonstrate that the Met35(red) --> Met35(ox) conversion prevents aggregation by reducing both hydrophobic and electrostatic association and that the A beta(1-40)Met35(red), A beta(1-40)Met35(ox), A beta(1-42)Met35(red), and A beta(1-42)Met35(ox) peptides may associate differently, through specific, sharp changes in structure during the initial stages of aggregation.     

Investigation of the neighboring residue effects on protein chemical shifts

In this study, we report nearest neighbor residue effects statistically determined from a chemical shift database. For an amino acid sequence XYZ, we define two correction factors, Delta((X)Y)n,s and Delta(Y(Z))n,s, representing the effects on Y's chemical shifts from the preceding residue (X) and the following residue (Z), respectively, where X, Y, and Z are any of the 20 naturally occurring amino acids, n stands for (1)H(N), (15)N, (1)H(alpha), (13)C(alpha), (13)C(beta), and (13)C' nuclei, and s represents the three secondary structural types beta-strand, random coil, and alpha-helix. A total of approximately 14400 Delta((X)Y)n,s and Delta(Y(Z))n,s, representing nearly all combinations of X, Y, Z, n, and s, have been quantitatively determined. Our approach overcomes the limits of earlier experimental methods using short model peptides, and the resulting correction factors have important applications such as chemical shift prediction for the folded proteins. More importantly, we have found, for the first time, a linear correlation between the Delta((X)Y)n,s (n = (15)N) and the (13)C(alpha) chemical shifts of the preceding residue X. Since (13)C(alpha) chemical shifts of the 20 amino acids, which span a wide range of 40-70 ppm, are largely dominated by one property, the electron density of the side chain, the correlation indicates that the same property is responsible for the effect on the following residue. The influence of the secondary structure on both the chemical shifts and the nearest neighbor residue effect are also investigated.     

Structure determination of a peptide model of the repeated helical domain in Samia cynthia ricini silk fibroin before spinning by a combination of advanced solid-state NMR methods

Fibrous proteins unlike globular proteins, contain repetitive amino acid sequences, giving rise to very regular secondary protein structures. Silk fibroin from a wild silkworm, Samia cynthia ricini, consists of about 100 repeats of alternating polyalanine (poly-Ala) regions of 12-13 residues in length and Gly-rich regions. In this paper, the precise structure of the model peptide, GGAGGGYGGDGG(A)(12)GGAGDGYGAG, which is a typical repeated sequence of the silk fibroin, was determined using a combination of three kinds of solid-state NMR studies; a quantitative use of (13)C CP/MAS NMR chemical shift with conformation-dependent (13)C chemical shift contour plots, 2D spin diffusion (13)C solid-state NMR under off magic angle spinning and rotational echo double resonance. The structure of the model peptide corresponding to the silk fibroin structure before spinning was determined. The torsion angles of the central Ala residue, Ala(19), in the poly-Ala region were determined to be (phi, psi) = (-59 degrees, -48 degrees ) which are values typically associated with alpha-helical structures. However, the torsion angles of the Gly(25) residue adjacent to the C-terminal side of the poly-Ala chain were determined to be (phi, psi) = (-66 degrees, -22 degrees ) and those of Gly(12) and Ala(13) residues at the N-terminal of the poly-Ala chain to be (phi, psi) = (-70 degrees, -30 degrees ). In addition, REDOR experiments indicate that the torsion angles of the two C-terminal Ala residues, Ala(23) and Ala(24), are (phi, psi) = (-66 degrees, -22 degrees ) and those of N-terminal two Ala residues, Ala(13) and Ala(14) are (phi, psi) = (-70 degrees, -30 degrees ). Thus, the local structure of N-terminal and C-terminal residues, and also the neighboring residues of alpha-helical poly-Ala chain in the model peptide is a more strongly wound structure than found in typical alpha-helix structures.     

Spin dynamics in the modulation frame: application to homonuclear recoupling in magic angle spinning solid-state NMR

We introduce a family of solid-state NMR pulse sequences that generalizes the concept of second averaging in the modulation frame and therefore provides a new approach to perform magic angle spinning dipolar recoupling experiments. Here, we focus on two particular recoupling mechanisms-cosine modulated rotary resonance (CMpRR) and cosine modulated recoupling with isotropic chemical shift reintroduction (COMICS). The first technique, CMpRR, is based on a cosine modulation of the rf phase and yields broadband double-quantum (DQ) (13)C recoupling using >70 kHz omega(1,C)/2pi rf field for the spinning frequency omega(r)/2=10-30 kHz and (1)H Larmor frequency omega(0,H)/2pi up to 900 MHz. Importantly, for p>or=5, CMpRR recouples efficiently in the absence of (1)H decoupling. Extension to lower p values (3.5相似文献   

Complete 1H and 13C NMR assignments of three new polyhydroxylated sterols from the South China Sea gorgonian Subergorgia suberosa

Three new polyhydroxylated sterols, 3beta,6alpha,11,20beta,24-pentahydroxy- 9,11-seco-5alpha-24-ethylcholest-7,28-diene-9-one (1), 3-(1',2'-ethandiol)-24- methylcholest-8(9),22E-diene-3beta,5alpha,6alpha,7alpha,11alpha-pentaol (2), 24-methylcholest-7,22 E-diene-3beta,5alpha,6beta,25-tetraol (3) together with five known sterols, were isolated from the EtOH/CH2Cl2 extract of the South China Sea gorgonian Subergorgia suberosa. The complete assignments of the 1H and 13C NMR chemical shifts for these new compounds were achieved by means of 1D and 2D NMR techniques, including HSQC, HMBC, 1H--1H COSY, and NOESY spectra.     

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 (+).     

Cucurbitacins from Cayaponia racemosa: isolation and total assignment of 1H and 13C NMR spectra

Two new cucurbitane-type triterpenoids, 2beta,3beta,16alpha,20(R),25-pentahydroxy-9-methyl-19-norlanost-5-en-7,22-dione and 2beta,3beta,16alpha,20(R),25-pentahydroxy-9-methyl-19-norlanost-5-en-7,11,22-trione, were isolated from fruits of Cayaponia racemosa. The total (1)H and (13)C chemical shift assignment of these two closely related compounds is described, making use of one- and two-dimensional NMR techniques.     

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.     

Ramachandran-type plots for glycosidic linkages: Examples from molecular dynamic simulations using the Glycam06 force field

The goals of this article are to (1) provide further validation of the Glycam06 force field, specifically for its use in implicit solvent molecular dynamic (MD) simulations, and (2) to present the extension of G.N. Ramachandran's idea of plotting amino acid phi and psi angles to the glycosidic phi, psi, and omega angles formed between carbohydrates. As in traditional Ramachandran plots, these carbohydrate Ramachandran-type (carb-Rama) plots reveal the coupling between the glycosidic angles by displaying the allowed and disallowed conformational space. Considering two-bond glycosidic linkages, there are 18 possible conformational regions that can be defined by (alpha, phi, psi) and (beta, phi, psi), whereas for three-bond linkages, there are 54 possible regions that can be defined by (alpha, phi, psi, omega) and (beta, phi, psi, omega). Illustrating these ideas are molecular dynamic simulations on an implicitly hydrated oligosaccharide (700 ns) and its eight constituent disaccharides (50 ns/disaccharide). For each linkage, we compare and contrast the oligosaccharide and respective disaccharide carb-Rama plots, validate the simulations and the Glycam06 force field through comparison to experimental data, and discuss the general trends observed in the plots.