Analysis of nitroxide spin label motion in a protein-protein complex using multiple frequency EPR spectroscopy

Electron paramagnetic resonance spectroscopy of nitroxide-labeled calmodulin

Calmodulin (CaM) is a highly conserved calcium-binding protein consisting of two homologous domains, each of which contains two EF-hands, that is known to bind well over 300 proteins and peptides. In most cases the (Ca(2+))(4-)form of CaM leads to the activation of a key regulatory enzyme or protein in a myriad of biological processes. Using the nitroxide spin-labeling reagent, 3-(2-iodoacetamido)-2,2,5,5-tetramethyl-1-pyrrolidinyl oxyl, bovine brain CaM was modified at 2-3 methionines with retention of activity as judged by the activation of cyclic nucleotide phosphodiesterase. X-band electron paramagnetic resonance (EPR) spectroscopy was used to measure the spectral changes upon addition of Ca(2+) to the apo-form of spin-labeled protein. A significant loss of spectral intensity, arising primarily from reductions in the heights of the low, intermediate, and high field peaks, accompanied Ca(2+) binding. The midpoint of the Ca(2+)-mediated transition determined by EPR occurred at a higher Ca(2+) concentration than that measured with circular dichroic spectroscopy and enzyme activation. Recent data have indicated that the transition from the apo-state of CaM to the fully saturated form, [(Ca(2+))(4-)CaM], contains a compact intermediate corresponding to [(Ca(2+))(2-)CaM], and the present results suggest that the spin probes are reporting on Ca(2+) binding to the last two sites in the N-terminal domain, i.e. for the [(Ca(2+))(2)-CaM] → [(Ca(2+))(4-)CaM] transition in which the compact structure becomes more extended. EPR of CaM, spin-labeled at methionines, offers a different approach for studying Ca(2+)-mediated conformational changes and may emerge as a useful technique for monitoring interactions with target proteins.

Continuous wave W- and D-band EPR spectroscopy offer "sweet-spots" for characterizing conformational changes and dynamics in intrinsically disordered proteins

Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful tool for characterizing conformational sampling and dynamics in biological macromolecules. Here we demonstrate that nitroxide spectra collected at frequencies higher than X-band (∼9.5 GHz) have sensitivity to the timescale of motion sampled by highly dynamic intrinsically disordered proteins (IDPs). The 68 amino acid protein IA3, was spin-labeled at two distinct sites and a comparison of X-band, Q-band (35 GHz) and W-band (95 GHz) spectra are shown for this protein as it undergoes the helical transition chemically induced by tri-fluoroethanol. Experimental spectra at W-band showed pronounced line shape dispersion corresponding to a change in correlation time from ∼0.3 ns (unstructured) to ∼0.6 ns (α-helical) as indicated by comparison with simulations. Experimental and simulated spectra at X- and Q-bands showed minimal dispersion over this range, illustrating the utility of SDSL EPR at higher frequencies for characterizing structural transitions and dynamics in IDPs.

Information-driven modeling of large macromolecular assemblies using NMR data

Availability of high-resolution atomic structures is one of the prerequisites for a mechanistic understanding of biomolecular function. This atomic information can, however, be difficult to acquire for interesting systems such as high molecular weight and multi-subunit complexes. For these, low-resolution and/or sparse data from a variety of sources including NMR are often available to define the interaction between the subunits. To make best use of all the available information and shed light on these challenging systems, integrative computational tools are required that can judiciously combine and accurately translate the sparse experimental data into structural information. In this Perspective we discuss NMR techniques and data sources available for the modeling of large and multi-subunit complexes. Recent developments are illustrated by particularly challenging application examples taken from the literature. Within this context, we also position our data-driven docking approach, HADDOCK, which can integrate a variety of information sources to drive the modeling of biomolecular complexes. It is the synergy between experimentation and computational modeling that will provides us with detailed views on the machinery of life and lead to a mechanistic understanding of biomolecular function.

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.

Multifrequency electron spin resonance study of the dynamics of spin labeled T4 lysozyme

An extensive set of electron spin resonance spectra was obtained over a wide range of frequencies (9, 95, 170, and 240 GHz) and temperatures (2 to 32 degrees C) to explore the dynamic modes of nitroxide-labeled T4 lysozyme in solution. A commonly used nitroxide side chain (R1), or a methylated analogue with hindered internal motion (R2), was substituted for the native side chain at solvent-exposed helical sites, 72 or 131. The spectra at all four frequencies were simultaneously fit with the slowly relaxing local structure (SRLS) model. Good fits were achieved at all the temperatures. Two principle dynamic modes are included in the SRLS model, the global tumbling of the protein and the internal motion consisting of backbone fluctuations and side chain isomerizations. Three distinct spectral components were required for R1 and two for R2 to account for the spectra at all temperatures. One is a highly ordered and slow motional component, which is observed in the spectra of both R1 and R2; it may correspond to conformers stabilized by interaction with the protein surface. The fraction of this component decreases with increasing temperature and is more populated in the R2 spectra, possibly arising from stronger interaction of the nitroxide ring with the protein surface due to the additional methyl group. The other two components of R1 and the second component of R2 are characterized by fast anisotropic diffusion and relatively low ordering, most likely corresponding to conformers having little or no interactions with nearby residues. Ficoll of different concentrations was added to increase the solution viscosity, thereby slowing down the global tumbling of the protein. A significant effect of Ficoll on the internal motion of an immobilized component was apparent in R2 but not in R1. The ability of such multifrequency studies to separate the effects of faster internal modes of motion from slower overall motions is clearly demonstrated, and its utility in future studies is considered.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

Electron paramagnetic resonance spectra simulation directly from molecular dynamics trajectories of a liquid crystal with a doped paramagnetic spin probe

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.

A novel approach to the simulation of nitroxide spin label EPR spectra from a single truncated dynamical trajectory

A simple effective method for calculation of EPR spectra from a single truncated dynamical trajectory of spin probe orientations is reported. It is shown that an accurate simulation can be achieved from the small initial fraction of a dynamical trajectory until the point when the autocorrelation function of re-orientational motion of spin label has relaxed. This substantially reduces the amount of time for spectra simulation compared to previous approaches, which require multiple full length trajectories (normally of several microseconds) to achieve the desired resolution of EPR spectra. Our method is applicable to trajectories generated from both Brownian dynamics and molecular dynamics (MD) calculations. Simulations of EPR spectra from Brownian dynamical trajectories under a variety of motional conditions including bi-modal dynamics with different hopping rates between the modes are compared to those performed by conventional method. Since the relatively short timescales of spin label motions are realistically accessible by modern MD computational methods, our approach, for the first time, opens the prospect of the simulation of EPR spectra entirely from MD trajectories of real proteins structures.