Journal of the American Chemical Society, Vol.129, No.18, 5898-5903, 2007
Correlation of protein structure and dynamics to scalar couplings across hydrogen bonds
NMR-observable scalar couplings across hydrogen bonds in nucleic acids and proteins present a quantitative measure for the geometry and by the implicit experimental time averaging dynamics of hydrogen bonds. We have carried out in-depth molecular dynamics (MD) simulations with various force fields on three proteins: ubiquitin, the GB1 domain of protein G, and the SMN Tudor domain, for which experimental (h3)J(NC') scalar couplings of backbone hydrogen bonds and various high-resolution X-ray structures are available. Theoretical average values for (h3)J(NC') were calculated from the snapshots of these MD simulations either by density functional theory or by a geometric parametrization (Barfield, M. J. Am. Chem. Soc. 2002, 124, 4158-4168). No significant difference was found between the two methods. The results indicate that time-averaging using explicit water solvation in the MD simulations improves significantly the agreement between experimental and theoretical values for the lower resolution structures ubiquitin (1.8 angstrom), Tudor domain (1.8 angstrom), and protein G (2.1 angstrom). Only marginal improvement is found for the high-resolution structure (1.1 angstrom) of protein G. Hence, experimental (h3)J(NC') values are compatible with a static, high-resolution structural model. The MD averaging of the low-resolution structures moves the averages of the r(HO) distance and the H center dot center dot center dot OC angle theta closer to their respective values in the high-resolution structures, thereby improving the agreement using experimental (h3)J(NC') data. In contrast, MD averaging with implicit water models deteriorates the agreement with experiment for all proteins. The differing behavior can be explained by an artifactual lengthening of H-bonds caused by the implicit water models.