Quantitative structure–property relationships (QSPRs) on the basis of constitutional, topological, geometrical, and electrostatic descriptors are developed for 2454 13 C α NMR chemical shifts of 21 structure-known, high-quality monomeric proteins. In this procedure, heuristic approach is employed to perform variable-selection for obtaining few independent and significant descriptors. Coupled with various machine learning methods, including MLR, PLS, LSSVM, RF, and GP, these selected variables are then used to create both linear and nonlinear statistical models with the experimentally determined 13 C α NMR chemical shifts of proteins. In addition, the secondary structural effect and environmental influence on protein chemical shifts are also investigated in detail through structural survey and quantum-mechanical calculations. We demonstrate that (i) relationship between 13 C α NMR chemical shifts and local structural features is, to some extent, nonlinear, and (ii) the 13 C α chemical shift values are not only determined by corresponding side-chain conformations, but also affected from the arrangement and configuration of spatially vicinal residues.