The Synergetic Effects of Combining Structural Biology and EPR Spectroscopy on Membrane Proteins

Protein structures as provided by structural biology such as X-ray crystallography, cryo-electron microscopy and NMR spectroscopy are key elements to understand the function of a protein on the molecular level. Nonetheless, they might be error-prone due to crystallization artifacts or, in particular in case of membrane-imbedded proteins, a mostly artificial environment. In this review, we will introduce different EPR spectroscopy methods as powerful tools to complement and validate structural data gaining insights in the dynamics of proteins and protein complexes such that functional cycles can be derived. We will highlight the use of EPR spectroscopy on membrane-embedded proteins and protein complexes ranging from receptors to secondary active transporters as structural information is still limited in this field and the lipid environment is a particular challenge.


Introduction
Structural biology aims for studying molecular structures and dynamics of biological macromolecules, in particular proteins and nucleic acids, to gain a deeper comprehension of how alterations in their structures affect their function. Here, mostly X-ray crystallography and cryo-electron microscopy provide structural information of the protein of interest at atomic resolution for a broad molecular weight range [1][2][3]. Based on those structures, complex interfaces can be analyzed, structure-based drugs designed and various mechanisms predicted. However, for different targets, in particular membrane proteins, structural biology is still challenging and often, rare conformations are difficult to trap. In addition, for the overall understanding of the proteins' function a dynamic picture is essential. Time-resolved crystallography has been successfully conducted on proteins that are biologically active in the crystalline state with recent developments on X-ray free-electron lasers allowing for femtosecond time resolution on macromolecular structures, in both single crystals and solution. However, more frequently information on rare conformations and protein dynamics is gained by nuclear magnetic resonance (NMR) spectroscopy, Förster resonance energy transfer (FRET) spectroscopy and electron paramagnetic resonance (EPR) spectroscopy as well as combinations thereof. Furthermore, NMR spectroscopy allows for high-resolution structure determination. Based on the nuclear Overhauser effect (NOE) intra-and intermolecular distances below 1 nm can be assigned. In addition, NMR spectroscopy provides information about changes in the protein conformation and protein folding/unfolding in native-like conditions. However, the timescale of this technique is in the nano-to millisecond time range and may impede the study of short-lived intermediates. Further challenges are its limitation by the size and molecular weight of the investigated system and the low sensitivity, which is based on the magnetic moment of the nucleus. In comparison, FRET spectroscopy permits measurements of proteins in solution and is not limited by the size or the molecular weight of the investigated system. As a fluorescence-based method FRET spectroscopy exhibits very high sensitivity and allows for time-resolved experiments especially in single-molecule approaches. A particular challenge of this method is the uncertainty of the orientation factor, one of the key parameters of FRET spectroscopy, allowing rather the accurate determination of distance changes than exact distance values in a distance range of 1-8 nm. In addition, the usage of two different FRET labels complicates the required labeling strategy in particular for multimeric proteins. Simultaneously, the molecular size of these labels may cause structural perturbations and prevent native folding of the labeled protein.
Here, EPR spectroscopy, especially continuous wave (cw) and pulsed EPR, can complement the above discussed methods and beyond that contribute novel structural and dynamic information. Its inherent strength is substantiated in its various applications to proteins, notably to membrane proteins solubilized in detergent or embedded in the lipid bilayer and protein complexes in native-like conditions. Thus, this technique is not limited by the size or the complexity of the investigated system (for further reviews see e.g., [4][5][6]). Cw EPR spectroscopy probes the dynamics of a nitroxide side chain, the polarity of its environment and the structural properties of proteins such as interspin distances in the range of 0.8-1.8 nm in a time regime ranging from picoseconds to seconds. The potential of pulsed EPR spectroscopy, in particular pulsed electron double resonance (PELDOR) or referred to as double electron-electron resonance (DEER), lies in the accurate determination of interspin distances in the range of 2-8 nm including the membrane thickness of approximately 4 nm and covering the most important distance regime necessary for structural investigations.
This review summarizes the versatile aspects of EPR spectroscopy and the synergistic effects resulting from the combination of EPR spectroscopy and structural biology. Furthermore, this paper highlights the latest progress with special emphasis on applications to different classes of membrane proteins.

Site-Directed Spin Labeling
A prerequisite for the utilization of EPR spectroscopy is the introduction of paramagnetic centers. For this purpose several paramagnetic species can be chosen e.g., radical centers, metal ions, Cu 2+ or FeS clusters or nitroxide spin labels [7,8]. The choice of the spin label can be critical as it should be short enough to reduce its intrinsic flexibility but long enough to allow the native folding of the spin-labeled protein. Due to its specificity to cysteine residues, chemical stability and only small structural perturbations the (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl)methanethiosulfonate (MTS) spin label is the most commonly used EPR spin probe in proteins [9]. The spin label is bound to the sulfhydryl group of a cysteine residue by the formation of a disulfide bond. By the formation of the disulfide bond created spin label side chains are frequently abbreviated as R1. Cysteines chosen for labeling are either endogenous amino acids or, more often, via site-directed mutagenesis engineered cysteine residues at defined positions. The latter requires not only the insertion of the cysteine but also the exchange of endogenous cysteine residues to achieve a defined system. Therefore, alanine or serine residues mostly replace endogenous cysteines. An alternative approach for spin labeling is the click-chemistry which utilizes a highly selective carbon-heteroatom bond under mild conditions [10]. Initially, the highly exergonic character of the copper-catalyzed azide-alkyne 1,3-dipolar cycloaddition (CuAAC), also called copper(I)-catalyzed Huisgen-Sharpless-Meldal alkyneazide cycloaddition has been used for spin labeling of DNA on solid-phase by the Sigurdsson group [11] and in solution by the Steinhoff group [12]. Although this approach has been advanced for the successful spin labeling of proteins by means of applying azide-and alkyne-functionalized spin labels, this review focuses on the application of EPR spectroscopy in combination with the commonly used MTS spin label. Additional functional approaches are required to confirm the functionality of the mutated, spin-labeled variant for all spin labeling approaches.

Spin Label Side Chain Mobility
In general, cw EPR measurements performed at different temperatures yield diverse information about the spin label dynamics, the interaction with its local environment, the position of the spin label with respect to the membrane boundaries and interspin distances in the range of 0.8-1.8 nm.
Cw EPR experiments at room temperature (RT) and the adjacent analysis gain information about the motion of the nitroxide side chain and allow for conclusions about the spin label's microenvironment as different environments effect the spin label's motion differently. The time scale of EPR spectroscopy is in the order of nanoseconds, such that it is suitable for a broad range of fundamental biological processes. Extensive studies of the relation of the EPR spectra and the protein structure are reported in literature [13,14]. In Section 3.2 we introduce an EPR study on the activation mechanisms in a functional chimera of the extracellular domain of amine-gated Erwinia chrysanthemi ligand-gated ion channel and the transmembrane domain of Gloeobacter violaceus ligand-gated ion channel. This investigation is based on the application of cw EPR experiments at room temperature.
The intrinsic spin label motion is characterized by its rotational correlation time, which depends on the spin label's intrinsic flexibility and its length. However, the measured effective rotational correlation time of a spin label bound to a protein is a result of the spin label's interaction with the primary, secondary, tertiary and eventually quaternary structure of the protein. Additionally, the protein itself tumbles in solution, which may add to the overall rotational correlation time. Although, the rotational correlation time for proteins, protein complexes or proteins inserted into detergent micelles or liposomes with a total molecular weight above 200 kDa exceeds the sensitivity of the EPR time scale and thus is not resolved in the room temperature cw EPR spectra. Consequently, the term 'mobility' can be mostly interpreted to describe the interaction of the nitroxide side chain with its local environment. Weak interactions and thus a high degree of spin label mobility are defined as only marginal or no restrictions of the motional freedom of the spin probe by its close proximity. This situation applies to helix surfaces or loop regions as depicted in Figure 1. As a consequence the linewidth and the apparent hyperfine splitting of the EPR spectrum, reflected in the difference between the low and high field peak, are narrow. In contrast, the linewidth and the apparent hyperfine splitting increase for a restricted motional freedom of the spin label. This is the case for strong interaction of the nitroxide side chain with neighboring side chains or backbone atoms as found in tertiary contacts or at buried sites. However, in an equilibrium between two conformations with different dynamics or in the presence of two spin label rotamers with different degrees of motion, the detected cw EPR spectrum is the sum of the respective cw EPR spectra of each contribution. These multiple components in the spectra can be deduced by analyzing the cw EPR spectrum and may provide additional valuable information about the protein or protein complex.

Spin Label Side Chain Mobility
In general, cw EPR measurements performed at different temperatures yield diverse information about the spin label dynamics, the interaction with its local environment, the position of the spin label with respect to the membrane boundaries and interspin distances in the range of 0.8-1.8 nm. Cw EPR experiments at room temperature (RT) and the adjacent analysis gain information about the motion of the nitroxide side chain and allow for conclusions about the spin label's microenvironment as different environments effect the spin label's motion differently. The time scale of EPR spectroscopy is in the order of nanoseconds, such that it is suitable for a broad range of fundamental biological processes. Extensive studies of the relation of the EPR spectra and the protein structure are reported in literature [13,14]. In Section 3.2 we introduce an EPR study on the activation mechanisms in a functional chimera of the extracellular domain of amine-gated Erwinia chrysanthemi ligand-gated ion channel and the transmembrane domain of Gloeobacter violaceus ligand-gated ion channel. This investigation is based on the application of cw EPR experiments at room temperature.
The intrinsic spin label motion is characterized by its rotational correlation time, which depends on the spin label's intrinsic flexibility and its length. However, the measured effective rotational correlation time of a spin label bound to a protein is a result of the spin label's interaction with the primary, secondary, tertiary and eventually quaternary structure of the protein. Additionally, the protein itself tumbles in solution, which may add to the overall rotational correlation time. Although, the rotational correlation time for proteins, protein complexes or proteins inserted into detergent micelles or liposomes with a total molecular weight above 200 kDa exceeds the sensitivity of the EPR time scale and thus is not resolved in the room temperature cw EPR spectra. Consequently, the term 'mobility' can be mostly interpreted to describe the interaction of the nitroxide side chain with its local environment. Weak interactions and thus a high degree of spin label mobility are defined as only marginal or no restrictions of the motional freedom of the spin probe by its close proximity. This situation applies to helix surfaces or loop regions as depicted in Figure 1. As a consequence the linewidth and the apparent hyperfine splitting of the EPR spectrum, reflected in the difference between the low and high field peak, are narrow. In contrast, the linewidth and the apparent hyperfine splitting increase for a restricted motional freedom of the spin label. This is the case for strong interaction of the nitroxide side chain with neighboring side chains or backbone atoms as found in tertiary contacts or at buried sites. However, in an equilibrium between two conformations with different dynamics or in the presence of two spin label rotamers with different degrees of motion, the detected cw EPR spectrum is the sum of the respective cw EPR spectra of each contribution. These multiple components in the spectra can be deduced by analyzing the cw EPR spectrum and may provide additional valuable information about the protein or protein complex. Overview of the cw EPR methods at room temperature. (A) Representative cw EPR spectra at room temperature for unbound spin label and spin label with minor interaction (*) and strong interaction (#) with the local environment; (B) Cartoon representation of KtrB (PDB 4J7C) embedded in the membrane bilayer with the indicated paramagnetic quenchers (oxygen, NiEDDA and CrOx) and potential rotamers of the spin label side chain located in the loop region (*) and buried in the protein (#).

Figure 1.
Overview of the cw EPR methods at room temperature. (A) Representative cw EPR spectra at room temperature for unbound spin label and spin label with minor interaction (*) and strong interaction (#) with the local environment; (B) Cartoon representation of KtrB (PDB 4J7C) embedded in the membrane bilayer with the indicated paramagnetic quenchers (oxygen, NiEDDA and CrOx) and potential rotamers of the spin label side chain located in the loop region (*) and buried in the protein (#).
The semi-empirical mobility parameters ∆H −1 , the inverse linewidth of the central EPR peak, and the inverse second moment <H 2 > −1 , the spectral breadth, allow a qualitative analysis of room temperature EPR spectra [14,15]. They refer to the motional freedom of the nitroxide side chain and therefore correlate to its interaction with the local environment and accordingly to the proximate structure of the investigated protein. As a result of the interaction between the nitroxide side chain and its environment, the motion of the spin label is anisotropic [16,17]. Both semi-empirical mobility parameters ∆H −1 and <H 2 > −1 are affected by the degree of averaging the anisotropic g and the hyperfine tensors. Thereby, the anisotropic g reflects the orientation dependency of the magnetic moment of the electron with respect to the molecule or on the orientation of the molecule with respect to the magnetic field direction, respectively.
In the case of multiple components in the cw EPR spectrum the inverse linewidth ∆H −1 is most affected by the mobile component and is sensitive to changes in the mobile component rather than to changes in the immobile component. On the contrary, the second moment <H 2 > −1 is dominated by the immobile component. Plotting the semi-empirical mobility parameters as a function of a series of amino acid residues discloses a periodic variation of the nitroxide side chain mobility and permits the identification of secondary structure elements such as protein surfaces, tertiary contacts or buried sites and thus the assignment of α-helices, β-strands or unstructured regions.
A quantitative analysis of room temperature cw EPR spectra requires the simulation of the EPR spectra on the basis of dynamic models [18][19][20] or molecular dynamics (MD) simulations [21,22]. Both approaches provide additional information for local interactions of the nitroxide side chain and the dynamics of the investigated system. In the latter case the analysis of cw EPR spectra supports the optimization of structural models as provided by structural biology.

Accessibility of the Spin Label Side Chain
A more advanced methodology compared to the beforehand-described cw EPR measurements is the determination of the accessibility of the nitroxide side chain towards paramagnetic quenchers. These cw EPR experiments carried out at room temperature provide additional information about the positioning of the EPR spin probe with respect to membrane/water/protein boundaries [23][24][25][26]. For instance, the analysis of accessibility experiments in combination with side chain mobility and interspin distance measurements has verified a beforehand proposed model for the HAMP domain of the NpHtrII/NpSRII complex from Natronomonas pharaonic (cf. Section 3.1).
The term accessibility is reflected by the Heisenberg exchange frequency W ex of the spin label with a paramagnetic reagent in solution. To determine the Heisenberg exchange frequency of a spin label with different paramagnetic quenchers initially the saturation behavior of a specific nitroxide side chain is determined. Here, saturation describes the state where the rate of excitation of the spin system exceeds the ability of the spin system to relax to equilibrium. In this context, the relaxation time is given by the spin-lattice relaxation time T 1 . Plotting the amplitude of the central EPR peak as a function of the square root of the incident microwave yields the so-called saturation curve. With increasing microwave power the amplitude of the EPR signal increases to its maximum followed by a decrease caused by saturation. Fitting of the saturation curve provides the parameter P 1/2 which is defined as the microwave power at which the amplitude of the EPR signal is half of its maximum in the absence of saturation and correlates to the relaxation rates of the nitroxide side chain. The presence of paramagnetic quenchers during accessibility measurements leads to frequent collisions with the spin label and causes Heisenberg exchange. This in turn results in alterations of the relaxation time of the spin label and consequently alterations in the values for P 1/2 . The difference of the P 1/2 values in the presence and in the absence of the paramagnetic reagent allows the determination of the Heisenberg exchange frequency W ex . A more detailed discussion of accessibility measurements can be found in [4]. For the analysis of the insertion of a membrane protein inside a lipid bilayer different paramagnetic quenchers with different permeabilities for different environments as visualized in Figure 1B are used. Possible paramagnetic reagents are polar metal complexes (Ni(II) complexes or chromium oxalate (CrOx)), molecular oxygen or organic radicals. Commonly utilized polar metal reagents such as NiEDDA and CrOx are not membrane permeable and therefore accumulate in the aqueous phase providing insights about the accessibility of the lipid/water interface. However, the application of CrOx potentially impedes the analysis due to its negative charge and the possibility of Coulomb interactions with charged membrane surfaces or protein segments, such that the use of NiEDDA is preferred. Molecular oxygen, on the other hand, preferably remains in the lipid bilayer with its maximum in the center of the membrane [26]. By determining the Heisenberg exchange frequencies W ex for a series of spin labeled residues secondary structure elements and the location of the spin label side chains with respect to the membrane interfaces can be assigned.
Throughout the last years further EPR approaches have been introduced, which allow the direct determination of the Heisenberg exchange frequency W ex and/or the relaxation time T 1 such as multiquantum EPR [27] or saturation recovery EPR [28,29].

Polarity of the Spin Label Side Chain Environment
Low temperature cw EPR measurements (T = 160 K) are applied to achieve information about the polarity of the spin label side chain environment and to utilize the dipolar interaction between nitroxide side chains to determine interspin distances [30]. The determination of the polarity of the spin label environment is based on the hyperfine coupling which is directly obvious in the cw EPR spectrum in the difference between the low and high field peak. In more detail, the strength of the hyperfine coupling is described by the hyperfine splitting A zz . In general, a high accessibility of water molecules to the spin label side chain causes a strong hyperfine coupling and consequently a high A zz value (cf. Figure 2A). In contrast, a low polarity of the microenvironment results in a weak hyperfine coupling and therefore small A zz values. This is associated with a low accessibility of the EPR spin probe to water molecules (cf. Figure 2B). Thereby, the determination of A zz values identifies spin labeled residues with high, low or no accessibility for water molecules and thus their orientation with respect to the lipid-water interface. Precise values for the hyperfine splitting A zz are obtained by fitting the cw EPR spectra with the EPR spectra simulation program DipFit [31] (cf. Figure 2C). Plotting the values for A zz as a function of a series of amino acid residues reveals additional structural and topological information and evidence for protein folding or existing hydrophobic barriers. and therefore accumulate in the aqueous phase providing insights about the accessibility of the lipid/water interface. However, the application of CrOx potentially impedes the analysis due to its negative charge and the possibility of Coulomb interactions with charged membrane surfaces or protein segments, such that the use of NiEDDA is preferred. Molecular oxygen, on the other hand, preferably remains in the lipid bilayer with its maximum in the center of the membrane [26]. By determining the Heisenberg exchange frequencies Wex for a series of spin labeled residues secondary structure elements and the location of the spin label side chains with respect to the membrane interfaces can be assigned. Throughout the last years further EPR approaches have been introduced, which allow the direct determination of the Heisenberg exchange frequency Wex and/or the relaxation time T1 such as multiquantum EPR [27] or saturation recovery EPR [28,29].

Polarity of the Spin Label Side Chain Environment
Low temperature cw EPR measurements (T = 160 K) are applied to achieve information about the polarity of the spin label side chain environment and to utilize the dipolar interaction between nitroxide side chains to determine interspin distances [30]. The determination of the polarity of the spin label environment is based on the hyperfine coupling which is directly obvious in the cw EPR spectrum in the difference between the low and high field peak. In more detail, the strength of the hyperfine coupling is described by the hyperfine splitting Azz. In general, a high accessibility of water molecules to the spin label side chain causes a strong hyperfine coupling and consequently a high Azz value (cf. Figure 2A). In contrast, a low polarity of the microenvironment results in a weak hyperfine coupling and therefore small Azz values. This is associated with a low accessibility of the EPR spin probe to water molecules (cf. Figure 2B). Thereby, the determination of Azz values identifies spin labeled residues with high, low or no accessibility for water molecules and thus their orientation with respect to the lipid-water interface. Precise values for the hyperfine splitting Azz are obtained by fitting the cw EPR spectra with the EPR spectra simulation program DipFit [31] (cf. Figure 2C). Plotting the values for Azz as a function of a series of amino acid residues reveals additional structural and topological information and evidence for protein folding or existing hydrophobic barriers.  [31] with equal g-and A-tensor elements and variable interspin distances. Twice the hyperfine coupling (2Azz) is indicated. All plots are normalized by spin number.

Interspin Distance Measurements
Today, the determination of interspin distances is the main application of EPR spectroscopy on proteins. Instead of one, two or more spin probes are inserted into the studied system to evaluate distance distributions and their changes at varying conditions. By this, existing protein structures can be validated and unknown structures predicted. Distance changes disclose conformational The external electric field component E x is caused by the polar micro-environment which shift the π-electron density along the N-O bond towards the N-atom represented by increased orbitals; (B) In a non-polar environment the π-electron density is shifted towards the O-atom; (C) Simulated low temperature (T = 160 K) cw EPR powder spectra corresponding to interspin distances of 1.3 nm (black line), 1.5 nm (red line), 1.8 nm (green line) and 2.0 nm (blue line), respectively. Spectra are obtained by using DipFit [31] with equal g-and A-tensor elements and variable interspin distances. Twice the hyperfine coupling (2A zz ) is indicated. All plots are normalized by spin number.

Interspin Distance Measurements
Today, the determination of interspin distances is the main application of EPR spectroscopy on proteins. Instead of one, two or more spin probes are inserted into the studied system to evaluate distance distributions and their changes at varying conditions. By this, existing protein structures can be validated and unknown structures predicted. Distance changes disclose conformational rearrangements or complex formations. Cw and pulsed EPR spectroscopy at low temperature are the methods of choice to determine interspin distances between two spin label side chains in a distance range of 0.8-6 nm [32][33][34]. The studies presented in Sections 3.3 and 3.4 on the ABC exporter TM287/288 from Thermotoga maritima and on the secondary active transporter Glt Ph from Pyrococcus horikoshii provide excellent examples for the application of interspin distance measurements. Both cases visualize the potential for identifying the molecular triggers that lead to conformational changes, here from an open to a closed NBD dimer in TM287/288 and a possible transport mechanism for the trimeric Glt Ph , respectively.
As discussed before each spin label possesses a specific length and an intrinsic flexibility, which'causes the existence of a distance distribution instead of a defined interspin distance. In addition, protein motion occurs over time scales ranging from femtoseconds (vibrations of individual bonds) to nano-and microseconds and longer (collective motions of groups of bonded and non-bonded atoms). In frozen solutions below the so-called glass transition temperature of proteins, at approximately 200 K, the dynamical properties of proteins change and the collective motions of atom groups and backbone motions are disabled. Thus, only harmonic vibrations take place. At this point, the distance distributions and the relative orientations of the nitroxide side chains are homogeneously distributed and assumed to reflect the situation present at room temperature.
Cw EPR measurements at 160 K allow the determination of interspin distances in the range of 0.8-1.8 nm based on the dipolar interaction of two spin label side chains [35]. Here, a dipolar interaction results in a considerable line broadening of the cw EPR spectra compared to the cw EPR spectra of single labeled proteins or of double labeled proteins with interspin distances larger than 1.8 nm as depicted in Figure 2C. The dipolar broadening increases with decreasing interspin distances. Fitting of the dipolar broadened EPR spectra is carried out for instance with the EPR spectra simulation program DipFit [31], which allows a quantitative analysis of EPR distance measurements. Here, accurate values for the interspin distances are derived from the precise line shape analysis of the broadened EPR spectrum with respect to beforehand-defined parameters of the respective non-broadened EPR spectrum. The best-fit parameters are computated taking a Gaussian distribution of interspin distances and an arbitrary fraction of single labeled protein under consideration. The fraction of single labeled protein can be estimated based on the spin label efficiency. Therefore, the respective experimental cw EPR spectrum at RT is double integrated and compared with the double integral of a MTS spin label standard of a known concentration to determine the concentration of the spin label in the protein sample. A subsequent comparison with the known protein concentration allows the rough calculation of the spin label efficiency. More qualitative approaches of dipolar interactions employ the correlation between the line width of the EPR spectra and the interspin distances [36]. For this purpose the spectral line amplitude ratios or the spectral second moment values are calculated, which identifies spin probes in close or distant proximity. This is commonly executed for a series of double-labeled protein [37] as this approach does not allow for the determination of the exact interspin distance values. Interspin distance values below 0.8 nm are not reliable due to the increased partial overlap of the nitrogen π-orbitals of the two interacting spin label side chains [38].
The currently most prominent applications for distance determinations are pulsed EPR methods, which allow to enlarge the measurable distance range of interspin distances up to 8 nm [33]. These measurements are performed at 50 K, which is a compromise designed to maximize T 2 (spin-spin relaxation time) and minimize T 1 (spin-lattice relaxation time) permitting experiments with a long persistent spin echo and a high repetition rate. Pulsed electron-electron double resonance (PELDOR) or also called double electron-electron resonance (DEER) experiments resolve long range distances by commonly utilizing a 4-pulse sequence at X-or Q-band (9.4 GHz and 34 GHz, respectively) as shown in Figure 3A [33,39]. These experiments are based on the selective excitation of two interacting spins fractions A and B using two different frequencies, the observer frequency ν obs and the pump frequency ν pump . The pump frequency ν pump is positioned at the maximum of the echo-detected nitroxide EPR absorption spectrum and the observer frequency ν obs at the low field local maximum resulting commonly in a 65 MHz and 50 MHz frequency offset for X-and Q-band, respectively. The excitation of the pump spin fraction B results in a resonance frequency shift of the dipolar coupled spins A. The Hahn-echo sequence, consisting of the first two pulses with the time delay of τ 1 at the observer frequency ν obs , creates the spin echo, which is also called Hahn-echo, after another time delay τ 1 . The π-pulse after time 2τ 1 + τ 2 leads to the refocused and monitored echo signal at 2τ 1 + 2τ 2 . Applying a π-pulse at the pump frequency ν pump selectively excites the spin fraction B which causes a resonance frequency shift of the dipolar coupled spin fraction A. For the experimental DEER signal V(t) the echo amplitude is determined in dependence of time t, corresponding to the pump pulse position [40]. Different t values result in a modulation of the detected echo signal with the dipolar frequency ω DD as depicted in Figure 3B. Furthermore, the observed experimental signal V(t) comprises the desired intramolecular contribution F(t) modulated with the frequency ω DD and an unmodulated background contribution B(t) (cf. Figure 3B, dotted line) resulting from homogeneous distributed molecules and therefore homogeneous distributed interspin distances within the sample. The modulated contribution shown in Figure 3C, also called form factor or background-corrected dipolar evolution data, has to be separated from the background by division of V(t) by B(t). Conventional, the final distance distribution is obtained by Tikhonov regularization [41,42] as utilized in the program DeerAnalysis [43] (cf. Figure 3D). In addition to the mean interspin distances the interspin distance distribution provides information about the distribution of the occupied conformations of the macromolecule and the flexibility of the spin label side chain reflected in the line width. To evaluate the quality and the reliability of the determined distance distribution a validation should be performed as implemented in the validation tool of DeerAnalysis. In order to compare the gained mean interspin distances and distance distributions with existing structural models conformational searching approaches such as the Multiscale Modeling of Macromolecules (MMM) software proved to be insightful [44]. Furthermore, the elastic network model (ENM), as implemented in MMM, has developed into an important tool in the construction of structural models. Here, the utilization of distance distribution constraints obtained by pulsed EPR measurements and an initial structure render new models possible. For instance cAMP-bound and unbound models of hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels have been successfully generated based on interspin distance constraints and the cAMP-bound HCN crystal structure as the starting structure [45]. excitation of the pump spin fraction B results in a resonance frequency shift of the dipolar coupled spins A. The Hahn-echo sequence, consisting of the first two pulses with the time delay of τ1 at the observer frequency νobs, creates the spin echo, which is also called Hahn-echo, after another time delay τ1. The π-pulse after time 2τ1 + τ2 leads to the refocused and monitored echo signal at 2τ1 + 2τ2. Applying a π-pulse at the pump frequency νpump selectively excites the spin fraction B which causes a resonance frequency shift of the dipolar coupled spin fraction A. For the experimental DEER signal V(t) the echo amplitude is determined in dependence of time t, corresponding to the pump pulse position [40]. Different t values result in a modulation of the detected echo signal with the dipolar frequency ωDD as depicted in Figure 3B. Furthermore, the observed experimental signal V(t) comprises the desired intramolecular contribution F(t) modulated with the frequency ωDD and an unmodulated background contribution B(t) (cf. Figure 3B, dotted line) resulting from homogeneous distributed molecules and therefore homogeneous distributed interspin distances within the sample. The modulated contribution shown in Figure 3C, also called form factor or background-corrected dipolar evolution data, has to be separated from the background by division of V(t) by B(t). Conventional, the final distance distribution is obtained by Tikhonov regularization [41,42] as utilized in the program DeerAnalysis [43] (cf. Figure 3D). In addition to the mean interspin distances the interspin distance distribution provides information about the distribution of the occupied conformations of the macromolecule and the flexibility of the spin label side chain reflected in the line width. To evaluate the quality and the reliability of the determined distance distribution a validation should be performed as implemented in the validation tool of DeerAnalysis. In order to compare the gained mean interspin distances and distance distributions with existing structural models conformational searching approaches such as the Multiscale Modeling of Macromolecules (MMM) software proved to be insightful [44]. Furthermore, the elastic network model (ENM), as implemented in MMM, has developed into an important tool in the construction of structural models.
Here, the utilization of distance distribution constraints obtained by pulsed EPR measurements and an initial structure render new models possible. For instance cAMP-bound and unbound models of hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels have been successfully generated based on interspin distance constraints and the cAMP-bound HCN crystal structure as the starting structure [45]. The π/2 pulse and π pulse with a time delay of τ1 at the observer frequency νobs generates a spin echo after time delay τ1. The third π-pulse after 2τ1 + τ2 results in the refocused and monitored echo signal at 2τ1 + 2τ2. The π-pulse at the pump frequency νpump selectively excites spin fraction B causing a resonance frequency shift of dipolar coupled spin fraction A; (B) Low temperature (T = 50 K) pulsed EPR spectrum recorded at Xband with indicated background (grey, dotted line) (C) background corrected dipolar evolution data F(t); (D) Distance distributions P(r) obtained by Tikhonov regularization as utilized in the program DeerAnalysis [43]. The π/2 pulse and π pulse with a time delay of τ 1 at the observer frequency ν obs generates a spin echo after time delay τ 1 . The third π-pulse after 2τ 1 + τ 2 results in the refocused and monitored echo signal at 2τ 1 + 2τ 2 . The π-pulse at the pump frequency ν pump selectively excites spin fraction B causing a resonance frequency shift of dipolar coupled spin fraction A; (B) Low temperature (T = 50 K) pulsed EPR spectrum recorded at X-band with indicated background (grey, dotted line) (C) background corrected dipolar evolution data F(t); (D) Distance distributions P(r) obtained by Tikhonov regularization as utilized in the program DeerAnalysis [43].

EPR Spectroscopy Applied to Different Classes of Membrane Proteins
In the following sections we will review different examples for the elegant use of EPR spectroscopy methods on membrane proteins and their combination with structural biology approaches. We focus on four different classes, namely receptors, ion channels, ABC transporters and secondary active transporters, since each of these classes has its own challenges. For each class we highlight one particular paper, which in our opinion represents an excellent example for the use of one or more of the above-mentioned EPR methods to answer defined research questions. Simultaneously, on the basis of these papers we emphasize the complementary and novel information derived from EPR spectroscopy. For instance, EPR spectroscopy has provided information about an equilibrium between two conformations distinguishable by their dynamics and the accessibilities for water, where NMR spectroscopy could only capture the compact conformation (cf. Section 3.1). Furthermore, with the application of EPR experiments technical challenges related to crystallization in the presence of 3-amino propanol and under acidic conditions have been overcome monitoring amine-and pH-induced conformational changes as demonstrated in Section 3.2. Upon the addition of different nucleotides in combination with inhibitors and functional mutations specific conformations can be trapped as performed in Section 3.3. In the last section we show how EPR spectroscopy can be applied to propose a transport mechanism.

Receptors
Receptors are transmembrane proteins, which bind external signaling molecules (ligands) or respond directly to light or pressure and subsequently transmit a signal into a cell or tissue. Membrane receptors are divided into ionotropic receptors, which upon activation open intrinsic ion channels, metabotropic or G protein-coupled receptors, and enzyme-linked e.g., kinase-linked receptors. Upon signal perception different intracellular cascades are initiated which consequently affect the cell's function. Since receptors are extracellular accessible they serve as major drug targets: Agonists activate the receptor signaling while antagonists bind to a receptor and block it. The development of new drugs is driven by structural and mechanistic information as gained by combinations of mostly X-ray crystallography and EPR spectroscopy. A prominent example studied for decades is rhodopsin, which has been investigated by the whole range of EPR methods discussed above (see review [46]). However, also recent publications on the β 2 -adrenergic receptor [47] and its interaction with heterotrimeric G proteins [48] highlight the synergistic effects of combining structural biology and EPR spectroscopy to understand the dynamic processes of receptor signaling. Light-activated ionotropic receptors are key players in optogenetics. Here, genetically modified cells e.g., neurons expressing light-sensitive receptors are controlled by light, which enables a high spatial and temporal resolution of activation. A function-specific engineering of receptors again requires the combined use of structural approaches like X-ray crystallography and NMR with EPR spectroscopy, which provides dynamic details and allows the prediction of non-resolved structures. Recent examples are two studies on Channelrhodopsin-2 [49,50], which identified helix movements upon light activation.
Intensively studied examples for light-activated metabotropic receptors are sensory rhodopsins I and II. In the haloarchaeon Halobacterium salinarum light sensing is executed by two sensory rhodopsins, SRI and SRII. According to the wavelength of the incident light these phototaxis receptors send signals to their associated halobacterial transducer proteins (HtrI and HtrII) and provide the cell with the ability to respond to light as an attractant or repellent [51]. For light with a wavelength below 500 nm sensory rhodopsin II (SRII) initiates a strong repellent response to avoid photooxidative stress in the presence of oxygen. After light activation of SRII the signal is transferred to the tightly bound transducer protein HtrII which in turn triggers a phosphorylation cascade modulating the cell's flagellar motors including the histidine kinase CheA and the response regulators CheY and CheB [52,53]. Sensory rhodopsins are transmembrane proteins comprising seven helices and a retinal cofactor bound to a lysine residue on helix G via a protonated Schiff base [54,55]. Light excitation triggers the isomerization of the retinal chromophore from the all-trans to a 3-cis configuration which is the starting event for the signal transfer. The overall structure of the transducer is similar to that of chemoreceptors. The transducer protein HtrII contains a transmembrane and a cytoplasmic signaling domain, linked by a HAMP5 domain (conserved in histidine kinases, adenylyl cyclases, methyl-accepting chemotaxis proteins, and phosphatases [56]). The transmembrane N-terminal domain is comprised of two transmembrane helices (TM1 and TM2) and incorporates part of the receptor-binding site. The HAMP domain in turn is composed of two helical amphipathic sequences (AS-1 and AS-2) separated by a connector sequence and functions as a ubiquitous signal transduction module [56,57]. The succeeding rod-shaped cytoplasmic domain includes the sub-domains for signal transfer and for adaptation. In membranes, the phototransducers constitute a 2:2 complex with their cognate photoreceptors [58,59].
In the elaborate characterization of the phototransducer from Natronomonas pharaonis (NpHtrII) in complex with sensory rhodopsin II (NpSRII) one breakthrough has been the publication of the solution NMR structure of the HAMP domain from Archaeoglobus fulgidus [60]. The NMR structure displays a homodimeric, four helical, parallel coiled coil with an unusual interhelical packing. The assumption that the parallel four-helical bundle structure is likely to represent at least one possible conformation of the HAMP domain is supported by the crystal structures for the concatenated HAMP domains from the Pseudomonas aeruginosa soluble receptor Aer2 [61] revealing the same general structural characteristics. Based on the first NMR structure, a new model has been proposed in which the HAMP domain alternates between two conformations, one conformation is captured by the NMR structure, as shown in Figure 4A, whereas the second conformation resembles a canonical coiled coil via concerted helix rotations [60]. In the EPR study selected here this proposed model is confirmed by characterizing the structural properties of the HAMP domain of the NpHtrII/NpSRII complex from Natronomonas pharaonic [62]. Therefore, mobility and accessibility values for a series of spin-labeled variants at low and high salt conditions (cf. Figure 4B,C) at room temperature and interspin distances at 160 K and 50 K for cw and pulsed EPR experiments, respectively, have been determined. Positions for spin labeling have been chosen based on a sequence alignment of the two HAMP domains from Af1503 (A. fulgidus) and from NpHtrII (N. pharaonic) as there has been no structure available of the HAMP domain from NpHtrII. Concentrating on cw power saturation EPR experiments this paper provides an excellent example for its use and visualizes its potential for identifying residues with high or low accessibility for water and verifying a conformational equilibrium between two conformations distinguishable by their dynamics and the accessibilities for water.
The analysis of the mobility values has provided first evidence for the presence of a dynamic HAMP domain (dHAMP) at low salt concentrations and a compact HAMP domain (cHAMP) at high salt concentrations (for details see [62]). This assumption is supported by accessibility measurements which differentiate between nitroxide side chains exposed to buffer, oriented towards the lipid phase or buried in the protein interior as discussed before. At low salt concentrations the accessibility values for oxygen relate to spin labels exposed to water (thus are comparably low) and are almost constant which is obvious in Figure 4B. This accessibility relates to a dynamic, water-exposed HAMP domain (dHAMP). Whereas at high salt concentrations the accessibility values decrease even further and a periodic pattern is obvious starting at position 95, indicating a compact HAMP domain (cHAMP) structure with spin-labeled residues buried inside the protein at high salt condition with a well-defined secondary and tertiary interactions. The accessibility values for NiEDDA confirm this conclusion ( Figure 4C) as at low salt concentrations and 3 mM NiEDDA relatively high Heisenberg exchange rates are determined indicative for a high water accessibility. At high salt concentrations as well as at low salt concentrations but high NiEDDA concentration (20 mM) low Heisenberg exchange rates with a periodic pattern are observed which are in agreement with the cHAMP structure. While expected for the high salt condition based on the oxygen accessibility measurements, this result is kind of surprising at low salt concentrations. However, the authors show by additional cw EPR measurements at room temperature that with 20 mM NiEDDA the mobile component, which agrees to the dHAMP structure, is almost completely suppressed such that also under this condition the HAMP domain to 90% adopts its compact structure. The authors thus demonstrate by the discussed accessibility measurements that the side chains in the cHAMP conformation are restricted in their mobility and poorly accessible from the bulk water. This conformation resembles the four helical bundle structure obtained by solution NMR for the HAMP domain of A. fulgidus. In the dHAMP conformation the spin label dynamics and its accessibility for water are very high, which suggests a rather unfolded secondary structure. Thus, changes in the dynamics of the HAMP domains of these systems can be hypothesized as the language of signaling. A detailed review about the progress in the EPR study of sensory rhodopsin-transducer complex can be found elsewhere [63]. However, the discussed paper highlights the necessity of the combined approach of structural biology and EPR spectroscopy and its potential to complement missing information even in the absence of a corresponding structure and to enhance our knowledge about structural properties and their relation to the function of proteins.
Crystals 2017, 7, 117 10 of 28 mobility and poorly accessible from the bulk water. This conformation resembles the four helical bundle structure obtained by solution NMR for the HAMP domain of A. fulgidus. In the dHAMP conformation the spin label dynamics and its accessibility for water are very high, which suggests a rather unfolded secondary structure. Thus, changes in the dynamics of the HAMP domains of these systems can be hypothesized as the language of signaling. A detailed review about the progress in the EPR study of sensory rhodopsin-transducer complex can be found elsewhere [63]. However, the discussed paper highlights the necessity of the combined approach of structural biology and EPR spectroscopy and its potential to complement missing information even in the absence of a corresponding structure and to enhance our knowledge about structural properties and their relation to the function of proteins.

Ion Channels
Ion channels are pore-forming membrane proteins that allow the more or less selective translocation of cations and anions along an electrochemical gradient. Ion channels are found in all living cells and are responsible for a plethora of tasks, including pH homeostasis, cell growth, maintenance of osmolarity and cell volume, movement and, probably most importantly, electrical signaling. As diverse as their function is the regulation of ion channels. Besides being part of the aforementioned receptors, ion channels group into voltage-gated ion channels, intracellular ligandgated ion channels, cyclic nucleotide-gated ion channels, temperature-gated ion channels, and mechanosensitive and volume-regulated ion channels. Understanding the mechanisms of activation, deactivation and channel gating are the key questions which are investigated. In addition, the molecular basis of ion selectivity and possibly rectification are in focus. The prime example of an ion channel is KcsA, the first ion channel ever crystallized [64]. Besides nowadays various structures a plethora of EPR studies has been performed with regard to conformational changes in the selectivity filter [65,66], the cytoplasmic helical bundle [67,68] and the communication between both features [69]. Further excellent examples for the combined use of structural biology and EPR spectroscopy are present for all classes of ion channels, e.g., for the intracellular Mg 2+ -gated Mg 2+ channel CorA [70][71][72], the nucleotide-and Na + -gated K + channel KtrAB [73], the mechanosensitive channel MscL [74,75],

Ion Channels
Ion channels are pore-forming membrane proteins that allow the more or less selective translocation of cations and anions along an electrochemical gradient. Ion channels are found in all living cells and are responsible for a plethora of tasks, including pH homeostasis, cell growth, maintenance of osmolarity and cell volume, movement and, probably most importantly, electrical signaling. As diverse as their function is the regulation of ion channels. Besides being part of the aforementioned receptors, ion channels group into voltage-gated ion channels, intracellular ligand-gated ion channels, cyclic nucleotide-gated ion channels, temperature-gated ion channels, and mechanosensitive and volume-regulated ion channels. Understanding the mechanisms of activation, deactivation and channel gating are the key questions which are investigated. In addition, the molecular basis of ion selectivity and possibly rectification are in focus. The prime example of an ion channel is KcsA, the first ion channel ever crystallized [64]. Besides nowadays various structures a plethora of EPR studies has been performed with regard to conformational changes in the selectivity filter [65,66], the cytoplasmic helical bundle [67,68] and the communication between both features [69]. Further excellent examples for the combined use of structural biology and EPR spectroscopy are present for all classes of ion channels, e.g., for the intracellular Mg 2+ -gated Mg 2+ channel CorA [70][71][72], the nucleotide-and Na + -gated K + channel KtrAB [73], the mechanosensitive channel MscL [74,75], the cyclic AMP-gated HCN channel [45,76], the voltage-gated Na + channel NavMs [77] and the pentameric receptor channels GLIC [78]. In all cases EPR measurements have been applied to validate existing protein structures, determine the structure of missing regions or distorted elements in structures, or add valuable information of the dynamics of the studied systems to provide a model of activation and regulation. A broader overview about EPR studies on ion channels can be found in a recent review [79]. A particular challenge for the investigation of ion channels by EPR spectroscopy is their mostly homo-oligomeric nature, which requires the interpretation of various distance distributions even if only one spin label is inserted per protomer [80]. A great advantage over most structural studies is the possibility to perform EPR measurements on liposome-reconstituted proteins such that electrochemical gradients as needed for the open conformation can be applied [73].
Synaptic transmissions in the central and peripheral nervous systems are conducted by pentameric ligand-gated ion channels (pLGICs). Nicotinic acetylcholine receptors (nAChRs), the γ-aminobutyric acid receptors (GABAARs), glycine receptors (GlyRs), and serotonin receptors (5HT3Rs) are members of the family of pLGICs [81][82][83][84]. The dysfunction of the synaptic transmission is associated with a multitude of neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease or schizophrenia. Thus, pLGICs are suitable candidates for therapeutic investigations and approaches. Here, neurotransmitter binding initiates a rapid switch between the closed, open, and desensitized conformations in the pLGICs. The switch between these states occurs over a 6-nm distance along the length of the channel. The addition of drugs and modulators allosterically shifts the equilibrium toward the open and closed/desensitized conformation, respectively, leading to an enhanced or inhibited channel activity. However, the underlying gating and drug modulation at the molecular level is still under debate. A recent breakthrough has been the publications of several high resolution structures such as the glutamate-gated chloride channel [85], the human 3-GABAAR [86], the mouse 5-HT3 [87], and the prokaryotic homologues Gloeobacter violaceus ligand-gated ion channel (GLIC; pH gated) [88,89] and Erwinia chrysanthemi ligand-gated ion channel (ELIC; primary amine gated) [90]. Due to their high structural similarity to eukaryotic homologues, GLIC and ELIC are considered to be suitable structural surrogates to examine channel gating and its underlying structural changes and to identify conserved drug binding regions. First information about agonist and modulator binding sites but not about gating or modulation has been conducted from ELIC structures [91][92][93][94]. The here selected advanced approach of cw and pulsed EPR spectroscopy in combination with X-ray crystallography demonstrates impressively the advantages of such a combination and how EPR spectroscopy is able to compensate for technical challenges related to crystallization. The paper investigated the activation mechanisms in a functional chimera of the extracellular domain (ECD) of amine-gated Erwinia chrysanthemi ligand-gated ion channel and the transmembrane domain (TMD) of Gloeobacter violaceus ligand-gated ion channel by applying cw and pulsed EPR spectroscopy [95]. However, this review focuses on the analyses of the cw EPR data proving that cw EPR studies at room temperature and the derived information about the spin label side chain dynamics yield valuable insights for the understanding of structural mechanism.
The chimera has been crystallized at pH 7.0 and in the absence of 3-AP (3-amino propanol) to capture the resting state (4.6-Å resolution). The overall arrangement of the TMD and the orientation of M2 in the ELIC-GLIC chimera are similar to GLIC at pH 7.0. A comparison of the chimera and the apo-ELIC-ECD highlighted a rotation of the chimera-ECD about the central symmetry axis such that the twist is most pronounced at the ECD-TMD interface. The chimera is in its closed-pore conformation. To achieve an enhanced understanding of gating further crystal structures in the presence of 3-AP and under acidic conditions are necessary. To overcome technical challenges related to crystals under these conditions cw and pulsed EPR experiments have been performed to monitor amine-and pH-induced conformational changes. It could be shown by pulsed EPR experiments at 50 K that structural rearrangements of loop C are not necessary for channel opening and gating and that the motions of loop F and loop C are not coupled as part of channel opening. Therefore, room temperature cw EPR measurements at representative positions in loops β1-β2, β4-β5, and β6-β7 in the ECD and in the M2-M3 linker and the C-terminal end of M4 in the TMD have been performed in the presence of 3-AP ( Figure 5, orange lines) and in the presence of protons ( Figure 5, green lines), respectively. These positions are chosen based on their corresponding positions in GLIC [96] and potential conformational changes.  First, mobility changes in the ECD upon activation are investigated. The cw EPR spectra for V82R1 (β4-β5 linker) display a slightly increased spin label mobility in the presence of 3-AP ( Figure  5A) and a more pronounced increase of the mobility in the presence of protons. In contrast, for L29R1 the addition of 3-AP or protons results in a decrease of spin label mobility with a more pronounced decrease at pH 3.0. Similarly, for both positions in the loop β6-β7 protons induce stronger changes in the line shape compared to 3-AP binding, although the trend is opposite for the positions M114 and F116. For variant M114R1 the cw EPR spectra exhibit a slightly decreased spin label mobility at pH 3.0 and upon addition of 3-AP, respectively. For variant F116R1 an increased spin label mobility is obvious at pH 3.0 and no significant change in the presence of 3-AP. Also in the TMD residues undergo larger conformational changes in the presence of protons compared to 3-AP. In addition to the mobility changes, variant K255R1 (M2-M3 loop) displays strong spectral broadening most likely due to dipolar coupling caused by the close proximity of the nitroxide side chains (Cβ-Cβ distance of 1.7 nm between adjacent subunits in the pentamer [96]) as shown in Figure 5B. In the presence of protons the spectral broadening decreased, which is in agreement with an outward movement of loop M2-M3 as reported in GLIC [96,97]. In comparison, the effect is much weaker upon 3-AP First, mobility changes in the ECD upon activation are investigated. The cw EPR spectra for V82R1 (β4-β5 linker) display a slightly increased spin label mobility in the presence of 3-AP ( Figure 5A) and a more pronounced increase of the mobility in the presence of protons. In contrast, for L29R1 the addition of 3-AP or protons results in a decrease of spin label mobility with a more pronounced decrease at pH 3.0. Similarly, for both positions in the loop β6-β7 protons induce stronger changes in the line shape compared to 3-AP binding, although the trend is opposite for the positions M114 and F116. For variant M114R1 the cw EPR spectra exhibit a slightly decreased spin label mobility at pH 3.0 and upon addition of 3-AP, respectively. For variant F116R1 an increased spin label mobility is obvious at pH 3.0 and no significant change in the presence of 3-AP. Also in the TMD residues undergo larger conformational changes in the presence of protons compared to 3-AP. In addition to the mobility changes, variant K255R1 (M2-M3 loop) displays strong spectral broadening most likely due to dipolar coupling caused by the close proximity of the nitroxide side chains (C β -C β distance of 1.7 nm between adjacent subunits in the pentamer [96]) as shown in Figure 5B. In the presence of protons the spectral broadening decreased, which is in agreement with an outward movement of loop M2-M3 as reported in GLIC [96,97]. In comparison, the effect is much weaker upon 3-AP addition. For variant F322R1, located at the C-terminal end of M4, the addition of 3-AP results in only minimal changes of the cw EPR spectrum. Whereas at pH 3.0 the immobile component becomes less pronounced. Further investigation of the ECD-TMD interface in the chimera and those at the equivalent positions in ELIC and GLIC provide evidence that for positions at the ECD-TMD interface the obtained differences between the chimera and the parent channels are consistent with a twist in the chimera ECD with respect to the TMD. While both 3-AP and protons are able to gate the ELIC-GLIC chimera, the degree of mobility changes upon gating depends on the type of ligand. In response to an acidification, the domain interfaces undergo large conformational rearrangements, consistent with conformational changes found in GLIC. Moreover, in response to 3-AP modest changes are observed in agreement with effects observed in ELIC. Thus, the ECD and TMD of the ELIC-GLIC chimera preserve the distinct conformational changes of the parent channels. By means of applying EPR spectroscopy in combination with structural information derived from crystal structures the authors suggested that the chimeric channel is gated not only by two different ligands but by two different mechanisms. These results enable the design of better allosteric modulators.

ABC Transporters
ATP-binding cassette transporters (ABC transporters) are a family of membrane proteins with representatives in all phyla of life from prokaryotes to humans. All ABC transporters consist of at least a pair of cytoplasmic nucleotide binding domains (NBDs), which are associated to two transmembrane domains (TMDs). Those domains can be arranged on a single polypeptide, as homo-or heterodimer with each one NBD and one TMD, with fused NBDs but single TMDs, or as single subunits. Despite those differences, the general overall topology is similar allowing a mechanism, in which ATP binding and hydrolysis, and ADP and P i release in the NBDs catalyze the transport of diverse substrates by driving distinct conformational changes of the TMDs. Nonetheless, one cannot define a uniform mechanism of transport because ABC transporters in particular vary in the architecture of their TMDs. Based on their structural and functional similarities ABC transporters are currently grouped into four types namely type 1 and type 2 importers, exporters and ECF transporters (type 3 importers). ABC importers are exclusively found in prokaryotes. Type 1 and 2 importers contain an additional, mostly soluble substrate binding domain (SBD). Type 1 importers facilitate the uptake of diverse nutrients including ions, sugars, amino acids, short peptides, and oligosaccharides; type 2 importers support the high-affinity uptake of metal chelates including heme and other iron-containing complexes, and cobalamin; ECF transporters take up vitamins with high affinity; and ABC exporters mediate the transport of a wide range of cargos such as lipids, peptides, ions and drugs across lipid bilayers [98].
The mechanistic diversity of ABC transporters has been summarized in a recent review [99]. The mechanisms as proposed are a result of combining structural, functional, biophysical and spectroscopic studies. For type 1 and type 2 importers each one common transport model has been proposed, which are supported by several EPR studies [100][101][102][103][104][105][106][107][108][109][110]. First EPR studies on the dynamics of ECF transporters were presented recently [111,112], while the final mechanism is still under debate. Most EPR studies with a particular focus on pulsed methods, however, are available for ABC exporters [113][114][115][116][117][118][119][120][121][122][123][124] but their data are discussed as controversial as the present structures such that no final, common model has been agreed on.
Mostly, the origin of the power stroke for the transport is still under debate [125] and currently two models are postulated. In the first model, called the ATP-switch model, ATP binding results in the closure of the NBDs and the transition from the inward-facing (IF) to the outward-facing (OF) state [126]. Upon ATP hydrolysis the NBD interaction is destabilized, the substrate is released and the transporter returns to its IF state with disengaged NBDs. In this model, ATP hydrolysis is not required for substrate release but rather resets the exporter for the next ATP binding and substrate transport reaction. A series of crystal structures of ABC exporters in the absence and in the presence of nucleotides and pulsed EPR investigations support the ATP-switch model with few exceptions [127,128]. In the absence of nucleotides all ABC exporter structures exhibit IF states with partially or fully separated NBDs [118,128,129]. In the presence of nucleotides the NBDs are closed and the TMDs are in an outward-occluded or outward-facing conformation [130,131]. However, the ATP-switch model is challenged by the finding that the asymmetric occlusion of one nucleotide is important to trigger the affinity switch and thus substrate release [132,133]. Though, no structural evidence for the asymmetric nucleotide occlusion is published. Furthermore, a recent pulsed EPR study on the heterodimeric ABC exporter BmrCD stated that BmrCD requires ATP hydrolysis for complete NBD closure and TMD reorientation [113], which contradicts the ATP-switch model. A striking structural feature of BmrCD and other heterodimeric ABC exporters is the so-called degenerate ATP binding site in one NBD, which binds ATP tightly but is catalytically impaired [134,135]. The authors hypothesized that the power stroke for heterodimeric exporters could be different than for homodimeric exporters.
In the following paragraphs a couple of highlights of a recent EPR study on TM287/288 will be emphasized [119], dealing with the molecular trigger leading from an open to a closed NBD dimer with two nucleotides sandwiched at the interface. The study is an exquisite example for the combined use of pulsed EPR with different nucleotides in combination with inhibitors and functional mutations. As BmrCD, TM287/288 from Thermotoga maritima is a heterodimeric ABC exporter comprising one catalytically active and one degenerate nucleotide binding site. The study based on available inward-facing crystal structures in the presence of AMP-PNP-Mg and in the apo state, respectively, and on a former EPR study [118,136]. Other than in homodimeric exporters, in TM287/288 the NBDs remain in contact in the IF conformation. Moreover, in the presence of AMP-PNP-Mg the NBDs are not fully closed but a transition to the OF state is supported [118]. In the chosen study six spin-labeled TM287/288 variants are created based on the cys-less TM287/288. Two selected spin-labeled pairs each are located in the NBDs to follow interspin distance changes, detected at 50 K, at the consensus and degenerate ATPase sites, in the intracellular part of the TMDs and in the extracellular part of the TMDs, respectively ( Figure 6A). The expected interspin distances for these variants are predicted by means of using the software MMM on the apo structure of TM287/288 and a homology model based on Sav1866. As mentioned before, this tool has proven to be very helpful to identify suitable residues for the intended measurements if structural information is available. The functionality of all spin-labeled variants has been tested by an ATPase activity test.
For all spin-labeled variants the experimentally determined interspin distances in the apo state are in agreement with the predicted interspin distances for the IF conformation. In the following experimental setup various nucleotides and inhibitors are applied to stabilize distinct conformations. Upon the assumptions of the ATP-switch model the following conformations are to be expected: ATP-Mg in the presence of vanadate (ATP-Vi-Mg) and ATPγS-Mg should only be slowly hydrolysable probably favoring the OF state or at least promoting a significant shift towards the OF state. AMP-PNP in the presence of Mg as well as ATP-EDTA (thus absence of magnesium) are non-hydrolysable, thus ideally the transporter should be trapped in an all OF state. ATP-Mg should enable the complete transport cycle allowing all conformations, while ADP-Mg mimics the post-hydrolysis state, in which the IF conformation should be favored. In fact, the addition of the different nucleotides has led to very different dipolar evolution traces for all spin pairs tested as can be clearly seen from the changing oscillations ( Figure 6A left). By fitting those oscillations, either with a two-Gaussian fit DEER analysis as shown here or by Tikhonov regularization as used elsewhere in the paper and introduced previously in this review, distance distributions are derived, which allow for the interpretation of conformational changes in different areas of the protein. In fact, the addition of ATP-Vi-Mg ( Figure 6B, magenta lines) results in the most pronounced change towards the OF conformation. Compared to the apo state, the interspin distance decreases for the spin-labeled variants in the NBDs and in intracellular TMD. In contrast, spin-labeled pairs in the extracellular region exhibit an increase of the mean distances ( Figure 6B, magenta lines). In the same way, the effect of the other nucleotides could be analyzed. Surprisingly, both in the presence of ADP-Mg ( Figure 6B, orange lines) and AMP-PNP-Mg ( Figure 6B, purple lines), respectively, the IF state is stabilized. While this is expected for ADP-Mg, an OF state is hypothesized for AMP-PNP-Mg in agreement with the ATP-switch model. The presence of both, ATP-Mg ( Figure 6B, red lines) and ATP-γS ( Figure 6B, black lines), respectively, leads to a significant shift in the mean distances with a broad distance distribution indicating that the transporter has undergone a complete transport cycle sampling different conformations. With respect to the nucleotides and combinations discussed so far, one would hypothesize that ATP hydrolysis is in fact needed for the conformational change from the IF to the OF state as at least residual ATPase activity is given for all combinations leading to an IF to OF shift. Only the additional sampling with ATP-EDTA ( Figure 6B, green lines) refutes this conclusion as a significant shift towards the OF state is observed (though not 100%) without any ATPase activity. Thus, a conformational change is possible without ATP hydrolysis. broad distance distribution indicating that the transporter has undergone a complete transport cycle sampling different conformations. With respect to the nucleotides and combinations discussed so far, one would hypothesize that ATP hydrolysis is in fact needed for the conformational change from the IF to the OF state as at least residual ATPase activity is given for all combinations leading to an IF to OF shift. Only the additional sampling with ATP-EDTA ( Figure 6B, green lines) refutes this conclusion as a significant shift towards the OF state is observed (though not 100%) without any ATPase activity. Thus, a conformational change is possible without ATP hydrolysis.  To further confirm the ATP-switch model similar pulsed EPR measurements are applied to the catalytically inactive E-to-Q mutant, in which glutamate 517 in TM288 is mutated to a glutamine, which reduces the ATPase activity to almost zero and thus should allow to further stabilize the OF state with the corresponding nucleotides added. Such stabilization in fact is seen for ATP-Mg, ATP-γS and ATP-EDTA ( Figure 6C) supporting the hypothesis that ATP hydrolysis is not required for the IF to OF switch. Surprisingly, AMP-PNP-Mg remains incapable in switching the transporter to the OF conformation ( Figure 6C, purple lines), which is interpreted with the proteins' inability to sandwich two AMP-PNP as needed for the switch due to the degenerated site in the NBD of TM287. By ranking the nucleotides according to their ability to stabilize the OF state ( Figure 6D,E) ATP-Vi-Mg is identified as the best. In combination with the E-to-Q mutant this nucleotide probably has the highest potential to support the crystallization of TM287/288 in its OF state. In this respect EPR spectroscopy may help to further improve structural studies.
Finally, the authors extend their study to the heterodimeric BmrCD and the homodimeric MsbA and confirm the ATP-switch model with these two alternatives. However, it is very interesting to note that the energy landscapes for the different transporters are very different as can be derived from the distance distributions. Thus, EPR analyses allow for the interpretation of transporter dynamics instead of only providing structural snap shots. In summary, the discussed study has led to the suggestion of a comprehensive mechanism for heterodimeric ABC exporters, which likely applies to both TM287/288 and BmrCD as well as to other even human homologs and summarizes both differences and similarities with homodimeric exporters.

Secondary Active Proteins
Secondary active transporters form the largest and most diverse group of membrane proteins while transporting substrates across the membrane. In contrast to primary active transport, secondary active transport is not directly coupled to the hydrolysis of ATP. Instead, the transport of one substrate is coupled to the symport or antiport of another substrate such that at least one substrate is translocated along its electrochemical gradient, which can fuel the uphill transport of the second substrate. Thus, the transport of molecules across the membrane results from an electrochemical gradient which beforehand is generated by pumping ions in or out of the cell. A general concept for secondary active transport is the alternating access mechanism, which includes an outward-facing, an inward-facing and a substrate-occluded state. This concept ensures the strict coupling of substrate transport and hinders the dissipation of the electrochemical gradient. Based on their structural folds secondary transporters are grouped in different classes, with each of them having a more or less different underlying transport mechanism. Based on the structural information available the different mechanisms have been discussed recently by Forrest et al. [137]. In recent years, various EPR studies have complemented this structural knowledge, of which we will name only a few examples. As member of the major facilitator family (MFS) the H + -coupled multidrug antiporter LmrP has been studied extensively and is shown to function in the proposed rocker-switch mechanism [138]. In addition, the lipid dependency of the conformational-energy landscape has been evaluated [139]. The SLC6A family, which is better known as LeuT family, has been exemplified based on the Na + -coupled amino acid symporter LeuT [140,141]), the Na + /hydantoin symporter Mhp1 [142] and the Na + /proline symporter PutP [143][144][145][146]. All studies are in agreement with the suggested combination of the rocker-switch and gating-type mechanisms in order to achieve alternating access. However, they also highlight the different coupling schemes for different transporters in this family. The small multidrug resistance exporter family (SMR), which supports proton/substrate antiport, is represented by EmrE. Here in particular a study unraveling the protonation-dependent conformational dynamics has contributed significantly to the understanding of the underlying transport mechanism [147]. A combined effort of NMR, EPR and MD simulations shed light on global structural changes in the Cl − /H + exchanger family CLC, which are necessary to form the previously unknown outward-facing open conformation [148]. Two studies [149,150] investigated the conformational dynamics of the aspartate/Na + symporter Glt Ph , which is a member of the excitatory amino acid transporter family (EAAT), also referred to as SLC1. In the following section we will focus on the study by Hänelt et al. [149] to exemplify the combined use of cw and pulsed EPR at 160 K and 50 K, respectively, for the determination of distance distributions and dynamics and to highlight the conformational dependence of a protein on its environment.
Glt Ph from the thermophilic archaeon Pyrococcus horikoshii serves as a model protein for EAATs. In the central nervous system EAATs transport extracellular glutamate to the inside of the nerve cells and thereby remove excessive extracellular fluid and facilitate efficient synaptic communication between the neurons. Dysfunction of the EAATs may cause various human diseases, such as Alzheimer's disease, epilepsy and strokes [151,152]. In humans, there are five subtypes of glutamate transporters (EAAT1-5) [153]. EAATs couple the transport of glutamate or aspartate to the co-and counter-transport of different cations such as sodium ions, potassium ions or protons across the membrane. The type and the amount of the cations vary [154][155][156][157]. Glt Ph couples the transport of one aspartate to the co-transport of 3 Na + [155]. Structurally, Glt Ph forms a homotrimeric protein, in which each protomer is composed of two domains, the trimerization domain and the transporting domain. Moreover, characteristic structural features of each subunit are the two helical hairpins called HP1 and HP2, which probably act as molecular gates closing the substrate binding site and translocation path [158]. The trimerization domain is formed by four transmembrane helices (TM1, TM2, TM4 and TM5) and functions as a rigid scaffold [158,159]. Whereas the four transmembrane helices TM3, TM6, TM7 and TM8 constitute the transporting domain which binds the substrates and is proposed to undergo large elevator-like movements spanning the membrane diameter to expose the helical hairpins to the intracellular and extracellular space, respectively (cf. Figure 7A). These movements are suggested by various crystal structures covering outward-(o) and inward-(i) oriented states as well as intermediate states (PDB 1XFH, 2NWL, 3KBC and 3V8G, respectively).
To follow structural changes of Glt Ph upon substrate binding by EPR spectroscopy single and double labeled variants have been designed with spin-labeled residues in the trimerization domain and in the transporting domain ( Figure 7A). Subsequently, distance distribution of Glt Ph in detergent solution and reconstituted in liposomes are determined. Generally, single labeled variants have been applied to pulsed EPR experiments at low temperature (T = 50 K) with the apo-protein ( Figure 7B, black lines), the protein in the presence of saturating concentrations of Na + alone ( Figure 7B, red lines), or the presence of saturating concentrations of both Na + and aspartate ( Figure 7B, blue lines) to determine inter-protomer interspin distances. Therefore, two sets of single labeled variants have been selected: T166R1 and V176R1 in the trimerization domain and S300R1, S331R1 and A364R1 in the transporting domain, of which we will focus on variants V176R1 and S331R1 exemplarily.
For V176R1 in detergent solution one main interspin distance around 3.9 nm and a less populated interspin distance at 3.1 nm have been identified both with a very narrow distance distribution at all conditions ( Figure 7B, upper panel). A comparison of the experimental interspin distances of T176R1 with the ones of the rotamer library analysis gives evidence that the two different experimental distances result from two rotamers of the nitroxide side chain rather than from two different conformations of the protein. Thus, it has been concluded that the trimerization domain forms a stable scaffold.
In contrast, significant differences in the distance distributions are found for single labeled variant S331R1 located in the transporting domain upon substrate addition. Based on the crystal structures the mean interspin distance as derived by a rotamer library analysis for variant S331R1 is expected to increase when the transporting domain moves from an inward-to an outward-facing orientation. The experimentally determined distance distributions however display very broad distribution widths at all conditions indicating the conformational heterogeneity of the transporting domains ( Figure 7B, upper panel). The addition of substrates shifts the distance distribution towards longer distances and thus to the outward-facing conformation, nonetheless at any condition several conformations are populated which leads to hypothesize a flat energy landscape with a relatively small difference in free energy between the conformations. This hypothesis is supported by the observation that already the insertion of a spin label at different positions (S300R1 and A364R1, respectively) populates different conformations.
To analyze the potential effect of the membrane environment on the protein conformations and their relative population EPR experiments for selected liposome-reconstituted variants are conducted ( Figure 7B, lower panel). For variant V176R1 similar interspin distances are determined in the absence and in the presence of the substrates in detergent solution as well as reconstituted in liposomes. In contrast, different interspin distances are revealed in proteoliposomes compared to detergent solution for variant S331R1. The effects on the interspin distances upon addition of the substrates are qualitatively similar in detergent solution and reconstituted into liposomes. However, the extreme long and short distances are less pronounced and the population of intermediate distances increases. Furthermore, the relative distance changes are smaller in proteoliposomes compared to those in detergent solution indicating that the extreme conformations as found in the crystal structures are less occupied in the native membrane environment. To analyze the potential effect of the membrane environment on the protein conformations and their relative population EPR experiments for selected liposome-reconstituted variants are conducted ( Figure 7B, lower panel). For variant V176R1 similar interspin distances are determined in the absence and in the presence of the substrates in detergent solution as well as reconstituted in liposomes. In contrast, different interspin distances are revealed in proteoliposomes compared to detergent solution for variant S331R1. The effects on the interspin distances upon addition of the substrates are qualitatively similar in detergent solution and reconstituted into liposomes. However, the extreme long and short distances are less pronounced and the population of intermediate distances increases. Furthermore, the relative distance changes are smaller in proteoliposomes compared to those in detergent solution indicating that the extreme conformations as found in the crystal structures are less occupied in the native membrane environment. . The residues used for cysteine mutagenesis and spin labeling are color-coded: purple, Thr166; orange, Val176; blue, Ser300; yellow, Ser331; red, Ala364; cyan, Lys55; (B) Normalized background-corrected dipolar evolution data F(t); tick marks are separated by 0.1 (left) and distance distributions obtained by Tikhonov regularization and by a rotamer library analysis of the inward-facing (i, dotted lines) and outward-facing (o, dashed lines) crystal structures (right). DEER data of GltPh variants recorded at X band (9.4 GHz and 50 K; for mutant V176R1) or Q band (34 GHz and 50 K; for mutant S331R1) in detergent solution (top) and reconstituted in liposomes (bottom). Data are shown for the apo protein (black lines) and for the protein in the presence of Na + (NaCl; red lines) or Na + and aspartate (NaCl, Asp; blue lines); (C) Lowtemperature (T = 160 K) cw EPR spectra measured for the detergent-solubilized (top) or membranereconstituted (bottom) GltPh variant K55R1A364R1. Adapted from [149]. and in the inward-facing orientation (3KBC) (right). The residues used for cysteine mutagenesis and spin labeling are color-coded: purple, Thr166; orange, Val176; blue, Ser300; yellow, Ser331; red, Ala364; cyan, Lys55; (B) Normalized background-corrected dipolar evolution data F(t); tick marks are separated by 0.1 (left) and distance distributions obtained by Tikhonov regularization and by a rotamer library analysis of the inward-facing (i, dotted lines) and outward-facing (o, dashed lines) crystal structures (right). DEER data of Glt Ph variants recorded at X band (9.4 GHz and 50 K; for mutant V176R1) or Q band (34 GHz and 50 K; for mutant S331R1) in detergent solution (top) and reconstituted in liposomes (bottom). Data are shown for the apo protein (black lines) and for the protein in the presence of Na + (NaCl; red lines) or Na + and aspartate (NaCl, Asp; blue lines); (C) Low-temperature (T = 160 K) cw EPR spectra measured for the detergent-solubilized (top) or membrane-reconstituted (bottom) Glt Ph variant K55R1A364R1. Adapted from [149].
In addition to the pulsed EPR measurements low temperature cw EPR measurement have been performed to determine shorter intraprotomer distances based on dipolar broadening. For this, residues K55 in the trimerization domain and A364 in the transporting domain are mutated to cysteines and subsequently spin labeled. A previously published inward-facing crystal structure was gained by crosslinking these two residues [158]. In fact, Glt Ph in detergent solution shows broadened powder spectra at all conditions tested with an increased broadening upon substrate addition ( Figure 7C, upper spectra). Thus, an inward-facing conformation is favored upon substrate addition, which is in agreement with the pulsed EPR measurements on single labeled residue A364R1. Additionally, the cw EPR measurements very clearly display the conformational changes upon reconstitution into liposomes. In comparison to the detergent-solubilized sample the apo protein shows very strong spin-spin interaction indicating distances in the 1-nm range. The addition of substrates decreases the spectral broadening, indicating an increase in average interspin distances ( Figure 7C, lower spectra). Thus, substrate addition induces opposite distance changes compared to K55R1A364R1 in detergent solution, highlighting the role of the environment in the stabilization of particular conformations.
Based on the results of this cw and pulsed EPR study and existing crystal structures [158,160-162] a possible transport mechanism has been proposed in which the trimerization domain remains rigid during the transport process forming a stable core of trimeric Glt Ph while the transporting domains are constantly shuttling across the membrane. However, this shuttling likely is only a prerequisite for alternately exposing the binding site. Finally, local rearrangements like those shown for HP2 [158,163] may trigger the alternating accessibility of the binding sites by coordinately opening and closing the inward-and outward-facing gates in dependence of the presence of substrate and coupling ions.

New Developments
Due to the importance of biological membranes for membrane-protein structure, stability and function [164][165][166] EPR investigations of membrane proteins are commonly performed either in detergent micelles or in liposomes, though detergent micelles can be poor membrane mimics. For membrane proteins which are sensitive to their lipid surroundings [167,168] the potential of proteoliposomes lies in the possibility to monitor their function, to study structural details and to follow conformational transitions in a physiologically relevant environment. A disadvantage however is that membrane proteins can insert in two different orientations into the liposomes creating heterogeneity in the accessibility for e.g., ligands. Due to its orientation a membrane transporter might not respond to the addition of nucleotides or ligands to the proteoliposomes [101]. Hence, this heterogeneity might result in subpopulations contributing to the distance distribution in pulsed EPR experiments and therefore impedes the analysis and interpretation of the data. Over the last years the reconstitution of membrane proteins into nanodiscs evolved as a powerful alternative [140,169]. Nanodiscs are self-assembling systems comprising a lipid bilayer stabilized by modified apolipoproteins called membrane scaffold proteins (MSPs) [170] which form a belt-like antiparallel dimer [171,172]. The size of the nanodiscs, and so the area of the lipid bilayer, depends on the length of the MSP proteins. Next to the size of the nanodiscs the lipid composition is fully adaptable to the protein under investigation. Thus, nanodiscs are able to provide adequate lipid environments, support protein oligomeric states and allow access to both the intra-and extracellular part of the studied protein [173,174]. For instance the maltose transporter (MalFGK2) from E. coli [169] and the proton-powered multidrug transporter LmrP from Lactococcus lactis [140] have been successfully reconstituted in nanodiscs. For the latter study pulsed EPR measurements revealed concurrently two conformations of LmrP in detergent micelles and reconstituted in nanodiscs, respectively. However, the apparent pK value of the transition between the inward-and the outward-facing conformations is shifted by approximately two pH units for the reconstituted protein. This shift likely is caused by a proposed interaction between the lipid head groups and a conserved motif of charged residues stabilizing the inward-facing conformation. Thus, lipid interactions in more general may play a crucial role in modulating the proteins' activity which promotes the comprehension of membrane-inserted proteins in lipid bilayers.
However, the nanodisc approach is strictly restricted to membrane proteins whose function and structural dynamics are not based on membrane potentials. Therefore, before reconstitution into nanodiscs or liposomes for EPR investigations a thoroughly survey of the system under investigation and a good understanding of its physical properties are essential for selecting the best possible experimental conditions.
A very recent and promising new development in the field of EPR spectroscopy is the application of distance measurements in whole cells (for now Escherichia coli) and in native membranes [175,176]. The function of membrane proteins and their underlying mechanism may be effected by their surrounding environment emphasizing the necessity for investigations under native conditions. First pulsed EPR measurements have been performed on the cobalamin transporter BtuB acting as a model system for outer membrane transporter [175]. Crystal structures of BtuB in its apo, vitamin B 12 -bound [177] and TonB-bound [178] states are available. Despite this information the mechanism by which the substrate is released is still not known. A particular challenge working with in vivo systems is the spin labeling of reactive cysteines. Proteins in the outer membrane are suitable model proteins as they are generally free of reactive cysteines, thereby minimizing potential sources of non-specific signals [179,180]. For the first in vivo studies BtuB has been overexpressed and spin labeled in whole cells and outer membranes. Interspin distance measurements between the spin labeled BtuB and the spin labeled cobalamin reveal comparable interspin distances of measurements done in whole cells, outer membranes and synthetic vesicles demonstrating the successful application of in vivo pulsed EPR spectroscopy. However, further applications to inner membrane proteins or eukaryotic membranes might be challenging as higher background signals from endogenous cysteine-containing proteins are to be expected and reducing conditions are found inside the cytoplasm such that cysteine labeling by disulfide bonds is significantly reduced. An alternative coupling strategy to overcome this problem is the click chemistry in combination with unnatural amino acids using azide-and alkyne-functionalized spin labels. For this purpose the copper-catalyzed azide-alkyne 1,3-dipolar cycloaddition (CuAAC) or the copper-free, strain-promoted azide-alkyne cycloaddition (SPAAC) [181] can be implemented for spin labeling. A spin labeling protocol has been recently developed where unnatural amino acids are spin labeled with synthesized azide-or alkyne-functionalized spin labels [182]. The authors performed CuAAC spin labeling in living E. coli cells as a proof of principle for in cell spin labeling based on the click chemistry and first pulsed ERR experiments proofing the applicability of this approach for further EPR studies.
To further improve EPR studies and to render a more comprehensive view of membrane proteins and how alterations in their structures affect their function it is essential to extend the upper limit and accuracy of interspin distances that can be determined in membrane protein complexes. These method developments rely on an improved signal-to-noise ratio and increased modulation depth. In this context, a 7-pulse CP-PELDOR sequence has been introduced based on the Carr-Purcell (CP) decoupling scheme on the observer spin, a specific series of observer spin refocusing pulses [183,184]. This method allows a substantially longer observation time window extending the reliable and accurate interspin distance determination. The 7-pulse CP-PELDOR has been successfully applied to the trimeric betaine transporter BetP providing reliable interspin distance values around 6 nm which are not detectable with the commonly used 4-pulse PELDOR sequence [184]. These CP-PELDOR experiments result in a better determination of the intermolecular background contribution and an improved signal-to-noise ratio. Further developments in this direction might lead to a further increased sensitivity of distance measurements allowing the application of pulsed EPR spectroscopy on single molecules in the future.