Toward the fourth dimension of membrane protein structure: insight into dynamics from spin-labeling EPR spectroscopy

Protein Modeling with DEER Spectroscopy

Double electron-electron resonance (DEER) combined with site-directed spin labeling can provide distance distributions between selected protein residues to investigate protein structure and conformational heterogeneity. The utilization of the full quantitative information contained in DEER data requires effective protein and spin label modeling methods. Here, we review the application of DEER data to protein modeling. First, we discuss the significance of spin label modeling for accurate extraction of protein structural information and review the most popular label modeling methods. Next, we review several important aspects of protein modeling with DEER, including site selection, how DEER restraints are applied, common artifacts, and the unique potential of DEER data for modeling structural ensembles and conformational landscapes. Finally, we discuss common applications of protein modeling with DEER data and provide an outlook.

Enhanced Binding Site Identification in Protein-Ligand Complexes with a Combined Blind Docking and Dipolar Electron Paramagnetic Resonance Approach

Understanding protein-drug complex structures is crucial for elucidating therapeutic mechanisms and side effects. Blind docking facilitates site identification but is hindered by computational complexity and imprecise scoring, causing ambiguity. Dipolar electron paramagnetic resonance (EPR) provides spin-spin distances but struggles to determine relative positions within complexes. We present a novel approach combining GPU-accelerated blind docking with EPR distance constraints to enhance binding site detection. Our algorithm uses a single EPR distance distribution to filter and validate docking results. Ligand poses from blind docking are clustered, filtered by expected distances, and refined through focused docking. To illustrate our approach, we investigated human serum albumin binding with porphyrin-based photosensitizers used in photodynamic therapy. Combining docking and EPR, we identified possible binding sites, demonstrating that EPR data significantly reduce possible configurations and provide experimentally validated information. This strategy produces a detailed map of photoligand binding sites, revealing that binding may occur away from standard albumin sites and often involves multiple locations. Furthermore, it overcomes key limitations of fluorescence-based methods, which are prone to misinterpretation in albumin studies due to non one-to-one donor-acceptor relationships. By resolving ambiguities in both blind docking and EPR, our framework provides a versatile platform for investigating EPR-active ligands.

Characterization of conformationally heterogeneous proteins by electron paramagnetic resonance spectroscopy

The Anfinsen paradigm of representing a protein by a single conformer is challenged by the uncertainty predictions that come with AlphaFold models, which suggest a greater extent of disorder. Characterization of such conformation heterogeneity requires experimental approaches that do not depend on long-range order. Site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is such an approach. The double electron-electron resonance (DEER) technique can access site-pair distance distributions in the 15-100 Å range, directly informing on ensemble width. SDSL-EPR can be applied in cellular environments, and recent work indicates that protein disorder is even more pervasive than predicted by AlphaFold. This suggests that the Anfinsen paradigm should be replaced by an ensemble paradigm.

Localized Reconstruction of Multimodal Distance Distribution from DEER Data of Biopolymers

Pulsed Dipolar ESR Spectroscopy (PDS) is a uniquely powerful technique to characterize the structural property of intrinsically disordered proteins (IDPs) and polymers and the conformational evolution of IDPs and polymers, e.g. during assembly, by offering the probability distribution of segment end-to-end distances. However, it is challenging to determine distance distribution P(r) of IDPs by PDS because of the uncertain and broad shape information that is intrinsic to the distance distribution of IDPs. We demonstrate here that the Srivastava-Freed Singular Value Decomposition (SF-SVD) point-wise mathematical inversion method along with wavelet denoising (WavPDS) can aid in obtaining reliable shapes for the distance distribution, P(r), for IDPs. We show that broad regions of P(r) as well as mixed narrow and broad features within the captured distance distribution range can be effectively resolved and differentiated without a priori knowledge. The advantage of SF-SVD and WavPDS is that the methods are transparent, requiring no adjustable parameters, the processing of the magnitude for the probability distribution is performed separately for each distance increment, and the outcome of the analysis is independent of the user's judgement. We demonstrate the performance and present the application of WavPDS and SF-SVD on model ruler molecules, model polyethylene glycol polymers with end-to-end spin labeling, and IDPs with pairwise labeling spanning different segments of the protein tau to generate the transparent solutions to the P(r)'s including their uncertainties and error analysis.

Mapping Protein Conformational Landscapes from Crystallographic Drug Fragment Screens

Proteins are dynamic macromolecules. Knowledge of a protein's thermally accessible conformations is critical to determining important transitions and designing therapeutics. Accessible conformations are highly constrained by a protein's structure such that concerted structural changes due to external perturbations likely track intrinsic conformational transitions. These transitions can be thought of as paths through a conformational landscape. Crystallographic drug fragment screens are high-throughput perturbation experiments, in which thousands of crystals of a drug target are soaked with small-molecule drug precursors (fragments) and examined for fragment binding, mapping potential drug binding sites on the target protein. Here, we describe an open-source Python package, COnformational LAndscape Visualization (COLAV), to infer conformational landscapes from such large-scale crystallographic perturbation studies. We apply COLAV to drug fragment screens of two medically important systems: protein tyrosine phosphatase 1B (PTP1B), which regulates insulin signaling, and the SARS CoV-2 Main Protease (MPro). With enough fragment-bound structures, we find that such drug screens enable detailed mapping of proteins' conformational landscapes.

Conformational cycle of a protease-containing ABC transporter in lipid nanodiscs reveals the mechanism of cargo-protein coupling

Protease-containing ABC transporters (PCATs) couple the energy of ATP hydrolysis to the processing and export of diverse cargo proteins across cell membranes to mediate antimicrobial resistance and quorum sensing. Here, we combine biochemical analysis, single particle cryoEM, and DEER spectroscopy in lipid bilayers along with computational analysis to illuminate the structural and energetic underpinnings of coupled cargo protein export. Our integrated investigation uncovers competitive interplay between nucleotides and cargo protein binding that ensures the latter's orderly processing and subsequent transport. The energetics of cryoEM structures in lipid bilayers are congruent with the inferred mechanism from ATP turnover analysis and reveal a snapshot of a high-energy outward-facing conformation that provides an exit pathway into the lipid bilayer and/or the extracellular side. DEER investigation of the core ABC transporter suggests evolutionary tuning of the energetic landscape to fulfill the function of substrate processing prior to export.

Mapping protein conformational landscapes from crystallographic drug fragment screens

Proteins are dynamic macromolecules. Knowledge of a protein's thermally accessible conformations is critical to determining important transitions and designing therapeutics. Accessible conformations are highly constrained by a protein's structure such that concerted structural changes due to external perturbations likely track intrinsic conformational transitions. These transitions can be thought of as paths through a conformational landscape. Crystallographic drug fragment screens are high-throughput perturbation experiments, in which thousands of crystals of a drug target are soaked with small-molecule drug precursors (fragments) and examined for fragment binding, mapping potential drug binding sites on the target protein. Here, we describe an open-source Python package, COLAV (COnformational LAndscape Visualization), to infer conformational landscapes from such large-scale crystallographic perturbation studies. We apply COLAV to drug fragment screens of two medically important systems: protein tyrosine phosphatase 1B (PTP-1B), which regulates insulin signaling, and the SARS CoV-2 Main Protease (MPro). With enough fragment-bound structures, we find that such drug screens also enable detailed mapping of proteins' conformational landscapes.

Designing surface exposed sites on Bacillus subtilis lipase A for spin-labeling and hydration studies

Spin-labeling with electron paramagnetic resonance spectroscopy (EPR) is a facile method for interrogating macromolecular flexibility, conformational changes, accessibility, and hydration. Within we present a computationally based approach for the rational selection of reporter sites in Bacillus subtilis lipase A (BSLA) for substitution to cysteine residues with subsequent modification with a spin-label that are expected to not significantly perturb the wild-type structure, dynamics, or enzymatic function. Experimental circular dichroism spectroscopy, Michaelis-Menten kinetic parameters and EPR spectroscopy data validate the success of this approach to computationally select reporter sites for future magnetic resonance investigations of hydration and hydration changes induced by polymer conjugation, tethering, immobilization, or amino acid substitution in BSLA. Analysis of molecular dynamic simulations of the impact of substitutions on the secondary structure agree well with experimental findings. We propose that this computationally guided approach for choosing spin-labeled EPR reporter sites, which evaluates relative surface accessibility coupled with hydrogen bonding occupancy of amino acids to the catalytic pocket via atomistic simulations, should be readily transferable to other macromolecular systems of interest including selecting sites for paramagnetic relaxation enhancement NMR studies, other spin-labeling EPR studies or any method requiring a tagging method where it is desirable to not alter enzyme stability or activity.

Integrating Electron Paramagnetic Resonance Spectroscopy and Computational Modeling to Measure Protein Structure and Dynamics

Electron paramagnetic resonance (EPR) has become a powerful probe of conformational heterogeneity and dynamics of biomolecules. In this Review, we discuss different computational modeling techniques that enrich the interpretation of EPR measurements of dynamics or distance restraints. A variety of spin labels are surveyed to provide a background for the discussion of modeling tools. Molecular dynamics (MD) simulations of models containing spin labels provide dynamical properties of biomolecules and their labels. These simulations can be used to predict EPR spectra, sample stable conformations and sample rotameric preferences of label sidechains. For molecular motions longer than milliseconds, enhanced sampling strategies and de novo prediction software incorporating or validated by EPR measurements are able to efficiently refine or predict protein conformations, respectively. To sample large-amplitude conformational transition, a coarse-grained or an atomistic weighted ensemble (WE) strategy can be guided with EPR insights. Looking forward, we anticipate an integrative strategy for efficient sampling of alternate conformations by de novo predictions, followed by validations by systematic EPR measurements and MD simulations. Continuous pathways between alternate states can be further sampled by WE-MD including all intermediate states.

HIV-1 Vpu protein forms stable oligomers in aqueous solution via its transmembrane domain self-association

We report our findings on the assembly of the HIV-1 protein Vpu into soluble oligomers. Vpu is a key HIV-1 protein. It has been considered exclusively a single-pass membrane protein. Previous observations show that this protein forms stable oligomers in aqueous solution, but details about these oligomers still remain obscure. This is an interesting and rather unique observation, as the number of proteins transitioning between soluble and membrane embedded states is limited. In this study we made use of protein engineering, size exclusion chromatography, cryoEM and electron paramagnetic resonance (EPR) spectroscopy to better elucidate the nature of the soluble oligomers. We found that Vpu oligomerizes via its N-terminal transmembrane domain (TM). CryoEM suggests that the oligomeric state most likely is a hexamer/heptamer equilibrium. Both cryoEM and EPR suggest that, within the oligomer, the distal C-terminal region of Vpu is highly flexible. Our observations are consistent with both the concept of specific interactions among TM helices or the core of the oligomers being stabilized by hydrophobic forces. While this study does not resolve all of the questions about Vpu oligomers or their functional role in HIV-1 it provides new fundamental information about the size and nature of the oligomeric interactions.

Eukaryotic Cell Membranes: Structure, Composition, Research Methods and Computational Modelling

This paper deals with the problems encountered in the study of eukaryotic cell membranes. A discussion on the structure and composition of membranes, lateral heterogeneity of membranes, lipid raft formation, and involvement of actin and cytoskeleton networks in the maintenance of membrane structure is included. Modern methods for the study of membranes and their constituent domains are discussed. Various simplified models of biomembranes and lipid rafts are presented. Computer modelling is considered as one of the most important methods. This is stated that from the study of the plasma membrane structure, it is desirable to proceed to the diverse membranes of all organelles of the cell. The qualitative composition and molar content of individual classes of polar lipids, free sterols and proteins in each of these membranes must be considered. A program to create an open access electronic database including results obtained from the membrane modelling of individual cell organelles and the key sites of the membranes, as well as models of individual molecules composing the membranes, has been proposed.

HIV-1 Vpu protein forms stable oligomers in aqueous solution via its transmembrane domain self-association

We report our findings on the assembly of the HIV-1 protein Vpu into soluble oligomers. Vpu is a key to HIV-1 protein. It has been considered exclusively a single-pass membrane protein. However, we revealed that this protein forms stable oligomers in aqueous solution, which is an interesting and rather unique observation, as the number of proteins transitioning between soluble and membrane embedded states is limited. Therefore, we undertook a study to characterize these oligomers by utilizing protein engineering, size exclusion chromatography, cryoEM and electron paramagnetic resonance (EPR) spectroscopy. We found that Vpu oligomerizes via its N-terminal transmembrane domain (TM). CryoEM analyses suggest that the oligomeric state most likely is a hexamer or hexamer-to-heptamer equilibrium. Both cryoEM and EPR suggest that, within the oligomer, the distant C-terminal region of Vpu is highly flexible. To the best of our knowledge, this is the first comprehensive study on soluble Vpu. We propose that these oligomers are stabilized via possibly hydrophobic interactions between Vpu TMs. Our findings contribute valuable information about this protein properties and about protein supramolecular complexes formation. The acquired knowledge could be further used in protein engineering, and could also help to uncover possible physiological function of these Vpu oligomers.

Perspective on the Effect of Membrane Mimetics on Dynamic Properties of Integral Membrane Proteins

Integral membrane proteins are embedded into cell membranes by spanning the width of the lipid bilayer. They play an essential role in important biological functions for the survival of living organisms. Their functions include the transportation of ions and molecules across the cell membrane and initiating signaling pathways. The dynamic behavior of integral membrane proteins is very important for their function. Due to the complex behavior of integral membrane proteins in the cell membrane, studying their structural dynamics using biophysical approaches is challenging. Here, we concisely discuss challenges and recent advances in technical and methodological aspects of biophysical approaches for gleaning dynamic properties of integral membrane proteins to answer pertinent biological questions associated with these proteins.

"Store-bought is fine": Sensitivity considerations using shaped pulses for DEER measurements on Cu(II) labels

The narrow excitation bandwidth of monochromic pulses is a sensitivity limitation for pulsed dipolar spectroscopy on Cu(II)-based measurements. In response, frequency-swept pulses with large excitation bandwidths have been adopted to probe a greater range of the EPR spectrum. However, much of the work utilizing frequency-swept pulses in Cu(II) distance measurements has been carried out on home-built spectrometers and equipment. Herein, we carry out systematic Cu(II) based distance measurements to demonstrate the capability of chirp pulses on commercial instrumentation. More importantly we delineate sensitivity considerations under acquisition schemes that are necessary for robust distance measurements using Cu(II) labels for proteins. We show that a 200 MHz sweeping bandwidth chirp pulse can improve the sensitivity of long-range distance measurements by factors of three to four. The sensitivity of short-range distances only increases slightly due to special considerations for the chirp pulse duration relative to the period length of the modulated dipolar signal. Enhancements in sensitivity also dramatically reduce measurement collection times enabling rapid collection of orientationally averaged Cu(II) distance measurements in under two hours.

Enzymatic Spin-Labeling of Protein N- and C-Termini for Electron Paramagnetic Resonance Spectroscopy

Electron paramagnetic resonance (EPR) spectroscopy is a powerful tool for investigating the structure and dynamics of proteins. The introduction of paramagnetic moieties at specific positions in a protein enables precise measurement of local structure and dynamics. This technique, termed site-directed spin-labeling, has traditionally been performed using cysteine-reactive radical-containing probes. However, large proteins are more likely to contain multiple cysteine residues and cysteine labeling at specific sites may be infeasible or impede function. To address this concern, we applied three peptide-ligating enzymes (sortase, asparaginyl endopeptidase, and inteins) for nitroxide labeling of N- and C-termini of select monomeric and dimeric proteins. Continuous wave and pulsed EPR (double electron electron resonance) experiments reveal specific attachment of nitroxide probes to ether N-termini (OaAEP1) or C-termini (sortase and intein) across three test proteins (CheY, CheA, and iLOV), thereby enabling a straightforward, highly specific, and general method for protein labeling. Importantly, the linker length (3, 5, and 9 residues for OaAEP1, intein, and sortase reactions, respectively) between the probe and the target protein has a large impact on the utility of distance measurements by pulsed EPR, with longer linkers leading to broader distributions. As these methods are only dependent on accessible N- and C-termini, we anticipate application to a wide range of protein targets for biomolecular EPR spectroscopy.

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, PepT_So 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.

Unveiling the orientation and dynamics of enzymes in unstructured artificial compartments of metal-organic frameworks (MOFs)

Confining enzymes in well-shaped MOF compartments is a promising approach to mimic the cellular environment of enzymes and determine enzyme structure-function relationship therein. Under the cellular crowding, however, enzymes can also be confined in unstructured spaces that are close to the shapes/outlines of the enzyme. Therefore, for a better understanding of enzymes in their physiological environments, it is necessary to study enzymes in these unstructured spaces. However, practically it is challenging to create compartments that are close to the outline of an enzyme and probe enzyme structural information therein. Here, for proof-of-principle, we confined a model enzyme, lysozyme, in the crystal defects of a MOF via co-crystallization, where lysozyme served as the nuclei for MOF crystal scaffolds to grow on so that unstructured spaces close to the outline of lysozyme are created, and determined enzyme relative orientation and dynamics. This effort is important for understanding enzymes in near-native environments and guiding the rational design of biocatalysts that mimic how nature confines enzymes.

Designer Nanodiscs to Probe and Reprogram Membrane Biology in Synapses

Signal transduction at the synapse is mediated by a variety of protein-lipid interactions, which are vital for the spatial and temporal regulation of synaptic vesicle biogenesis, neurotransmitter release, and postsynaptic receptor activation. Therefore, our understanding of synaptic transmission cannot be completed until the elucidation of these critical protein-lipid interactions. On this front, recent advances in nanodiscs have vastly expanded our ability to probe and reprogram membrane biology in synapses. Here, we summarize the progress of the nanodisc toolbox and discuss future directions in this exciting field.

Comparative evaluation of spin-label modeling methods for protein structural studies

Site-directed spin-labeling electron paramagnetic resonance spectroscopy is a powerful technique for the investigation of protein structure and dynamics. Accurate spin-label modeling methods are essential to make full quantitative use of site-directed spin-labeling electron paramagnetic resonance data for protein modeling and model validation. Using a set of double electron-electron resonance data from seven different site pairs on maltodextrin/maltose-binding protein under two different conditions using five different spin labels, we compare the ability of two widely used spin-label modeling methods, based on accessible volume sampling and rotamer libraries, to predict experimental distance distributions. We present a spin-label modeling approach inspired by canonical side-chain modeling methods and compare modeling accuracy with the established methods.

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.

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.

Integrative structural dynamics probing of the conformational heterogeneity in synaptosomal-associated protein 25

SNAP-25 (synaptosomal-associated protein of 25 kDa) is a prototypical intrinsically disordered protein (IDP) that is unstructured by itself but forms coiled-coil helices in the SNARE complex. With high conformational heterogeneity, detailed structural dynamics of unbound SNAP-25 remain elusive. Here, we report an integrative method to probe the structural dynamics of SNAP-25 by combining replica-exchange discrete molecular dynamics (rxDMD) simulations and label-based experiments at ensemble and single-molecule levels. The rxDMD simulations systematically characterize the coil-to-molten globular transition and reconstruct structural ensemble consistent with prior ensemble experiments. Label-based experiments using Förster resonance energy transfer and double electron-electron resonance further probe the conformational dynamics of SNAP-25. Agreements between simulations and experiments under both ensemble and single-molecule conditions allow us to assign specific helix-coil transitions in SNAP-25 that occur in submillisecond timescales and potentially play a vital role in forming the SNARE complex. We expect that this integrative approach may help further our understanding of IDPs.

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

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

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

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

Lipid Membrane Mimetics in Functional and Structural Studies of Integral Membrane Proteins

Integral membrane proteins (IMPs) fulfill important physiological functions by providing cell-environment, cell-cell and virus-host communication; nutrients intake; export of toxic compounds out of cells; and more. However, some IMPs have obliterated functions due to polypeptide mutations, modifications in membrane properties and/or other environmental factors-resulting in damaged binding to ligands and the adoption of non-physiological conformations that prevent the protein from returning to its physiological state. Thus, elucidating IMPs' mechanisms of function and malfunction at the molecular level is important for enhancing our understanding of cell and organism physiology. This understanding also helps pharmaceutical developments for restoring or inhibiting protein activity. To this end, in vitro studies provide invaluable information about IMPs' structure and the relation between structural dynamics and function. Typically, these studies are conducted on transferred from native membranes to membrane-mimicking nano-platforms (membrane mimetics) purified IMPs. Here, we review the most widely used membrane mimetics in structural and functional studies of IMPs. These membrane mimetics are detergents, liposomes, bicelles, nanodiscs/Lipodisqs, amphipols, and lipidic cubic phases. We also discuss the protocols for IMPs reconstitution in membrane mimetics as well as the applicability of these membrane mimetic-IMP complexes in studies via a variety of biochemical, biophysical, and structural biology techniques.

Single-Molecule FRET of Membrane Transport Proteins

Uncovering the structure and function of biomolecules is a fundamental goal in structural biology. Membrane-embedded transport proteins are ubiquitous in all kingdoms of life. Despite structural flexibility, their mechanisms are typically studied by ensemble biochemical methods or by static high-resolution structures, which complicate a detailed understanding of their dynamics. Here, we review the recent progress of single molecule Förster Resonance Energy Transfer (smFRET) in determining mechanisms and timescales of substrate transport across membranes. These studies do not only demonstrate the versatility and suitability of state-of-the-art smFRET tools for studying membrane transport proteins but they also highlight the importance of membrane mimicking environments in preserving the function of these proteins. The current achievements advance our understanding of transport mechanisms and have the potential to facilitate future progress in drug design.

Assessing the Role of Lipids in the Molecular Mechanism of Membrane Proteins

Membrane proteins have evolved to work optimally within the complex environment of the biological membrane. Consequently, interactions with surrounding lipids are part of their molecular mechanism. Yet, the identification of lipid-protein interactions and the assessment of their molecular role is an experimental challenge. Recently, biophysical approaches have emerged that are compatible with the study of membrane proteins in an environment closer to the biological membrane. These novel approaches revealed specific mechanisms of regulation of membrane protein function. Lipids have been shown to play a role in oligomerization, conformational transitions or allosteric coupling. In this review, we summarize the recent biophysical approaches, or combination thereof, that allow to decipher the role of lipid-protein interactions in the mechanism of membrane proteins.

Methodology for rigorous modeling of protein conformational changes by Rosetta using DEER distance restraints

We describe an approach for integrating distance restraints from Double Electron-Electron Resonance (DEER) spectroscopy into Rosetta with the purpose of modeling alternative protein conformations from an initial experimental structure. Fundamental to this approach is a multilateration algorithm that harnesses sets of interconnected spin label pairs to identify optimal rotamer ensembles at each residue that fit the DEER decay in the time domain. Benchmarked relative to data analysis packages, the algorithm yields comparable distance distributions with the advantage that fitting the DEER decay and rotamer ensemble optimization are coupled. We demonstrate this approach by modeling the protonation-dependent transition of the multidrug transporter PfMATE to an inward facing conformation with a deviation to the experimental structure of less than 2Å Cα RMSD. By decreasing spin label rotamer entropy, this approach engenders more accurate Rosetta models that are also more closely clustered, thus setting the stage for more robust modeling of protein conformational changes.

dHis-troying Barriers: Deuteration Provides a Pathway to Increase Sensitivity and Accessible Distances for Cu2+ Labels

Recently, site-directed Cu²⁺ labeling has emerged as an incisive biophysical tool to directly report on distance constraints that pertain to the structure, conformational transitions, and dynamics of proteins and nucleic acids. However, short phase memory times inherent to the Cu²⁺ labels limit measurable distances to 4-5 nm. In this work we systematically examine different methods to dampen electron-nuclear and electron-electron coupled interactions to decrease rapid relaxation. We show that using Cu²⁺ spin concentrations up to ca. 800 μM has an invariant effect on relaxation and that increasing the cryoprotectant concentration reduces contributions of solvent protons to relaxation. On the other hand, the deuteration of protein and solvent dramatically increases the duration of the dipolar modulated signal by over 6-fold to 32 μs. Based on this increase in signal longevity, distances up to 9 nm and beyond can potentially be measured with Cu²⁺ labels.

AlphaFold2 predicts the inward-facing conformation of the multidrug transporter LmrP

As part of the CASP competition, the protein structure prediction algorithm AlphaFold2 generated multiple models of the proton/drug antiporter LmrP. Previous distance restraints from double electron-electron resonance spectroscopy, a technique which reports distance distributions between spin labels attached to proteins, suggest that one of the lower-ranked models may have captured a conformation that has so far eluded experimental structure determination.

Psychomotor impairments and therapeutic implications revealed by a mutation associated with infantile Parkinsonism-Dystonia

Parkinson disease (PD) is a progressive, neurodegenerative disorder affecting over 6.1 million people worldwide. Although the cause of PD remains unclear, studies of highly penetrant mutations identified in early-onset familial parkinsonism have contributed to our understanding of the molecular mechanisms underlying disease pathology. Dopamine (DA) transporter (DAT) deficiency syndrome (DTDS) is a distinct type of infantile parkinsonism-dystonia that shares key clinical features with PD, including motor deficits (progressive bradykinesia, tremor, hypomimia) and altered DA neurotransmission. Here, we define structural, functional, and behavioral consequences of a Cys substitution at R445 in human DAT (hDAT R445C), identified in a patient with DTDS. We found that this R445 substitution disrupts a phylogenetically conserved intracellular (IC) network of interactions that compromise the hDAT IC gate. This is demonstrated by both Rosetta molecular modeling and fine-grained simulations using hDAT R445C, as well as EPR analysis and X-ray crystallography of the bacterial homolog leucine transporter. Notably, the disruption of this IC network of interactions supported a channel-like intermediate of hDAT and compromised hDAT function. We demonstrate that Drosophila melanogaster expressing hDAT R445C show impaired hDAT activity, which is associated with DA dysfunction in isolated brains and with abnormal behaviors monitored at high-speed time resolution. We show that hDAT R445C Drosophila exhibit motor deficits, lack of motor coordination (i.e. flight coordination) and phenotypic heterogeneity in these behaviors that is typically associated with DTDS and PD. These behaviors are linked with altered dopaminergic signaling stemming from loss of DA neurons and decreased DA availability. We rescued flight coordination with chloroquine, a lysosomal inhibitor that enhanced DAT expression in a heterologous expression system. Together, these studies shed some light on how a DTDS-linked DAT mutation underlies DA dysfunction and, possibly, clinical phenotypes shared by DTDS and PD.

Probing Structural Dynamics of Membrane Proteins Using Electron Paramagnetic Resonance Spectroscopic Techniques

Membrane proteins are essential for the survival of living organisms. They are involved in important biological functions including transportation of ions and molecules across the cell membrane and triggering the signaling pathways. They are targets of more than half of the modern medical drugs. Despite their biological significance, information about the structural dynamics of membrane proteins is lagging when compared to that of globular proteins. The major challenges with these systems are low expression yields and lack of appropriate solubilizing medium required for biophysical techniques. Electron paramagnetic resonance (EPR) spectroscopy coupled with site directed spin labeling (SDSL) is a rapidly growing powerful biophysical technique that can be used to obtain pertinent structural and dynamic information on membrane proteins. In this brief review, we will focus on the overview of the widely used EPR approaches and their emerging applications to answer structural and conformational dynamics related questions on important membrane protein systems.

PELDOR/DEER: An Electron Paramagnetic Resonance Method to Study Membrane Proteins in Lipid Bilayers

Every membrane protein is involved in close interactions with the lipid environment of cellular membranes. The annular lipids, that are in direct contact with the polypeptide, can in principle be seen as an integral part of its structure, akin to the first hydration shell of soluble proteins. It is therefore desirable to investigate the structure of membrane proteins and especially their conformational flexibility under conditions that are as close as possible to their native state. This can be achieved by reconstituting the protein into proteoliposomes, nanodiscs, or bicelles. In recent years, PELDOR/DEER spectroscopy has proved to be a very useful method to study the structure and function of membrane proteins in such artificial membrane environments. The technique complements both X-ray crystallography and cryo-EM and can be used in combination with virtually any artificial membrane environment and under certain circumstances even in native membranes. Of the above-mentioned membrane mimics, bicelles are currently the least often used for PELDOR studies, although they offer some advantages, especially their ease of use. Here, we provide a step-by-step protocol for studying a bicelle reconstituted membrane protein with PELDOR/DEER spectroscopy.

Prokaryotic Solute/Sodium Symporters: Versatile Functions and Mechanisms of a Transporter Family

The solute/sodium symporter family (SSS family; TC 2.A.21; SLC5) consists of integral membrane proteins that use an existing sodium gradient to drive the uphill transport of various solutes, such as sugars, amino acids, vitamins, or ions across the membrane. This large family has representatives in all three kingdoms of life. The human sodium/iodide symporter (NIS) and the sodium/glucose transporter (SGLT1) are involved in diseases such as iodide transport defect or glucose-galactose malabsorption. Moreover, the bacterial sodium/proline symporter PutP and the sodium/sialic acid symporter SiaT play important roles in bacteria–host interactions. This review focuses on the physiological significance and structural and functional features of prokaryotic members of the SSS family. Special emphasis will be given to the roles and properties of proteins containing an SSS family domain fused to domains typically found in bacterial sensor kinases.

Conserved binding site in the N-lobe of prokaryotic MATE transporters suggests a role for Na+ in ion-coupled drug efflux

In both prokaryotes and eukaryotes, multidrug and toxic-compound extrusion (MATE) transporters catalyze the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are secondary-active antiporters, i.e., their drug-efflux activity is coupled to, and powered by, the uptake of ions down a preexisting transmembrane electrochemical gradient. Key aspects of this mechanism, however, remain to be delineated, such as its ion specificity and stoichiometry. We previously revealed the existence of a Na⁺-binding site in a MATE transporter from Pyroccocus furiosus (PfMATE) and hypothesized that this site might be broadly conserved among prokaryotic MATEs. Here, we evaluate this hypothesis by analyzing VcmN and ClbM, which along with PfMATE are the only three prokaryotic MATEs whose molecular structures have been determined at atomic resolution, i.e. better than 3 Å. Reinterpretation of existing crystallographic data and molecular dynamics simulations indeed reveal an occupied Na⁺-binding site in the N-terminal lobe of both structures, analogous to that identified in PfMATE. We likewise find this site to be strongly selective against K⁺, suggesting it is mechanistically significant. Consistent with these computational results, DEER spectroscopy measurements for multiple doubly-spin-labeled VcmN constructs demonstrate Na⁺-dependent changes in protein conformation. The existence of this binding site in three MATE orthologs implicates Na⁺ in the ion-coupled drug-efflux mechanisms of this class of transporters. These results also imply that observations of H⁺-dependent activity likely stem either from a site elsewhere in the structure, or from H⁺ displacing Na⁺ under certain laboratory conditions, as has been noted for other Na⁺-driven transport systems.

Hydrogen-deuterium exchange mass spectrometry captures distinct dynamics upon substrate and inhibitor binding to a transporter

Proton-coupled transporters use transmembrane proton gradients to power active transport of nutrients inside the cell. High-resolution structures often fail to capture the coupling between proton and ligand binding, and conformational changes associated with transport. We combine HDX-MS with mutagenesis and MD simulations to dissect the molecular mechanism of the prototypical transporter XylE. We show that protonation of a conserved aspartate triggers conformational transition from outward-facing to inward-facing state. This transition only occurs in the presence of substrate xylose, while the inhibitor glucose locks the transporter in the outward-facing state. MD simulations corroborate the experiments by showing that only the combination of protonation and xylose binding, and not glucose, sets up the transporter for conformational switch. Overall, we demonstrate the unique ability of HDX-MS to distinguish between the conformational dynamics of inhibitor and substrate binding, and show that a specific allosteric coupling between substrate binding and protonation is a key step to initiate transport.

XylE is a bacterial xylose transporter and homologue of human glucose transporters GLUTs 1-4. HDX-MS, mutagenesis and MD simulations suggest that protonation of a conserved aspartate triggers conformational transition from outward- to inward facing state only in the presence of substrate xylose. In contrast, inhibitor glucose locks the transporter in the outward facing state.

Viewing rare conformations of the β2 adrenergic receptor with pressure-resolved DEER spectroscopy

The β₂ adrenergic receptor (β₂AR) is an archetypal G protein coupled receptor (GPCR). One structural signature of GPCR activation is a large-scale movement (ca. 6 to 14 Å) of transmembrane helix 6 (TM6) to a conformation which binds and activates a cognate G protein. The β₂AR exhibits a low level of agonist-independent G protein activation. The structural origin of this basal activity and its suppression by inverse agonists is unknown but could involve a unique receptor conformation that promotes G protein activation. Alternatively, a conformational selection model proposes that a minor population of the canonical active receptor conformation exists in equilibrium with inactive forms, thus giving rise to basal activity of the ligand-free receptor. Previous spin-labeling and fluorescence resonance energy transfer experiments designed to monitor the positional distribution of TM6 did not detect the presence of the active conformation of ligand-free β₂AR. Here we employ spin-labeling and pressure-resolved double electron-electron resonance spectroscopy to reveal the presence of a minor population of unliganded receptor, with the signature outward TM6 displacement, in equilibrium with inactive conformations. Binding of inverse agonists suppresses this population. These results provide direct structural evidence in favor of a conformational selection model for basal activity in β₂AR and provide a mechanism for inverse agonism. In addition, they emphasize 1) the importance of minor populations in GPCR catalytic function; 2) the use of spin-labeling and variable-pressure electron paramagnetic resonance to reveal them in a membrane protein; and 3) the quantitative evaluation of their thermodynamic properties relative to the inactive forms, including free energy, partial molar volume, and compressibility.

Protein Conformational Dynamics upon Association with the Surfaces of Lipid Membranes and Engineered Nanoparticles: Insights from Electron Paramagnetic Resonance Spectroscopy

Detailed study of conformational rearrangements and dynamics of proteins is central to our understanding of their physiological functions and the loss of function. This review outlines the applications of the electron paramagnetic resonance (EPR) technique to study the structural aspects of proteins transitioning from a solution environment to the states in which they are associated with the surfaces of biological membranes or engineered nanoobjects. In the former case these structural transitions generally underlie functional protein states. The latter case is mostly relevant to the application of protein immobilization in biotechnological industries, developing methods for protein purification, etc. Therefore, evaluating the stability of the protein functional state is particularly important. EPR spectroscopy in the form of continuous-wave EPR or pulse EPR distance measurements in conjunction with protein spin labeling provides highly versatile and sensitive tools to characterize the changes in protein local dynamics as well as large conformational rearrangements. The technique can be widely utilized in studies of both protein-membrane and engineered nanoobject-protein complexes.

Capturing Peptide-GPCR Interactions and Their Dynamics

Many biological functions of peptides are mediated through G protein-coupled receptors (GPCRs). Upon ligand binding, GPCRs undergo conformational changes that facilitate the binding and activation of multiple effectors. GPCRs regulate nearly all physiological processes and are a favorite pharmacological target. In particular, drugs are sought after that elicit the recruitment of selected effectors only (biased ligands). Understanding how ligands bind to GPCRs and which conformational changes they induce is a fundamental step toward the development of more efficient and specific drugs. Moreover, it is emerging that the dynamic of the ligand–receptor interaction contributes to the specificity of both ligand recognition and effector recruitment, an aspect that is missing in structural snapshots from crystallography. We describe here biochemical and biophysical techniques to address ligand–receptor interactions in their structural and dynamic aspects, which include mutagenesis, crosslinking, spectroscopic techniques, and mass-spectrometry profiling. With a main focus on peptide receptors, we present methods to unveil the ligand–receptor contact interface and methods that address conformational changes both in the ligand and the GPCR. The presented studies highlight a wide structural heterogeneity among peptide receptors, reveal distinct structural changes occurring during ligand binding and a surprisingly high dynamics of the ligand–GPCR complexes.

Thermodynamic Basis for Conformational Coupling in an ATP-Binding Cassette Exporter

ATP-binding cassette (ABC) transporters constitute one of the largest protein superfamilies, and they mediate the transport of diverse substrates across the membrane. The molecular mechanism for transducing the energy from ATP binding and hydrolysis into the conformational changes remains elusive. Here, we determined the thermodynamics underlying the ATP-induced global conformational switching for the ABC exporter TmrAB using temperature-resolved pulsed electron-electron double resonance (PELDOR or DEER) spectroscopy. We show that a strong entropy-enthalpy compensation mechanism enables the closure of the nucleotide-binding domains (NBDs) over a wide temperature range. This is mechanically coupled with an outward opening of the transmembrane domains (TMDs) accompanied by an entropy gain. The conserved catalytic glutamate plays a key role in the overall energetics. Our results reveal the thermodynamic basis for the chemomechanical energy coupling in an ABC exporter and present a new strategy to explore the energetics of similar membrane protein complexes.

Electron Paramagnetic Resonance as a Tool for Studying Membrane Proteins

Membrane proteins possess a variety of functions essential to the survival of organisms. However, due to their inherent hydrophobic nature, it is extremely difficult to probe the structure and dynamic properties of membrane proteins using traditional biophysical techniques, particularly in their native environments. Electron paramagnetic resonance (EPR) spectroscopy in combination with site-directed spin labeling (SDSL) is a very powerful and rapidly growing biophysical technique to study pertinent structural and dynamic properties of membrane proteins with no size restrictions. In this review, we will briefly discuss the most commonly used EPR techniques and their recent applications for answering structure and conformational dynamics related questions of important membrane protein systems.

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.

Sequence and structural determinants of ligand-dependent alternating access of a MATE transporter

Multidrug and toxic compound extrusion (MATE) transporters are ubiquitous ion-coupled antiporters that extrude structurally and chemically dissimilar cytotoxic compounds and have been implicated in conferring multidrug resistance. Here, we integrate double electron-electron resonance (DEER) with functional assays and site-directed mutagenesis of conserved residues to illuminate principles of ligand-dependent alternating access of PfMATE, a proton-coupled MATE from the hyperthermophilic archaeon Pyrococcus furiosus Pairs of spin labels monitoring the two sides of the transporter reconstituted into nanodiscs reveal large-amplitude movement of helices that alter the orientation of a putative substrate binding cavity. We found that acidic pH favors formation of an inward-facing (IF) conformation, whereas elevated pH (>7) and the substrate rhodamine 6G stabilizes an outward-facing (OF) conformation. The lipid-dependent PfMATE isomerization between OF and IF conformation is driven by protonation of a previously unidentified intracellular glutamate residue that is critical for drug resistance. Our results can be framed in a mechanistic model of transport that addresses central aspects of ligand coupling and alternating access.

Rapid Simulation of Unprocessed DEER Decay Data for Protein Fold Prediction

Despite advances in sampling and scoring strategies, Monte Carlo modeling methods still struggle to accurately predict de novo the structures of large proteins, membrane proteins, or proteins of complex topologies. Previous approaches have addressed these shortcomings by leveraging sparse distance data gathered using site-directed spin labeling and electron paramagnetic resonance spectroscopy to improve protein structure prediction and refinement outcomes. However, existing computational implementations entail compromises between coarse-grained models of the spin label that lower the resolution and explicit models that lead to resource-intense simulations. These methods are further limited by their reliance on distance distributions, which are calculated from a primary refocused echo decay signal and contain uncertainties that may require manual refinement. Here, we addressed these challenges by developing RosettaDEER, a scoring method within the Rosetta software suite capable of simulating double electron-electron resonance spectroscopy decay traces and distance distributions between spin labels fast enough to fold proteins de novo. We demonstrate that the accuracy of resulting distance distributions match or exceed those generated by more computationally intensive methods. Moreover, decay traces generated from these distributions recapitulate intermolecular background coupling parameters even when the time window of data collection is truncated. As a result, RosettaDEER can discriminate between poorly folded and native-like models by using decay traces that cannot be accurately converted into distance distributions using regularized fitting approaches. Finally, using two challenging test cases, we demonstrate that RosettaDEER leverages these experimental data for protein fold prediction more effectively than previous methods. These benchmarking results confirm that RosettaDEER can effectively leverage sparse experimental data for a wide array of modeling applications built into the Rosetta software suite.

Liquid-phase electron microscopy imaging of cellular and biomolecular systems

Liquid-phase electron microscopy, a new method for real-time nanoscopic imaging in liquid, makes it possible to study cells or biomolecules with a singular combination of spatial and temporal resolution. We review the state of the art in biological research in this growing and promising field.

Mechanism of allosteric modulation of P-glycoprotein by transport substrates and inhibitors

The ATP-binding cassette subfamily B member 1 (ABCB1) multidrug transporter P-glycoprotein plays a central role in clearance of xenobiotics in humans and is implicated in cancer resistance to chemotherapy. We used double electron electron resonance spectroscopy to uncover the basis of stimulation of P-glycoprotein adenosine 5'-triphosphate (ATP) hydrolysis by multiple substrates and illuminate how substrates and inhibitors differentially affect its transport function. Our results reveal that substrate-induced acceleration of ATP hydrolysis correlates with stabilization of a high-energy, post-ATP hydrolysis state characterized by structurally asymmetric nucleotide-binding sites. By contrast, this state is destabilized in the substrate-free cycle and by high-affinity inhibitors in favor of structurally symmetric nucleotide binding sites. Together with previous data, our findings lead to a general model of substrate and inhibitor coupling to P-glycoprotein.