Comparison of MD simulations and NMR experiments for hen lysozyme. Analysis of local fluctuations, cooperative motions, and global changes.

Three 1000 ps molecular dynamics simulations of hen lysozyme have been compared with a range of experimental NMR parameters in order to gain insight into the dynamical properties of the protein and to assess the significance of the motional events observed in the simulations. The simulations, one in vacuum and two in water, were used to estimate interproton distances (for comparison with NOE data), 3JHN alpha and 3J alpha beta coupling constants and 1H-15N order parameters. Comparison of these values with experimental data, particularly NOEs, enabled force field-induced changes to the structure during the simulations to be recognized. It has been shown, however, that these changes can be largely eliminated by slight modifications to the force field. Using a simulation performed in water with this modified force field, it has been found that 1H-15N order parameters calculated for side chain groups in particular correlate well with experimental values and reflect the substantial dependence of these motional properties on the environment, particularly surface exposure, in which the side chain is found. In this case, the simulation then provides models for the motional processes giving rise to the observed experimental data. The results indicate that the order parameter values reflect primarily the number of torsion angles about which rotameric interchange occurs. In addition to local motions, the two different domains of lysozyme have been found to behave differently in the simulations. Possible implications of these differences for the interpretation of unfolding simulations and experimental observations of folding intermediates for lysozyme are discussed.