Double Electron-electron Resonance Research Articles

Cross-validation of distance measurements in proteins by PELDOR/DEER and single-molecule FRET

Pulsed electron-electron double resonance spectroscopy (PELDOR/DEER) and single-molecule Förster resonance energy transfer spectroscopy (smFRET) are frequently used to determine conformational changes, structural heterogeneity, and inter probe distances in biological macromolecules. They provide qualitative information that facilitates mechanistic understanding of biochemical processes and quantitative data for structural modelling. To provide a comprehensive comparison of the accuracy of PELDOR/DEER and smFRET, we use a library of double cysteine variants of four proteins that undergo large-scale conformational changes upon ligand binding. With either method, we use established standard experimental protocols and data analysis routines to determine inter-probe distances in the presence and absence of ligands. The results are compared to distance predictions from structural models. Despite an overall satisfying and similar distance accuracy, some inconsistencies are identified, which we attribute to the use of cryoprotectants for PELDOR/DEER and label-protein interactions for smFRET. This large-scale cross-validation of PELDOR/DEER and smFRET highlights the strengths, weaknesses, and synergies of these two important and complementary tools in integrative structural biology.

The mechanism of damage recognition by apurinic/apyrimidinic endonuclease Nfo from Escherichia coli

Apurinic/apyrimidinic (AP) endonuclease Nfo from Escherichia coli recognises AP sites in DNA and catalyses phosphodiester bond cleavage on the 5′ side of AP sites and some damaged or undamaged nucleotides. Here, the mechanism of target nucleotide recognition by Nfo was analysed by pulsed electron–electron double resonance (PELDOR, also known as DEER) spectroscopy and pre–steady-state kinetic analysis with Förster resonance energy transfer detection of DNA conformational changes during DNA binding. The efficiency of endonucleolytic cleavage of a target nucleotide in model DNA substrates was ranked as (2R,3S)-2-(hydroxymethyl)-3-hydroxytetrahydrofuran [F-site] > 5,6-dihydro-2′-deoxyuridine > α-anomer of 2′-deoxyadenosine >2′-deoxyuridine > undamaged DNA. Real-time conformational changes of DNA during interaction with Nfo revealed an increase of distances between duplex ends during the formation of the initial enzyme–substrate complex. The use of rigid-linker spin-labelled DNA duplexes in DEER measurements indicated that double-helix bending and unwinding by the target nucleotide itself is one of the key factors responsible for indiscriminate recognition of a target nucleotide by Nfo. The results for the first time show that AP endonucleases from different structural families utilise a common strategy of damage recognition, which globally may be integrated with the mechanism of searching for specific sites in DNA by other enzymes.

Solution-State Inter-Copper Distribution of Redox Partner-Linked Copper Nitrite Reductases: A Pulsed Electron-Electron Double Resonance Spectroscopy Study.

Copper nitrite reductases (CuNiRs) catalyze the reduction of nitrite to form nitric oxide. In recent years, new classes of redox partner linked CuNiRs have been isolated and characterized by crystallographic techniques. Solution-state biophysical studies have shed light on the complex catalytic mechanisms of these enzymes and implied that protein dynamics may play a role in CuNiR catalysis. To investigate the structural, dynamical, and functional relationship of these CuNiRs, we have used protein reverse engineering and pulsed electron–electron double resonance (PELDOR) spectroscopy to determine their solution-state inter-copper distributions. Data show the multidimensional conformational landscape of this family of enzymes and the role of tethering in catalysis. The importance of combining high-resolution crystallographic techniques and low-resolution solution-state approaches in determining the structures and mechanisms of metalloenzymes is emphasized by our approach.

The effect of spin polarization on double electron-electron resonance (DEER) spectroscopy.

Double electron-electron resonance (DEER) spectroscopy measures the distribution of distances between two electron spins in the nanometer range, often on doubly spin-labeled proteins, via the modulation of a refocused spin echo by the dipolar interaction between the spins. DEER is commonly conducted under conditions where the polarization of the spins is small. Here, we examine the DEER signal under conditions of high spin polarization, thermally obtainable at low temperatures and high magnetic fields, and show that the signal acquires a polarization-dependent out-of-phase component both for the intramolecular and intermolecular contributions. For the latter, this corresponds to a phase shift of the spin echo that is linear in the pump pulse position. We derive a compact analytical form of this phase shift and show experimental measurements using monoradical and biradical nitroxides at several fields and temperatures. The effect highlights a novel aspect of the fundamental spin physics underlying DEER spectroscopy.

Quantitative Agreement between Conformational Substates of Holo Calcium-Loaded Calmodulin Detected by Double Electron–Electron Resonance EPR and Predicted by Molecular Dynamics Simulations

Calcium-loaded calmodulin (CaM/4Ca2+) comprises two domains that undergo rigid body reorientation from a predominantly extended conformation to a compact one upon binding target peptides. A recent replica-exchange molecular dynamics (MD) simulation on holo CaM/4Ca2+ suggested the existence of distinct structural clusters (substates) along the path from extended to compact conformers in the absence of substrates. Here, we experimentally demonstrate the existence of CaM/4Ca2+ substates trapped in local minima by three freezing/annealing regimes (slow, 40 s; intermediate, 1.5 s; fast, 0.5 ms) using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure interdomain distances between nitroxide spin-labels positioned at A17C and A128C in the N- and C-terminal domains, respectively. The DEER echo curves were directly fit to population-optimized P(r) pairwise distance distributions calculated from the coordinates of the MD clusters and compact crystal structure. DEER data on fully deuterated CaM/4Ca2+ were acquired at multiple values of the second echo period (10-35 μs) and analyzed globally to eliminate instrumental and overfitting artifacts and ensure accurate populations, peak positions, and widths. The DEER data for all three freezing regimes are quantitatively accounted for within experimental error by 5-6 distinct conformers comprising a predominantly populated extended form (60-75%) and progressively more compact states whose populations decrease as the degree of compactness increases. The shortest interdomain separation is found in the compact crystal structure, which has an occupancy of 4-6%. Thus, CaM/4Ca2+ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

OOP-ESEEM Spectroscopy: Accuracies of Distances of Spin-Correlated Radical Pairs in Biomolecules.

In addition to the commonly used electron–electron double resonance (ELDOR) technique, there are several other electron paramagnetic resonance (EPR) methods by which structure information can be obtained by exploiting the dipolar coupling between two radicals based on its characteristic r −3 dependence. In this contribution, we explore the potential of out-of-phase-electron-spin echo envelope modulation (OOP-ESEEM) spectroscopy to collect accurate distance information in photo-sensitive (bio) molecules. Although the method has already been applied to spin-correlated radical pairs in several classes of light-active proteins, the accuracy of the information obtained has not yet been extensively evaluated. To do this in a system-independent fashion, OOP-ESEEM time traces simulated with different values of the dipolar and exchange couplings were generated and analyzed in a best-possible way. Excellent agreement between calculated and numerically fitted values over a wide range of distances (between 15 and 45 Å) was obtained. Furthermore, the limitations of the method and the dependence on various experimental parameters could be evaluated.

Mechanism of Electron Spin Decoherence in a Partially Deuterated Glassy Matrix.

Long electron spin coherence lifetimes are essential for applications in quantum information science and electron paramagnetic resonance, for instance, for nanoscale distance measurements in biomolecular systems using double electron-electron resonance. We experimentally investigate the decoherence dynamics under the Hahn echo sequence of the organic radical d18-TEMPO in a variably deuterated frozen water:glycerol matrix. The coherence time (phase memory time) TM scales with proton concentration as [1H]-0.65. For selectively deuterated matrices, decoherence is accelerated in the presence of proton clustering, that is, substantial short-range density in the proton-proton radial distribution functions (<3 Å). Simulations using molecular dynamics and many-body spin quantum dynamics show excellent agreement with experiment and show that geminal proton pairs such as CH2 and OH2 groups are major decoherence drivers. This provides a predictive tool for designing molecular systems with long electron spin coherence times.

19F Electron-Nuclear Double Resonance Reveals Interaction between Redox-Active Tyrosines across the α/β Interface of E. coli Ribonucleotide Reductase.

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to deoxyribonucleotides, thereby playing a key role in DNA replication and repair. Escherichia coli class Ia RNR is an α2β2 enzyme complex that uses a reversible multistep radical transfer (RT) over 32 Å across its two subunits, α and β, to initiate, using its metallo-cofactor in β2, nucleotide reduction in α2. Each step is proposed to involve a distinct proton-coupled electron-transfer (PCET) process. An unresolved step is the RT involving Y356(β) and Y731(α) across the α/β interface. Using 2,3,5-F3Y122-β2 with 3,5-F2Y731-α2, GDP (substrate) and TTP (allosteric effector), a Y356• intermediate was trapped and its identity was verified by 263 GHz electron paramagnetic resonance (EPR) and 34 GHz pulse electron–electron double resonance spectroscopies. 94 GHz 19F electron-nuclear double resonance spectroscopy allowed measuring the interspin distances between Y356• and the 19F nuclei of 3,5-F2Y731 in this RNR mutant. Similar experiments with the double mutant E52Q/F3Y122-β2 were carried out for comparison to the recently published cryo-EM structure of a holo RNR complex. For both mutant combinations, the distance measurements reveal two conformations of 3,5-F2Y731. Remarkably, one conformation is consistent with 3,5-F2Y731 within the H-bond distance to Y356•, whereas the second one is consistent with the conformation observed in the cryo-EM structure. The observations unexpectedly suggest the possibility of a colinear PCET, in which electron and proton are transferred from the same donor to the same acceptor between Y356 and Y731. The results highlight the important role of state-of-the-art EPR spectroscopy to decipher this mechanism.

DEER Data Analysis Software: A Comparative Guide.

Pulsed dipolar electron paramagnetic resonance (PDEPR) spectroscopy experiments measure the dipolar coupling, and therefore nanometer-scale distances and distance distributions, between paramagnetic centers. Of the family of PDEPR experiments, the most commonly used pulsed sequence is four-pulse double electron resonance (DEER, also known as PELDOR). There are several ways to analyze DEER data to extract distance distributions, and this may appear overwhelming at first. This work compares and reviews six of the packages, and a brief getting started guide for each is provided.

Modeling of protein conformational changes with Rosetta guided by limited experimental data

Conformational changes are an essential component of functional cycles of many proteins, but their characterization often requires an integrative structural biology approach. Here, we introduce and benchmark ConfChangeMover (CCM), a new method built into the widely used macromolecular modeling suite Rosetta that is tailored to model conformational changes in proteins using sparse experimental data. CCM can rotate and translate secondary structural elements and modify their backbone dihedral angles in regions of interest. We benchmarked CCM on soluble and membrane proteins with simulated Cα-Cα distance restraints and sparse experimental double electron-electron resonance (DEER) restraints, respectively. In both benchmarks, CCM outperformed state-of-the-art Rosetta methods, showing that it can model a diverse array of conformational changes. In addition, the Rosetta framework allows a wide variety of experimental data to be integrated with CCM, thus extending its capability beyond DEER restraints. This method will contribute to the biophysical characterization of protein dynamics.

Functional Dynamics of Hsp70/CHIP/E2‐Ub Complex in Ubiquitination of Substrate Proteins

Heat shock protein (Hsp) 70 is the main mediator of protein quality control in the cytosol of eukaryotic cells. One of the major functions of Hsp70 is mediate the ubiquitination of client proteins which subsequently leads to proteasomal degradation. This process occurs via a multiprotein complex composed of Hsp70, carboxy terminus of Hsc70 interacting protein (CHIP), ubiquitin‐conjugating enzyme (E2), and ubiquitin (Ub). This complex is thought to be flexible, however, knowledge of the degree of dynamics are still lacking and detailed experimental evidence is needed. Here we present the models derived from Electron Paramagnetic Resonance (EPR) data and small‐angle X‐ray scattering (SAXS) ensembles for the complexes of Hsp70/CHIP/E2‐Ub with ATP or ADP. The model of ATP bound Hsp70/CHIP/E2‐Ub complex has an average distance of ~120 Å between the substrate‐binding domain (SBD) and the E2‐ubiquitin active site. The same distance of the multiprotein complex is found to be decreased to ~60 Å when Hsp70 is bound to ADP. Since the ADP bound conformation is considered as the active conformation that tightly binds client proteins, the alteration of the average distance between the Hsp70 SBD and the E2‐Ub active site implies that the Lys residue of the substrate might need to be positioned within appropriate distance to span the gap in order to efficiently promote ubiquitination. Here, we used molecular dynamic simulations, double electron‐electron resonance (DEER) EPR spectroscopy and in‐vitroubiquitination assays with designed model substrate proteins that place lysine residues at variable positions. We found the dynamics of the Hsp70/CHIP/E2‐Ubiquitin complex is playing a major role in the ubiquitination of various client proteins and these dynamics are necessary to overcome the distance barrier between E2‐Ub active site and the Hsp70 SBD. Our data also suggests that the Hsp70/CHIP complex is capable of ubiquitinating a client protein substrate when lysine residues are placed at several disparate positions. In sum, our data reveals the functional dynamics of the Hsp70/CHIP/E2‐Ub complex which enable ubiquitination of a client protein.

Compactness regularization in the analysis of dipolar EPR spectroscopy data

Dipolar electron paramagnetic resonance (EPR) experiments, such as double electron–electron resonance (DEER), measure distributions of nanometer-scale distances between paramagnetic centers, which are valuable for structural characterization of proteins and other macromolecular systems. One challenge in the least-squares fitting analysis of dipolar EPR data is the separation of the inter-molecular contribution (background) and the intra-molecular contribution. For noisy experimental traces of insufficient length, this separation is not unique, leading to identifiability problems for the background model parameters and the long-distance region of the intra-molecular distance distribution. Here, we introduce a regularization approach that mitigates this by including an additional penalty term in the objective function that is proportional to the variance of the distance distribution and thereby penalizes non-compact distributions. We examine the reliability of this approach statistically on a large set of synthetic data and illustrate it with an experimental example. The results show that the introduction of compactness can improve identifiability.

DEER experiments reveal fundamental differences between calmodulin complexes with IQ and MARCKS peptides in solution

Calmodulin (CaM) is a calcium-binding protein that regulates the function of many proteins by indirectly conferring Ca2+ sensitivity, and it undergoes a large conformational change on partners' binding. We compared the solution binding mode of the target peptides MARCKS and IQ by double electron-electron resonance (DEER) distance measurements and paramagnetic NMR. We combined nitroxide and Gd(III) spin labels, including specific substitution of one of the Ca2+ ions in the CaM mutant N60D by a Gd(III) ion. The binding of MARCKS to holo-CaM resulted neither in a closed conformation nor in a unique relative orientation between the two CaM domains, in contrast with the crystal structure. Binding of IQ to holo-CaM did generate a closed conformation. Using elastic network modeling and 12 distance restraints obtained from multiple holo-CaM/IQ DEER data, we derived a model of the solution structure, which is in reasonable agreement with the crystal structure.

Calculation of pulsed EPR DEER signal for two coupled Gd3+ ions by dipolar interaction using rotating frames

A method, based on the use of rotating frames, to either consider only the rotating frame of the observer spin, which is the main rotating frame, or to consider only the rotating frame of the pump spin, used only during the application of the pump pulse, has been exploited to calculate accurately the four-pulse EPR (Electron Paramagnetic Resonance) DEER (Double Electron-Electron Resonance) signal for two Gd3+ ions, coupled by the dipolar interaction in a biradical in a biological system, including the zero-field splitting of each Gd3+ ion. It has been illustrated here by calculating the four-pulse DEER signals at W-band for three values of the dipolar constant: ωdd=0.63MHz, ωdd=0.49MHz and ωdd=0.39MHz, corresponding to the range of typical Gd3+−Gd3+ distance (r) values, r=4.36 nm, 4.73 nm, 5.10 nm, respectively, in a sample of Gd-DOTA.

Enhancement of the creation yield of NV ensembles in a chemically vapour deposited diamond

In this work we investigate the properties of negatively charged nitrogen-vacancy (NV−) centres created during single crystal diamond growth by chemical vapour deposition (CVD) on [113]-oriented substrates and with N2O as a dopant gas. The use of spin echo and double electron-electron resonance (DEER) allows us to assess NV− ratio with respect to substitutional nitrogen impurities (Ns0) incorporated during growth, a critical figure of merit for quantum technologies. We demonstrate that, at moderate growth temperatures (800 °C), dense NV− ensembles of several hundreds of ppb (800 ppb for 50 ppm of added N2O) and with exceptionally high NV−/ Ns0 ratios of up to 25% can be achieved. This NV− creation yield is higher by at least an order of magnitude to that typically obtained on standard [100]-grown diamonds and comparable to the best values reported for electron-irradiated diamonds. The material obtained here thus advantageously combines a high NV− density, high creation yield, long coherence times of several tens of μs together with a partial preferential orientation. These are highly desirable requirements for diamond-based quantum sensors that may spur new developments with this crystalline orientation leading to improved performance and sensitivity.

Neural networks in pulsed dipolar spectroscopy: A practical guide

This is a methodological guide to the use of deep neural networks in the processing of pulsed dipolar spectroscopy (PDS) data encountered in structural biology, organic photovoltaics, photosynthesis research, and other domains featuring long-lived radical pairs and paramagnetic metal ions. PDS uses distance dependence of magnetic dipolar interactions; measuring a single well-defined distance is straightforward, but extracting distance distributions is a hard and mathematically ill-posed problem requiring careful regularisation and background fitting. Neural networks do this exceptionally well, but their “robust black box” reputation hides the complexity of their design and training – particularly when the training dataset is effectively infinite. The objective of this paper is to give insight into training against simulated databases, to discuss network architecture choices, to describe options for handling DEER (double electron-electron resonance) and RIDME (relaxation-induced dipolar modulation enhancement) experiments, and to provide a practical data processing flowchart.

Human CEACAM1 N-domain dimerization is independent from glycan modifications

Carcinoembryonic cellular adhesion molecules (CEACAMs) serve diverse roles in cell signaling, proliferation, and survival and are made up of one or several immunoglobulin (Ig)-like ectodomains glycosylated invivo. The physiological oligomeric state and how it contributes to protein function are central to understanding CEACAMs. Two putative dimer conformations involving different CEACAM1 N-terminal Ig-like domain (CCM1) protein faces (ABED and GFCC'C″) were identified from crystal structures. GFCC'C″ was identified as the dominant CCM1 solution dimer, but ambiguity regarding the effect of glycosylation on dimer formation calls its physiological relevance into question. We present the first crystal structure of minimally glycosylated CCM1 in the GFCC'C″ dimer conformation and characterization in solution by continuous-wave and double electron-electron resonance electron paramagnetic resonance spectroscopy. Our results suggest the GFCC'C″ dimer is dominant in solution with different levels of glycosylation, and structural conservation and co-evolved residues support that the GFCC'C″ dimer is conserved across CEACAMs.

Neutralizing antibodies induced in immunized macaques recognize the CD4-binding site on an occluded-open HIV-1 envelope trimer

Broadly-neutralizing antibodies (bNAbs) against HIV-1 Env can protect from infection. We characterize Ab1303 and Ab1573, heterologously-neutralizing CD4-binding site (CD4bs) antibodies, isolated from sequentially-immunized macaques. Ab1303/Ab1573 binding is observed only when Env trimers are not constrained in the closed, prefusion conformation. Fab-Env cryo-EM structures show that both antibodies recognize the CD4bs on Env trimer with an ‘occluded-open’ conformation between closed, as targeted by bNAbs, and fully-open, as recognized by CD4. The occluded-open Env trimer conformation includes outwardly-rotated gp120 subunits, but unlike CD4-bound Envs, does not exhibit V1V2 displacement, 4-stranded gp120 bridging sheet, or co-receptor binding site exposure. Inter-protomer distances within trimers measured by double electron-electron resonance spectroscopy suggest an equilibrium between occluded-open and closed Env conformations, consistent with Ab1303/Ab1573 binding stabilizing an existing conformation. Studies of Ab1303/Ab1573 demonstrate that CD4bs neutralizing antibodies that bind open Env trimers can be raised by immunization, thereby informing immunogen design and antibody therapeutic efforts.

In situ distance measurements in a membrane transporter using maleimide functionalized orthogonal spin labels and 5-pulse electron-electron double resonance spectroscopy

Spectroscopic investigation of membrane proteins in their native environment is a challenging task. Earlier we demonstrated the feasibility to measure precise distances within outer membrane proteins in E. coli and native membranes using methanethiosulfonate (MTS) functionalized labels combined with pulsed electron-electron double resonance spectroscopy. Here we show the application of maleimide functionalized Gd(III), nitroxide, and trityl labels for in situ distance measurement using the cobalamin transporter BtuB. These labels enabled distance measurements for BtuB in E. coli and native outer membranes and in the membranes maleimide-Gd-DOTA also is effective. Further, we show that the observable dipolar evolution time can be significantly prolonged in the native environments using the Carr-Purcell 5-pulse electron double resonance sequence. For a nitroxide-nitroxide pair, application of sech/tanh inversion pulses substantially suppressed the 4-pulse artifact at the Q-band frequency. In the case of a nitroxide-trityl pair, Gaussian pump pulses of varying amplitude are sufficient to suppress the artifact to the typical noise level. The feasibility of a range of bioresistant spin labels and the 5-pulse electron double resonance offers promising tools for investigating heterooligomeric membrane protein complexes in their native environment.

Accessing protein conformational distributions in mammalian cells with fast, bio-orthogonal spin-labeling and DEER spectroscopy

The ability to measure protein structure and dynamics directly within the cellular environment is fundamental to understanding the molecular mechanisms of protein function and dysfunction. Site-directed spin-labeling, in combination with double electron-electron resonance (DEER) spectroscopy, is a powerful method for determining both the structural states and the conformational equilibria of biomacromolecules. The application of DEER to expressed proteins within cells, however, has been hindered by issues of sensitivity, toxicity, specificity and extent of labeling, and spin-label stability.