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

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 β₁-adrenergic receptor (β₁AR) 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 β₁AR with CHS and sodium cholate reveal the following: CHS binds specifically to the β₁AR 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 β₁AR receptor was observed by DEER, sodium cholate induces specific β₁AR 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 β₁AR dimer in the presence of excess cholate results in dimer dissociation with species occupancies as predicted from the individual K_D 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 (Asp60^TM2:Arg289^TM7) 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.

The Spns lipid transporters use ion gradients to drive substrate transport, including bioactive sphingolipids. Here, Dastvan et al. investigated how binding of protons powers the conformational changes that enable broad transport by a bacterial Spns.

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.