Double Electron-electron Resonance Spectroscopy Research Articles

The Trajectory of Damaged-Base Eversion into the Active Site of Apurinic/Apyrimidinic Endonuclease APE1 Regulates This Enzyme's Substrate Specificity.

Apurinic/apyrimidinic endonuclease 1 (APE1) is responsible for the hydrolysis of the phosphodiester bond on the 5' side of an apurinic/apyrimidinic site during base excision repair. Moreover, in DNA, this enzyme can recognize nucleotides containing such damaged bases as 5,6-dihydro-2'-deoxyuridine (DHU), 2'-deoxyuridine (dU), alpha-2'-deoxyadenosine (αA), and 1,N6-ethenoadenosine (εA). Previously, by pulsed electron-electron double resonance spectroscopy and pre-steady-state kinetic analysis, we have revealed multistep DNA rearrangements during the formation of the catalytic complex. In the present study, the modeling of the eversion trajectory of nucleotides with various damaged bases was performed by directed molecular dynamics simulations. It was found that each damaged base at the beginning of the eversion interacts with protein loop Val196-Arg201, which should be moved to enable further nucleotide eversion. This movement involves a shift in loop Val196-Arg201 away from loop Asn253-Thr257 and requires the disruption of contacts between these loops. The Glu260Ala substitution facilitates the separation of the two loops. Moreover, conformational changes in the Asn253-Thr257 loop should occur in the second half of the lesion eversion trajectory. All these perturbations within the protein globule tend to reduce steric interactions of each damaged base with the protein during the eversion of the nucleotide from DNA and movement to the active site. These perturbations are important determinants of substrate specificity of endonuclease APE1.

Bayesian Probabilistic Inference of Nonparametric Distance Distributions in DEER Spectroscopy.

Double electron-electron resonance (DEER) spectroscopy measures distance distributions between spin labels in proteins, yielding important structural and energetic information about conformational landscapes. Analysis of an experimental DEER signal in terms of a distance distribution is a nontrivial task due to the ill-posed nature of the underlying mathematical inversion problem. This work introduces a Bayesian probabilistic inference approach to analyze DEER data, assuming a nonparametric distance distribution with a Tikhonov smoothness prior. The method uses Markov Chain Monte Carlo sampling with a compositional Gibbs sampler to determine a posterior probability distribution over the entire parameter space, including the distance distribution, given an experimental data set. This posterior contains all of the information available from the data, including a full quantification of the uncertainty about the model parameters. The corresponding uncertainty about the distance distribution is visually captured via an ensemble of posterior predictive distributions. Several examples are presented to illustrate the method. Compared with bootstrapping, it performs faster and provides slightly larger uncertainty intervals.

In-Cell DEER Spectroscopy of Nanodisc-Delivered Membrane Proteins in Living Cell Membranes.

Membrane proteins are integral to numerous cellular processes, yet their conformational dynamics in native environments remains difficult to study. This study introduces a nanodelivery method using nanodiscs to transport spin-labeled membrane proteins into the membranes of living cells, enabling direct in-cell double electron-electron resonance (DEER) spectroscopy measurements. We investigated the membrane protein BsYetJ, incorporating spin labels at key positions to monitor conformational changes. Our findings demonstrate successful delivery and high-quality DEER data for BsYetJ in both Gram-negative E. coli and Gram-positive B. subtilis membranes. The delivered BsYetJ retains its ability to transport calcium ions. DEER analysis reveals distinct conformational states of BsYetJ in different membrane environments, highlighting the influence of lipid composition on the protein structure. This nanodelivery method overcomes traditional limitations, enabling the study of membrane proteins in more physiologically relevant conditions.

Alternating access of a bacterial homolog of neurotransmitter: sodium symporters determined from AlphaFold2 ensembles and DEER spectroscopy

Neurotransmitter:sodium symporters (NSSs) play critical roles in neural signaling by regulating neurotransmitter uptake into cells powered by sodium electrochemical gradients. Bacterial NSSs orthologs, including MhsT from Bacillus halodurans, have emerged as model systems to understand the structural motifs of alternating access in NSSs and the extent of conservation of these motifs across the family. Here, we apply a computational/experimental methodology to illuminate the conformational landscape of MhsT alternating access. Capitalizing on our recently developed method, Sampling Protein Ensembles and Conformational Heterogeneity with AlphaFold2 (SPEACH_AF), we derived clusters of MhsT models spanning the transition from inward-facing to outward-facing conformations. Systematic application of double electron-electron resonance (DEER) spectroscopy revealed ligand-dependent movements of multiple structural motifs that underpin MhsT's conformational cycle. Remarkably, comparative DEER analysis in detergent micelles and lipid nanodiscs highlights the profound effect of the environment on the energetics of conformational changes. Through experimentally derived selection of collective variables, we present a model of ion and substrate-powered transport by MhsT consistent with the conformational cycle derived from DEER. Our findings not only advance the understanding of MhsT's function but also uncover motifs of conformational dynamics conserved within the broader context of the NSS family and within the LeuT-fold class of transporters. Importantly, our methodological blueprint introduces an approach that can be applied across a diverse spectrum of transporters to describe their conformational landscapes.

Pulsed Dipolar EPR for Self-Limited Complexes of Oligonucleotides Studies.

Pulsed electron-electron double resonance (PELDOR) spectroscopy is a powerful method for determining nucleic acid (NA) structure and conformational dynamics. PELDOR with molecular dynamics (MD) simulations opens up unique possibilities for defining the conformational ensembles of flexible, three-dimensional, self-assembled complexes of NA. Understanding the diversity and structure of these complexes is vital for uncovering matrix and regulative biological processes in the human body and artificially influencing them for therapeutic purposes. To explore the reliability of PELDOR and MD simulations, we site-specifically attached nitroxide spin labels to oligonucleotides, which form self-assembled complexes between NA chains and exhibit significant conformational flexibility. The DNA complexes assembled from a pair of oligonucleotides with different linker sizes showed excellent agreement between the distance distributions obtained from PELDOR and calculated from MD simulations, both for the mean inter-spin distance and the distance distribution width. These results prove that PELDOR with MD simulations has significant potential for studying the structure and dynamics of conformational flexible complexes of NA.

Partial wrapping of single-stranded DNA by replication protein A and modulation through phosphorylation.

Single-stranded DNA (ssDNA) intermediates which emerge during DNA metabolic processes are shielded by replication protein A (RPA). RPA binds to ssDNA and acts as a gatekeeper to direct the ssDNA towards downstream DNA metabolic pathways with exceptional specificity. Understanding the mechanistic basis for such RPA-dependent functional specificity requires knowledge of the structural conformation of ssDNA when RPA-bound. Previous studies suggested a stretching of ssDNA by RPA. However, structural investigations uncovered a partial wrapping of ssDNA around RPA. Therefore, to reconcile the models, in this study, we measured the end-to-end distances of free ssDNA and RPA-ssDNA complexes using single-molecule FRET and double electron-electron resonance (DEER) spectroscopy and found only a small systematic increase in the end-to-end distance of ssDNA upon RPA binding. This change does not align with a linear stretching model but rather supports partial wrapping of ssDNA around the contour of DNA binding domains of RPA. Furthermore, we reveal how phosphorylation at the key Ser-384 site in the RPA70 subunit provides access to the wrapped ssDNA by remodeling the DNA-binding domains. These findings establish a precise structural model for RPA-bound ssDNA, providing valuable insights into how RPA facilitates the remodeling of ssDNA for subsequent downstream processes.

Conformations of influenza A M2 protein in DOPC/DOPS and E. coli native lipids and proteins

We compared the conformations of the transmembrane domain (TMD) of influenza A M2 (IM2) protein reconstituted in 1,2-dioleoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phospho-L-serine (DOPC/DOPS) bilayers to those in isolated Escherichia coli (E. coli) membranes, having preserved its native proteins and lipids. IM2 is a single-pass transmembrane protein known to assemble into a homo-tetrameric proton channel. To represent this channel, we made a construct containing the IM2’s TMD region flanked by the juxtamembrane residues. The single cysteine substitution, L43C, of leucine located in the bilayer polar region was paramagnetically tagged with a methanethiosulfonate nitroxide label for the electron spin resonance (ESR) study. For this particular residue, we probed the conformations of the spin-labeled IM2 reconstituted in DOPC/DOPS and isolated E. coli membranes using continuous-wave ESR and double electron-electron resonance (DEER) spectroscopy. The total protein-to-lipid molar ratio spanned the range from 1:230 to 1:10,400. The continuous-wave ESR spectra corresponded to very slow spin-label motion in both environments. In all cases, the DEER data were reconstructed into distance distributions with well-resolved peaks at 1.68 and 2.37 nm in distance and amplitude ratios of 1.41 ± 0.2 and 2:1, respectively. This suggests four nitroxide spin labels located at the corners of a square, indicative of an axially symmetric tetramer. The distance modeling of DEER data with molecular modeling software applied to the NMR molecular structures (PDB: 2L0J) confirmed the symmetry and closed state of the C-terminal exit pore of the IM2 TMD tetramer in agreement with the model. Thus, we can conclude that, under conditions of pH 7.4 used in this study, IM2 TMD has similar conformations in model lipid bilayers and membranes made of native E. coli lipids and proteins of comparable thickness and fluidity, notwithstanding the complexity of the E. coli membranes caused by their lipid diversity and the abundance of integral and peripheral membrane proteins.

Structural Flexibility of Tau in Its Interaction with Microtubules as Viewed by Site-Directed Spin Labeling EPR Spectroscopy.

Tau is a microtubule-associated protein that belongs to the Intrinsically Disordered Proteins (IDPs) family. IDPs or Intrinsically Disordered Regions (IDRs) play key roles in protein interaction networks and their dysfunctions are often related to severe diseases. Defined by their lack of stable secondary and tertiary structures in physiological conditions while being functional, these proteins use their inherent structural flexibility to adapt to and interact with various binding partners. Knowledges on the structural dynamics of IDPs and their different conformers are crucial to finely decipher fundamental biological processes controlled by mechanisms such as conformational adaptations or switches, induced fit, or conformational selection events. Different mechanisms of binding have been proposed: among them, the so-called folding-upon-binding in which the IDP adopts a certain conformation upon interacting with a partner protein, or the formation of a "fuzzy" complex in which the IDP partly keeps its dynamical character at the surface of its partner. The dynamical nature and physicochemical properties of unbound as well as bound IDPs make this class of proteins particularly difficult to characterize by classical bio-structural techniques and require specific approaches for the fine description of their inherent dynamics.Among other techniques, Site-Directed Spin Labeling combined with Electron Paramagnetic Resonance (SDSL-EPR) spectroscopy has gained much interest in this last decade for the study of IDPs. SDSL-EPR consists in grafting a paramagnetic label (mainly a nitroxide radical) at selected site(s) of the macromolecule under interest followed by its observation using and/or combining different EPR strategies. These nitroxide spin labels detected by continuous wave (cw) EPR spectroscopy are used as perfect reporters or "spy spins" of their local environment, being able to reveal structural transitions, folding/unfolding events, etc. Another approach is based on the measurement of inter-label distance distributions in the 1.5-8.0nm range using pulsed dipolar EPR experiments, such as Double Electron-Electron Resonance (DEER) spectroscopy. The technique is then particularly well suited to study the behavior of Tau in its interaction with its physiological partner: microtubules (MTs). In this chapter we provide a detailed experimental protocol for the labeling of Tau protein and its EPR study while interacting with preformed (Paclitaxel-stabilized) MTs, or using Tau as MT inducer. We show how the choice of nitroxide label can be crucial to obtain functional information on Tau/tubulin complexes.

Simultaneous Localization of Two High Affinity Divalent Metal Ion Binding Sites in the Tetracycline RNA Aptamer with Mn2+-Based Pulsed Dipolar EPR Spectroscopy.

Mg2+ ions play an essential part in stabilizing the tertiary structure of nucleic acids. While the importance of these ions is well documented, their localization and elucidation of their role in the structure and dynamics of nucleic acids are often challenging. In this work, pulsed electron-electron double resonance spectroscopy (PELDOR, also known as DEER) was used to localize two high affinity divalent metal ion binding sites in the tetracycline RNA aptamer with high accuracy. For this purpose, the aptamer was labeled at different positions with a semirigid nitroxide spin label and diamagnetic Mg2+ was replaced with paramagnetic Mn2+, which did not alter the folding process or ligand binding. Out of the several divalent metal ion binding sites that are known from the crystal structure, two binding sites with high affinity were detected: one that is located at the ligand binding center and another at the J1/2 junction of the RNA.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter

Here we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize ATP- and substrate-bound inward-facing (IF) and occluded (OC) conformational states of the heterodimeric ATP binding cassette (ABC) multidrug exporter BmrCD in lipid nanodiscs. Supported by DEER analysis, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain to the nucleotide binding domain. The structures uncover asymmetric substrate and Mg2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrate that multiple lipid molecules differentially bind the IF versus the OC conformation thus establishing that lipid interactions modulate BmrCD energy landscape. Our findings are framed in a model that highlights the role of asymmetric conformations in the ATP-coupled transport with general implications to the mechanism of ABC transporters.

Design of stimulus-responsive two-state hinge proteins.

In nature, proteins that switch between two conformations in response to environmental stimuli structurally transduce biochemical information in a manner analogous to how transistors control information flow in computing devices. Designing proteins with two distinct but fully structured conformations is a challenge for protein design as it requires sculpting an energy landscape with two distinct minima. Here we describe the design of "hinge" proteins that populate one designed state in the absence of ligand and a second designed state in the presence of ligand. X-ray crystallography, electron microscopy, double electron-electron resonance spectroscopy, and binding measurements demonstrate that despite the significant structural differences the two states are designed with atomic level accuracy and that the conformational and binding equilibria are closely coupled.

Inner Amino Acid Contacts Are Key Factors of Multistage Structural Rearrangements of DNA and Affect Substrate Specificity of Apurinic/Apyrimidinic Endonuclease APE1.

Apurinic/apyrimidinic endonuclease 1 (APE1) is one of the most important enzymes in base excision repair. Studies on this enzyme have been conducted for a long time, but some aspects of its activity remain poorly understood. One such question concerns the mechanism of damaged-nucleotide recognition by the enzyme, and the answer could shed light on substrate specificity control in all enzymes of this class. In the present study, by pulsed electron-electron double resonance (DEER, also known as PELDOR) spectroscopy and pre-steady-state kinetic analysis along with wild-type (WT) APE1 from Danio rerio (zAPE1) or three mutants (carrying substitution N253G, A254G, or E260A), we aimed to elucidate the molecular events in the process of damage recognition. The data revealed that the zAPE1 mutant E260A has much higher activity toward DNA substrates containing 5,6-dihydro-2'-deoxyuridine (DHU), 2'-deoxyuridine (dU), alpha-2'-deoxyadenosine (αA), or 1,N6-ethenoadenosine (εA). Examination of conformational changes in DNA clearly revealed multistep DNA rearrangements during the formation of the catalytic complex. These structural rearrangements of DNA are directly associated with the capacity of damaged DNA for enzyme-induced bending and unwinding, which are required for eversion of the damaged nucleotide from the DNA duplex and for its placement into the active site of the enzyme. Taken together, the results experimentally prove the factors that control substrate specificity of the AP endonuclease zAPE1.

Ultrafast Bioorthogonal Spin-Labeling and Distance Measurements in Mammalian Cells Using Small, Genetically Encoded Tetrazine Amino Acids.

Site-directed spin-labeling (SDSL)─in combination with double electron-electron resonance (DEER) spectroscopy─has emerged as a powerful technique for determining both the structural states and the conformational equilibria of biomacromolecules. DEER combined with in situ SDSL in live cells is challenging since current bioorthogonal labeling approaches are too slow to allow for complete labeling with low concentrations of spin label prior to loss of signal from cellular reduction. Here, we overcome this limitation by genetically encoding a novel family of small, tetrazine-bearing noncanonical amino acids (Tet-v4.0) at multiple sites in proteins expressed in Escherichia coli and in human HEK293T cells. We achieved specific and quantitative spin-labeling of Tet-v4.0-containing proteins by developing a series of strained trans-cyclooctene (sTCO)-functionalized nitroxides─including a gem-diethyl-substituted nitroxide with enhanced stability in cells─with rate constants that can exceed 106 M-1 s-1. The remarkable speed of the Tet-v4.0/sTCO reaction allowed efficient spin-labeling of proteins in live cells within minutes, requiring only sub-micromolar concentrations of sTCO-nitroxide. DEER recorded from intact cells revealed distance distributions in good agreement with those measured from proteins purified and labeled in vitro. Furthermore, DEER was able to resolve the maltose-dependent conformational change of Tet-v4.0-incorporated and spin-labeled MBP in vitro and support assignment of the conformational state of an MBP mutant within HEK293T cells. We anticipate the exceptional reaction rates of this system, combined with the relatively short and rigid side chains of the resulting spin labels, will enable structure/function studies of proteins directly in cells, without any requirements for protein purification.

Darobactin B Stabilises a Lateral-Closed Conformation of the BAM Complex in E. coli Cells.

The β-barrel assembly machinery (BAM complex) is essential for outer membrane protein (OMP) folding in Gram-negative bacteria, and represents a promising antimicrobial target. Several conformational states of BAM have been reported, but all have been obtained under conditions which lack the unique features and complexity of the outer membrane (OM). Here, we use Pulsed Electron-Electron Double Resonance (PELDOR, or DEER) spectroscopy distance measurements to interrogate the conformational ensemble of the BAM complex in E. coli cells. We show that BAM adopts a broad ensemble of conformations in the OM, while in the presence of the antibiotic darobactin B (DAR-B), BAM's conformational equilibrium shifts to a restricted ensemble consistent with the lateral closed state. Our in-cell PELDOR findings are supported by new cryoEM structures of BAM in the presence and absence of DAR-B. This work demonstrates the utility of PELDOR to map conformational changes in BAM within its native cellular environment.

Darobactin B Stabilises a Lateral-Closed Conformation of the BAM Complex in E. coli Cells.

The β-barrel assembly machinery (BAM complex) is essential for outer membrane protein (OMP) folding in Gram-negative bacteria, and represents a promising antimicrobial target. Several conformational states of BAM have been reported, but all have been obtained under conditions which lack the unique features and complexity of the outer membrane (OM). Here, we use Pulsed Electron-Electron Double Resonance (PELDOR, or DEER) spectroscopy distance measurements to interrogate the conformational ensemble of the BAM complex in E. coli cells. We show that BAM adopts a broad ensemble of conformations in the OM, while in the presence of the antibiotic darobactin B (DAR-B), BAM's conformational equilibrium shifts to a restricted ensemble consistent with the lateral closed state. Our in-cell PELDOR findings are supported by new cryoEM structures of BAM in the presence and absence of DAR-B. This work demonstrates the utility of PELDOR to map conformational changes in BAM within its native cellular environment.

Conducting a three-pulse DEER experiment without dead time: A review

Double electron-electron resonance spectroscopy (DEER, also known as PELDOR) is used to study spin-spin dipolar interactions between spin labels, at the nanoscale range of distances. The DEER effect is obtained as a signal generated by echo-forming microwave (mw) pulses with an additional mw pump pulse applied at a different frequency. It is important to carry out measurements without artefacts induced by overlap of the pulses in the time scale. Such an experiment without the dead-time effect is achieved using the 4-pulse (4p) DEER method. The analysis of the literature performed here shows however that the 3-pulse (3p) DEER can also be free of the dead time problem, for which there are two possibilities. The first occurs using a specially designed bimodal resonator, for which the two frequencies are completely decoupled. The second possibility, which can be implemented for any commercial spectrometer, involves the signal correction based on an additional “blank” measurement with the pump pulse applied outside the EPR resonance. A detailed comparison of the 3p and 4p DEER data obtained previously by Milov et al. [Appl. Magn. Reson. 41 (2011) 59–67] shows that 3p and 4p approaches give similar results. The advantages of the 3p DEER techniques are discussed.

Quantitative analysis of sterol-modulated monomer–dimer equilibrium of the β1-adrenergic receptor by DEER spectroscopy

G protein-coupled receptors (GPCR) activate numerous intracellular signaling pathways. The oligomerization properties of GPCRs, and hence their cellular functions, may be modulated by various components within the cell membrane (such as the presence of cholesterol). Modulation may occur directly via specific interaction with the GPCR or indirectly by affecting the physical properties of the membrane. Here, we use pulsed Q-band double electron-electron resonance (DEER) spectroscopy to probe distances between R1 nitroxide spin labels attached to Cys163 and Cys344 of the β1-adrenergic receptor (β1AR) in n-dodecyl-β-D-maltoside micelles upon titration with two soluble cholesterol analogs, cholesteryl hemisuccinate (CHS) and sodium cholate. The former, like cholesterol, inserts itself into the lipid membrane, parallel to the phospholipid chains; the latter is aligned parallel to the surface of membranes. Global quantitative analysis of DEER echo curves upon titration of spin-labeled β1AR with CHS and sodium cholate reveal the following: CHS binds specifically to the β1AR monomer at a site close to the Cys163-R1 spin label with an equilibrium dissociation constant [Formula: see text]~1.4± 0.4mM. While no direct binding of sodium cholate to the β1AR receptor was observed by DEER, sodium cholate induces specific β1AR dimerization ([Formula: see text]~35± 6mM and a Hill coefficient n~ 2.5±0.4) with intersubunit contacts between transmembrane helices 1 and 2 and helix 8. Analysis of the DEER data obtained upon the addition of CHS to the β1AR dimer in the presence of excess cholate results in dimer dissociation with species occupancies as predicted from the individual KD values.

A structural insight into the PpIX degradation activity of the translocator membrane protein TSPO.

The 18-kDa translocator protein (TSPO) is an evolutionarily well-conserved membrane protein, found primarily in the outer membrane of mitochondria among various species. TSPO binds with high affinity to a variety of drugs, for which it serves as biomarkers and therapeutic targets. Protoporphyrin IX (PpIX) is an endogenous ligand of TSPO and plays a vital role in organisms as a precursor to other critical compounds. The interaction between TSPO and PpIX has been studied for years, but the molecular details of the interaction have yet to be revealed. It has remained challenging to resolve their detailed structure and functional mechanism in a functionally active state in a membrane environment. Here we report our studies on the structure-function relationship of BcTSPO, a bacterial homologue of TSPO, in lipid nanodiscs. Spin-label double electron-electron resonance (DEER) spectroscopy was utilized to investigate the conformation of BcTSPO in nanodiscs with or without PpIX. The results show that the conformation of BcTSPO in lipid nanodiscs is significantly different from the known crystal structures determined in detergent micelles. In addition to the structural conformation of TSPO in nanodiscs by DEER, we have investigated the functional activity of BcTSPO using fluorescence and absorption spectroscopy. We used fluorescence spectroscopy to monitor TSPO-mediated PpIX degradation, and absorption spectroscopy to observe the resulting products. Fluorescence spectroscopy reveals that wild-type BcTSPO, when embedded in nanodiscs or detergent micelles, can catalyze PpIX degradation with a similar reaction rate. When single mutations are applied, we observed large changes in the catalytic behaviors concerning the PpIX degradation, providing new insights into the TSPO-mediated porphyrin catabolism in a light-dependent manner.

Gd3+-Trityl-Nitroxide Triple Labeling and Distance Measurements in the Heterooligomeric Cobalamin Transport Complex in the Native Lipid Bilayers.

Increased efforts are being made for observing proteins in their native environments. Pulsed electron-electron double resonance spectroscopy (PELDOR, also known as DEER) is a powerful tool for this purpose. Conventionally, PELDOR employs an identical spin pair, which limits the output to a single distance for monomeric samples. Here, we show that the Gd3+-trityl-nitroxide (NO) three-spin system is a versatile tool to study heterooligomeric membrane protein complexes, even within their native membrane. This allowed for an independent determination of four different distances (Gd3+-trityl, Gd3+-NO, trityl-NO, and Gd3+-Gd3+) within the same sample. We demonstrate the feasibility of this approach by observing sequential ligand binding and the dynamics of complex formation in the cobalamin transport system involving four components (cobalamin, BtuB, TonB, and BtuF). Our results reveal that TonB binding alone is sufficient to release cobalamin from BtuB in the native asymmetric bilayers. This approach provides a potential tool for the structural and quantitative analysis of dynamic protein-protein interactions in oligomeric complexes, even within their native surroundings.

Reconciling membrane protein simulations with experimental DEER spectroscopy data.

Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepTSo and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.