Calculation of chemical shift anisotropy in proteins

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 N-15, C-13' and H-1 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.