Predicting protein backbone chemical shifts from Cα coordinates: extracting high resolution experimental observables from low resolution models.

Abstract

Given the demonstrated utility of coarse-grained modeling and simulations approaches in studying protein structure and dynamics, developing methods that allow experimental observables to be directly recovered from coarse-grained models is of great importance. In this work, we develop one such method that enables protein backbone chemical shifts (1HN, 1Hα, 13Cα, 13C, 13Cβ, and 15N) to be predicted from Cα coordinates. We show that our Cα-based method, LARMORCα, predicts backbone chemical shifts with comparable accuracy to some all-atom approaches. More importantly, we demonstrate that LARMORCα predicted chemical shifts are able to resolve native structure from decoy pools that contain both native and non-native models, and so it is sensitive to protein structure. As an application, we use LARMORCα to characterize the transient state of the fast-folding protein gpW using recently published NMR relaxation dispersion derived backbone chemical shifts. The model we obtain is consistent with the previously proposed model based on independent analysis of the chemical shift dispersion pattern of the transient state. We anticipate that LARMORCα will find utility as a tool that enables important protein conformational substates to be identified by “parsing” trajectories and ensembles generated using coarse-grained modeling and simulations.