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

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.

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.

Exploring protein structural ensembles: Integration of sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling methods

Under physiological conditions, proteins continuously undergo structural fluctuations on different timescales. Some conformations are only sparsely populated, but still play a key role in protein function. Thus, meaningful structure-function frameworks must include structural ensembles rather than only the most populated protein conformations. To detail protein plasticity, modern structural biology combines complementary experimental and computational approaches. In this review, we survey available computational approaches that integrate sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling techniques to derive all-atom structural models of rare protein conformations. We also propose strategies to increase the reliability and improve efficiency using deep learning approaches, thus advancing the field of integrative structural biology.

Impact of Cellular Crowding on Protein Structural Dynamics Investigated by EPR Spectroscopy

The study of how the intracellular medium influences protein structural dynamics and protein-protein interactions is a captivating area of research for scientists aiming to comprehend biomolecules in their native environment. As the cellular environment can hardly be reproduced in vitro, direct investigation of biomolecules within cells has attracted growing interest in the past two decades. Among magnetic resonances, site-directed spin labeling coupled to electron paramagnetic resonance spectroscopy (SDSL-EPR) has emerged as a powerful tool for studying the structural properties of biomolecules directly in cells. Since the first in-cell EPR experiment was reported in 2010, substantial progress has been made, and this Review provides a detailed overview of the developments and applications of this spectroscopic technique. The strategies available for preparing a cellular sample and the EPR methods that can be applied to cells will be discussed. The array of spin labels available, along with their strengths and weaknesses in cellular contexts, will also be described. Several examples will illustrate how in-cell EPR can be applied to different biological systems and how the cellular environment affects the structural and dynamic properties of different proteins. Lastly, the Review will focus on the future developments expected to expand the capabilities of this promising technique.

Deconvoluting Monomer- and Dimer-Specific Distance Distributions between Spin Labels in a Monomer/Dimer Mixture Using T1-Edited DEER EPR Spectroscopy

Double electron-electron resonance (DEER) EPR is a powerful tool in structural biology, providing distances between pairs of spin labels. When the sample consists of a mixture of oligomeric species (e.g., monomer and dimer), the question arises as to how to assign the peaks in the DEER-derived probability distance distribution to the individual species. Here, we propose incorporating an EPR longitudinal electron relaxation (T1) inversion recovery experiment within a DEER pulse sequence to resolve this problem. The apparent T1 between dipolar coupled electron spins measured from the inversion recovery time (τinv) dependence of the peak intensities in the T1-edited DEER-derived probability P(r) distance distribution will be affected by the number of nitroxide labels attached to the biomolecule of interest, for example, two for a monomer and four for a dimer. We show that global fitting of all the T1-edited DEER echo curves, recorded over a range of τinv values, permits the deconvolution of distances between spin labels originating from monomeric (longer T1) and dimeric (shorter T1) species. This is especially useful when the trapping of spin labels in different conformational states during freezing gives rise to complex P(r) distance distributions. The utility of this approach is demonstrated for two systems, the β1 adrenergic receptor and a construct of the huntingtin exon-1 protein fused to the immunoglobulin domain of protein G, both of which exist in a monomer-dimer equilibrium.

Monitoring the conformational ensemble and lipid environment of a mechanosensitive channel under cyclodextrin-induced membrane tension

Membrane forces shift the equilibria of mechanosensitive channels enabling them to convert mechanical cues into electrical signals. Molecular tools to stabilize and methods to capture their highly dynamic states are lacking. Cyclodextrins can mimic tension through the sequestering of lipids from membranes. Here we probe the conformational ensemble of MscS by EPR spectroscopy, the lipid environment with NMR, and function with electrophysiology under cyclodextrin-induced tension. We show the extent of MscS activation depends on the cyclodextrin-to-lipid ratio, and that lipids are depleted slower when MscS is present. This has implications in MscS' activation kinetics when distinct membrane scaffolds such as nanodiscs or liposomes are used. We find MscS transits from closed to sub-conducting state(s) before it desensitizes, due to the lack of lipid availability in its vicinity required for closure. Our approach allows for monitoring tension-sensitive states in membrane proteins and screening molecules capable of inducing molecular tension in bilayers.

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, directing the ssDNA towards downstream DNA metabolic pathways with exceptional specificity. Understanding the mechanistic basis for such RPA-dependent specificity requires a comprehensive understanding of the structural conformation of ssDNA when bound to RPA. 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.

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.

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.

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.4 mM. 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 ± 6 mM 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.

Proton-driven alternating access in a spinster lipid transporter

Spinster (Spns) lipid transporters are critical for transporting sphingosine-1-phosphate (S1P) across cellular membranes. In humans, Spns2 functions as the main S1P transporter in endothelial cells, making it a potential drug target for modulating S1P signaling. Here, we employed an integrated approach in lipid membranes to identify unknown conformational states of a bacterial Spns from Hyphomonas neptunium (HnSpns) and to define its proton- and substrate-coupled conformational dynamics. Our systematic study reveals conserved residues critical for protonation steps and their regulation, and how sequential protonation of these proton switches coordinates the conformational transitions in the context of a noncanonical ligand-dependent alternating access. A conserved periplasmic salt bridge (Asp60TM2:Arg289TM7) keeps the transporter in a closed conformation, while proton-dependent conformational dynamics are significantly enhanced on the periplasmic side, providing a pathway for ligand exchange.

Integrated AlphaFold2 and DEER investigation of the conformational dynamics of a pH-dependent APC antiporter

The Amino Acid-Polyamine-Organocation (APC) transporter GadC contributes to the survival of pathogenic bacteria under extreme acid stress by exchanging extracellular glutamate for intracellular γ-aminobutyric acid (GABA). Its structure, determined in an inward-facing conformation at alkaline pH, consists of the canonical LeuT-fold with a conserved five-helix inverted repeat, thereby resembling functionally divergent transporters such as the serotonin transporter SERT and the glucose-sodium symporter SGLT1. However, despite this structural similarity, it is unclear if the conformational dynamics of antiporters such as GadC follow the blueprint of these or other LeuT-fold transporters. Here, we used double electron-electron resonance (DEER) spectroscopy to monitor the conformational dynamics of GadC in lipid bilayers in response to acidification and substrate binding. To guide experimental design and facilitate the interpretation of the DEER data, we generated an ensemble of structural models in multiple conformations using a recently introduced modification of AlphaFold2 . Our experimental results reveal acid-induced conformational changes that dislodge the Cterminus from the permeation pathway coupled with rearrangement of helices that enables isomerization between inward- and outward-facing states. The substrate glutamate, but not GABA, modulates the dynamics of an extracellular thin gate without shifting the equilibrium between inward- and outward-facing conformations. In addition to introducing an integrated methodology for probing transporter conformational dynamics, the congruence of the DEER data with patterns of structural rearrangements deduced from ensembles of AlphaFold2 models illuminates the conformational cycle of GadC underpinning transport and exposes yet another example of the divergence between the dynamics of different families in the LeuT-fold.

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/4Ca²⁺) 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/4Ca²⁺ 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/4Ca²⁺ 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/4Ca²⁺ 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/4Ca²⁺ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

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.

Time-resolved DEER EPR and solid-state NMR afford kinetic and structural elucidation of substrate binding to Ca2+-ligated calmodulin

Complex formation between calmodulin and target proteins underlies numerous calcium signaling processes in biology, yet structural and mechanistic details, which entail major conformational changes in both calmodulin and its substrates, have been unclear. We show that a combination of time-resolved electron paramagnetic and NMR measurements can elucidate the molecular mechanism, at the quantitative kinetic and structural levels, of the binding pathway of a peptide substrate from skeletal muscle myosin light-chain kinase to calcium-loaded calmodulin. The mechanism involves coupled folding and binding and comprises a bifurcated process, with rapid, direct complex formation when the peptide interacts first with the C-terminal domain of calmodulin or a slower, two-step complex formation when the peptide interacts initially with the N-terminal domain.

Recent advances in rapid mixing and freeze quenching have opened the path for time-resolved electron paramagnetic resonance (EPR)-based double electron-electron resonance (DEER) and solid-state NMR of protein–substrate interactions. DEER, in conjunction with phase memory time filtering to quantitatively extract species populations, permits monitoring time-dependent probability distance distributions between pairs of spin labels, while solid-state NMR provides quantitative residue-specific information on the appearance of structural order and the development of intermolecular contacts between substrate and protein. Here, we demonstrate the power of these combined approaches to unravel the kinetic and structural pathways in the binding of the intrinsically disordered peptide substrate (M13) derived from myosin light-chain kinase to the universal eukaryotic calcium regulator, calmodulin. Global kinetic analysis of the data reveals coupled folding and binding of the peptide associated with large spatial rearrangements of the two domains of calmodulin. The initial binding events involve a bifurcating pathway in which the M13 peptide associates via either its N- or C-terminal regions with the C- or N-terminal domains, respectively, of calmodulin/4Ca²⁺ to yield two extended “encounter” complexes, states A and A*, without conformational ordering of M13. State A is immediately converted to the final compact complex, state C, on a timescale τ ≤ 600 μs. State A*, however, only reaches the final complex via a collapsed intermediate B (τ ∼ 1.5 to 2.5 ms), in which the peptide is only partially ordered and not all intermolecular contacts are formed. State B then undergoes a relatively slow (τ ∼ 7 to 18 ms) conformational rearrangement to state C.

Asymmetric drug binding in an ATP-loaded inward-facing state of an ABC transporter

Substrate efflux by ATP-binding cassette (ABC) transporters, which play a major role in multidrug resistance, entails the ATP-powered interconversion between transporter intermediates. Despite recent progress in structure elucidation, a number of intermediates have yet to be visualized and mechanistically interpreted. Here, we combine cryogenic-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy and molecular dynamics simulations to profile a previously unobserved intermediate of BmrCD, a heterodimeric multidrug ABC exporter from Bacillus subtilis. In our cryo-EM structure, ATP-bound BmrCD adopts an inward-facing architecture featuring two molecules of the substrate Hoechst-33342 in a striking asymmetric head-to-tail arrangement. Deletion of the extracellular domain capping the substrate-binding chamber or mutation of Hoechst-coordinating residues abrogates cooperative stimulation of ATP hydrolysis. Together, our findings support a mechanistic role for symmetry mismatch between the nucleotide binding and the transmembrane domains in the conformational cycle of ABC transporters and is of notable importance for rational design of molecules for targeted ABC transporter inhibition.