Search

We report a general approach for the simulation of the electron paramagnetic resonance (EPR) spectra of spin labels attached to peptides and proteins directly from replica-exchange molecular dynamics (REMD) trajectories. Conventional MD trajectories are generally inadequate for the prediction of EPR line shapes since the label can become trapped in one or more of a set of rotameric states, thus preventing both conformational sampling and accurate estimates of the exchange rates between different rotamers. The advantage of using REMD is its ability to provide both efficient conformational sampling and kinetic information for spin-label dynamics. Our approach is illustrated with spin-labeled peptide. This approach to REMD-EPR simulation paves the way for the wider application of MD modeling to the simulation and interpretation of EPR spectra of spin-labeled molecules.

We report simulation of EPR spectra directly and entirely from trajectories generated from molecular dynamics simulations. Results are reported for a model 3beta-DOXYL-5alpha-cholestane spin probe in a coarse-grained solvent representing a 5CB nematic host. The results are in excellent agreement with the experimental spectra. The calculated order parameters associated with the paramagnetic probe show strong correlation with the order parameter of 5CB mesogens and are in agreement with those reported in the literature. Simulation of EPR spectra entirely from molecular dynamics of real structures provides direct correlation between molecular motions and the features observed in the spectra, allowing unambiguous interpretation of the spectra. This method opens the possibility for "computer engineering" of spin-labeled materials with the desired properties, such as spin-labeled proteins, prior to experiment.