Benchmark Test and Guidelines for DEER/PELDOR Experiments on Nitroxide-Labeled Biomolecules.

Distance distribution information obtained by pulsed dipolar EPR spectroscopy provides an important contribution to many studies in structural biology. Increasingly, such information is used in integrative structural modeling, where it delivers unique restraints on the width of conformational ensembles. In order to ensure reliability of the structural models and of biological conclusions, we herein define quality standards for sample preparation and characterization, for measurements of distributed dipole-dipole couplings between paramagnetic labels, for conversion of the primary time-domain data into distance distributions, for interpreting these distributions, and for reporting results. These guidelines are substantiated by a multi-laboratory benchmark study and by analysis of data sets with known distance distribution ground truth. The study and the guidelines focus on proteins labeled with nitroxides and on double electron-electron resonance (DEER aka PELDOR) measurements and provide suggestions on how to proceed analogously in other cases.

Confidence Analysis of DEER Data and Its Structural Interpretation with Ensemble-Biased Metadynamics.

Given its ability to measure multicomponent distance distributions between electron-spin probes, double electron-electron resonance (DEER) spectroscopy has become a leading technique to assess the structural dynamics of biomolecules. However, methodologies to evaluate the statistical error of these distributions are not standard, often hampering a rigorous interpretation of the experimental results. Distance distributions are often determined from the experimental DEER data through a mathematical method known as Tikhonov regularization, but this approach makes rigorous error estimates difficult. Here, we build upon an alternative, model-based approach in which the distance probability distribution is represented as a sum of Gaussian components, and use propagation of errors to calculate an associated confidence band. Our approach considers all sources of uncertainty, including the experimental noise, the uncertainty in the fitted background signal, and the limited time span of the data collection. The resulting confidence band reveals the most and least reliable features of the probability distribution, thereby informing the structural interpretation of DEER experiments. To facilitate this interpretation, we also generalize the molecular simulation method known as ensemble-biased metadynamics (EBMetaD). This method, originally designed to generate maximal-entropy structural ensembles consistent with one or more probability distributions, now also accounts for the uncertainty in those target distributions exactly as dictated by their confidence bands. After careful benchmarks, we demonstrate the proposed techniques using DEER results from spin-labeled T4 lysozyme.

Long Distance Measurements up to 160†Šin the GroEL Tetradecamer Using Q-Band DEER EPR Spectroscopy.

Current distance measurements between spin-labels on multimeric protonated proteins using double electron-electron resonance (DEER) EPR spectroscopy are generally limited to the 15-60 Šrange. Here we show how DEER experiments can be extended to dipolar evolution times of ca. 80 µs, permitting distances up to 170 Što be accessed in multimeric proteins. The method relies on sparse spin-labeling, supplemented by deuteration of protein and solvent, to minimize the deleterious impact of multispin effects and substantially increase the apparent spin-label phase memory relaxation time, complemented by high sensitivity afforded by measurements at Q-band. We demonstrate the approach using the tetradecameric molecular machine GroEL as an example. Two engineered surface-exposed mutants, R268C and E315C, are used to measure pairwise distance distributions with mean values ranging from 20 to 100 Šand from 30 to 160 Å, respectively, both within and between the two heptameric rings of GroEL. The measured distance distributions are consistent with the known crystal structure of apo GroEL. The methodology presented here should significantly expand the use of DEER for the structural characterization of conformational changes in higher order oligomers.

Structural arrangement of the intracellular Ca2+ binding domains of the cardiac Na+/Ca2+ exchanger (NCX1.1): effects of Ca2+ binding.

The cardiac Na(+)/Ca(2+) exchanger (NCX1.1) serves as the primary means of Ca(2+) extrusion across the plasma membrane of cardiomyocytes after the rise in intracellular Ca(2+) during contraction. The exchanger is regulated by binding of Ca(2+) to its intracellular domain, which contains two structurally homologous Ca(2+) binding domains denoted as CBD1 and CBD2. NMR and x-ray crystallographic studies have provided structures for the isolated CBD1 and CBD2 domains and have shown how Ca(2+) binding affects their structures and motional dynamics. However, structural information on the entire Ca(2+) binding domain, denoted CBD12, and how binding of Ca(2+) alters its structure and dynamics is more limited. Site-directed spin labeling has been employed in this work to address these questions. Electron paramagnetic resonance measurements on singly labeled constructs of CBD12 have identified the regions that undergo changes in dynamics as a result of Ca(2+) binding. Double electron-electron resonance (DEER) measurements on doubly labeled constructs of CBD12 have shown that the ß-sandwich regions of the CBD1 and CBD2 domains are largely insensitive to Ca(2+) binding and that these two domains are widely separated at their N and C termini. Interdomain distances measured by DEER have been employed to construct structural models for CBD12 in the presence and absence of Ca(2+). These models show that there is not a major change in the relative orientation of the two Ca(2+) binding domains as a result of Ca(2+) binding in the NCX1.1 isoform. Additional measurements have shown that there are significant changes in the dynamics of the F-G loop region of CBD2 that merit further characterization with regard to their possible involvement in regulation of NCX1.1 activity.

Competition between homodimerization and cholesterol binding to the C99 domain of the amyloid precursor protein.

The 99-residue transmembrane C-terminal domain (C99, also known as ß-CTF) of the amyloid precursor protein (APP) is the product of the ß-secretase cleavage of the full-length APP and is the substrate for γ-secretase cleavage. The latter cleavage releases the amyloid-ß polypeptides that are closely associated with Alzheimer's disease. C99 is thought to form homodimers; however, the free energy in favor of dimerization has not previously been quantitated. It was also recently documented that cholesterol forms a 1:1 complex with monomeric C99 in bicelles. Here, the affinities for both homodimerization and cholesterol binding to C99 were measured in bilayered lipid vesicles using both electron paramagnetic resonance (EPR) and Förster resonance energy transfer (FRET) methods. Homodimerization and cholesterol binding were seen to be competitive processes that center on the transmembrane G700XXXG704XXXG708 glycine-zipper motif and adjacent Gly709. On one hand, the observed Kd for cholesterol binding (Kd = 2.7 ± 0.3 mol %) is on the low end of the physiological cholesterol concentration range in mammalian cell membranes. On the other hand, the observed K(d) for homodimerization (K(d) = 0.47 ± 0.15 mol %) likely exceeds the physiological concentration range for C99. These results suggest that the 1:1 cholesterol/C99 complex will be more highly populated than C99 homodimers under most physiological conditions. These observations are of relevance for understanding the γ-secretase cleavage of C99.

Determination of structural models of the complex between the cytoplasmic domain of erythrocyte band 3 and ankyrin-R repeats 13-24.

The adaptor protein ankyrin-R interacts via its membrane binding domain with the cytoplasmic domain of the anion exchange protein (AE1) and via its spectrin binding domain with the spectrin-based membrane skeleton in human erythrocytes. This set of interactions provides a bridge between the lipid bilayer and the membrane skeleton, thereby stabilizing the membrane. Crystal structures for the dimeric cytoplasmic domain of AE1 (cdb3) and for a 12-ankyrin repeat segment (repeats 13-24) from the membrane binding domain of ankyrin-R (AnkD34) have been reported. However, structural data on how these proteins assemble to form a stable complex have not been reported. In the current studies, site-directed spin labeling, in combination with electron paramagnetic resonance (EPR) and double electron-electron resonance, has been utilized to map the binding interfaces of the two proteins in the complex and to obtain inter-protein distance constraints. These data have been utilized to construct a family of structural models that are consistent with the full range of experimental data. These models indicate that an extensive area on the peripheral domain of cdb3 binds to ankyrin repeats 18-20 on the top loop surface of AnkD34 primarily through hydrophobic interactions. This is a previously uncharacterized surface for binding of cdb3 to AnkD34. Because a second dimer of cdb3 is known to bind to ankyrin repeats 7-12 of the membrane binding domain of ankyrin-R, the current models have significant implications regarding the structural nature of a tetrameric form of AE1 that is hypothesized to be involved in binding to full-length ankyrin-R in the erythrocyte membrane.

Protein functional dynamics from the rigorous global analysis of DEER data: Conditions, components, and conformations.

The potential of spin labeling to reveal the dynamic dimension of macromolecules has been recognized since the dawn of the methodology in the 1960s. However, it was the development of pulsed electron paramagnetic resonance spectroscopy to detect dipolar coupling between spin labels and the availability of turnkey instrumentation in the 21st century that realized the full promise of spin labeling. Double electron-electron resonance (DEER) spectroscopy has seen widespread applications to channels, transporters, and receptors. In these studies, distance distributions between pairs of spin labels obtained under different biochemical conditions report the conformational states of macromolecules, illuminating the key movements underlying biological function. These experimental studies have spurred the development of methods for the rigorous analysis of DEER spectroscopic data along with methods for integrating these distributions into structural models. In this tutorial, we describe a model-based approach to obtaining a minimum set of components of the distance distribution that correspond to functionally relevant protein conformations with a set of fractional amplitudes that define the equilibrium between these conformations. Importantly, we review and elaborate on the error analysis reflecting the uncertainty in the various parameters, a critical step in rigorous structural interpretation of the spectroscopic data.

Conformational equilibria of bulged sites in duplex DNA studied by EPR spectroscopy.

Conformational flexibility in nucleic acids provides a basis for complex structures, binding, and signaling. One-base bulges directly neighboring single-base mismatches in nucleic acids can be present in a minimum of two distinct conformations, complicating the examination of the thermodynamics by calorimetry or UV-monitored melting techniques. To provide additional information about such structures, we demonstrate how electron paramagnetic resonance (EPR) active spin-labeled base analogues, base-specifically incorporated into the DNA, are monitors of the superposition of different bulge-mismatch conformations. EPR spectra provide information about the dynamic environments of the probe. This information is cast in terms of "dynamic signatures" that have an underlying basis in structural variations. By examining the changes in the equilibrium of the different states across a range of temperatures, the enthalpy and entropy of the interconversion among possible conformations can be determined. The DNA constructs with a single bulge neighboring a single-base mismatch ("bulge-mismatches") may be approximately modeled as an equilibrium between two possible conformations. This structural information provides insight into the local composition of the bulge-mismatch sequences. Experiments on the bulge-mismatches show that basepairing across the helix can be understood in terms of purine and pyrimidine interactions, rather than specific bases. Measurements of the enthalpy and entropy of formation for the bulge-mismatches by differential scanning calorimetry and UV-monitored melting confirm that the formation of bulge-mismatches is in fact more complicated than a simple two-state process, consistent with the base-specific spectral data that bulge-mismatches exist in multiple conformations in the premelting temperature region. We find that the calculations with the nearest-neighbor (NN) model for the two likely conformations do not correlate well with the populations of structures and thermodynamic parameters inferred from the base-specific EPR dynamics probe. We report that the base-specific spin probes are able to identify a bistable, temperature dependent, switching between conformations for a particular complex bulged construct.

CW EPR and DEER Methods to Determine BCL-2 Family Protein Structure and Interactions: Application of Site-Directed Spin Labeling to BAK Apoptotic Pores.

The continuous wave (CW) and pulse electron paramagnetic resonance (EPR) methods enable the measurement of distances between spin-labeled residues in biopolymers including proteins, providing structural information. Here we describe the CW EPR deconvolution/convolution method and the four-pulse double electron-electron resonance (DEER) approach for distance determination, which were applied to elucidate the organization of the BAK apoptotic pores formed in the lipid bilayers.

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

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

Assembly of Bak homodimers into higher order homooligomers in the mitochondrial apoptotic pore.

In mitochondrial apoptosis, Bak is activated by death signals to form pores of unknown structure on the mitochondrial outer membrane via homooligomerization. Cytochrome c and other apoptotic factors are released from the intermembrane space through these pores, initiating downstream apoptosis events. Using chemical crosslinking and double electron electron resonance (DEER)-derived distance measurements between specific structural elements in Bak, here we clarify how the Bak pore is assembled. We propose that previously described BH3-in-groove homodimers (BGH) are juxtaposed via the 'α3/α5' interface, in which the C-termini of helices α3 and α5 are in close proximity between two neighboring Bak homodimers. This interface is observed concomitantly with the well-known 'α6:α6' interface. We also mapped the contacts between Bak homodimers and the lipid bilayer based on EPR spectroscopy topology studies. Our results suggest a model for the lipidic Bak pore, whereby the mitochondrial targeting C-terminal helix does not change topology to accommodate the lining of the pore lumen by BGH.

A Straightforward Approach to the Analysis of Double Electron-Electron Resonance Data.

Double electron-electron resonance (DEER) is now widely utilized to measure distance distributions in the 20-70Å range. DEER is frequently applied to biological systems that have multiple conformational states leading to complex distance distributions. These complex distributions raise issues regarding the best approach to analyze DEER data. A widely used method utilizes a priori background correction followed by Tikhonov regularization. Unfortunately, the underlying assumptions of this approach can impact the analysis. In this chapter, a method of analyzing DEER data is presented that is ideally suited to obtain these complex distance distributions. The approach allows the fitting of raw experimental data without a priori background correction as well as the rigorous determination of uncertainties for all fitting parameters. This same methodological approach can be used for the simultaneous or global analysis of multiple DEER data sets using variable ratios of a common set of components, thus allowing direct correlation of distance components with functionally relevant conformational and biochemical states. Examples are given throughout to highlight this robust fitting approach.

Formulation of Zeeman modulation as a signal filter.

The Bloch equation containing a Zeeman modulation field is solved analytically by treating the Zeeman modulation frequency as a perturbation. The absorption and dispersion signals at both 0 degrees and 90 degrees modulation phase are obtained. The solutions are valid to first order in the modulation frequency, but are otherwise valid for any value of modulation amplitude or microwave amplitude. A first order treatment of modulation frequency is shown to be a valid approximation over a wide range of typical experimental EPR conditions. The solutions derived from the Bloch equation suggest that the effect of over-modulation on first and second harmonic EPR spectra can be formulated as a mathematical filter that smoothes and broadens the under-modulated signal. The only adjustable filter parameter is a width that is equivalent to the applied peak-to-peak modulation amplitude. The true spin-spin and spin-lattice relaxation rates are completely determined from the under-modulated spectrum. The filters derived from the analytic solutions of the Bloch equation in the linear limit of modulation frequency are tested against numerical solutions of the Bloch equation that are valid for any modulation frequency to show their applicability. The filters are further tested using experimental EPR spectra. Experimental under-modulated spectra are mathematically filtered and compared with the experimental over-modulated spectra. The application of modulation filters to STEPR spectra is explored and limitations are discussed.

CW dipolar broadening EPR spectroscopy and mechanically aligned bilayers used to measure distance and relative orientation between two TOAC spin labels on an antimicrobial peptide.

An EPR membrane alignment technique was applied to measure distance and relative orientations between two spin labels on a protein oriented along the surface of the membrane. Previously we demonstrated an EPR membrane alignment technique for measuring distances and relative orientations between two spin labels using a dual TOAC-labeled integral transmembrane peptide (M2δ segment of Acetylcholine receptor) as a test system. In this study we further utilized this technique and successfully measured the distance and relative orientations between two spin labels on a membrane peripheral peptide (antimicrobial peptide magainin-2). The TOAC-labeled magainin-2 peptides were mechanically aligned using DMPC lipids on a planar quartz support, and CW-EPR spectra were recorded at specific orientations. Global analysis in combination with rigorous spectral simulation was used to simultaneously analyze data from two different sample orientations for both single- and double-labeled peptides. We measured an internitroxide distance of 15.3Šfrom a dual TOAC-labeled magainin-2 peptide at positions 8 and 14 that closely matches with the 13.3Šdistance obtained from a model of the labeled magainin peptide. In addition, the angles determining the relative orientations of the two nitroxides have been determined, and the results compare favorably with molecular modeling. This study demonstrates the utility of the technique for proteins oriented along the surface of the membrane in addition to the previous results for proteins situated within the membrane bilayer.

Automated structure refinement for a protein heterodimer complex using limited EPR spectroscopic data and a rigid-body docking algorithm: a three-dimensional model for an ankyrin-CDB3 complex.

We report here specialized functions incorporated recently in the rigid-body docking software toolkit TagDock to utilize electron paramagnetic resonance derived (EPR-derived) interresidue distance measurements and spin-label accessibility data. The TagDock package extensions include a custom methanethiosulfonate spin label rotamer library to enable explicit, all-atom spin-label side-chain modeling and scripts to evaluate spin-label surface accessibility. These software enhancements enable us to better utilize the biophysical data routinely available from various spin-labeling experiments. To illustrate the power and utility of these tools, we report the refinement of an ankyrin:CDB3 complex model that exhibits much improved agreement with the EPR distance measurements, compared to model structures published previously.

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

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

The global analysis of DEER data.

Double Electron-Electron Resonance (DEER) has emerged as a powerful technique for measuring long range distances and distance distributions between paramagnetic centers in biomolecules. This information can then be used to characterize functionally relevant structural and dynamic properties of biological molecules and their macromolecular assemblies. Approaches have been developed for analyzing experimental data from standard four-pulse DEER experiments to extract distance distributions. However, these methods typically use an a priori baseline correction to account for background signals. In the current work an approach is described for direct fitting of the DEER signal using a model for the distance distribution which permits a rigorous error analysis of the fitting parameters. Moreover, this approach does not require a priori background correction of the experimental data and can take into account excluded volume effects on the background signal when necessary. The global analysis of multiple DEER data sets is also demonstrated. Global analysis has the potential to provide new capabilities for extracting distance distributions and additional structural parameters in a wide range of studies.

The amyloid precursor protein has a flexible transmembrane domain and binds cholesterol.

C99 is the transmembrane carboxyl-terminal domain of the amyloid precursor protein that is cleaved by γ-secretase to release the amyloid-ß polypeptides, which are associated with Alzheimer's disease. Nuclear magnetic resonance and electron paramagnetic resonance spectroscopy show that the extracellular amino terminus of C99 includes a surface-embedded "N-helix" followed by a short "N-loop" connecting to the transmembrane domain (TMD). The TMD is a flexibly curved α helix, making it well suited for processive cleavage by γ-secretase. Titration of C99 reveals a binding site for cholesterol, providing mechanistic insight into how cholesterol promotes amyloidogenesis. Membrane-buried GXXXG motifs (G, Gly; X, any amino acid), which have an established role in oligomerization, were also shown to play a key role in cholesterol binding. The structure and cholesterol binding properties of C99 may aid in the design of Alzheimer's therapeutics.

Distance measurements on a dual-labeled TOAC AChR M2δ peptide in mechanically aligned DMPC bilayers via dipolar broadening CW-EPR spectroscopy.

A membrane alignment technique has been used to measure the distance between two TOAC nitroxide spin labels on the membrane-spanning M2δ, peptide of the nicotinic acetylcholine receptor (AChR), via CW-EPR spectroscopy. The TOAC-labeled M2δ peptides were mechanically aligned using DMPC lipids on a planar quartz support, and CW-EPR spectra were recorded at specific orientations. Global analysis in combination with rigorous spectral simulation was used to simultaneously analyze data from two different sample orientations for both single- and double-labeled peptides. We measured an internitroxide distance of 14.6 Šfrom a dual TOAC-labeled AChR M2δ peptide at positions 7 and 13 that closely matches with the 14.5 Šdistance obtained from a model of the labeled AChR M2δ peptide. In addition, the angles determining the relative orientation of the two nitroxides have been determined, and the results compare favorably with molecular modeling. The global analysis of the data from the aligned samples gives much more precise estimates of the parameters defining the geometry of the two labels than can be obtained from a randomly dispersed sample.

Structure of the cytoplasmic domain of erythrocyte band 3 hereditary spherocytosis variant P327R: band 3 Tuscaloosa.

Previous studies have shown that a single P327R point mutation in the cytoplasmic domain of band 3 (cdb3) protein, known as band 3 Tuscaloosa, leads to a reduction in protein 4.2 content of the erythrocyte membrane and hemolytic anemia. Recent studies have shown that this point mutation does not dissociate the cdb3 dimer, nor does it lead to large-scale rearrangement of the protein structure (Bustos, S. P., and Reithmeier, R. A. F. (2006) Biochemistry 45, 1026-1034). To better define the structural changes in cdb3 that lead to the hemolytic anemia phenotype, site-directed spin labeling (SDSL), in combination with continuous wave electron paramagnetic resonance (EPR) and pulsed double electron-electron resonance (DEER) spectroscopies, has been employed in this study to compare the structure of the R327 variant with wild type P327 cdb3. It is confirmed that the P327R mutation does not dissociate the cdb3 dimer, nor does it change the spatial orientation of the two peripheral domains relative to the dimer interface. However, it does affect the packing of the C-terminal end of helix 10 of the dimerization arms in a subpopulation of cdb3 dimers, it leads to spectral changes at some residues in beta-strand 11 and in the N-terminal end of helix10, and it produces measurable spectral changes at other residues that are near the mutation site. The data indicate that the structural changes are subtle and are localized to one surface of the cdb3 dimer. The spectroscopic description of structural features of the P327R variant provides important clues about the location of one potential protein 4.2 binding surface on cdb3 as well as new insight into the structural basis of the membrane destabilization.