Conformational dynamics of a crystalline protein from microsecond-scale molecular dynamics simulations and diffuse X-ray scattering

X-ray diffraction from protein crystals includes both sharply peaked Bragg reflections and diffuse intensity between the peaks. The information in Bragg scattering is limited to what is available in the mean electron density. The diffuse scattering arises from correlations in the electron density variations and therefore contains information about collective motions in proteins. Previous studies using molecular-dynamics (MD) simulations to model diffuse scattering have been hindered by insufficient sampling of the conformational ensemble. To overcome this issue, we have performed a 1.1-μs MD simulation of crystalline staphylococcal nuclease, providing 100-fold more sampling than previous studies. This simulation enables reproducible calculations of the diffuse intensity and predicts functionally important motions, including transitions among at least eight metastable states with different active-site geometries. The total diffuse intensity calculated using the MD model is highly correlated with the experimental data. In particular, there is excellent agreement for the isotropic component of the diffuse intensity, and substantial but weaker agreement for the anisotropic component. Decomposition of the MD model into protein and solvent components indicates that protein–solvent interactions contribute substantially to the overall diffuse intensity. We conclude that diffuse scattering can be used to validate predictions from MD simulations and can provide information to improve MD models of protein motions.Proteins explore many conformations while carrying out their functions in biological systems (1–3). X-ray crystallography is the dominant source of information about protein structure; however, crystal structure models usually consist of just a single major conformation and at most a small portion of the model as alternate conformations. Crystal structures therefore are missing many details about the underlying conformational ensemble (4).Proteins assembled in crystalline arrays, like proteins in solution, exhibit rich conformational diversity (4) and often can perform their native functions (5). Many methods have emerged for using Bragg data to model conformational diversity in protein crystals (6–17). The development of these methods has been important as conformational diversity can lead to inaccuracies in protein structure models (9, 18–20). A key limitation of using the Bragg data, however, is that different models of conformational diversity can yield the same mean electron density.Whereas the Bragg scattering only contains information about the mean electron density, diffuse scattering (diffraction resulting in intensity between the Bragg peaks) is sensitive to spatial correlations in electron density variations (21–28) and therefore contains information about the way that atomic positions vary together in protein crystals. Because models that yield the same mean electron density can yield different correlations in electron density variations, diffuse scattering provides a means to increase the accuracy of crystallography for determining protein conformational variations (29). Peter Moore (30) and Mark Wilson (31) have argued that diffuse scattering should be used to test models of conformational diversity in X-ray crystallography.Several pioneering studies used diffuse scattering to reveal insights into correlated motions in proteins (17, 30, 32–49). Some of these studies used diffuse scattering to experimentally validate predictions of correlated motions from molecular-dynamics (MD) simulations (35–37, 40, 42–44). These studies revealed important insights but were limited by inadequate sampling of the conformational ensemble, leading to lack of convergence of the diffuse scattering calculations (35). Microsecond-scale simulations of staphylococcal nuclease were predicted to be adequate for convergence of diffuse scattering calculations (42). Modern simulation algorithms and computer hardware now enable microsecond or longer MD simulations of protein crystals (50).Here, we present calculations of diffuse X-ray scattering using a 1.1-μs MD simulation of crystalline staphylococcal nuclease. The results demonstrate that we have overcome the past limitation of inadequate sampling. We chose staphylococcal nuclease because the experiments of Wall et al. (49) still represent the only complete, high-quality, 3D diffuse scattering data set from a protein crystal. The calculated diffuse intensity is very similar using two independent halves of the trajectory; the results therefore are reproducible and can be meaningfully compared with the experimental data. The MD simulation provides a rich picture of conformational diversity in the energy landscape of a protein crystal, consisting of at least eight metastable states. Like previous MD studies of crystalline staphylococcal nuclease (42–44), the agreement of the simulation with the total experimental diffuse intensity is excellent, supporting the use of MD simulations to model diffuse scattering data. Unlike previous MD studies, we separately compared the more finely structured, anisotropic component of the diffuse intensity with experimental data. The agreement is substantial but weaker than for the isotropic component, indicating there are inaccuracies in the MD models. Our results therefore point toward using diffuse scattering to improve MD models of protein motions.