Fluctuations and Correlations in Crystalline Protein Dynamics: A Simulation Analysis of Staphylococcal Nuclease

Understanding collective motions in protein crystals is likely to furnish insight into functional protein dynamics and will improve models for refinement against diffraction data. Here, four 10 ns molecular dynamics simulations of crystalline Staphylococcal nuclease are reported and analyzed in terms of fluctuations and correlations in atomic motion. The simulation-derived fluctuations strongly correlate with, but are slightly higher than, the values derived from the experimental B-factors. Approximately 70% of the atomic fluctuations are due to internal protein motion. For 65% of the protein atoms the internal fluctuations converge on the nanosecond timescale. Convergence is much slower for the elements of the interatomic displacement correlation matrix — of these, >80% converge within 1 ns for interatomic distances ≲6 Å, but only 10% for separations ≃12 Å. Those collective motions that converged on the nanosecond timescale involve mostly correlations within the β-barrel or between α-helices of the protein. The R-factor with the experimental x-ray diffuse scattering for the crystal, which is determined by the displacement variance-covariance matrix, decreases to 8% after 10 ns simulation. Both the number of converged correlation matrix elements and the R-factor depend logarithmically on time, consistent with a model in which the number of energy minima sampled depends exponentially on the maximum energy barrier crossed. The logarithmic dependence is also extrapolated to predict a convergence time for the whole variance-covariance matrix of ∼1 μs.