Quantitative Resolution of Monomer-Dimer Populations by Inversion Modulated DEER EPR Spectroscopy

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

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

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.

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.

NMR methods for exploring 'dark' states in ligand binding and protein-protein interactions

A survey, primarily based on work in the authors' laboratory during the last 10 years, is provided of recent developments in NMR studies of exchange processes involving protein-ligand and protein-protein interactions. We start with a brief overview of the theoretical background of Dark state Exchange Saturation Transfer (DEST) and lifetime line-broadening (ΔR₂) NMR methodology. Some limitations of the DEST/ΔR₂ methodology in applications to molecular systems with intermediate molecular weights are discussed, along with the means of overcoming these limitations with the help of closely related exchange NMR techniques, such as the measurements of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion, exchange-induced chemical shifts or rapidly-relaxing components of relaxation decays. Some theoretical underpinnings of the quantitative description of global dynamics of proteins on the surface of very high molecular weight particles (nanoparticles) are discussed. Subsequently, several applications of DEST/ΔR₂ methodology are described from a methodological perspective with an emphasis on providing examples of how kinetic and relaxation parameters for exchanging systems can be reliably extracted from NMR data for each particular model of exchange. Among exchanging systems that are not associated with high molecular weight species, we describe several exchange NMR-based studies that focus on kinetic modelling of transient pre-nucleation oligomerization of huntingtin peptides that precedes aggregation and fibril formation.

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.

Bicelles Rich in both Sphingolipids and Cholesterol and Their Use in Studies of Membrane Proteins

How the distinctive lipid composition of mammalian plasma membranes impacts membrane protein structure is largely unexplored, partly because of the dearth of isotropic model membrane systems that contain abundant sphingolipids and cholesterol. This gap is addressed by showing that sphingomyelin and cholesterol-rich (SCOR) lipid mixtures with phosphatidylcholine can be cosolubilized by n-dodecyl-β-melibioside to form bicelles. Small-angle X-ray and neutron scattering, as well as cryo-electron microscopy, demonstrate that these assemblies are stable over a wide range of conditions and exhibit the bilayered-disc morphology of ideal bicelles even at low lipid-to-detergent mole ratios. SCOR bicelles are shown to be compatible with a wide array of experimental techniques, as applied to the transmembrane human amyloid precursor C99 protein in this medium. These studies reveal an equilibrium between low-order oligomer structures that differ significantly from previous experimental structures of C99, providing an example of how ordered membranes alter membrane protein structure.

Submillisecond Freezing Permits Cryoprotectant-Free EPR Double Electron-Electron Resonance Spectroscopy

Double electron-electron resonance (DEER) EPR spectroscopy is a powerful method for obtaining distance distributions between pairs of engineered nitroxide spin-labels in proteins and other biological macromolecules. These measurements require the use of cryogenic temperatures (77 K or less) to prolong the phase memory relaxation time (T_m ) sufficiently to enable detection of a DEER echo curve. Generally, a cryoprotectant such as glycerol is added to protein samples to facilitate glass formation and avoid protein clustering (which can result in a large decrease in T_m ) during relatively slow flash freezing in liquid N₂ . However, cryoprotectants are osmolytes and can influence protein folding/unfolding equilibria, as well as species populations in weak multimeric systems. Here we show that submillisecond rapid freezing, achieved by high velocity spraying of the sample onto a rapidly spinning, liquid nitrogen cooled copper disc obviates the requirement for cryoprotectants and permits high quality DEER data to be obtained in absence of glycerol. We demonstrate this approach on five different protein systems: protein A, the metastable drkN SH3 domain, urea-unfolded drkN SH3, HIV-1 reverse transcriptase, and the transmembrane domain of HIV-1 gp41 in lipid bicelles.

Spatial domain organization in the HIV-1 reverse transcriptase p66 homodimer precursor probed by double electron-electron resonance EPR

HIV type I (HIV-1) reverse transcriptase (RT) catalyzes the conversion of viral RNA into DNA, initiating the chain of events leading to integration of proviral DNA into the host genome. RT is expressed as a single polypeptide chain within the Gag-Pol polyprotein, and either prior to or following excision by HIV-1 protease forms a 66 kDa chain (p66) homodimer precursor. Further proteolytic attack by HIV-1 protease cleaves the ribonuclease H (RNase H) domain of a single subunit to yield the mature p66/p51 heterodimer. Here, we probe the spatial domain organization within the p66 homodimer using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure a large number of intra- and intersubunit distances between spin labels attached to surface-engineered cysteines. The DEER-derived distances are fully consistent with the structural subunit asymmetry found in the mature p66/p51 heterodimer in which catalytic activity resides in the p66 subunit, while the p51 subunit purely serves as a structural scaffold. Furthermore, the p66 homodimer precursor undergoes a conformational change involving the thumb, palm, and finger domains in one of the subunits (corresponding to the p66 subunit in the mature p66/p51 heterodimer) from a closed to a partially open state upon addition of a nonnucleoside inhibitor. The relative orientation of the domains was modeled by simulated annealing driven by the DEER-derived distances. Finally, the RNase H domain that is cleaved to generate p51 in the mature p66/p51 heterodimer is present in 2 major conformers. One conformer is fully solvent accessible thereby accounting for the observation that only a single subunit of the p66 homodimer precursor is susceptible to HIV-1 protease.

Probing initial transient oligomerization events facilitating Huntingtin fibril nucleation at atomic resolution by relaxation-based NMR

The N-terminal region of the huntingtin protein, encoded by exon-1, comprises an amphiphilic domain (htt^NT), a polyglutamine (Q _n ) tract, and a proline-rich sequence. Polyglutamine expansion results in an aggregation-prone protein responsible for Huntington's disease. Here, we study the earliest events involved in oligomerization of a minimalistic construct, htt^NTQ₇, which remains largely monomeric over a sufficiently long period of time to permit detailed quantitative NMR analysis of the kinetics and structure of sparsely populated [Formula: see text] oligomeric states, yet still eventually forms fibrils. Global fitting of concentration-dependent relaxation dispersion, transverse relaxation in the rotating frame, and exchange-induced chemical shift data reveals a bifurcated assembly mechanism in which the NMR observable monomeric species either self-associates to form a productive dimer (τ_ex ∼ 30 μs, K _diss ∼ 0.1 M) that goes on to form a tetramer ([Formula: see text] μs; K _diss ∼ 22 μM), or exchanges with a "nonproductive" dimer that does not oligomerize further (τ_ex ∼ 400 μs; K _diss ∼ 0.3 M). The excited state backbone chemical shifts are indicative of a contiguous helix (residues 3-17) in the productive dimer/tetramer, with only partial helical character in the nonproductive dimer. A structural model of the productive dimer/tetramer was obtained by simulated annealing driven by intermolecular paramagnetic relaxation enhancement data. The tetramer comprises a D ₂ symmetric dimer of dimers with largely hydrophobic packing between the helical subunits. The structural model, validated by EPR distance measurements, illuminates the role of the htt^NT domain in the earliest stages of prenucleation and oligomerization, before fibril formation.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.

Interaction of Huntingtin Exon-1 Peptides with Lipid-Based Micellar Nanoparticles Probed by Solution NMR and Q-Band Pulsed EPR

Lipid-based micellar nanoparticles promote aggregation of huntingtin exon-1 peptides. Here we characterize the interaction of two such peptides, htt^NTQ ₇ and htt^NTQ ₁₀ comprising the N-terminal amphiphilic domain of huntingtin followed by 7 and 10 glutamine repeats, respectively, with 8 nm lipid micelles using NMR chemical exchange saturation transfer (CEST), circular dichroism and pulsed Q-band EPR. Exchange between free and micelle-bound htt^NTQ  _n peptides occurs on the millisecond time scale with a K_D ∼ 0.5-1 mM. Upon binding micelles, residues 1-15 adopt a helical conformation. Oxidation of Met⁷ to a sulfoxide reduces the binding affinity for micelles ∼3-4-fold and increases the length of the helix by a further two residues. A structure of the bound monomer unit is calculated from the backbone chemical shifts of the micelle-bound state obtained from CEST. Pulsed Q-band EPR shows that a monomer-dimer equilibrium exists on the surface of the micelles and that the two helices of the dimer adopt a parallel orientation, thereby bringing two disordered polyQ tails into close proximity which may promote aggregation upon dissociation from the micelle surface.

Probing Conformational States of the Finger and Thumb Subdomains of HIV-1 Reverse Transcriptase Using Double Electron-Electron Resonance Electron Paramagnetic Resonance Spectroscopy

The configurational space sampled by the finger and thumb subdomains of the p66 subunit of HIV-1 reverse transcriptase was investigated by Q-band double electron-electron resonance pulsed electron paramagnetic resonance spectroscopy, a method for determining long-range distances between pairs of nitroxide spin-labels introduced via surface-engineered cysteine residues. Four constructs were examined, each containing two spin-labels in the p66 subunit, one in the finger subdomain and the other in the thumb subdomain. In the unliganded state, open and closed configurations for the finger and thumb subdomains are observed with the distribution between these states modulated by the spin-labels and associated mutations, in contrast to crystallographic data in which the unliganded state crystallizes in the closed conformation. Upon addition of double-stranded DNA, all constructs adopt open conformations consistent with previous crystallographic data in which the position of the thumb and finger subdomains is determined by contacts with the bound oligonucleotide duplex (DNA or DNA/RNA). Likewise, binary complexes with five different non-nucleoside reverse transcriptase inhibitors are in open or partially open conformations, indicating that binding of the inhibitor to the palm subdomain indirectly restricts the conformational space sampled by the finger and thumb subdomains.

Conformational Heterogeneity in the Activation Mechanism of Bax

Bax is known for its pro-apoptotic role within the mitochondrial pathway of apoptosis. However, the mechanism for transitioning Bax from cytosolic to membrane-bound oligomer remains elusive. Previous nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) studies defined monomeric Bax as conformationally homogeneous. Yet it has recently been proposed that monomeric Bax exists in equilibrium with a minor state that is distinctly different from its NMR structure. Here, we revisited the structural analysis of Bax using methods uniquely suited for unveiling "invisible" states of proteins, namely, NMR paramagnetic relaxation enhancements and EPR double electron-electron resonance (DEER). Additionally we examined the effect of glycerol, the co-solvent of choice in DEER studies, on the structure of Bax using NMR chemical-shift perturbations and residual dipolar couplings. Based on our combined NMR and EPR results, Bax is a conformationally homogeneous protein prior to its activation.