Calculation of electron paramagnetic resonance spectra from Brownian dynamics trajectories: application to nitroxide side chains in proteins.

We present a method to simulate electron paramagnetic resonance spectra of spin-labeled proteins that explicitly includes the protein structure in the vicinity of the attached spin label. The method is applied to a spin-labeled polyleucine alpha-helix trimer. From short (6 ns) stochastic dynamics simulations of this trimer, an effective potential energy function is calculated. Interaction with secondary and tertiary structures determine the reorientational motion of the spin label side chains. After reduction to a single particle problem, long stochastic dynamic trajectories (700 ns) of the spin label side-chain reorientation are calculated from which the Lamor frequency trajectory and subsequently the electron paramagnetic resonance spectrum is determined. The simulated spectra agree well with experimental electron paramagnetic resonance spectra of bacteriorhodopsin mutants with spin labels in similar secondary and tertiary environments as in the polyleucine.     

Simulation of nitroxide electron paramagnetic resonance spectra from brownian trajectories and molecular dynamics simulations

A simulated continuous wave electron paramagnetic resonance spectrum of a nitroxide spin label can be obtained from the Fourier transform of a free induction decay. It has been previously shown that the free induction decay can be calculated by solving the time-dependent stochastic Liouville equation for a set of Brownian trajectories defining the rotational dynamics of the label. In this work, a quaternion-based Monte Carlo algorithm has been developed to generate Brownian trajectories describing the global rotational diffusion of a spin-labeled protein. Also, molecular dynamics simulations of two spin-labeled mutants of T4 lysozyme, T4L F153R1, and T4L K65R1 have been used to generate trajectories describing the internal dynamics of the protein and the local dynamics of the spin-label side chain. Trajectories from the molecular dynamics simulations combined with trajectories describing the global rotational diffusion of the protein are used to account for all of the dynamics of a spin-labeled protein. Spectra calculated from these combined trajectories correspond well to the experimental spectra for the buried site T4L F153R1 and the helix surface site T4L K65R1. This work provides a framework to further explore the modeling of the dynamics of the spin-label side chain in the wide variety of labeling environments encountered in site-directed spin labeling studies.     

Toward the fourth dimension of membrane protein structure: insight into dynamics from spin-labeling EPR spectroscopy

Trapping membrane proteins in the confines of a crystal lattice obscures dynamic modes essential for interconversion between multiple conformations in the functional cycle. Moreover, lattice forces could conspire with detergent solubilization to stabilize a minor conformer in an ensemble thus confounding mechanistic interpretation. Spin labeling in conjunction with electron paramagnetic resonance (EPR) spectroscopy offers an exquisite window into membrane protein dynamics in the native-like environment of a lipid bilayer. Systematic application of spin labeling and EPR identifies sequence-specific secondary structures, defines their topology and their packing in the tertiary fold. Long range distance measurements (60 ?-80 ?) between pairs of spin labels enable quantitative analysis of equilibrium dynamics and triggered conformational changes. This review highlights the contribution of spin labeling to bridging structure and mechanism. Efforts to develop methods for determining structures from EPR restraints and to increase sensitivity and throughput promise to expand spin labeling applications in membrane protein structural biology.     

Probing iso-1-cytochrome c structure by site-directed spin labeling and electron paramagnetic resonance techniques

A cysteine-specific methanethiosulfonate spin label was introduced into yeast iso-1-cytochrome c at three different positions. The modified forms of cytochrome c included: the wild-type protein labeled at naturally occurring C102, and two mutated proteins, S47C and L85C, labeled at positions 47 and 85, respectively (both S47C and L85C derived from the protein in which C102 had been replaced by threonine). All three spin-labeled protein derivatives were characterized using electron paramagnetic resonance (EPR) techniques. The continuous wave (CW) EPR spectrum of spin label attached to L85C differed from those recorded for spin label attached to C102 or S47C, indicating that spin label at position 85 was more immobilized and exhibited more complex tumbling than spin label at two other positions. The temperature dependence of the CW EPR spectra and CW EPR power saturation revealed further differences of spin-labeled L85C. The results were discussed in terms of application of the site-directed spin labeling technique in probing the local dynamic structure of iso-1-cytochrome c.     

Recent advances and applications of site-directed spin labeling

Site-directed spin labeling has become a popular biophysical tool for the characterization of protein structure, dynamics and conformational change. This method is well suited and widely used to study small soluble proteins, membrane proteins and large protein complexes. Recent advances in site-directed spin labeling methodology have occurred in two areas. The first involves an understanding of the conformations and local dynamics of the spin-labeled sidechain, including the features of proteins that influence electron paramagnetic resonance lineshape. The second advance is the application of pulse techniques to determine long-range distances and distance distributions in proteins. During the past two years, these technical developments have been used to address several important problems concerning the molecular function of proteins.     

Identifying conformational changes with site-directed spin labeling

Site-direct spin labeling combined with electron paramagnetic resonance (EPR) spectroscopy is a powerful tool for detecting structural changes in proteins. This review provides examples that illustrate strategies for interpreting the data in terms of specific rearrangements in secondary and tertiary structure. The changes in the mobility and solvent accessibility of the spin label side chains, and in the distances between spin labels, report (i) rigid body motions of alpha-helices and beta-strands (ii) relative movements of domains and (iii) changes in secondary structure. Such events can be monitored in the millisecond time-scale, making it possible to follow structural changes during function. There is no upper limit to the size of proteins that can be investigated, and only 50-100 picomoles of protein are required. These features make site-directed spin labeling an attractive approach for the study of structure and dynamics in a wide range of systems.     

NMR resonance assignment of selectively labeled proteins by the use of paramagnetic ligands

Selective isotopic labeling of larger proteins greatly simplifies protein NMR spectra and reduces signal overlap, but selectively labeled proteins cannot be easily assigned since the sequential assignment method is not applicable. Here we describe a strategy for resonance assignment in selectively labeled proteins. Our approach involves a spin-labeled analog of a ligand of which the three-dimensional structure in complex with the target protein is known. Other methods for introduction of the spin label are possible. The paramagnetic center causes faster relaxation of all neighboring nuclei in a distance-dependent manner. Measurement of this effect allows to deduce distances between isotopically labeled residues and the paramagnetic center which can be used for resonance assignment. The method is demonstrated for the catalytic domain of Abl kinase in complex with the inhibitor, STI571.     

Integrated analysis of the conformation of a protein-linked spin label by crystallography,EPR and NMR spectroscopy

Long-range structural information derived from paramagnetic relaxation enhancement observed in the presence of a paramagnetic nitroxide radical is highly useful for structural characterization of globular, modular and intrinsically disordered proteins, as well as protein–protein and protein-DNA complexes. Here we characterized the conformation of a spin-label attached to the homodimeric protein CylR2 using a combination of X-ray crystallography, electron paramagnetic resonance (EPR) and NMR spectroscopy. Close agreement was found between the conformation of the spin label observed in the crystal structure with interspin distances measured by EPR and signal broadening in NMR spectra, suggesting that the conformation seen in the crystal structure is also preferred in solution. In contrast, conformations of the spin label observed in crystal structures of T4 lysozyme are not in agreement with the paramagnetic relaxation enhancement observed for spin-labeled CylR2 in solution. Our data demonstrate that accurate positioning of the paramagnetic center is essential for high-resolution structure determination.     

Association of alpha-synuclein and mutants with lipid membranes: spin-label ESR and polarized IR

Alpha-synuclein is a presynaptic protein, the A53T and A30P mutants of which are linked independently to early-onset familial Parkinson's disease. The association of wild-type alpha-synuclein with lipid membranes was characterized previously by electron spin resonance (ESR) spectroscopy with spin-labeled lipids [Ramakrishnan, M., Jensen, P. H., and Marsh, D. (2003) Biochemistry 42, 12919-12926]. Here, we study the interaction of the A53T and A30P alpha-synuclein mutants and a truncated form that lacks the acidic C-terminal domain with phosphatidylglycerol bilayer membranes, using anionic phospholipid spin labels. The strength of the interaction with phosphatidylglycerol membranes lies in the order: wild type approximately truncated > A53T > A30P > fibrils approximately 0, and only the truncated form interacts with phosphatidylcholine membranes. The selectivity of the interaction of the mutant alpha-synucleins with different spin-labeled lipid species is reduced considerably, relative to the wild-type protein, whereas that of the truncated protein is increased. Polarized infrared (IR) spectroscopy is used to study the interactions of the wild-type and truncated proteins with aligned lipid membranes and additionally to characterize the fibrillar form. Wild-type alpha-synuclein is natively unfolded in solution and acquires secondary structure upon binding to membranes containing phosphatidylglycerol. Up to 30-40% of the amide I band intensity of the membrane-bound wild-type and truncated proteins is attributable to beta-sheet structure, at the surface densities used for IR spectroscopy. The remainder is alpha-helix and residual unordered structure. Fibrillar alpha-synuclein contains 62% antiparallel beta-sheet and is oriented on the substrate surface but does not interact with deposited lipid membranes. The beta-sheet secondary-structural elements of the wild-type and truncated proteins are partially oriented on the surface of membranes with which they interact.     

Low microwave-amplitude ESR spectroscopy: measuring spin-relaxation interactions of moderately immobilized spin labels in proteins

Electron spin resonance (ESR) spectroscopy in combination with site-directed spin labeling (SDSL) is a powerful tool for determining protein structure, dynamics and interactions. We report here a method for determining interactions between spin labels and paramagnetic relaxation agents, which is performed under subsaturating conditions. The low microwave-field amplitude employed (h(1)<0.36 G) only requires standard, commercially available ESR equipment. The effect of relaxation enhancement on the spin-spin-relaxation time, T(2e), is measured by this method, and compared to classical progressive power saturation performed on a free spin label, (1-oxyl-2,2,5,5-tetramethyl-Delta(3)-pyrroline-3-methyl)methanethiosulfonate (MTSL), and a spin-labeled protein (Thermomyces lanuginosa lipase, TLL-I252C), employing the water-soluble relaxation agent chromium(III) oxalate (Crox) in concentrations between 0-10 mM. The low-amplitude theory showed excellent agreement with that of classical power saturation in quantifying Crox-induced relaxation enhancement. Low-amplitude measurements were then performed using a standard resonator, with Crox, on 11 spin-labeled TLL mutants displaying rotational correlation times in the motional narrowing regime. All spin-labeled proteins exhibited significant changes in T(2e). We postulate that this novel method is especially suitable for studying moderately immobilized spin labels, such as those positioned at exposed sites in a protein. This method should prove useful for research groups with access to any ESR instrumentation.     

Membrane orientation and position of the C2 domain from cPLA2 by site-directed spin labeling

The C2 domain is a ubiquitous Ca(2+)-binding motif that triggers the membrane docking of many key signaling proteins during intracellular Ca(2+) signals. Site-directed spin labeling was carried out on the C2 domain of cytosolic phospholipase A(2) in order to determine the depth of penetration and orientation of the domain at the membrane interface. Membrane depth parameters, Phi, were obtained by EPR spectroscopy for a series of selectively spin-labeled C2 domain cysteine mutants, and for spin-labeled lipids and spin-labeled bacteriorhodopsin cysteine mutants. Values of Phi were combined with several other constraints, including the solution NMR structure, to generate a model for the position of the C2 domain at the membrane interface. This modeling yielded an empirical expression for Phi, which for the first time defines its behavior from the bulk aqueous phase to the center of the lipid bilayer. In this model, the backbones of both the first and third Ca(2+)-binding loops are inserted approximately 10 A into the bilayer, with residues inserted as deep as 15 A. The backbone of the second Ca(2+)-binding loop is positioned near the lipid phosphate, and the two beta-sheets of the C2 domain are oriented so that the individual strands make angles of 30-45 degrees with respect to the bilayer surface. Upon membrane docking, spin labels in the Ca(2+)-binding loops exhibit decreases in local motion, suggesting either changes in tertiary contacts due to protein conformational changes and/or interactions with lipid.     

Motion of spin-labeled side chains in T4 lysozyme: effect of side chain structure

Previous studies have shown that the mobility of nitroxide side chains in a protein, inferred from the electron paramagnetic resonance (EPR) spectra, can be used to classify particular sites as helix surface sites, tertiary contact sites, buried sites, or loop sites. In addition, the sequence dependence of mobility can identify regular secondary structure. However, in the most widely used side chain, an apparent interaction of the nitroxide ring with the protein at some helix surface sites gives rise to EPR spectra degenerate with those at tertiary contact sites. In the present study, we use selected sites in T4 lysozyme to evaluate novel nitroxide side chains designed to resolve this degeneracy. The results indicate that the reagent 3-(methanesulfonylthiomethyl)-2,2, 5,5-tetramethylpyrrolidin-1-yloxy reacts with cysteine to give a nitroxide side chain that has a high contrast in mobility between helix surface and tertiary contact sites, effectively resolving the degeneracy. The reagent 3-(iodomercuriomethyl)-2,2,5,5-tetramethyl-2, 5-dihydro-1H-pyrrol-1-yloxy reacts with cysteine to provide a mercury-linked nitroxide that also shows reduced interaction with the protein at most helix surface sites. Thus, these new side chains may be the preferred choices for structure determination using site-directed spin labeling.     

Comparison of electron paramagnetic resonance methods to determine distances between spin labels on human carbonic anhydrase II

Four doubly spin-labeled variants of human carbonic anhydrase II and corresponding singly labeled variants were prepared by site-directed spin labeling. The distances between the spin labels were obtained from continuous-wave electron paramagnetic resonance spectra by analysis of the relative intensity of the half-field transition, Fourier deconvolution of line-shape broadening, and computer simulation of line-shape changes. Distances also were determined by four-pulse double electron-electron resonance. For each variant, at least two methods were applicable and reasonable agreement between methods was obtained. Distances ranged from 7 to 24 A. The doubly spin-labeled samples contained some singly labeled protein due to incomplete labeling. The sensitivity of each of the distance determination methods to the non-interacting component was compared.     

Monomeric and dimeric conformation of the vinculin tail five-helix bundle in solution studied by EPR spectroscopy

The cytoskeletal adaptor protein vinculin plays an important role in the control of cell adhesion and migration, linking the actin cytoskeleton to adhesion receptor complexes in cell adhesion sites. The conformation of the vinculin tail dimer, which is crucial for protein function, was analyzed using site-directed spin labeling in electron paramagnetic resonance spectroscopy. Interspin distances for a set of six singly and four doubly spin-labeled mutants of the tail domain of vinculin were determined and used as constraints for modeling of the vinculin tail dimer. A comparison of the results obtained by molecular dynamic simulations and a rotamer library approach reveals that the crystal structure of the vinculin tail monomer is essentially preserved in aqueous solution. The orientation of monomers within the dimer observed previously by x-ray crystallography agrees with the solution electron paramagnetic resonance data. Furthermore, the distance between positions 1033 is shown to increase by >3 nm upon interaction of the vinculin tail domain with F-actin.     

Effect of membrane protein on lipid bilayer structure: a spin-label electron spin resonance study of vesicular stomatitis virus.

Spin-label electron spin resonance (ESR) methods have been used to study the structure of the envelope of vesicular stomatitis virus (VSV). The data indicate that the lipid is organized in a bilayer structure. Proteolytic digestion of the glycoproteins which are the spike-like projections on the outer surface of the virus particle increases the fluidity of the lipid bilayer. Since the lipid composition of the virion reflects the composition of the host plasma membrane and the protein composition is determined by the viral genome, VSV was grown in both MDBK and BHK21-F cells to determine the effect of a change in lipid composition on the structure of the lipid bilayer of VSV. The lipid bilayer of the virion was found to be more rigid when derived from MDBK cells than from BHK21-F cells. Studies comparing spin-labeled intact cells and cell membrane fractions suggest that upon labeling the whole cell the spin label probes the plasma membrane. Comparison of spin-labeled VSV particles and their host cells indicates that the lipid bilayer of the plasma membrane is considerably more fluid than that of the virion. These results are discussed in terms of the effect of membrane-associated protein on the structure of the lipid bilayer.     

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.     

Oxygen permeation profile in lipid membranes: comparison with transmembrane polarity profile

Permeation of oxygen into membranes is relevant not only to physiological function, but also to depth determinations in membranes by site-directed spin labeling. Spin-lattice (T(1)) relaxation enhancements by air or molecular oxygen were determined for phosphatidylcholines spin labeled at positions (n = 4-14, 16) of the sn-2 chain in fluid membranes of dimyristoyl phosphatidylcholine, by using nonlinear continuous-wave electron paramagnetic resonance (EPR). Both progressive saturation and out-of-phase continuous-wave EPR measurements yield similar oxygen permeation profiles. With pure oxygen, the T(2)-relaxation enhancements determined from homogeneous linewidths of the linear EPR spectra are equal to the T(1)-relaxation enhancements determined by nonlinear EPR. This confirms that both relaxation enhancements occur by Heisenberg exchange, which requires direct contact between oxygen and spin label. Oxygen concentrates in the hydrophobic interior of phospholipid bilayer membranes with a sigmoidal permeation profile that is the inverse of the polarity profile established earlier for these spin-labeled lipids. The shape of the oxygen permeation profile in fluid lipid membranes is controlled partly by the penetration of water, via the transmembrane polarity profile. At the protein interface of the KcsA ion channel, the oxygen profile is more diffuse than that in fluid lipid bilayers.     

Characterization of the solution structure of human serum albumin loaded with a metal porphyrin and fatty acids

The structure of human serum albumin loaded with a metal porphyrin and fatty acids in solution is characterized by orientation-selective double electron-electron resonance (DEER) spectroscopy. Human serum albumin, spin-labeled fatty acids, and Cu(II) protoporphyrin IX—a hemin analog—form a fully self-assembled system that allows obtaining distances and mutual orientations between the paramagnetic guest molecules. We report a simplified analysis for the orientation-selective DEER data which can be applied when the orientation selection of one spin in the spin pair dominates the orientation selection of the other spin. The dipolar spectra reveal a dominant distance of 3.85 nm and a dominant orientation of the spin-spin vectors between Cu(II) protoporphyrin IX and 16-doxyl stearic acid, the electron paramagnetic resonance reporter group of the latter being located near the entry points to the fatty acid binding sites. This observation is in contrast to crystallographic data that suggest an asymmetric distribution of the entry points in the protein and hence the occurrence of various distances. In conjunction with the findings of a recent DEER study, the obtained data are indicative of a symmetric distribution of the binding site entries on the protein's surface. The overall anisotropic shape of the protein is reflected by one spin-spin vector orientation dominating the DEER data.     

Membrane hydrocarbon thickness modulates the dynamics of a membrane transport protein

Nitroxide spin labels were incorporated into selected sites within the β-barrel of the bacterial outer-membrane transport protein BtuB by site-directed mutagenesis, followed by chemical modification with a methanethiosufonate spin label. The electron paramagnetic resonance lineshapes of the spin-labeled side chain (R1) from these sites are highly variable, and have spectral parameters that reflect secondary structure and local steric constraints. In addition, these lineshape parameters correlate with crystallographic structure factors for Cα carbons, suggesting that the motion of the spin label is modulated by both the local modes of motion of the spin label and the local dynamics of the protein backbone. Experiments performed as a function of lipid composition and sample temperature indicate that nitroxide spin labels on the exterior surface of BtuB, which face the membrane hydrocarbon, are not strongly influenced by the phase state of the bulk lipids. However, these spectra are modulated by membrane hydrocarbon thickness. Specifically, the values of the scaled mobility parameter for the R1 lineshapes are inversely proportional to the hydrocarbon thickness. These data suggest that protein dynamics and structure in BtuB are directly coupled to membrane hydrophobic thickness.     

Pulsed EPR Determination of Water Accessibility to Spin-Labeled Amino Acid Residues in LHCIIb

Membrane proteins reside in a structured environment in which some of their residues are accessible to water, some are in contact with alkyl chains of lipid molecules, and some are buried in the protein. Water accessibility of residues may change during folding or function-related structural dynamics. Several techniques based on the combination of pulsed electron paramagnetic resonance (EPR) with site-directed spin labeling can be used to quantify such water accessibility. Accessibility parameters for different residues in major plant light-harvesting complex IIb are determined by electron spin echo envelope modulation spectroscopy in the presence of deuterated water, deuterium contrast in transversal relaxation rates, analysis of longitudinal relaxation rates, and line shape analysis of electron-spin-echo-detected EPR spectra as well as by the conventional techniques of measuring the maximum hyperfine splitting and progressive saturation in continuous-wave EPR. Systematic comparison of these parameters allows for a more detailed characterization of the environment of the spin-labeled residues. These techniques are applicable independently of protein size and require ∼10-20 nmol of singly spin-labeled protein per sample. For a residue close to the N-terminus, in a domain unresolved in the existing x-ray structures of light-harvesting complex IIb, all methods indicate high water accessibility.