Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

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.     

Utilization of 13C-labeled amino acids to probe the α-helical local secondary structure of a membrane peptide using electron spin echo envelope modulation (ESEEM) spectroscopy

Electron spin echo envelope modulation (ESEEM) spectroscopy in combination with site-directed spin labeling (SDSL) has been established as a valuable biophysical technique to provide site-specific local secondary structure of membrane proteins. This pulsed electron paramagnetic resonance (EPR) method can successfully distinguish between α-helices, β-sheets, and 3₁₀-helices by strategically using ²H-labeled amino acids and SDSL. In this study, we have explored the use of ¹³C-labeled residues as the NMR active nuclei for this approach for the first time. ¹³C-labeled d₅-valine (Val) or ¹³C-labeled d₆-leucine (Leu) were substituted at a specific Val or Leu residue (i), and a nitroxide spin label was positioned 2 or 3 residues away (denoted i-2 and i-3) on the acetylcholine receptor M2δ (AChR M2δ) in a lipid bilayer. The ¹³C ESEEM peaks in the FT frequency domain data were observed for the i-3 samples, and no ¹³C peaks were observed in the i-2 samples. The resulting spectra were indicative of the α-helical local secondary structure of AChR M2δ in bicelles. This study provides more versatility and alternative options when using this ESEEM approach to study the more challenging recombinant membrane protein secondary structures.     

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.     

A new spin on protein dynamics

Site-directed spin labeling is a general method for investigating structure and conformational switching in soluble and membrane proteins. It will also be an important tool for exploring protein backbone dynamics. A semi-empirical analysis of nitroxide sidechain dynamics in spin-labeled proteins reveals contributions from fluctuations in backbone dihedral angles and rigid-body (collective) motions of alpha helices. Quantitative analysis of sidechain dynamics is sometimes possible, and contributions from backbone modes can be expressed in terms of relative order parameters and rates. Dynamic sequences identified by site-directed spin labeling correlate with functional domains, and so nitroxide scanning could provide an efficient strategy for identifying such domains in high-molecular weight proteins, supramolecular complexes and membrane proteins.     

Site-directed spin labeling measurements of nanometer distances in nucleic acids using a sequence-independent nitroxide probe

In site-directed spin labeling (SDSL), local structural and dynamic information is obtained via electron paramagnetic resonance (EPR) spectroscopy of a stable nitroxide radical attached site-specifically to a macromolecule. Analysis of electron spin dipolar interactions between pairs of nitroxides yields the inter-nitroxide distance, which provides quantitative structural information. The development of pulse EPR methods has enabled such distance measurements up to 70Å in bio-molecules, thus opening up the possibility of SDSL global structural mapping. This study evaluates SDSL distance measurement using a nitroxide (designated as R5) that can be attached, in an efficient and cost-effective manner, to a phosphorothioate backbone position at arbitrary DNA or RNA sequences. R5 pairs were attached to selected positions of a dodecamer DNA duplex with a known NMR structure, and eight distances, ranging from 20 to 40Å, were measured using double electron-electron resonance (DEER). The measured distances correlated strongly (R² = 0.98) with the predicted values calculated based on a search of sterically allowable R5 conformations in the NMR structure, thus demonstrating accurate distance measurements using R5. Furthermore, distance measurement in a 42 kD DNA was demonstrated. The results establish R5 as a sequence-independent probe for global structural mapping of DNA and DNA–protein complexes.     

Solutes alter the conformation of the ligand binding loops in outer membrane transporters

The binding and recognition of ligands by bacterial outer membrane transport proteins is mediated in part by interactions made through their extracellular loops. Here, site-directed spin labeling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy were used to examine the effect of stabilizing solutes on the extracellular loops in BtuB, the vitamin B12 transporter, and FecA, the ferric citrate transporter. EPR spectra from the extracellular loops of FecA and BtuB arise from dynamic backbone segments, and distance measurements made by double electron-electron resonance indicate that the second extracellular loop in BtuB samples a wide range of conformations. These conformations are dramatically restricted upon substrate binding. In addition, the EPR spectra from nitroxide labels attached to the extracellular loops in BtuB and FecA are highly sensitive to solutes, and at every site examined the motion of the label is significantly reduced in the presence of stabilizing osmolytes, such as polyethylene glycols. For the second extracellular loop in BtuB, the solute-induced structural changes are small, but they are sufficient to bring spin-labeled side chains into tertiary contact with other portions of the protein. The spectroscopic changes seen by SDSL suggest that high concentrations of stabilizing solutes, such as those used to generate membrane protein crystals, result in a more compact and ordered state of the protein than is seen under more physiological conditions.     

High-Resolution NMR Reveals Secondary Structure and Folding of Amino Acid Transporter from Outer Chloroplast Membrane

Solving high-resolution structures for membrane proteins continues to be a daunting challenge in the structural biology community. In this study we report our high-resolution NMR results for a transmembrane protein, outer envelope protein of molar mass 16 kDa (OEP16), an amino acid transporter from the outer membrane of chloroplasts. Three-dimensional, high-resolution NMR experiments on the ¹³C, ¹⁵N, ²H-triply-labeled protein were used to assign protein backbone resonances and to obtain secondary structure information. The results yield over 95% assignment of N, H_N, CO, C_α, and C_β chemical shifts, which is essential for obtaining a high resolution structure from NMR data. Chemical shift analysis from the assignment data reveals experimental evidence for the first time on the location of the secondary structure elements on a per residue basis. In addition T _1Z and T₂ relaxation experiments were performed in order to better understand the protein dynamics. Arginine titration experiments yield an insight into the amino acid residues responsible for protein transporter function. The results provide the necessary basis for high-resolution structural determination of this important plant membrane protein.     

Structure and dynamics of a membrane protein in micelles from three solution NMR experiments

Three solution NMR experiments on a uniformly ¹⁵N labeled membrane protein in micelles provide sufficient information to describe the structure, topology, and dynamics of its helices, as well as additional information that characterizes the principal features of residues in terminal and inter-helical loop regions. The backbone amide resonances are assigned with an HMQC-NOESY experiment and the backbone dynamics are characterized by a ¹H-¹⁵N heteronuclear NOE experiment, which clearly distinguishes between the structured helical residues and the more mobile residues in the terminal and interhelical loop regions of the protein. The structure and topology of the helices are described by Dipolar waves and PISA wheels derived from experimental measurements of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). The results show that the membrane-bound form of Pf1 coat protein has a 20-residue trans-membrane hydrophobic helix with an orientation that differs by about 90° from that of an 8-residue amphipathic helix. This combination of three-experiments that yields Dipolar waves and PISA wheels has the potential to contribute to high-throughput structural characterizations of membrane proteins.     

Quantification of protein backbone hydrogen-deuterium exchange rates by solid state NMR spectroscopy

We present the quantification of backbone amide hydrogen-deuterium exchange rates (HDX) for immobilized proteins. The experiments make use of the deuterium isotope effect on the amide nitrogen chemical shift, as well as on proton dilution by deuteration. We find that backbone amides in the microcrystalline α-spectrin SH3 domain exchange rather slowly with the solvent (with exchange rates negligible within the individual ¹⁵N–T ₁ timescales). We observed chemical exchange for 6 residues with HDX exchange rates in the range from 0.2 to 5 s⁻¹. Backbone amide ¹⁵N longitudinal relaxation times that we determined previously are not significantly affected for most residues, yielding no systematic artifacts upon quantification of backbone dynamics (Chevelkov et al. 2008b). Significant exchange was observed for the backbone amides of R21, S36 and K60, as well as for the sidechain amides of N38, N35 and for W41ε. These residues could not be fit in our previous motional analysis, demonstrating that amide proton chemical exchange needs to be considered in the analysis of protein dynamics in the solid-state, in case D₂O is employed as a solvent for sample preparation. Due to the intrinsically long ¹⁵N relaxation times in the solid-state, the approach proposed here can expand the range of accessible HDX rates in the intermediate regime that is not accessible so far with exchange quench and MEXICO type experiments.     

Antioxidant reactivity toward nitroxide probes anchored into human serum albumin. A new model for studying antioxidant repairing capacity of protein radicals

A new strategy to evaluate accessibility of antioxidants to radical proteins has been developed using nitroxide prefluorescent probes anchored into human serum albumin (HSA). Binding association constants for the nitroxide probes C₃₄₃T and QT with HSA were 5 × 10⁴ and 9 × 10⁴ M⁻¹, respectively. Rate constants for the nitroxide reduction by antioxidants in HSA were determined finding k_HSA/k_buffer ratio of 0.8, 1.9, and 0.075 for ascorbic acid, Trolox, and caffeic acid, respectively, for the nitroxide C₃₄₃T reduction.     

Solid-state H and N NMR studies of side-chain and backbone dynamics of phospholamban in lipid bilayers: Investigation of the N27A mutation

Phospholamban (PLB) is an integral membrane protein regulating Ca²⁺ transport through inhibitory interaction with sarco(endo)plasmic reticulum calcium ATPase (SERCA). The Asn27 to Ala (N27A) mutation of PLB has been shown to function as a superinhibitor of the affinity of SERCA for Ca²⁺ and of cardiac contractility in vivo. The effects of this N27A mutation on the side-chain and backbone dynamics of PLB were investigated with ²H and ¹⁵N solid-state NMR spectroscopy in phospholipid multilamellar vesicles (MLVs). ²H and ¹⁵N NMR spectra indicate that the N27A mutation does not significantly change the side-chain or backbone dynamics of the transmembrane and cytoplasmic domains when compared to wild-type PLB. However, dynamic changes are observed for the hinge region, in which greater mobility is observed for the CD₃-labeled Ala24 N27A-PLB. The increased dynamics in the hinge region of PLB upon N27A mutation may allow the cytoplasmic helix to more easily interact with the Ca²⁺-ATPase; thus, showing increased inhibition of Ca²⁺-ATPase.     

Nanometer distance measurements in RNA using site-directed spin labeling

The method of site-directed spin labeling (SDSL) utilizes a stable nitroxide radical to obtain structural and dynamic information on biomolecules. Measuring dipolar interactions between pairs of nitroxides yields internitroxide distances, from which quantitative structural information can be derived. This study evaluates SDSL distance measurements in RNA using a nitroxide probe, designated as R5, which is attached in an efficient and cost-effective manner to backbone phosphorothioate sites that are chemically substituted in arbitrary sequences. It is shown that R5 does not perturb the global structure of the A-form RNA helix. Six sets of internitroxide distances, ranging from 20 to 50 A, were measured on an RNA duplex with a known X-ray crystal structure. The measured distances strongly correlate (R(2) = 0.97) with those predicted using an efficient algorithm for determining the expected internitroxide distances from the parent RNA structure. The results enable future studies of global RNA structures for which high-resolution structural data are absent.     

Combined solid state and solution NMR studies of α,ɛ-<Superscript>15</Superscript>N labeled bovine rhodopsin

Rhodopsin is the visual pigment of the vertebrate rod photoreceptor cell and is the only member of the G protein coupled receptor family for which a crystal structure is available. Towards the study of dynamics in rhodopsin, we report NMR-spectroscopic investigations of α,ɛ-¹⁵N-tryptophan labeled rhodopsin in detergent micelles and reconstituted in phospholipids. Using a combination of solid state ¹³C,¹⁵N-REDOR and HETCOR experiments of all possible ¹³C′ _i-1 carbonyl/¹⁵N _i -tryptophan isotope labeled amide pairs, and H/D exchange ¹H,¹⁵N-HSQC experiments conducted in solution, we assigned chemical shifts to all five rhodopsin tryptophan backbone ¹⁵N nuclei and partially to their bound protons. ¹H,¹⁵N chemical shift assignment was achieved for indole side chains of Trp35^1.30 and Trp175^4.65. ¹⁵N chemical shifts were found to be similar when comparing those obtained in the native like reconstituted lipid environment and those obtained in detergent micelles for all tryptophans except Trp175^4.65 at the membrane interface. The results suggest that the integrated solution and solid state NMR approach presented provides highly complementary information in the study of structure and dynamics of large membrane proteins like rhodopsin.     

Optimization of amino acid type-specific <Superscript>13</Superscript>C and <Superscript>15</Superscript>N labeling for the backbone assignment of membrane proteins by solution- and solid-state NMR with the UPLABEL algorithm

We present a computational method for finding optimal labeling patterns for the backbone assignment of membrane proteins and other large proteins that cannot be assigned by conventional strategies. Following the approach of Kainosho and Tsuji (Biochemistry 21:6273–6279 (1982)), types of amino acids are labeled with ¹³C or/and ¹⁵N such that cross peaks between ¹³CO(i – 1) and ¹⁵NH(i) result only for pairs of sequentially adjacent amino acids of which the first is labeled with ¹³C and the second with ¹⁵N. In this way, unambiguous sequence-specific assignments can be obtained for unique pairs of amino acids that occur exactly once in the sequence of the protein. To be practical, it is crucial to limit the number of differently labeled protein samples that have to be prepared while obtaining an optimal extent of labeled unique amino acid pairs. Our computer algorithm UPLABEL for optimal unique pair labeling, implemented in the program CYANA and in a standalone program, and also available through a web portal, uses combinatorial optimization to find for a given amino acid sequence labeling patterns that maximize the number of unique pair assignments with a minimal number of differently labeled protein samples. Various auxiliary conditions, including labeled amino acid availability and price, previously known partial assignments, and sequence regions of particular interest can be taken into account when determining optimal amino acid type-specific labeling patterns. The method is illustrated for the assignment of the human G-protein coupled receptor bradykinin B2 (B₂R) and applied as a starting point for the backbone assignment of the membrane protein proteorhodopsin.     

1H, 13C,and 15N backbone resonance assignments of the 37 kDa voltage-gated Ca2+ channel β4 subunit core SH3-GK domains

The β subunit of the voltage-gated Ca²⁺ channel (α₁, α₂δ, and β subunits) is a member of the MAGUK family of proteins and plays an essential role in regulating Ca²⁺ channel trafficking and gating. It also serves as a central interaction partner for various Ca²⁺ channel regulatory proteins. We report here the nearly complete ¹H, ¹³C, and ¹⁵N backbone resonance assignments of the 37 kDa core SH3-GK domains of the β4 subunit. This is the first report of solution assignments for β subunits, and as such will lay the foundation for future investigations of interaction site mapping, functional dynamics, and protein complex structure determination.     

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.     

Protein global fold determination using site-directed spin and isotope labeling

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.     

Modeling a spin-labeled fusion peptide in a membrane: implications for the interpretation of EPR experiments

Site-directed spin-labeling and electron paramagnetic resonance are powerful tools for studying structure and conformational dynamics of proteins, especially in membranes. The position of the spin label is used as an indicator of the position of the site to which it is attached. The interpretation of these experiments is based on the assumptions that the spin label does not affect the peptide configuration and that it has a fixed orientation and distance with respect to the protein backbone. Here, the validity of these assumptions is examined through implicit membrane molecular dynamics simulations of the influenza hemagglutinin fusion peptide that has been labeled with methanethiosulfonate spin label. We find that the methanethiosulfonate spin label can occasionally induce peptide orientations that differ from those adopted by the wild-type peptide. Furthermore, the spin-label resides, on average, several Angstroms deeper in the membrane than the corresponding backbone C(alpha)-atom even at sites pointing toward the solvent. The nitroxide spin label exhibits flexibility and adopts various configurations depending on the surrounding residues.     

Quantification of superoxide radical production in thylakoid membrane using cyclic hydroxylamines

Applicability of two lipophilic cyclic hydroxylamines (CHAs), CM-H and TMT-H, and two hydrophilic CHAs, CAT1-H and DCP-H, for detection of superoxide anion radical (O₂^∙−) produced by the thylakoid photosynthetic electron transfer chain (PETC) of higher plants under illumination has been studied. ESR spectrometry was applied for detection of the nitroxide radical originating due to CHAs oxidation by O₂^∙−. CHAs and corresponding nitroxide radicals were shown to be involved in side reactions with PETC which could cause miscalculation of O₂^∙− production rate. Lipophilic CM-H was oxidized by PETC components, reducing the oxidized donor of Photosystem I, P₇₀₀⁺, while at the same concentration another lipophilic CHA, TMT-H, did not reduce P₇₀₀⁺. The nitroxide radical was able to accept electrons from components of the photosynthetic chain. Electrostatic interaction of stable cation CAT1-H with the membrane surface was suggested. Water-soluble superoxide dismutase (SOD) was added in order to suppress the reaction of CHA with O₂^∙− outside the membrane. SOD almost completely inhibited light-induced accumulation of DCP^∙, nitroxide radical derivative of hydrophilic DCP-H, in contrast to TMT^∙ accumulation. Based on the results showing that change in the thylakoid lumen pH and volume had minor effect on TMT^∙ accumulation, the reaction of TMT-H with O₂^∙− in the lumen was excluded. Addition of TMT-H to thylakoid suspension in the presence of SOD resulted in the increase in light-induced O₂ uptake rate, that argued in favor of TMT-H ability to detect O₂^∙− produced within the membrane core. Thus, hydrophilic DCP-H and lipophilic TMT-H were shown to be usable for detection of O₂^∙− produced outside and within thylakoid membranes.