Structural origin of weakly ordered nitroxide motion in spin-labeled proteins

Structural and energetic determinants of adhesive binding specificity in type I cadherins

Type I cadherin cell-adhesion proteins are similar in sequence and structure and yet are different enough to mediate highly specific cell-cell recognition phenomena. It has previously been shown that small differences in the homophilic and heterophilic binding affinities of different type I family members can account for the differential cell-sorting behavior. Here we use a combination of X-ray crystallography, analytical ultracentrifugation, surface plasmon resonance and double electron-electron resonance (DEER) electron paramagnetic resonance spectroscopy to identify the molecular determinants of type I cadherin dimerization affinities. Small changes in sequence are found to produce subtle structural and dynamical changes that impact relative affinities, in part via electrostatic and hydrophobic interactions, and in part through entropic effects because of increased conformational heterogeneity in the bound states as revealed by DEER distance mapping in the dimers. These findings highlight the remarkable ability of evolution to exploit a wide range of molecular properties to produce closely related members of the same protein family that have affinity differences finely tuned to mediate their biological roles.

CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface to prepare molecular simulation systems and input files to facilitate the usage of common and advanced simulation techniques. Since it is originally developed in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to setup a broad range of simulations including free energy calculation and large-scale coarse-grained representation. Here, we describe functionalities that have recently been integrated into CHARMM-GUI PDB Manipulator, such as ligand force field generation, incorporation of methanethiosulfonate spin labels and chemical modifiers, and substitution of amino acids with unnatural amino acids. These new features are expected to be useful in advanced biomolecular modeling and simulation of proteins.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.

High-resolution crystal structure of spin labelled (T21R1) azurin from Pseudomonas aeruginosa: a challenging structural benchmark for in silico spin labelling algorithms

Background

EPR-based distance measurements between spin labels in proteins have become a valuable tool in structural biology. The direct translation of the experimental distances into structural information is however often impaired by the intrinsic flexibility of the spin labelled side chains. Different algorithms exist that predict the approximate conformation of the spin label either by using pre-computed rotamer libraries of the labelled side chain (rotamer approach) or by simply determining its accessible volume (accessible volume approach). Surprisingly, comparisons with many experimental distances have shown that both approaches deliver the same distance prediction accuracy of about 3 Å.

Results

Here, instead of comparing predicted and experimental distances, we test the ability of both approaches to predict the actual conformations of spin labels found in a new high-resolution crystal structure of spin labelled azurin (T21R1). Inside the crystal, the label is found in two very different environments which serve as a challenging test for the in silico approaches.

Conclusions

Our results illustrate why simple and more sophisticated programs lead to the same prediciton error. Thus, a more precise treatment of the complete environment of the label and also its interactions with the environment will be needed to increase the accuracy of in silico spin labelling algorithms.

EPR assessment of protein sites for incorporation of Gd(III) MRI contrast labels

We have engineered apolipoprotein A-I (apoA-I), a major protein constituent of high-density lipoprotein (HDL), to contain DOTA-chelated Gd(III) as an MRI contrast agent for the purpose of imaging reconstituted HDL (rHDL) biodistribution, metabolism and regulation in vivo. This protein contrast agent was obtained by attaching the thiol-reactive Gd[MTS-ADO3A] label at Cys residues replaced at four distinct positions (52, 55, 76 and 80) in apoA-I. MRI of infused mice previously showed that the Gd-labeled apoA-I migrates to both the liver and the kidney, the organs responsible for HDL catabolism; however, the contrast properties of apoA-I are superior when the ADO3A moiety is located at position 55, compared with the protein labeled at positions 52, 76 or 80. It is shown here that continuous wave X-band (9 GHz) electron paramagnetic resonance (EPR) spectroscopy is capable of detecting differences in the Gd(III) signal when comparing the labeled protein in the lipid-free with the rHDL state. Furthermore, the values of NMR relaxivity obtained for labeled variants in both the lipid-free and rHDL states correlate to the product of the X-band Gd(III) spectral width and the collision frequency between a nitroxide spin label and a polar relaxation agent. Consistent with its superior relaxivity measured by NMR, the rHDL-associated apoA-I containing the Gd[MTS-ADO3A] probe attached to position 55 displays favorable dynamic and water accessibility properties as determined by X-band EPR. While room temperature EPR requires >1 m m Gd(III)-labeled and only >10 µ m nitroxide-labeled protein to resolve the spectrum, the volume requirement is exceptionally low (~5 µl). Thus, X-band EPR provides a practical assessment for the suitability of imaging candidates containing the site-directed ADO3A contrast probe.

Conformational selection and adaptation to ligand binding in T4 lysozyme cavity mutants

The studies presented here explore the relationship between protein packing and molecular flexibility using ligand-binding cavity mutants of T4 lysozyme. Although previously reported crystal structures of the mutants investigated show single conformations that are similar to the WT protein, site-directed spin labeling in solution reveals additional conformational substates in equilibrium exchange with a WT-like population. Remarkably, binding of ligands, including the general anesthetic halothane shifts the population to the WT-like state, consistent with a conformational selection model of ligand binding, but structural adaptation to the ligand is also apparent in one mutant. Distance mapping with double electron-electron resonance spectroscopy and the absence of ligand binding suggest that the new substates induced by the cavity-creating mutations represent alternate packing modes in which the protein fills or partially fills the cavity with side chains, including the spin label in one case; external ligands compete with the side chains for the cavity space, stabilizing the WT conformation. The results have implications for mechanisms of anesthesia, the response of proteins to hydrostatic pressure, and protein engineering.

RosettaEPR: rotamer library for spin label structure and dynamics

An increasingly used parameter in structural biology is the measurement of distances between spin labels bound to a protein. One limitation to these measurements is the unknown position of the spin label relative to the protein backbone. To overcome this drawback, we introduce a rotamer library of the methanethiosulfonate spin label (MTSSL) into the protein modeling program Rosetta. Spin label rotamers were derived from conformations observed in crystal structures of spin labeled T4 lysozyme and previously published molecular dynamics simulations. Rosetta’s ability to accurately recover spin label conformations and EPR measured distance distributions was evaluated against 19 experimentally determined MTSSL labeled structures of T4 lysozyme and the membrane protein LeuT and 73 distance distributions from T4 lysozyme and the membrane protein MsbA. For a site in the core of T4 lysozyme, the correct spin label conformation (Χ₁ and Χ₂) is recovered in 99.8% of trials. In surface positions 53% of the trajectories agree with crystallized conformations in Χ₁ and Χ₂. This level of recovery is on par with Rosetta performance for the 20 natural amino acids. In addition, Rosetta predicts the distance between two spin labels with a mean error of 4.4 Å. The width of the experimental distance distribution, which reflects the flexibility of the two spin labels, is predicted with a mean error of 1.3 Å. RosettaEPR makes full-atom spin label modeling available to a wide scientific community in conjunction with the powerful suite of modeling methods within Rosetta.

Technological advances in site-directed spin labeling of proteins

Molecular flexibility over a wide time range is of central importance to the function of many proteins, both soluble and membrane. Revealing the modes of flexibility, their amplitudes, and time scales under physiological conditions is the challenge for spectroscopic methods, one of which is site-directed spin labeling EPR (SDSL-EPR). Here we provide an overview of some recent technological advances in SDSL-EPR related to investigation of structure, structural heterogeneity, and dynamics of proteins. These include new classes of spin labels, advances in measurement of long range distances and distance distributions, methods for identifying backbone and conformational fluctuations, and new strategies for determining the kinetics of protein motion.

Structure and dynamics of an imidazoline nitroxide side chain with strongly hindered internal motion in proteins

A disulfide-linked imidazoline nitroxide side chain (V1) has a similar and highly constrained internal motion at diverse topological sites in a protein, unlike that for the disulfide-linked pyrroline nitroxide side chain (R1) widely used in site directed spin labeling EPR. Crystal structures of V1 at two positions in a helix of T4 Lysozyme and quantum mechanical calculations suggest the source of the constraints as intra-side chain interactions of the disulfide sulfur atoms with both the protein backbone and the 3-nitrogen in the imidazoline ring. These interactions apparently limit the conformation of the side chain to one of only three possible rotamers, two of which are observed in the crystal structure. An inter-spin distance measurement in frozen solution using double electron-electron resonance (DEER) gives a value essentially identical to that determined from the crystal structure of the protein containing two copies of V1, indicating that lattice forces do not dictate the rotamers observed. Collectively, the results suggest the possibility of predetermining a unique rotamer of V1 in helical structures. In general, the reduced rotameric space of V1 compared to R1 should simplify interpretation of inter-spin distance information in terms of protein structure, while the highly constrained internal motion is expected to extend the dynamic range for characterizing large amplitude nanosecond backbone fluctuations.

Conformational dynamics and distribution of nitroxide spin labels

Long-range distance measurements based on paramagnetic relaxation enhancement (PRE) in NMR, quantification of surface water dynamics near biomacromolecules by Overhauser dynamic nuclear polarization (DNP) and sensitivity enhancement by solid-state DNP all depend on introducing paramagnetic species into an otherwise diamagnetic NMR sample. The species can be introduced by site-directed spin labeling, which offers precise control for positioning the label in the sequence of a biopolymer. However, internal flexibility of the spin label gives rise to dynamic processes that potentially influence PRE and DNP behavior and leads to a spatial distribution of the electron spin even in solid samples. Internal dynamics of spin labels and their static conformational distributions have been studied mainly by electron paramagnetic resonance spectroscopy and molecular dynamics simulations, with a large body of results for the most widely applied methanethiosulfonate spin label MTSL. These results are critically discussed in a unifying picture based on rotameric states of the group that carries the spin label. Deficiencies in our current understanding of dynamics and conformations of spin labeled groups and of their influence on NMR observables are highlighted and directions for further research suggested.

Structural refinement from restrained-ensemble simulations based on EPR/DEER data: application to T4 lysozyme

DEER (double electron-electron resonance) is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of distance histograms between pairs of nitroxide spin-labels linked to a protein in a native-like solution environment. However, exploiting the huge amount of information provided by ESR/DEER histograms to refine structural models is extremely challenging. In this study, a restrained ensemble (RE) molecular dynamics (MD) simulation methodology is developed to address this issue. In RE simulation, the spin-spin distance distribution histograms calculated from a multiple-copy MD simulation are enforced, via a global ensemble-based energy restraint, to match those obtained from ESR/DEER experiments. The RE simulation is applied to 51 ESR/DEER distance histogram data from spin-labels inserted at 37 different positions in T4 lysozyme (T4L). The rotamer population distribution along the five dihedral angles connecting the nitroxide ring to the protein backbone is determined and shown to be consistent with available information from X-ray crystallography. For the purpose of structural refinement, the concept of a simplified nitroxide dummy spin-label is designed and parametrized on the basis of these all-atom RE simulations with explicit solvent. It is demonstrated that RE simulations with the dummy nitroxide spin-labels imposing the ESR/DEER experimental distance distribution data are able to systematically correct and refine a series of distorted T4L structures, while simple harmonic distance restraints are unsuccessful. This computationally efficient approach allows experimental restraints from DEER experiments to be incorporated into RE simulations for efficient structural refinement.

Mapping molecular flexibility of proteins with site-directed spin labeling: a case study of myoglobin

Site-directed spin labeling (SDSL) has potential for mapping protein flexibility under physiological conditions. The purpose of the present study was to explore this potential using 38 singly spin-labeled mutants of myoglobin distributed throughout the sequence. Correlation of the EPR spectra with protein structure provides new evidence that the site-dependent variation in line shape, and hence motion of the spin label, is due largely to differences in mobility of the helical backbone in the ns time range. Fluctuations between conformational substates, typically in the μs-ms time range, are slow on the EPR time scale, and the spectra provide a snapshot of conformational equilibria frozen in time as revealed by multiple components in the spectra. A recent study showed that osmolyte perturbation can positively identify conformational exchange as the origin of multicomponent spectra (López et al. (2009), Protein Sci. 18, 1637). In the present study, this new strategy is employed in combination with line shape analysis and pulsed-EPR interspin distance measurements to investigate the conformation and flexibility of myoglobin in three folded and partially folded states. The regions identified to be in conformational exchange in the three forms agree remarkably well with those assigned by NMR, but the faster time scale of EPR allows characterization of localized states not detected in NMR. Collectively, the results suggest that SDSL-EPR and osmolyte perturbation provide a facile means for mapping the amplitude of fast backbone fluctuations and for detecting sequences in slow conformational exchange in folded and partially folded protein sequences.

High-resolution structure of a protein spin-label in a solvent-exposed β-sheet and comparison with DEER spectroscopy

X-ray crystallography has been a useful tool in the development of site-directed spin labeling by resolving rotamers of the nitroxide spin-label side chain in a variety of α-helical environments. In this work, the crystal structure of a doubly spin-labeled N8C/K28C mutant of the B1 immunoglobulin-binding domain of protein G (GB1) was solved. The double mutant formed a domain-swapped dimer under crystallization conditions. Two rotameric states of the spin-label were resolved at the solvent-exposed α-helical site, at residue 28; these are in good agreement with rotamers previously reported for helical structures. The second site, at residue 8 on an interior β-strand, shows the presence of three distinct solvent-exposed side-chain rotamers. One of these rotamers is rarely observed within crystal structures of R1 sites and suggests that the H(α) and S(δ) hydrogen bond that is common to α-helical sites is absent at this interior β-strand residue. Variable temperature continuous wave (CW) experiments of the β-strand site showed two distinct components that were correlated to the rotameric states observed in crystallography. Interestingly, the CW data at room temperature could be fit without the use of an order parameter, which is consistent with the lack of the H(α) and S(δ) interaction. Additionally, double electron electron resonance (DEER) spectroscopy was performed on the GB1 double mutant in its monomeric form and yielded a most probable interspin distance of 25 ± 1 Å. In order to evaluate the accuracy of the measured DEER distance, the rotamers observed in the crystal structure of the domain-swapped GB1 dimer were modeled into a high-resolution structure of the wild type monomeric GB1. The distances generated in the resulting GB1 structural models match the most probable DEER distance within ~2 Å. The results are interesting as they indicate by direct experimental measurement that the rotameric states of R1 found in this crystal provide a very close match to the most probable distance measured by DEER.

Inter-spin distance determination using L-band (1-2 GHz) non-adiabatic rapid sweep electron paramagnetic resonance (NARS EPR)

Site-directed spin-labeling electron paramagnetic resonance (SDSL EPR) provides insight into the local structure and motion of a spin probe strategically attached to a molecule. When a second spin is introduced to the system, macromolecular information can be obtained through measurement of inter-spin distances either by continuous wave (CW) or pulsed electron double resonance (ELDOR) techniques. If both methodologies are considered, inter-spin distances of 8-80 Å can be experimentally determined. However, there exists a region at the upper limit of the conventional X-band (9.5 GHz) CW technique and the lower limit of the four-pulse double electron-electron resonance (DEER) experiment where neither method is particularly reliable. The work presented here utilizes L-band (1.9 GHz) in combination with non-adiabatic rapid sweep (NARS) EPR to address this opportunity by increasing the upper limit of the CW technique. Because L-band linewidths are three to seven times narrower than those at X-band, dipolar broadenings that are small relative to the X-band inhomogeneous linewidth become observable, but the signal loss, due to the frequency dependence of the Boltzmann factor, has made L-band especially challenging. NARS has been shown to increase sensitivity by a factor of five, and overcomes much of this loss, making L-band distance determination more feasible. Two different systems are presented, and distances of 18-30 Å have been experimentally determined at physiologically relevant temperatures. Measurements are in excellent agreement with a helical model and values determined by DEER.

Simulation vs. reality: a comparison of in silico distance predictions with DEER and FRET measurements

Site specific incorporation of molecular probes such as fluorescent- and nitroxide spin-labels into biomolecules, and subsequent analysis by Förster resonance energy transfer (FRET) and double electron-electron resonance (DEER) can elucidate the distance and distance-changes between the probes. However, the probes have an intrinsic conformational flexibility due to the linker by which they are conjugated to the biomolecule. This property minimizes the influence of the label side chain on the structure of the target molecule, but complicates the direct correlation of the experimental inter-label distances with the macromolecular structure or changes thereof. Simulation methods that account for the conformational flexibility and orientation of the probe(s) can be helpful in overcoming this problem. We performed distance measurements using FRET and DEER and explored different simulation techniques to predict inter-label distances using the Rpo4/7 stalk module of the M. jannaschii RNA polymerase. This is a suitable model system because it is rigid and a high-resolution X-ray structure is available. The conformations of the fluorescent labels and nitroxide spin labels on Rpo4/7 were modeled using in vacuo molecular dynamics simulations (MD) and a stochastic Monte Carlo sampling approach. For the nitroxide probes we also performed MD simulations with explicit water and carried out a rotamer library analysis. Our results show that the Monte Carlo simulations are in better agreement with experiments than the MD simulations and the rotamer library approach results in plausible distance predictions. Because the latter is the least computationally demanding of the methods we have explored, and is readily available to many researchers, it prevails as the method of choice for the interpretation of DEER distance distributions.

Simulating the dynamics and orientations of spin-labeled side chains in a protein-DNA complex

Site-directed spin labeling, wherein a nitroxide side chain is introduced into a protein at a selected mutant site, is increasingly employed to investigate biological systems by electron spin resonance (ESR) spectroscopy. An understanding of the packing and dynamics of the spin label is needed to extract the biologically relevant information about the macromolecule from ESR measurements. In this work, molecular dynamics (MD) simulations were performed on the spin-labeled restriction endonuclease, EcoRI in complex with DNA. Mutants of this homodimeric enzyme were previously constructed, and distance measurements were performed using the double electron electron resonance experiment. These correlated distance constraints have been leveraged with MD simulations to learn about side chain packing and preferred conformers of the spin label on sites in an α-helix and a β-strand. We found three dihedral angles of the spin label side chain to be most sensitive to the secondary structure where the spin label was located. Conformers sampled by the spin label differed between secondary structures as well. C(α)-C(α) distance distributions were constructed and used to extract details about the protein backbone mobility at the two spin labeled sites. These simulation studies enhance our understanding of the behavior of spin labels in proteins and thus expand the ability of ESR spectroscopy to contribute to knowledge of protein structure and dynamics.

MtsslWizard: In Silico Spin-Labeling and Generation of Distance Distributions in PyMOL

MtsslWizard is a computer program, which operates as a plugin for the PyMOL molecular graphics system. MtsslWizard estimates distances between spin labels on proteins quickly with user-configurable options through a simple graphical interface. In default mode, the program searches for ensembles of possible MTSSL conformations that do not clash with a static model of the protein. Once conformations are assigned, distance distributions between two or more ensembles are calculated, displayed, and can be exported to other software. The program's use is evaluated in a number of challenging test cases and its strengths and weaknesses evaluated. The benefits of the program are its accuracy and simplicity. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00723-012-0314-0) contains supplementary material, which is available to authorized users.

Identification and removal of nitroxide spin label contaminant: impact on PRE studies of α-helical membrane proteins in detergent

NMR paramagnetic relaxation enhancement (PRE) provides long-range distance constraints (~15-25 Å) that can be critical to determining overall protein topology, especially where long-range NOE information is unavailable such as in the case of larger proteins that require deuteration. However, several challenges currently limit the use of NMR PRE for α-helical membrane proteins. One challenge is the nonspecific association of the nitroxide spin label to the protein-detergent complex that can result in spurious PRE derived distance restraints. The effect of the nitroxide spin label contaminant is evaluated and quantified and a robust method for the removal of the contaminant is provided to advance the application of PRE restraints to membrane protein NMR structure determination.

DEER distance measurements on proteins

Distance distributions between paramagnetic centers in the range of 1.8 to 6 nm in membrane proteins and up to 10 nm in deuterated soluble proteins can be measured by the DEER technique. The number of paramagnetic centers and their relative orientation can be characterized. DEER does not require crystallization and is not limited with respect to the size of the protein or protein complex. Diamagnetic proteins are accessible by site-directed spin labeling. To characterize structure or structural changes, experimental protocols were optimized and techniques for artifact suppression were introduced. Data analysis programs were developed, and it was realized that interpretation of the distance distributions must take into account the conformational distribution of spin labels. First methods have appeared for deriving structural models from a small number of distance constraints. The present scope and limitations of the technique are illustrated.

The "beta-clasp" model of apolipoprotein A-I--a lipid-free solution structure determined by electron paramagnetic resonance spectroscopy

Apolipoprotein A-I (apoA-I) is the major protein component of high density lipoproteins (HDL) and plays a central role in cholesterol metabolism. The lipid-free/lipid-poor form of apoA-I is the preferred substrate for the ATP-binding cassette transporter A1 (ABCA1). The interaction of apoA-I with ABCA1 leads to the formation of cholesterol laden high density lipoprotein (HDL) particles, a key step in reverse cholesterol transport and the maintenance of cholesterol homeostasis. Knowledge of the structure of lipid-free apoA-I is essential to understanding its critical interaction with ABCA1 and the molecular mechanisms underlying HDL biogenesis. We therefore examined the structure of lipid-free apoA-I by electron paramagnetic resonance spectroscopy (EPR). Through site directed spin label EPR, we mapped the secondary structure of apoA-I and identified sites of spin coupling as residues 26, 44, 64, 167, 217 and 226. We capitalize on the fact that lipid-free apoA-I self-associates in an anti-parallel manner in solution. We employed these sites of spin coupling to define the central plane in the dimeric apoA-I complex. Applying both the constraints of dipolar coupling with the EPR-derived pattern of solvent accessibility, we assembled the secondary structure into a tertiary context, providing a solution structure for lipid-free apoA-I. This article is part of a Special Issue entitled Advances in High Density Lipoprotein Formation and Metabolism: A Tribute to John F. Oram (1945-2010).

Lipid and membrane mimetic environments modulate spin label side chain configuration in the outer membrane protein A

In the present work, the factors that determine EPR line shapes from spin labels at the protein-hydrocarbon interface of a β-barrel membrane protein are examined. The EPR spectra from hydrocarbon facing sites in the outer membrane protein A (OmpA) are highly dependent upon the detergent or lipid into which OmpA is reconstituted. In general, line shapes at these sites are correlated with the solvent accessibility in the supporting amphiphile. A notable exception is CHAPS, which yields rigid limit EPR line shapes for labels at every position along a transmembrane β-strand in OmpA. EPR line shapes from the surface of OmpA are not strongly influenced by steric interference with neighboring side chains, but are modulated by solutes that should interact with hydrophobic surfaces. These results suggest that differences in EPR spectra in different supporting environments are not the result of differences in protein dynamics but are a result of different configurations or rotameric states that are assumed by the label. This conclusion is supported by distance measurements across the OmpA β-barrel, which indicate that labels yielding more motionally restricted line shapes interact more closely with the protein surface. These results have implications for the use of spin-label-derived distance constraints in protein structure determination and demonstrate that spin labels on membrane proteins provide a highly sensitive probe for the environment surrounding a membrane protein.

Molecular origin of electron paramagnetic resonance line shapes on β-barrel membrane proteins: the local solvation environment modulates spin-label configuration

In this work, electron paramagnetic resonance (EPR) spectroscopy and X-ray crystallography were used to examine the origins of EPR line shapes from spin-labels at the protein-lipid interface on the β-barrel membrane protein BtuB. Two atomic-resolution structures were obtained for the methanethiosulfonate spin-label derivatized to cysteines on the membrane-facing surface of BtuB. At one of these sites, position 156, the label side chain resides in a pocket formed by neighboring residues; however, it extends from the protein surface and yields a single-component EPR spectrum in the crystal that results primarily from fast rotation about the fourth and fifth bonds linking the spin-label to the protein backbone. In lipid bilayers, site 156 yields a multicomponent spectrum resulting from different rotameric states of the labeled side chain. Moreover, changes in the lipid environment, such as variations in bilayer thickness, modulate the EPR spectrum by modulating label rotamer populations. At a second site, position 371, the labeled side chain interacts with a pocket on the protein surface, leading to a highly immobilized single-component EPR spectrum that is not sensitive to hydrocarbon thickness. This spectrum is similar to that seen at other sites that are deep in the hydrocarbon, such as position 170. This work indicates that the rotameric states of spin-labels on exposed hydrocarbon sites are sensitive to the environment at the protein-hydrocarbon interface, and that this environment may modulate weak interactions between the labeled side chain and the protein surface. In the case of BtuB, lipid acyl chain packing is not symmetric around the β-barrel, and EPR spectra from labeled hydrocarbon-facing sites in BtuB may reflect this asymmetry. In addition to facilitating the interpretation of EPR spectra of membrane proteins, these results have important implications for the use of long-range distance restraints in protein structure refinement that are obtained from spin-labels.

Structure and dynamics of a conformationally constrained nitroxide side chain and applications in EPR spectroscopy

A disulfide-linked nitroxide side chain (R1) is the most widely used spin label for determining protein topology, mapping structural changes, and characterizing nanosecond backbone motions by site-directed spin labeling. Although the internal motion of R1 and the number of preferred rotamers are limited, translating interspin distance measurements and spatial orientation information into structural constraints is challenging. Here, we introduce a highly constrained nitroxide side chain designated RX as an alternative to R1 for these applications. RX is formed by a facile cross-linking reaction of a bifunctional methanethiosulfonate reagent with pairs of cysteine residues at i and i + 3 or i and i + 4 in an α-helix, at i and i + 2 in a β-strand, or with cysteine residues in adjacent strands in a β-sheet. Analysis of EPR spectra, a crystal structure of RX in T4 lysozyme, and pulsed electron-electron double resonance (ELDOR) spectroscopy on an immobilized protein containing RX all reveal a highly constrained internal motion of the side chain. Consistent with the constrained geometry, interspin distance distributions between pairs of RX side chains are narrower than those from analogous R1 pairs. As an important consequence of the constrained internal motion of RX, spectral diffusion detected with ELDOR reveals microsecond internal motions of the protein. Collectively, the data suggest that the RX side chain will be useful for distance mapping by EPR spectroscopy, determining spatial orientation of helical segments in oriented specimens, and measuring structural fluctuations on the microsecond time scale.

A rigid disulfide-linked nitroxide side chain simplifies the quantitative analysis of PRE data

The measurement of (1)H transverse paramagnetic relaxation enhancement (PRE) has been used in biomolecular systems to determine long-range distance restraints and to visualize sparsely-populated transient states. The intrinsic flexibility of most nitroxide and metal-chelating paramagnetic spin-labels, however, complicates the quantitative interpretation of PREs due to delocalization of the paramagnetic center. Here, we present a novel, disulfide-linked nitroxide spin label, R1p, as an alternative to these flexible labels for PRE studies. When introduced at solvent-exposed α-helical positions in two model proteins, calmodulin (CaM) and T4 lysozyme (T4L), EPR measurements show that the R1p side chain exhibits dramatically reduced internal motion compared to the commonly used R1 spin label (generated by reacting cysteine with the spin labeling compound often referred to as MTSL). Further, only a single nitroxide position is necessary to account for the PREs arising from CaM S17R1p, while an ensemble comprising multiple conformations is necessary for those observed for CaM S17R1. Together, these observations suggest that the nitroxide adopts a single, fixed position when R1p is placed at solvent-exposed α-helical positions, greatly simplifying the interpretation of PRE data by removing the need to account for the intrinsic flexibility of the spin label.

Structural changes in bacteriorhodopsin during in vitro refolding from a partially denatured state

We report on the formation of the secondary and tertiary structure of bacteriorhodopsin during its in vitro refolding from an SDS-denatured state. We used the mobility of single spin labels in seven samples, attached at various locations to six of the seven helical segments to engineered cysteine residues, to follow coil-to-helix formation. Distance measurements obtained by spin dipolar quenching in six samples labeled at either the cytoplasmic or extracellular ends of pairs of helices revealed the time dependence of the recovery of the transmembrane helical bundle. The secondary structure in the majority of the helical segments refolds with a time constant of <100-140 ms. Recovery of the tertiary structure is achieved by sequential association of the helices and occurs in at least three distinct steps with time constants of 1), well below 1 s; 2), 3-4 s; and 3), 60-130 s (the latter depending on the helical pair). The slowest of these processes occurs in concert with recovery of the retinal chromophore.

Site-directed spin labeling electron paramagnetic resonance study of the ORF1 protein from a mouse L1 retrotransposon

Long interspersed nuclear element-1 is a highly abundant mammalian retrotransposon that comprises 17% of the human genome. L1 retrotransposition requires the protein encoded by open reading frame-1 (ORF1p), which binds single-stranded RNA with high affinity and functions as a nucleic acid chaperone. ORF1p has been shown to adopt a homo-trimeric, asymmetric dumbbell-shaped structure. However, its atomic-level structure and mechanism of RNA binding remains poorly understood. Here, we report the results of a site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) study of 27 residues within the RNA binding region of the full-length protein. The EPR data are compatible with the large RNA binding lobe of ORF1p containing a RNA recognition motif (RRM) domain and a carboxyl-terminal domain (CTD) that are predicted from crystallographic and NMR studies of smaller fragments of the protein. Interestingly, the EPR data indicate that residues in strands β3 and β4 of the RRM are structurally unstable, compatible with the previously observed sensitivity of this region to proteolysis. Affinity measurements and RNA-dependent EPR spectral changes map the RNA binding site on ORF1p to residues located in strands β3 and β4 of the RRM domain and to helix α1 of the CTD. Complementary in vivo studies also identify residues within the RRM domain that are required for retrotransposition. We propose that in the context of the full-length trimeric protein these distinct surfaces are positioned adjacent to one another providing a continuous surface that may interact with nucleic acids.

Interaction of a G protein with an activated receptor opens the interdomain interface in the alpha subunit

In G-protein signaling, an activated receptor catalyzes GDP/GTP exchange on the G(α) subunit of a heterotrimeric G protein. In an initial step, receptor interaction with G(α) acts to allosterically trigger GDP release from a binding site located between the nucleotide binding domain and a helical domain, but the molecular mechanism is unknown. In this study, site-directed spin labeling and double electron-electron resonance spectroscopy are employed to reveal a large-scale separation of the domains that provides a direct pathway for nucleotide escape. Cross-linking studies show that the domain separation is required for receptor enhancement of nucleotide exchange rates. The interdomain opening is coupled to receptor binding via the C-terminal helix of G(α), the extension of which is a high-affinity receptor binding element.

High-pressure EPR reveals conformational equilibria and volumetric properties of spin-labeled proteins

Identifying equilibrium conformational exchange and characterizing conformational substates is essential for elucidating mechanisms of function in proteins. Site-directed spin labeling has previously been employed to detect conformational changes triggered by some event, but verifying conformational exchange at equilibrium is more challenging. Conformational exchange (microsecond-millisecond) is slow on the EPR time scale, and this proves to be an advantage in directly revealing the presence of multiple substates as distinguishable components in the EPR spectrum, allowing the direct determination of equilibrium constants and free energy differences. However, rotameric exchange of the spin label side chain can also give rise to multiple components in the EPR spectrum. Using spin-labeled mutants of T4 lysozyme, it is shown that high-pressure EPR can be used to: (i) demonstrate equilibrium between spectrally resolved states, (ii) aid in distinguishing conformational from rotameric exchange as the origin of the resolved states, and (iii) determine the relative partial molar volume (ΔV°) and isothermal compressibility (Δβ(T)) of conformational substates in two-component equilibria from the pressure dependence of the equilibrium constant. These volumetric properties provide insight into the structure of the substates. Finally, the pressure dependence of internal side-chain motion is interpreted in terms of volume fluctuations on the nanosecond time scale, the magnitude of which may reflect local backbone flexibility.

Conformational analysis of a nitroxide side chain in an α-helix with density functional theory

In site directed spin labeling, a nitroxide side chain is introduced at a selected site in a protein; the most commonly used is a disulfide-linked side chain designated R1. The electron paramagnetic resonance (EPR) spectra of R1, and the interspin distance between pairs of R1 residues as determined by dipolar EPR spectroscopy, encode a wealth of information on the protein structure and dynamics. However, extracting this information requires structural and dynamical models of the R1 side chain, that is, the favored rotamers, the intraresidue interactions that stabilize them, and the internal modes of motion. X-ray crystal structures of R1 in proteins have revealed a set of preferred rotamers in the crystal lattice. To identify the intraresidue interactions that stabilize particular rotamers of R1 in the absence of interactions with nearby side chains in a helix, and to evaluate models for the internal motion of the side chain, quantum mechanical calculations were performed on a relevant fragment of R1 in a 10-residue α-helix. Relative rotamer energies were determined in the gas phase, and solvation energies were estimated from a continuum solvent model that includes both electrostatic and hydrophobic contributions. The results identified preferred rotamers that are in agreement with the X-ray crystallographic studies. The rotamers are apparently stabilized by intraresidue sulfur-backbone interactions, suggesting that the preferred rotamers may be the same at all solvent-exposed helix sites.

Crystal structure of a DNA containing the planar, phenoxazine-derived bi-functional spectroscopic probe C

Previously, we developed the deoxycytosine analog Ç (C-spin) as a bi-functional spectroscopic probe for the study of nucleic acid structure and dynamics using electron paramagnetic resonance (EPR) and fluorescence spectroscopy. To understand the effect of Ç on nucleic acid structure, we undertook a detailed crystallographic analysis. A 1.7 Å resolution crystal structure of Ç within a decamer duplex A-form DNA confirmed that Ç forms a non-perturbing base pair with deoxyguanosine, as designed. In the context of double-stranded DNA Ç adopted a planar conformation. In contrast, a crystal structure of the free spin-labeled base ç displayed a ∼ 20° bend at the oxazine linkage. Density function theory calculations revealed that the bent and planar conformations are close in energy and exhibit the same frequency for bending. These results indicate a small degree of flexibility around the oxazine linkage, which may be a consequence of the antiaromaticity of a 16-π electron ring system. Within DNA, the amplitude of the bending motion is restricted, presumably due to base-stacking interactions. This structural analysis shows that the Ç forms a planar, structurally non-perturbing base pair with G indicating it can be used with high confidence in EPR- or fluorescence-based structural and dynamics studies.

Rotamer libraries of spin labelled cysteines for protein studies

Studies of structure and dynamics of proteins using site-directed spin labelling rely on explicit modelling of spin label conformations. The large computational effort associated with such modelling with molecular dynamics (MD) simulations can be avoided by a rotamer library approach based on a coarse-grained representation of the conformational space of the spin label. We show here that libraries of about 200 rotamers, obtained by iterative projection of a long MD trajectory of the free spin label onto a set of canonical dihedral angles, provide a representation of the underlying trajectory adequate for EPR distance measurements. Rotamer analysis was performed on selected X-ray structures of spin labelled T4 lysozyme mutants to characterize the spin label rotamer ensemble on a single protein site. Furthermore, predictions based on the rotamer library approach are shown to be in nearly quantitative agreement with electron paramagnetic resonance (EPR) distance data on the Na(+)/H(+) antiporter NhaA and on the light-harvesting complex LHCII whose structures are known from independent cryo electron microscopy and X-ray studies, respectively. Suggestions for the selection of labelling sites in proteins are given, limitations of the approach discussed, and requirements for further development are outlined.

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.

Synaptotagmin 1 and SNAREs form a complex that is structurally heterogeneous

Synaptotagmin 1 (syt1) functions as a Ca(2+)-sensor for neuronal exocytosis. Here, site-directed spin labeling was used to examine the complex formed between a soluble fragment of syt1, which contains its two C2 domains, and the neuronal core soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) complex. Changes in electron paramagnetic resonance lineshape and accessibility for spin-labeled syt1 mutants indicate that in solution, the assembled core SNARE complex contacts syt1 in several regions. For the C2B domain, contact occurs in the polybasic face and sites opposite the Ca(2+)-binding loops. For the C2A domain, contact is seen with the SNARE complex in a region near loop 2. Double electron-electron resonance was used to estimate distances between the two C2 domains of syt1. These distances have broad distributions in solution, which do not significantly change when syt1 is fully associated with the core SNARE complex. The broad distance distributions indicate that syt1 is structurally heterogeneous when bound to the SNAREs and does not assume a well-defined structure. Simulated annealing using electron paramagnetic resonance-derived distance restraints produces a family of syt1 structures where the Ca(2+)-binding regions of each domain face in roughly opposite directions. The results suggest that when associated with the SNAREs, syt1 is configured to bind opposing bilayers, but that the syt1/SNARE complex samples multiple conformational states.

RosettaEPR: an integrated tool for protein structure determination from sparse EPR data

Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is often used for the structural characterization of proteins that elude other techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR). However, high-resolution structures are difficult to obtain due to uncertainty in the spin label location and sparseness of experimental data. Here, we introduce RosettaEPR, which has been designed to improve de novo high-resolution protein structure prediction using sparse SDSL-EPR distance data. The "motion-on-a-cone" spin label model is converted into a knowledge-based potential, which was implemented as a scoring term in Rosetta. RosettaEPR increased the fractions of correctly folded models ( [Formula: see text] <7.5Å) and models accurate at medium resolution ( [Formula: see text] <3.5Å) by 25%. The correlation of score and model quality increased from 0.42 when using no restraints to 0.51 when using bounded restraints and again to 0.62 when using RosettaEPR. This allowed for the selection of accurate models by score. After full-atom refinement, RosettaEPR yielded a 1.7Å model of T4-lysozyme, thus indicating that atomic detail models can be achieved by combining sparse EPR data with Rosetta. While these results indicate RosettaEPR's potential utility in high-resolution protein structure prediction, they are based on a single example. In order to affirm the method's general performance, it must be tested on a larger and more versatile dataset of proteins.

Electron paramagnetic resonance studies of functionally active, nitroxide spin-labeled peptide analogues of the C-terminus of a G-protein alpha subunit

The C-terminal tail of the transducin alpha subunit, Gtalpha(340-350), is known to bind and stabilize the active conformation of rhodopsin upon photoactivation (R*). Five spin-labeled analogues of Gtalpha(340-350) demonstrated native-like activity in their ability to bind and stabilize R*. The spin-label 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC) was employed at interior sites within the peptide, whereas a Proxyl (3-carboxyl-2,2,5,5-tetramethyl-pyrrolidinyloxy) spin-label was employed at the amino terminus of the peptide. Upon binding to R*, the electron paramagnetic resonance spectrum of TOAC(343)-Gtalpha(340-350) revealed greater immobilization of the nitroxide when compared to that of the N-terminally modified Proxyl-Gtalpha(340-350) analogue. A doubly labeled Proxyl/TOAC(348)-Gtalpha(340-350) was examined by DEER spectrocopy to determine the distribution of distances between the two nitroxides in the peptides when in solution and when bound to R*. TOAC and Proxyl spin-labels in this GPCR-G-protein alpha-peptide system provide unique biophysical probes that can be used to explore the structure and conformational changes at the rhodopsin-G-protein interface.

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.

Resolving Conformational and Rotameric Exchange in Spin-Labeled Proteins Using Saturation Recovery EPR

The function of many proteins involves equilibria between conformational substates, and to elucidate mechanisms of function it is essential to have experimental tools to detect the presence of conformational substates and to determine the time scale of exchange between them. Site-directed spin labeling (SDSL) has the potential to serve this purpose. In proteins containing a nitroxide side chain (R1), multicomponent electron paramagnetic resonance (EPR) spectra can arise either from equilibria involving different conformational substates or rotamers of R1. To employ SDSL to uniquely identify conformational equilibria, it is thus essential to distinguish between these origins of multicomponent spectra. Here we show that this is possible based on the time scale for exchange of the nitroxide between distinct environments that give rise to multicomponent EPR spectra; rotamer exchange for R1 lies in the ≈0.1-1 μs range, while conformational exchange is at least an order of magnitude slower. The time scales of exchange events are determined by saturation recovery EPR, and in favorable cases, the exchange rate constants between substates with lifetimes of approximately 1-70 μs can be estimated by the approach.

Site-directed spin labeling of a genetically encoded unnatural amino acid

The traditional site-directed spin labeling (SDSL) method, which utilizes cysteine residues and sulfhydryl-reactive nitroxide reagents, can be challenging for proteins that contain functionally important native cysteine residues or disulfide bonds. To make SDSL amenable to any protein, we introduce an orthogonal labeling strategy, i.e., one that does not rely on any of the functional groups found in the common 20 amino acids. In this method, the genetically encoded unnatural amino acid p-acetyl-L-phenylalanine (p-AcPhe) is reacted with a hydroxylamine reagent to generate a nitroxide side chain (K1). The utility of this scheme was demonstrated with seven mutants of T4 lysozyme, each containing a single p-AcPhe at a solvent-exposed helix site; the mutants were expressed in amounts qualitatively similar to the wild-type protein. In general, the EPR spectra of the resulting K1 mutants reflect higher nitroxide mobilities than the spectra of analogous mutants containing the more constrained disulfide-linked side chain (R1) commonly used in SDSL. Despite this increased flexibility, site dependence of the EPR spectra suggests that K1 will be a useful sensor of local structure and of conformational changes in solution. Distance measurements between pairs of K1 residues using double electron electron resonance (DEER) spectroscopy indicate that K1 will also be useful for distance mapping.

Osmolyte perturbation reveals conformational equilibria in spin-labeled proteins

Recent evidence suggests that proteins at equilibrium can exist in a manifold of conformational substates, and that these substates play important roles in protein function. Therefore, there is great interest in identifying regions in proteins that are in conformational exchange. Electron paramagnetic resonance spectra of spin-labeled proteins containing the nitroxide side chain (R1) often consist of two (or more) components that may arise from slow exchange between conformational substates (lifetimes > 100 ns). However, crystal structures of proteins containing R1 have shown that multicomponent spectra can also arise from equilibria between rotamers of the side chain itself. In this report, it is shown that these scenarios can be distinguished by the response of the system to solvent perturbation with stabilizing osmolytes such as sucrose. Thus, site-directed spin labeling (SDSL) emerges as a new tool to explore slow conformational exchange in proteins of arbitrary size, including membrane proteins in a native-like environment. Moreover, equilibrium between substates with even modest differences in conformation is revealed, and the simplicity of the method makes it suitable for facile screening of multiple proteins. Together with previously developed strategies for monitoring picosecond to millisecond backbone dynamics, the results presented here expand the timescale over which SDSL can be used to explore protein flexibility.