Microsecond time-scale conformational exchange in proteins: using long molecular dynamics trajectory to simulate NMR relaxation dispersion data

With the advent of ultra-long MD simulations it becomes possible to model microsecond time-scale protein dynamics and, in particular, the exchange broadening effects (R(ex)) as probed by NMR relaxation dispersion measurements. This new approach allows one to identify the exchanging species, including the elusive "excited states". It further helps to map out the exchange network, which is potentially far more complex than the commonly assumed 2- or 3-site schemes. Under fast exchange conditions, this method can be useful for separating the populations of exchanging species from their respective chemical shift differences, thus paving the way for structural analyses. In this study, recent millisecond-long MD trajectory of protein BPTI (Shaw et al. Science 2010, 330, 341) is employed to simulate the time variation of amide (15)N chemical shifts. The results are used to predict the exchange broadening of (15)N lines and, more generally, the outcome of the relaxation dispersion measurements using Carr-Purcell-Meiboom-Gill sequence. The simulated R(ex) effect stems from the fast (~10-100 μs) isomerization of the C14-C38 disulfide bond, in agreement with the prior experimental findings (Grey et al. J. Am. Chem. Soc. 2003, 125, 14324).