Rates and equilibrium constants of the ligand-induced conformational transition of an HCN ion channel protein domain determined by DEER spectroscopy

A high affinity switch for cAMP in the HCN pacemaker channels

Binding of cAMP to Hyperpolarization activated cyclic nucleotide gated (HCN) channels facilitates pore opening. It is unclear why the isolated cyclic nucleotide binding domain (CNBD) displays in vitro lower affinity for cAMP than the full-length channel in patch experiments. Here we show that HCN are endowed with an affinity switch for cAMP. Alpha helices D and E, downstream of the cyclic nucleotide binding domain (CNBD), bind to and stabilize the holo CNBD in a high affinity state. These helices increase by 30-fold cAMP efficacy and affinity measured in patch clamp and ITC, respectively. We further show that helices D and E regulate affinity by interacting with helix C of the CNBD, similarly to the regulatory protein TRIP8b. Our results uncover an intramolecular mechanism whereby changes in binding affinity, rather than changes in cAMP concentration, can modulate HCN channels, adding another layer to the complex regulation of their 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.

Time-Resolved Mn2+ -NO and NO-NO Distance Measurements Reveal That Catalytic Asymmetry Regulates Alternating Access in an ABC Transporter

ATP-binding cassette (ABC) transporters shuttle diverse substrates across biological membranes. Transport is often achieved through a transition between an inward-facing (IF) and an outward-facing (OF) conformation of the transmembrane domains (TMDs). Asymmetric nucleotide-binding sites (NBSs) are present among several ABC subfamilies and their functional role remains elusive. Here we addressed this question using concomitant NO-NO, Mn2+ -NO, and Mn2+ -Mn2+ pulsed electron-electron double-resonance spectroscopy of TmrAB in a time-resolved manner. This type-IV ABC transporter undergoes a reversible transition in the presence of ATP with a significantly faster forward transition. The impaired degenerate NBS stably binds Mn2+ -ATP, and Mn2+ is preferentially released at the active consensus NBS. ATP hydrolysis at the consensus NBS considerably accelerates the reverse transition. Both NBSs fully open during each conformational cycle and the degenerate NBS may regulate the kinetics of this process.

Enhanced sensitivity for pulse dipolar EPR spectroscopy using variable-time RIDME

Pulse dipolar EPR spectroscopy (PDS) measurements are an important complementary tool in structural biology and are increasingly applied to macromolecular assemblies implicated in human health and disease at physiological concentrations. This requires ever higher sensitivity, and recent advances have driven PDS measurements into the mid-nanomolar concentration regime, though optimization and acquisition of such measurements remains experimentally demanding and time expensive. One important consideration is that constant-time acquisition represents a hard limit for measurement sensitivity, depending on the maximum measured distance. Determining this distance a priori has been facilitated by machine-learning structure prediction (AlphaFold2 and RoseTTAFold) but is often confounded by non-representative behaviour in frozen solution that may mandate multiple rounds of optimization and acquisition. Herein, we endeavour to simultaneously enhance sensitivity and streamline PDS measurement optimization to one-step by benchmarking a variable-time acquisition RIDME experiment applied to CuII-nitroxide and CuII-CuII model systems. Results demonstrate marked sensitivity improvements of both 5- and 6-pulse variable-time RIDME of between 2- and 5-fold over the constant-time analogues.

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

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

Exploring protein conformations in vitro and in cell with EPR distance measurements

The electron-electron double resonance (DEER) method, which provides distance distributions between two spin labels, attached site specifically to biomolecules (proteins and nucleic acids), is currently a well-recognized biophysical tool in structural biology. The most commonly used spin labels are based on nitroxide stable radicals, conjugated to the proteins primarily via native or engineered cysteine residues. However, in recent years, new spin labels, along with different labeling chemistries, have been introduced, driven in part by the desire to study structural and dynamical properties of biomolecules in their native environment, the cell. This mini-review focuses on these new spin labels, which allow for DEER on orthogonal spin labels, and on the state of the art methods for in-cell DEER distance measurements.

Simple Cryoprotectant-Free Method to Advance Pulsed Dipolar ESR Spectroscopy for Capturing Protein Conformational Ensembles

Double electron-electron resonance (DEER) is a powerful technique for studying protein conformations. To preserve the room-temperature ensemble, proteins are usually shock-frozen in liquid nitrogen prior to DEER measurements. The use of cryoprotectant additives is, therefore, necessary to ensure the formation of a vitrified state. Here, we present a simple modification of the freezing process using a flexible fused silica microcapillary, which increases the freezing rates and thus enables DEER measurement without the use of cryoprotectants. The Bid protein, which is highly sensitive to cryoprotectant additives, is used as a model. We show that DEER with the simple modification can successfully reveal the cold denaturation of Bid, which was not possible with the conventional DEER preparations. The DEER result reveals the nature of Bid folding. Our method advances DEER for capturing the chemically and thermally induced conformational changes of a protein in a cryoprotectant-free medium.

Conformational Flexibility of the Protein Insertase BamA in the Native Asymmetric Bilayer Elucidated by ESR Spectroscopy

The β-barrel assembly machinery (BAM) consisting of the central β-barrel BamA and four other lipoproteins mediates the folding of the majority of the outer membrane proteins. BamA is placed in an asymmetric bilayer and its lateral gate is suggested to be the functional hotspot. Here we used in situ pulsed electron-electron double resonance spectroscopy to characterize BamA in the native outer membrane. In the detergent micelles, the data is consistent with mainly an inward-open conformation of BamA. The native membrane considerably enhanced the conformational heterogeneity. The lateral gate and the extracellular loop 3 exist in an equilibrium between different conformations. The outer membrane provides a favorable environment for occupying multiple conformational states independent of the lipoproteins. Our results reveal a highly dynamic behavior of the lateral gate and other key structural elements and provide direct evidence for the conformational modulation of a membrane protein in situ.

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.

Unraveling a ligand-induced twist of a homodimeric enzyme by pulsed electron-electron double resonance

Mechanistic insights into protein–ligand interactions can yield chemical tools for modulating protein function and enable their use for therapeutic purposes. For the homodimeric enzyme tRNA‐guanine transglycosylase (TGT), a putative virulence target of shigellosis, ligand binding has been shown by crystallography to transform the functional dimer geometry into an incompetent twisted one. However, crystallographic observation of both end states does neither verify the ligand‐induced transformation of one dimer into the other in solution nor does it shed light on the underlying transformation mechanism. We addressed these questions in an approach that combines site‐directed spin labeling (SDSL) with distance measurements based on pulsed electron–electron double resonance (PELDOR or DEER) spectroscopy. We observed an equilibrium between the functional and twisted dimer that depends on the type of ligand, with a pyranose‐substituted ligand being the most potent one in shifting the equilibrium toward the twisted dimer. Our experiments suggest a dissociation–association mechanism for the formation of the twisted dimer upon ligand binding.

cAMP binding to closed pacemaker ion channels is non-cooperative

Electrical activity in the brain and heart depends on rhythmic generation of action potentials by pacemaker ion channels (HCN) whose activity is regulated by cAMP binding1. Previous work has uncovered evidence for both positive and negative cooperativity in cAMP binding2,3, but such bulk measurements suffer from limited parameter resolution. Efforts to eliminate this ambiguity using single-molecule techniques have been hampered by the inability to directly monitor binding of individual ligand molecules to membrane receptors at physiological concentrations. Here we overcome these challenges using nanophotonic zero-mode waveguides4 to directly resolve binding dynamics of individual ligands to multimeric HCN1 and HCN2 ion channels. We show that cAMP binds independently to all four subunits when the pore is closed, despite a subsequent conformational isomerization to a flip state at each site. The different dynamics in binding and isomerization are likely to underlie physiologically distinct responses of each isoform to cAMP5 and provide direct validation of the ligand-induced flip-state model6-9. This approach for observing stepwise binding in multimeric proteins at physiologically relevant concentrations can directly probe binding allostery at single-molecule resolution in other intact membrane proteins and receptors.

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.

Spatiotemporal Resolution of Conformational Changes in Biomolecules by Combining Pulsed Electron-Electron Double Resonance Spectroscopy with Microsecond Freeze-Hyperquenching

The function of proteins is linked to their conformations that can be resolved with several high-resolution methods. However, only a few methods can provide the temporal order of intermediates and conformational changes, with each having its limitations. Here, we combine pulsed electron-electron double resonance spectroscopy with a microsecond freeze-hyperquenching setup to achieve spatiotemporal resolution in the angstrom range and lower microsecond time scale. We show that the conformational change of the C_α-helix in the cyclic nucleotide-binding domain of the Mesorhizobium loti potassium channel occurs within about 150 μs and can be resolved with angstrom precision. Thus, this approach holds great promise for obtaining 4D landscapes of conformational changes in biomolecules.

A general model to optimise CuII labelling efficiency of double-histidine motifs for pulse dipolar EPR applications

Electron paramagnetic resonance (EPR) distance measurements are making increasingly important contributions to studies of biomolecules underpinning health and disease by providing highly accurate and precise geometric constraints. Combining double-histidine (dH) motifs with CuII spin labels shows promise for further increasing the precision of distance measurements, and for investigating subtle conformational changes. However, non-covalent coordination-based spin labelling is vulnerable to low binding affinity. Dissociation constants of dH motifs for CuII-nitrilotriacetic acid were previously investigated via relaxation induced dipolar modulation enhancement (RIDME), and demonstrated the feasibility of exploiting the dH motif for EPR applications at sub-μM protein concentrations. Herein, the feasibility of using modulation depth quantitation in CuII-CuII RIDME to simultaneously estimate a pair of non-identical independent KD values in such a tetra-histidine model protein is addressed. Furthermore, we develop a general speciation model to optimise CuII labelling efficiency, depending upon pairs of identical or disparate KD values and total CuII label concentration. We find the dissociation constant estimates are in excellent agreement with previously determined values, and empirical modulation depths support the proposed model.

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.

Characterization of the ExoU activation mechanism using EPR and integrative modeling

ExoU, a type III secreted phospholipase effector of Pseudomonas aeruginosa, serves as a prototype to model large, dynamic, membrane-associated proteins. ExoU is synergistically activated by interactions with membrane lipids and ubiquitin. To dissect the activation mechanism, structural homology was used to identify an unstructured loop of approximately 20 residues in the ExoU amino acid sequence. Mutational analyses indicate the importance of specific loop amino acid residues in mediating catalytic activity. Engineered disulfide cross-links show that loop movement is required for activation. Site directed spin labeling EPR and DEER (double electron-electron resonance) studies of apo and holo states demonstrate local conformational changes at specific sites within the loop and a conformational shift of the loop during activation. These data are consistent with the formation of a substrate-binding pocket providing access to the catalytic site. DEER distance distributions were used as constraints in RosettaDEER to construct ensemble models of the loop in both apo and holo states, significantly extending the range for modeling a conformationally dynamic loop.

DeerLab: a comprehensive software package for analyzing dipolar electron paramagnetic resonance spectroscopy data

Dipolar EPR spectroscopy (DEER and other techniques) enables the structural characterization of macromolecular and biological systems by measurement of distance distributions between unpaired electrons on a nanometer scale. The inference of these distributions from the measured signals is challenging due to the ill-posed nature of the inverse problem. Existing analysis tools are scattered over several applications with specialized graphical user interfaces. This renders comparison, reproducibility, and method development difficult. To remedy this situation, we present DeerLab, an open-source software package for analyzing dipolar EPR data that is modular and implements a wide range of methods. We show that DeerLab can perform one-step analysis based on separable non-linear least squares, fit dipolar multi-pathway models to multi-pulse DEER data, run global analysis with non-parametric distributions, and use a bootstrapping approach to fully quantify the uncertainty in the analysis.

Pulse EPR in biological systems - Beyond the expert's courtyard

Application of EPR to biological systems includes many techniques and applications. In this short perspective, which dares to look into the future, I focus on pulse EPR, which is my field of expertise. Generally, pulse EPR techniques can be divided into two main groups: (1) hyperfine spectroscopy, which explores electron-nuclear interactions, and (2) pulse-dipolar (PD) EPR spectroscopy, which is based on electron-electron spin interactions. Here I focus on PD-EPR because it has a better chance of becoming a widely applied, easy-to-use table-top method to study the structural and dynamic aspects of bio-molecules. I will briefly introduce this technique, its current state of the art, the challenges it is facing, and finally I will describe futuristic scenarios of low-cost PD-EPR approaches that can cross the diffusion barrier from the core of experts to the bulk of the scientific community.

Multimerization of Homo sapiens TRPA1 ion channel cytoplasmic domains

The transient receptor potential Ankyrin-1 (TRPA1) ion channel is modulated by myriad noxious stimuli that interact with multiple regions of the channel, including cysteine-reactive natural extracts from onion and garlic which modify residues in the cytoplasmic domains. The way in which TRPA1 cytoplasmic domain modification is coupled to opening of the ion-conducting pore has yet to be elucidated. The cryo-EM structure of TRPA1 revealed a tetrameric C-terminal coiled-coil surrounded by N-terminal ankyrin repeat domains (ARDs), an architecture shared with the canonical transient receptor potential (TRPC) ion channel family. Similarly, structures of the TRP melastatin (TRPM) ion channel family also showed a C-terminal coiled-coil surrounded by N-terminal cytoplasmic domains. This conserved architecture may indicate a common gating mechanism by which modification of cytoplasmic domains can transduce conformational changes to open the ion-conducting pore. We developed an in vitro system in which N-terminal ARDs and C-terminal coiled-coil domains can be expressed in bacteria and maintain the ability to interact. We tested three gating regulators: temperature; the polyphosphate compound IP6; and the covalent modifier allyl isothiocyanate to determine whether they alter N- and C-terminal interactions. We found that none of the modifiers tested abolished ARD-coiled-coil interactions, though there was a significant reduction at 37˚C. We found that coiled-coils tetramerize in a concentration dependent manner, with monomers and trimers observed at lower concentrations. Our system provides a method for examining the mechanism of oligomerization of TRPA1 cytoplasmic domains as well as a system to study the transmission of conformational changes resulting from covalent modification.

Confidence Analysis of DEER Data and Its Structural Interpretation with Ensemble-Biased Metadynamics

Given its ability to measure multicomponent distance distributions between electron-spin probes, double electron-electron resonance (DEER) spectroscopy has become a leading technique to assess the structural dynamics of biomolecules. However, methodologies to evaluate the statistical error of these distributions are not standard, often hampering a rigorous interpretation of the experimental results. Distance distributions are often determined from the experimental DEER data through a mathematical method known as Tikhonov regularization, but this approach makes rigorous error estimates difficult. Here, we build upon an alternative, model-based approach in which the distance probability distribution is represented as a sum of Gaussian components, and use propagation of errors to calculate an associated confidence band. Our approach considers all sources of uncertainty, including the experimental noise, the uncertainty in the fitted background signal, and the limited time span of the data collection. The resulting confidence band reveals the most and least reliable features of the probability distribution, thereby informing the structural interpretation of DEER experiments. To facilitate this interpretation, we also generalize the molecular simulation method known as ensemble-biased metadynamics (EBMetaD). This method, originally designed to generate maximal-entropy structural ensembles consistent with one or more probability distributions, now also accounts for the uncertainty in those target distributions exactly as dictated by their confidence bands. After careful benchmarks, we demonstrate the proposed techniques using DEER results from spin-labeled T4 lysozyme.

Global analysis of complex PELDOR time traces

Pulsed electron-electron double resonance (PELDOR, alternatively called DEER for double electron-electron resonance) pulse sequences allow for the detection of echo decay curves that are modulated by dipole-dipole-coupling frequencies of interacting electron spins. With increasing distance between them, the echo decay needs to be monitored over a progressively extended time period. However, since the echo intensity typically falls off exponentially with increasing time, this might be problematic with respect to the minimum signal-to-noise ratio required for a sound data analysis. In this contribution we present the new PELDOR analysis tool GloPel (Global analysis of PELDOR data), an open-source Python-based application, that allows to extract improved-quality distance distributions from PELDOR data for which no ideal signal-to-noise ratio can be achieved for a very long observation window. By using Tikhonov regularization, GloPel allows for the simultaneous analysis of two time traces acquired for a sample in two different observation time windows, thus taking advantage of both, the typically high signal-to-noise ratio of the time trace acquired at early times of the echo decay, and the best possible background function fitted for the decay at later times, which is in most cases superimposed with considerable noise. In this way, short distances are not overseen in the higher noise of the longer time traces while long distances are not artificially shortened by limiting the observation time window of the experiment. Following our suggested data acquisition procedure, a significant reduction of the measurement time may also be achieved.

Mechanism for the inhibition of the cAMP dependence of HCN ion channels by the auxiliary subunit TRIP8b

TRIP8b, an accessory subunit of hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels, alters both the cell surface expression and cyclic nucleotide dependence of these channels. However, the mechanism by which TRIP8b exerts these dual effects is still poorly understood. In addition to binding to the carboxyl-terminal tripeptide of HCN channels, TRIP8b also binds directly to the cyclic nucleotide-binding domain (CNBD). That interaction, which requires a small central portion of TRIP8b termed TRIP8b_core, is both necessary and sufficient for reducing the cAMP-dependent regulation of HCN channels. Here, using fluorescence anisotropy, we report that TRIP8b binding to the CNBD of HCN2 channels decreases the apparent affinity of cAMP for the CNBD. We explored two possible mechanisms for this inhibition. A noncompetitive mechanism in which TRIP8b inhibits the conformational change of the CNBD associated with cAMP regulation and a competitive mechanism in which TRIP8b and cAMP compete for the same binding site. To test these two mechanisms, we used a combination of fluorescence anisotropy, biolayer interferometry, and double electron-electron resonance spectroscopy. Fitting these models to our fluorescence anisotropy binding data revealed that, surprisingly, the TRIP8b-dependent reduction of cAMP binding to the CNBD can largely be explained by partial competition between TRIP8b and cAMP. On the basis of these findings, we propose that TRIP8b competes with a portion of the cAMP-binding site or distorts the binding site by making interactions with the binding pocket, thus acting predominantly as a competitive antagonist that inhibits the cyclic-nucleotide dependence of HCN channels.