Project description:There is considerable interest in determining amide-(15)N chemical shift anisotropy (CSA) tensors from biomolecules and understanding their variation for structural and dynamics studies using solution and solid-state NMR spectroscopy and also by quantum chemical calculations. Due to the difficulties associated with the measurement of CSA tensors from membrane proteins, NMR-based structural studies heavily relied on the CSA tensors determined from model systems, typically single crystals of model peptides. In the present study, the principal components of backbone amide-(15)N CSA tensors have been determined using density functional theory for a 16.7 kDa membrane-bound paramagnetic heme containing protein, cytochrome-b(5) (cytb(5)). All the calculations were performed by taking residues within 5 Å distance from the backbone amide-(15)N nucleus of interest. The calculated amide-(15)N CSA spans agree less well with our solution NMR data determined for an effective internuclear distance r(N-H) = 1.023 Å and a constant angle β = 18° that the least shielded component (δ(11)) makes with the N-H bond. The variation of amide-(15)N CSA span obtained using quantum chemical calculations is found to be smaller than that obtained from solution NMR measurements, whereas the trends of the variations are found to be in close agreement. We believe that the results reported in this study will be useful in studying the structure and dynamics of membrane proteins and heme-containing proteins, and also membrane-bound protein-protein complexes such as cytochromes-b5-P450.
Project description:Chemical shift anisotropy (CSA) tensors are essential in the structural and dynamic studies of proteins using NMR spectroscopy. Results from relaxation studies in biomolecular solution and solid-state NMR experiments on aligned samples are routinely interpreted using well-characterized CSA tensors determined from model compounds. Since CSA tensors, particularly the (15)N CSA, highly depend on a number of parameters including secondary structure, electrostatic interaction, and the amino acid sequence, there is a need for accurately determined CSA tensors from proteins. In this study, we report the backbone amide-(15)N CSA tensors for a 16.7-kDa membrane-bound and paramagnetic-heme containing protein, rabbit Cytochrome b(5) (cytb(5)), determined using the (15)N CSA/(15)N-(1)H dipolar transverse cross-correlation rates. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions are -187.9, -166.0, and -161.1 ppm, respectively, with an overall average value of -171.7 ppm. While the average CSA value determined from this study is in good agreement with previous solution NMR experiments on small globular proteins, the CSA value determined for residues in helical conformation is slightly larger, which may be attributed to the paramagnetic effect from Fe(III) of the heme unit in cytb(5). However, like in previous solution NMR studies, the CSA values reported in this study are larger than the values measured from solid-state NMR experiments. We believe that the CSA parameters reported in this study will be useful in determining the structure, dynamics, and orientation of proteins, including membrane proteins, using NMR spectroscopy.
Project description:Individual peptide groups in proteins must exhibit some variation in the chemical shift anisotropy (CSA) of their constituent atoms, but not much is known about the extent or origins of this dispersion. Direct spectroscopic measurement of CSA remains technically challenging, and theoretical methods can help to overcome these limitations by estimating shielding tensors for arbitrary structures. Here we use an automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach to compute (15)N, (13)C' and (1)H chemical shift tensors for human ubiquitin and the GB1 and GB3 fragments of staphylococcal protein G. The average and range of variation of the anisotropies is in good agreement with experimental estimates from solid-state NMR, and the variation among residues is somewhat smaller than that estimated from solution-state measurements. Hydrogen-bond effects account for much of the variation, both between helix and sheet regions, and within elements of secondary structure, but other effects (including variations in torsion angles) may play a role as well.
Project description:It has been well realized that the dependence of chemical shift anisotropy (CSA) tensors on the amino acid sequence, secondary structure, dynamics, and electrostatic interactions can be utilized in the structural and dynamic studies of proteins by NMR spectroscopy. In addition, CSA tensors could also be utilized to measure the structural interactions between proteins in a protein-protein complex. To this end, we report the experimentally measured backbone amide-(15)N CSA tensors for a membrane-bound 16.7 kDa full-length rabbit cytochrome-b5 (cytb5), in complexation with a 55.8 kDa microsomal rabbit cytochrome P450 2B4 (cytP4502B4). The (15)N-CSAs, determined using the (15)N CSA/(15)N-(1)H dipolar coupling transverse cross-correlated rates, for free cytb5 are compared with those for the cytb5 bound to cytP4502B4. An overall increase in backbone amide-(15)N transverse cross-correlated rates for the cytb5 residues in the cytb5-cytP450 complex is observed as compared to the free cytb5 residues. Due to fast spin-spin relaxation (T2) and subsequent broadening of the signals in the complex, we could measure amide-(15)N CSAs only for 48 residues of cytb5 as compared to 84 residues of free cytb5. We observed a change in (15)N CSA for most residues of cytb5 in the complex, as compared to free cytb5, suggesting a dynamic interaction between the oppositely charged surfaces of anionic cytb5 and cationic cytP450. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions of cytb5 in the complex are -184.5, -146.8, and -146.2 ppm, respectively, with an overall average value of -165.5 ppm (excluding the values from residues in more flexible termini). The measured CSA value for residues in helical conformation is slightly larger as compared to previously reported values. This may be attributed to the paramagnetic effect from Fe(III) of the heme in cytb5, which is similar to our previously reported values for the free cytb5.
Project description:A new computer program, called SHIFTX2, is described which is capable of rapidly and accurately calculating diamagnetic (1)H, (13)C and (15)N chemical shifts from protein coordinate data. Compared to its predecessor (SHIFTX) and to other existing protein chemical shift prediction programs, SHIFTX2 is substantially more accurate (up to 26% better by correlation coefficient with an RMS error that is up to 3.3× smaller) than the next best performing program. It also provides significantly more coverage (up to 10% more), is significantly faster (up to 8.5×) and capable of calculating a wider variety of backbone and side chain chemical shifts (up to 6×) than many other shift predictors. In particular, SHIFTX2 is able to attain correlation coefficients between experimentally observed and predicted backbone chemical shifts of 0.9800 ((15)N), 0.9959 ((13)Cα), 0.9992 ((13)Cβ), 0.9676 ((13)C'), 0.9714 ((1)HN), 0.9744 ((1)Hα) and RMS errors of 1.1169, 0.4412, 0.5163, 0.5330, 0.1711, and 0.1231 ppm, respectively. The correlation between SHIFTX2's predicted and observed side chain chemical shifts is 0.9787 ((13)C) and 0.9482 ((1)H) with RMS errors of 0.9754 and 0.1723 ppm, respectively. SHIFTX2 is able to achieve such a high level of accuracy by using a large, high quality database of training proteins (>190), by utilizing advanced machine learning techniques, by incorporating many more features (χ(2) and χ(3) angles, solvent accessibility, H-bond geometry, pH, temperature), and by combining sequence-based with structure-based chemical shift prediction techniques. With this substantial improvement in accuracy we believe that SHIFTX2 will open the door to many long-anticipated applications of chemical shift prediction to protein structure determination, refinement and validation. SHIFTX2 is available both as a standalone program and as a web server ( http://www.shiftx2.ca ).
Project description:Residual dipolar coupling (RDC) and residual chemical shift anisotropy (RCSA) report on orientational properties of a dipolar bond vector and a chemical shift anisotropy principal axis system, respectively. They can be highly complementary in the analysis of backbone structure and dynamics in proteins as RCSAs generally include a report on vectors out of a peptide plane while RDCs usually report on in-plane vectors. Both RDC and RCSA average to zero in isotropic solutions and require partial orientation in a magnetic field to become observable. While the alignment and measurement of RDC has become routine, that of RCSA is less common. This is partly due to difficulties in providing a suitable isotopic reference spectrum for the measurement of the small chemical shift offsets coming from RCSA. Here we introduce a device (modified NMR tube) specifically designed for accurate measurement of reference and aligned spectra for RCSA measurements, but with a capacity for RDC measurements as well. Applications to both soluble and membrane anchored proteins are illustrated.
Project description:A method is introduced to measure residual chemical shift anisotropies conveniently and accurately in the mesophase of poly-?-(benzyl-l-glutamate). The alignment amplitude is substantially enhanced over common methods which greatly benefits measurements particularly on sp3 carbons. The approach offers significant improvements in data accuracy and utility for small molecule structure determination.
Project description:Chemical shift frequencies represent a time-average of all the conformational states populated by a protein. Thus, chemical shift prediction programs based on sequence and database analysis yield higher accuracy for rigid rather than flexible protein segments. Here we show that the prediction accuracy can be significantly improved by averaging over an ensemble of structures, predicted solely from amino acid sequence with the Rosetta program. This approach to chemical shift and structure prediction has the potential to be useful for guiding resonance assignments, especially in solid-state NMR structural studies of membrane proteins in proteoliposomes.
Project description:3D molecular structure determination is a challenge for organic compounds or natural products available in minute amounts. Proton/proton and proton/carbon correlations yield the constitution. J couplings and NOEs oftentimes supported by one-bond 1H,13C residual dipolar couplings (RDCs) or by 13C residual chemical shift anisotropies (RCSAs) provide the relative configuration. However, these RDCs or carbon RCSAs rely on 1% natural abundance of 13C preventing their use for compounds available only in quantities of a few 10's of µgs. By contrast, 1H RCSAs provide similar information on spatial orientation of structural moieties within a molecule, while using the abundant 1H spin. Herein, 1H RCSAs are accurately measured using constrained aligning gels or liquid crystals and applied to the 3D structural determination of molecules with varying complexities. Even more, deuterated alignment media allow the elucidation of the relative configuration of around 35 µg of a briarane compound isolated from Briareum asbestinum.
Project description:This study reports magnitudes and the orientation of the (13)C(alpha) chemical shift anisotropy (CSA) tensors of peptides obtained using quantum chemical calculations. The dependency of the CSA tensor parameters on the energy optimization of hydrogen atom positions and hydrogen bonding effects and the use of zwitterionic peptides in the calculations are examined. Our results indicate that the energy optimization of the hydrogen atom positions in crystal structures is necessary to obtain accurate CSA tensors. The inclusion of intermolecular effects such as hydrogen bonding in the calculations provided better agreement between the calculated and experimental values; however, the use of zwitterionic peptides in calculations, with or without the inclusion of hydrogen bonding, did not improve the results. In addition, our calculated values are in good agreement with tensor values obtained from solid-state NMR experiments on glycine-containing tripeptides. In the case of peptides containing an aromatic residue, calculations on an isolated peptide yielded more accurate isotropic shift values than the calculations on extended structures of the peptide. The calculations also suggested that the presence of an aromatic ring in the extended crystal peptide structure influences the magnitude of the delta(22) which the present level of ab initio calculations are unable to reproduce.