Calculation of ¹³C chemical shifts in RNA nucleosides: Structure-¹³C chemical shift relationships

Isotropic ¹³C chemical shifts of the ribose sugar in model RNA nucleosides are calculated using SCF and DFT-GIAO ab initio methods for different combinations of ribose sugar pucker, exocyclic torsion angle, and glycosidic torsion angle. Idealized conformations were obtained using structures that were fully optimized by ab initio DFT methods starting with averaged parameters from a collection of crystallographic data. Solid-state coordinates of accurate crystal or neutron diffraction structures were also examined directly without optimization. The resulting ¹³C chemical shifts for the two sets of calculations are then compared. The GIAO-DFT method overestimates the shifts by an average of 5 ppm while the GIAO-SCF underestimates the shifts by the same amount. However, in the majority of cases the errors appear to be systematic, as the slope of a plot of calculated vs experimental shifts is very close to unity, with minimal scatter. The values of the ¹³C NMR shifts of the ribose sugar are therefore sufficiently precise to allow for statistical separation of sugar puckering modes and exocyclic torsion angle conformers, based on the canonical equation model formulated in a previous paper.