A 4D approach for protein H-1 chemical shift prediction was explored. The 4th dimension is the molecular flexibility, mapped using molecular dynamics simulations. The chemical shifts were predicted with a principal component model based on atom coordinates from a database of 40 protein structures. When compared to the corresponding non-dynamic (3D) model, the 4th dimension improved prediction by 6-7%. The prediction method achieved RMS errors of 0.29 and 0.50 ppm for H alpha and HN shifts, respectively. However, for individual proteins the RMS errors were 0.17-0.34 and 0.34-0.65 ppm for the H alpha and HN shifts, respectively. X-ray structures gave better predictions than the corresponding NMR structures, indicating that chemical shifts contain invaluable information about local structures. The H-1 chemical shift prediction tool 4DSPOT is available from http://www.uku.fi/kemia/4dspot.