Conformational dynamics and distribution of nitroxide spin labels

Modeling of Cu(II)-based protein spin labels using rotamer libraries

The bifunctional spin label double-histidine copper-(II) capped with nitrilotriacetate [dHis-Cu(II)-NTA], used in conjunction with electron paramagnetic resonance (EPR) methods can provide high-resolution distance data for investigating protein structure and backbone conformational diversity. Quantitative utilization of this data is limited due to a lack of rapid and accurate dHis-Cu(II)-NTA modeling methods that can be used to translate experimental data into modeling restraints. Here, we develop two dHis-Cu(II)-NTA rotamer libraries using a set of recently published molecular dynamics simulations and a semi-empirical meta-dynamics-based conformational ensemble sampling tool for use with the recently developed chiLife bifunctional spin label modeling method. The accuracy of both the libraries and the modeling method are tested by comparing model predictions to experimentally determined distance distributions. We show that this method is accurate with absolute deviation between the predicted and experimental modes between 0.0-1.2 Å with an average of 0.6 Å over the test data used. In doing so, we also validate the generality of the chiLife bifunctional label modeling method. Taken together, the increased structural resolution and modeling accuracy of dHis-Cu(II)-NTA over other spin labels promise improvements in the accuracy and resolution of protein models by EPR.

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.

Differentiating between Label and Protein Conformers in Pulsed Dipolar EPR Spectroscopy with the dHis-Cu2+ (NTA) Motif

Pulsed dipolar EPR spectroscopy (PDS) in combination with site-directed spin labeling is a powerful tool in structural biology. However, the commonly used spin labels are conjugated to biomolecules via rather long and flexible linkers, which hampers the translation of distance distributions into biomolecular conformations. In contrast, the spin label copper(II)-nitrilotriacetic acid [Cu2+ (NTA)] bound to two histidines (dHis) is rigid and yields narrow distance distributions, which can be more easily translated into biomolecular conformations. Here, we use this label on the 71 kDa Yersinia outer protein O (YopO) to decipher whether a previously experimentally observed bimodal distance distribution is due to two conformations of the biomolecule or of the flexible spin labels. Two different PDS experiments, that is, pulsed electron-electron double resonance (PELDOR aka DEER) and relaxation-induced dipolar modulation enhancement (RIDME), yield unimodal distance distribution with the dHis-Cu2+ (NTA) motif; this result suggests that the α-helical backbone of YopO adopts a single conformation in frozen solution. In addition, we show that the Cu2+ (NTA) label preferentially binds to the target double histidine (dHis) sites even in the presence of 22 competing native histidine residues. Our results therefore suggest that the generation of a His-null background is not required for this spin labeling methodology. Together these results highlight the value of the dHis-Cu2+ (NTA) motif in PDS experiments.

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins

Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is an established tool for exploring protein structure and dynamics. Although nitroxide side chains attached to a single cysteine via a disulfide linkage are commonly employed in SDSL-EPR, their internal flexibility complicates applications to monitor slow internal motions in proteins and to structure determination by distance mapping. Moreover, the labile disulfide linkage prohibits the use of reducing agents often needed for protein stability. To enable the application of SDSL-EPR to the measurement of slow internal dynamics, new spin labels with hindered internal motion are desired. Here, we introduce a highly ordered nitroxide side chain, designated R9, attached at a single cysteine residue via a non-reducible thioether linkage. The reaction to introduce R9 is highly selective for solvent-exposed cysteine residues. Structures of R9 at two helical sites in T4 Lysozyme were determined by X-ray crystallography and the mobility in helical sequences was characterized by EPR spectral lineshape analysis, Saturation Transfer EPR, and Saturation Recovery EPR. In addition, interspin distance measurements between pairs of R9 residues are reported. Collectively, all data indicate that R9 will be useful for monitoring slow internal structural fluctuations, and applications to distance mapping via dipolar spectroscopy and relaxation enhancement methods are anticipated.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00723-023-01618-8.

A novel approach to modeling side chain ensembles of the bifunctional spin label RX

We introduce a novel approach to modeling side chain ensembles of bifunctional spin labels. This approach utilizes rotamer libraries to generate side chain conformational ensembles. Because the bifunctional label is constrained by two attachment sites, the label is split into two monofunctional rotamers which are first attached to their respective sites, then rejoined by a local optimization in dihedral space. We validate this method against a set of previously published experimental data using the bifunctional spin label, RX. This method is relatively fast and can readily be used for both experimental analysis and protein modeling, providing significant advantages over modeling bifunctional labels with molecular dynamics simulations. Use of bifunctional labels for site directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy dramatically reduces label mobility, which can significantly improve resolution of small changes in protein backbone structure and dynamics. Coupling the use of bifunctional labels with side chain modeling methods allows for improved quantitative application of experimental SDSL EPR data to protein modeling.

Disrupting Cu trafficking as a potential therapy for cancer

Copper ions play a crucial role in various cellular biological processes. However, these copper ions can also lead to toxicity when their concentration is not controlled by a sophisticated copper-trafficking system. Copper dys-homeostasis has been linked to a variety of diseases, including neurodegeneration and cancer. Therefore, manipulating Cu-trafficking to trigger selective cancer cell death may be a viable strategy with therapeutic benefit. By exploiting combined in silico and experimental strategies, we identified small peptides able to bind Atox1 and metal-binding domains 3-4 of ATP7B proteins. We found that these peptides reduced the proliferation of cancer cells owing to increased cellular copper ions concentration. These outcomes support the idea of harming copper trafficking as an opportunity for devising novel anti-cancer therapies.

Domain-Specific Phase Transitions in a Supramolecular Nanostructure

Understanding thermal phase behavior within nanomaterials can inform their rational design for medical technologies like drug delivery systems and vaccines, as well as for energy technologies and catalysis. This study resolves thermal phases of discrete domains within a supramolecular aramid amphiphile (AA) nanoribbon. Dynamics are characterized by X-band EPR spectroscopy of spin labels positioned at specific sites through the nanoribbon cross-section. The fitting of the electron paramagnetic resonance (EPR) line shapes reveals distinct conformational dynamics, with fastest dynamics at the surface water layer, intermediate dynamics within the flexible cationic head group domain, and slowest dynamics in the interior aramid domain. Measurement of conformational mobility as a function of temperature reveals first- and second-order phase transitions, with melting transitions observed in the surface and head group domains and a temperature-insensitive crystalline phase in the aramid domain. Arrhenius analysis yields activation energies of diffusion at each site. This work demonstrates that distinct thermal phase behaviors between adjacent nanodomains within a supramolecular nanostructure may be resolved and illustrates the utility of EPR spectroscopy for thermal phase characterization of nanostructures.

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.

GPCR structural characterization by NMR spectroscopy in solution

In the human proteome, 826 G-protein-coupled receptors (GPCRs) interact with extracellular stimuli to initiate cascades of intracellular signaling. Determining conformational dynamics and intermolecular interactions are key to understand GPCR function as a basis for drug design. X-ray crystallography and cryo-electron microscopy (cryo-EM) contribute molecular architectures of GPCRs and GPCR-signaling complexes. NMR spectroscopy is complementary by providing information on the dynamics of GPCR structures at physiological temperature. In this review, several NMR approaches in use to probe GPCR dynamics and intermolecular interactions are discussed. The topics include uniform stable-isotope labeling, amino acid residue-selective stable-isotope labeling, site-specific labeling by genetic engineering, the introduction of ¹⁹F-NMR probes, and the use of paramagnetic nitroxide spin labels. The unique information provided by NMR spectroscopy contributes to our understanding of GPCR biology and thus adds to the foundations for rational drug design.

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.

Large-Scale Conformational Changes of FhaC Provide Insights Into the Two-Partner Secretion Mechanism

The Two-Partner secretion pathway mediates protein transport across the outer membrane of Gram-negative bacteria. TpsB transporters belong to the Omp85 superfamily, whose members catalyze protein insertion into, or translocation across membranes without external energy sources. They are composed of a transmembrane β barrel preceded by two periplasmic POTRA domains that bind the incoming protein substrate. Here we used an integrative approach combining in vivo assays, mass spectrometry, nuclear magnetic resonance and electron paramagnetic resonance techniques suitable to detect minor states in heterogeneous populations, to explore transient conformers of the TpsB transporter FhaC. This revealed substantial, spontaneous conformational changes on a slow time scale, with parts of the POTRA2 domain approaching the lipid bilayer and the protein’s surface loops. Specifically, our data indicate that an amphipathic POTRA2 β hairpin can insert into the β barrel. We propose that these motions enlarge the channel and initiate substrate secretion. Our data propose a solution to the conundrum how TpsB transporters mediate protein secretion without the need for cofactors, by utilizing intrinsic protein dynamics.

A Chemical Biology Primer for NMR Spectroscopists

Among structural biology techniques, NMR spectroscopy offers unique capabilities that enable the atomic resolution studies of dynamic and heterogeneous biological systems under physiological and native conditions. Complex biological systems, however, often challenge NMR spectroscopists with their low sensitivity, crowded spectra or large linewidths that reflect their intricate interaction patterns and dynamics. While some of these challenges can be overcome with the development of new spectroscopic approaches, chemical biology can also offer elegant and efficient solutions at the sample preparation stage. In this tutorial, we aim to present several chemical biology tools that enable the preparation of selectively and segmentally labeled protein samples, as well as the introduction of site-specific spectroscopic probes and post-translational modifications. The four tools covered here, namely cysteine chemistry, inteins, native chemical ligation, and unnatural amino acid incorporation, have been developed and optimized in recent years to be more efficient and applicable to a wider range of proteins than ever before. We briefly introduce each tool, describe its advantages and disadvantages in the context of NMR experiments, and offer practical advice for sample preparation and analysis. We hope that this tutorial will introduce beginning researchers in the field to the possibilities chemical biology can offer to NMR spectroscopists, and that it will inspire new and exciting applications in the quest to understand protein function in health and disease.

Integrative ensemble modeling of proteins and their complexes with distance distribution restraints

Many proteins and protein complexes exhibit regions that are intrinsically disordered. Whereas an arsenal of techniques exists to characterize structured proteins or protein regions, characterization of the vast conformational space occupied by intrinsically disordered regions remains a challenging task due the ensemble-averaging nature of many techniques that provide mean value restraints. More representative information can be gained in the form of distribution restraints, such as EPR-derived distance distributions. Previously we developed the ensemble modeling tool MMM, where we partition the macromolecule into structured and unstructured domains and utilize an integrative structural approach with a focus on EPR-derived distance restraints. Here we present the successor program of MMM: MMMx. All the modeling functionality was ported to MMMx and is now accessed by a uniform script format, allowing to combine the different modules at will to modeling pipelines. During the conception of MMMx many of the tools were improved or updated. We discuss the general functionality of MMMx and its modules, and illustrate some of the modeling tools by application examples.

In-Cell Structural Biology by NMR: The Benefits of the Atomic Scale

In-cell structural biology aims at extracting structural information about proteins or nucleic acids in their native, cellular environment. This emerging field holds great promise and is already providing new facts and outlooks of interest at both fundamental and applied levels. NMR spectroscopy has important contributions on this stage: It brings information on a broad variety of nuclei at the atomic scale, which ensures its great versatility and uniqueness. Here, we detail the methods, the fundamental knowledge, and the applications in biomedical engineering related to in-cell structural biology by NMR. We finally propose a brief overview of the main other techniques in the field (EPR, smFRET, cryo-ET, etc.) to draw some advisable developments for in-cell NMR. In the era of large-scale screenings and deep learning, both accurate and qualitative experimental evidence are as essential as ever to understand the interior life of cells. In-cell structural biology by NMR spectroscopy can generate such a knowledge, and it does so at the atomic scale. This review is meant to deliver comprehensive but accessible information, with advanced technical details and reflections on the methods, the nature of the results, and the future of the field.

Is Cys(MTSL) the Best α-Amino Acid Residue to Electron Spin Labeling of Synthetically Accessible Peptide Molecules with Nitroxides?

Electron paramagnetic resonance spectroscopy, particularly its pulse technique double electron-electron resonance (DEER) (also termed PELDOR), is rapidly becoming an extremely useful tool for the experimental determination of side chain-to-side chain distances between free radicals in molecules fundamental for life, such as polypeptides. Among appropriate probes, the most popular are undoubtedly nitroxide electron spin labels. In this context, suitable biosynthetically derived, helical regions of proteins, along with synthetic peptides with amphiphilic properties and antibacterial activities, are the most extensively investigated compounds. A strict requirement for a precise distance measurement has been identified in a minimal dynamic flexibility of the two nitroxide-bearing α-amino acid side chains. To this end, in this study, we have experimentally compared in detail the side-chain mobility properties of the two currently most widely utilized residues, namely, Cys(MTSL) and 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC). In particular, two double-labeled, chemically synthesized 20-mer peptide molecules have been adopted as appropriate templates for our investigation on the determination of the model intramolecular separations. These double-Cys(MTSL) and double-TOAC compounds are both analogues of the almost completely rigid backbone peptide ruler which we have envisaged and 3D structurally analyzed as our original, unlabeled compound. Here, we have clearly found that the TOAC side-chain labels are largely more 3D structurally restricted than the MTSL labels. From this result, we conclude that the TOAC residue offers more precise information than the Cys(MTSL) residue on the side chain-to-side chain distance distribution in synthetically accessible peptide molecules.

A radical approach to radicals

Nitroxide radicals are characterized by a long-lived spin-unpaired electronic ground state and are strongly sensitive to their chemical surroundings. Combined with electron paramagnetic resonance spectroscopy, these electronic features have led to the widespread application of nitroxide derivatives as spin labels for use in studying protein structure and dynamics. Site-directed spin labelling requires the incorporation of nitroxides into the protein structure, leading to a new protein-ligand molecular model. However, in protein crystallographic refinement nitroxides are highly unusual molecules with an atypical chemical composition. Because macromolecular crystallography is almost entirely agnostic to chemical radicals, their structural information is generally less accurate or even erroneous. In this work, proteins that contain an example of a radical compound (Chemical Component Dictionary ID MTN) from the nitroxide family were re-refined by defining its ideal structural parameters based on quantum-chemical calculations. The refinement results show that this procedure improves the MTN ligand geometries, while at the same time retaining higher agreement with experimental data.

The Advantages of EPR Spectroscopy in Exploring Diamagnetic Metal Ion Binding and Transfer Mechanisms in Biological Systems

Electron paramagnetic resonance (EPR) spectroscopy has emerged as an ideal biophysical tool to study complex biological processes. EPR spectroscopy can follow minor conformational changes in various proteins as a function of ligand or protein binding or interactions with high resolution and sensitivity. Resolving cellular mechanisms, involving small ligand binding or metal ion transfer, is not trivial and cannot be studied using conventional biophysical tools. In recent years, our group has been using EPR spectroscopy to study the mechanism underlying copper ion transfer in eukaryotic and prokaryotic systems. This mini-review focuses on our achievements following copper metal coordination in the diamagnetic oxidation state, Cu(I), between biomolecules. We discuss the conformational changes induced in proteins upon Cu(I) binding, as well as the conformational changes induced in two proteins involved in Cu(I) transfer. We also consider how EPR spectroscopy, together with other biophysical and computational tools, can identify the Cu(I)-binding sites. This work describes the advantages of EPR spectroscopy for studying biological processes that involve small ligand binding and transfer between intracellular proteins.

A Comparison of Cysteine-Conjugated Nitroxide Spin Labels for Pulse Dipolar EPR Spectroscopy

The structure-function and materials paradigms drive research on the understanding of structures and structural heterogeneity of molecules and solids from materials science to structural biology. Functional insights into complex architectures are often gained from a suite of complementary physicochemical methods. In the context of biomacromolecular structures, the use of pulse dipolar electron paramagnetic resonance spectroscopy (PDS) has become increasingly popular. The main interest in PDS is providing long-range nanometre distance distributions that allow for identifying macromolecular topologies, validating structural models and conformational transitions as well as docking of quaternary complexes. Most commonly, cysteines are introduced into protein structures by site-directed mutagenesis and modified site-specifically to a spin-labelled side-chain such as a stable nitroxide radical. In this contribution, we investigate labelling by four different commercial labelling agents that react through different sulfur-specific reactions. Further, the distance distributions obtained are between spin-bearing moieties and need to be related to the protein structure via modelling approaches. Here, we compare two different approaches to modelling these distributions for all four side-chains. The results indicate that there are significant differences in the optimum labelling procedure. All four spin-labels show differences in the ease of labelling and purification. Further challenges arise from the different tether lengths and rotamers of spin-labelled side-chains; both influence the modelling and translation into structures. Our comparison indicates that the spin-label with the shortest tether in the spin-labelled side-group, (bis-(2,2,5,5-Tetramethyl-3-imidazoline-1-oxyl-4-yl) disulfide, may be underappreciated and could increase the resolution of structural studies by PDS if labelling conditions are optimised accordingly.

Resolving distance variations by single-molecule FRET and EPR spectroscopy using rotamer libraries

Förster resonance energy transfer (FRET) and electron paramagnetic resonance (EPR) spectroscopy are complementary techniques for quantifying distances in the nanometer range. Both approaches are commonly employed for probing the conformations and conformational changes of biological macromolecules based on site-directed fluorescent or paramagnetic labeling. FRET can be applied in solution at ambient temperature and thus provides direct access to dynamics, especially if used at the single-molecule level, whereas EPR requires immobilization or work at cryogenic temperatures but provides data that can be more reliably used to extract distance distributions. However, a combined analysis of the complementary data from the two techniques has been complicated by the lack of a common modeling framework. Here, we demonstrate a systematic analysis approach based on rotamer libraries for both FRET and EPR labels to predict distance distributions between two labels from a structural model. Dynamics of the fluorophores within these distance distributions are taken into account by diffusional averaging, which improves the agreement with experiment. Benchmarking this methodology with a series of surface-exposed pairs of sites in a structured protein domain reveals that the lowest resolved distance differences can be as small as ∼0.25 nm for both techniques, with quantitative agreement between experimental and simulated transfer efficiencies within a range of ±0.045. Rotamer library analysis thus establishes a coherent way of treating experimental data from EPR and FRET and provides a basis for integrative structural modeling, including studies of conformational distributions and dynamics of biological macromolecules using both techniques.

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.

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.

MMM: Integrative ensemble modeling and ensemble analysis

Proteins and their complexes can be heterogeneously disordered. In ensemble modeling of such systems with restraints from several experimental techniques the following problems arise: (a) integration of diverse restraints obtained on different samples under different conditions; (b) estimation of a realistic ensemble width; (c) sufficient sampling of conformational space; (d) representation of the ensemble by an interpretable number of conformers; (e) recognition of weak order with site resolution. Here, I introduce several tools that address these problems, focusing on utilization of distance distribution information for estimating ensemble width. The RigiFlex approach integrates such information with high-resolution structures of ordered domains and small-angle scattering data. The EnsembleFit module provides moderately sized ensembles by fitting conformer populations and discarding conformers with low population. EnsembleFit balances the loss in fit quality upon combining restraint subsets from different techniques. Pair correlation analysis for residues and local compaction analysis help in feature detection. The RigiFlex pipeline is tested on data simulated from the structure 70 kDa protein-RNA complex RsmE/RsmZ. It recovers this structure with ensemble width and difference from ground truth both being on the order of 4.2 Å. EnsembleFit reduces the ensemble of the proliferating-cell-nuclear-antigen-associated factor p15PAF from 4,939 to 75 conformers while maintaining good fit quality of restraints. Local compaction analysis for the PaaA2 antitoxin from E. coli O157 revealed correlations between compactness and enhanced residual dipolar couplings in the original NMR restraint set.

Secondary structure and transmembrane topology analysis of the N-terminal domain of the inner membrane protein EccE1 from M. smegmatis using site-directed spin labeling EPR

Protein EccE₁ is an essential component of the mycobacterial ESX-1 secretion system, which plays a crucial part in the process of virulence factors secretion, especially for pathogenic mycobacteria such as Mycobacterium tuberculosis. While EccE₁ was previously postulated to be the inner membrane pore-forming unit of a membrane complex through which substrates are transported, the structural properties of EccE₁ remains to be explored. In the present study, systematic Site-Directed Spin Labeling (SDSL) and Electron Paramagnetic Resonance (EPR) spectroscopic studies was carried out to reveal the secondary structure and transmembrane topology of the N-terminal Domain of EccE₁ protein (EccE₁-NTD) from M. smegmatis in detergent micelles. EPR-based mobility and accessibility analysis of the R1 side chain for 64 residue positions of EccE₁-NTD indicates that the transmembrane domain adopts two α-helices spanning Phe7-Cys30 and Leu36-Ile54. A tentative structural topology model of EccE₁-NTD embedded in membrane is also suggested based on EPR spectroscopic data in this study, which will provide further insights into this protein and the ESX secretion systems of mycobacteria.

Biophysical Characterization of Pro-apoptotic BimBH3 Peptides Reveals an Unexpected Capacity for Self-Association

Bcl-2 proteins orchestrate the mitochondrial pathway of apoptosis, pivotal for cell death. Yet, the structural details of the conformational changes of pro- and antiapoptotic proteins and their interactions remain unclear. Pulse dipolar spectroscopy (double electron-electron resonance [DEER], also known as PELDOR) in combination with spin-labeled apoptotic Bcl-2 proteins unveils conformational changes and interactions of each protein player via detection of intra- and inter-protein distances. Here, we present the synthesis and characterization of pro-apoptotic BimBH3 peptides of different lengths carrying cysteines for labeling with nitroxide or gadolinium spin probes. We show by DEER that the length of the peptides modulates their homo-interactions in the absence of other Bcl-2 proteins and solve by X-ray crystallography the structure of a BimBH3 tetramer, revealing the molecular details of the inter-peptide interactions. Finally, we prove that using orthogonal labels and three-channel DEER we can disentangle the Bim-Bim, Bcl-xL-Bcl-xL, and Bim-Bcl-xL interactions in a simplified interactome.

Conformational Dynamics in Extended RGD-Containing Peptides

RGD is a prolific example of a tripeptide used in biomaterials for cell adhesion, but the potency of free or surface-bound RGD tripeptide is orders-of-magnitude less than the RGD domain within natural proteins. We designed a set of peptides with varying lengths, composed of fragments of fibronectin protein whose central three residues are RGD, in order to vary their conformational behavior without changing the binding site's chemical environment. With these peptides, we measure the conformational dynamics and transient structure of the active site. Our studies reveal how flanking residues affect conformational behavior and integrin binding. We find that disorder of the binding site is important to the potency of RGD peptides and that transient hydrogen bonding near the RGD site affects both the energy landscape roughness of the peptides and peptide binding. This phenomenon is independent of longer-range folding interactions and helps explain why short binding sequences, including RGD itself, do not fully replicate the integrin-targeting properties of extracellular matrix proteins. Our studies reinforce that peptide binding is a holistic event and fragments larger than those directly involved in binding should be considered in the design of peptide epitopes for functional biomaterials.

Microsecond dynamics in proteins by two-dimensional ESR: Predictions

Two-dimensional electron-electron double resonance (2D-ELDOR) provides extensive insight into molecular motions. Recent developments permitting experiments at higher frequencies (95 GHz) provide molecular orientational resolution, enabling a clearer description of the nature of the motions. In this work, simulations are provided for the example of domain motions within proteins that are themselves slowly tumbling in solution. These show the nature of the exchange cross-peaks that are predicted to develop in real time from such domain motions. However, we find that the existing theoretical methods for computing 2D-ELDOR experiments over a wide motional range begin to fail seriously when applied to very slow motions characteristic of proteins in solution. One reason is the failure to obtain accurate eigenvectors and eigenvalues of the complex symmetric stochastic Liouville matrices describing the experiment when computed by the efficient Lanczos algorithm in the range of very slow motion. Another, perhaps more serious, issue is that these matrices are "non-normal," such that for the very slow motional range even rigorous diagonalization algorithms do not yield the correct eigenvalues and eigenvectors. We have employed algorithms that overcome both these issues and lead to valid 2D-ELDOR predictions even for motions approaching the rigid limit. They are utilized to describe the development of cross-peaks in 2D-ELDOR at 95 GHz for a particular case of domain motion.

Supramolecular Approach to Electron Paramagnetic Resonance Distance Measurement of Spin-Labeled Proteins

We demonstrate a host-guest molecular recognition approach to advance double electron-electron resonance (DEER) distance measurements of spin-labeled proteins. We synthesized an iodoacetamide derivative of 2,6-diazaadamantane nitroxide (DZD) spin label that could be doubly incorporated into T4 Lysozyme (T4L) by site-directed spin labeling with efficiency up to 50% per cysteine. The rigidity of the fused ring structure and absence of mobile methyl groups increase the spin echo dephasing time (T_m) at temperatures above 80 K. This enables DEER measurements of distances >4 nm in DZD-labeled T4L in glycerol/water at temperatures up to 150 K with increased sensitivity compared to that of a common spin label such as MTSL. Addition of β-cyclodextrin reduces the rotational correlation time of the label, slightly increases T_m, and most importantly, narrows (and slightly lengthens) the interspin distance distributions. The distance distributions are in good agreement with simulated distance distributions obtained by rotamer libraries. These results provide a foundation for developing supramolecular recognition to facilitate long-distance DEER measurements at near physiological temperatures.

2'-Alkynyl spin-labelling is a minimally perturbing tool for DNA structural analysis

The determination of distances between specific points in nucleic acids is essential to understanding their behaviour at the molecular level. The ability to measure distances of 2-10 nm is particularly important: deformations arising from protein binding commonly fall within this range, but the reliable measurement of such distances for a conformational ensemble remains a significant challenge. Using several techniques, we show that electron paramagnetic resonance (EPR) spectroscopy of oligonucleotides spin-labelled with triazole-appended nitroxides at the 2' position offers a robust and minimally perturbing tool for obtaining such measurements. For two nitroxides, we present results from EPR spectroscopy, X-ray crystal structures of B-form spin-labelled DNA duplexes, molecular dynamics simulations and nuclear magnetic resonance spectroscopy. These four methods are mutually supportive, and pinpoint the locations of the spin labels on the duplexes. In doing so, this work establishes 2'-alkynyl nitroxide spin-labelling as a minimally perturbing method for probing DNA conformation.

Modeling of spin–spin distance distributions for nitroxide labeled biomacromolecules

Combining CREST and MD simulations based on GFN-FF for the automated computation of distance distributions for nitroxide labeled (metallo-) proteins.

Structural heterogeneity of α-synuclein fibrils amplified from patient brain extracts

Parkinson’s disease (PD) and Multiple System Atrophy (MSA) are clinically distinctive diseases that feature a common neuropathological hallmark of aggregated α-synuclein. Little is known about how differences in α-synuclein aggregate structure affect disease phenotype. Here, we amplified α-synuclein aggregates from PD and MSA brain extracts and analyzed the conformational properties using fluorescent probes, NMR spectroscopy and electron paramagnetic resonance. We also generated and analyzed several in vitro α-synuclein polymorphs. We found that brain-derived α-synuclein fibrils were structurally different to all of the in vitro polymorphs analyzed. Importantly, there was a greater structural heterogeneity among α-synuclein fibrils from the PD brain compared to those from the MSA brain, possibly reflecting on the greater variability of disease phenotypes evident in PD. Our findings have significant ramifications for the use of non-brain-derived α-synuclein fibrils in PD and MSA studies, and raise important questions regarding the one disease-one strain hypothesis in the study of α-synucleinopathies.

Parkinson’s disease (PD) and Multiple System Atrophy (MSA) are characterized by the pathological accumulation of α-synuclein. Here the authors employ fluorescent probes, electron microscopy and NMR spectroscopy to study the properties of α-synuclein aggregates that were amplified from patient brain extracts and observe a greater structural diversity among PD patients compared to MSA patients.

gem-Diethyl Pyrroline Nitroxide Spin Labels: Synthesis, EPR Characterization, Rotamer Libraries and Biocompatibility

The availability of bioresistant spin labels is crucial for the optimization of site-directed spin labeling protocols for EPR structural studies of biomolecules in a cellular context. As labeling can affect proteins' fold and/or function, having the possibility to choose between different spin labels will increase the probability to produce spin-labeled functional proteins. Here, we report the synthesis and characterization of iodoacetamide- and maleimide-functionalized spin labels based on the gem-diethyl pyrroline structure. The two nitroxide labels are compared to conventional gem-dimethyl analogs by site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy, using two water soluble proteins: T4 lysozyme and Bid. To foster their use for structural studies, we also present rotamer libraries for these labels, compatible with the MMM software. Finally, we investigate the "true" biocompatibility of the gem-diethyl probes comparing the resistance towards chemical reduction of the NO group in ascorbate solutions and E. coli cytosol at different spin concentrations.

In situ observation of conformational dynamics and protein ligand-substrate interactions in outer-membrane proteins with DEER/PELDOR spectroscopy

Observing structure and conformational dynamics of membrane proteins at high-resolution in their native environments is challenging because of the lack of suitable techniques. We have developed an approach for high-precision distance measurements in the nanometer range for outer membrane proteins (OMPs) in intact E. coli and native membranes. OMPs in Gram-negative bacteria rarely have reactive cysteines. This enables in-situ labeling of engineered cysteines with a methanethiosulfonate functionalized nitroxide spin label (MTSL) with minimal background signals. Following overexpression of the target protein, spin labeling is performed with E. coli or isolated outer membranes (OM) under selective conditions. The interspin distances are measured in-situ using pulsed electron-electron double resonance spectroscopy (PELDOR or DEER). The residual background signals, which are problematic for in-situ structural biology, contributes specifically to the intermolecular part of the signal and can be selectively removed to extract the desired interspin distance distribution. The initial cloning stage can take 5–7 d and the subsequent protein expression, OM isolation, spin labeling, PELDOR experiment, and the data analysis typically take 4–5 d. The described protocol provides a general strategy for observing protein-ligand/substrate interactions, oligomerization, and conformational dynamics of OMPs in the native OM and intact E. coli.

EDITORIAL SUMMARY

This protocol describes how to label bacterial outer membrane proteins with spin labels to study conformational changes and their interaction with ligands and substrates in native membranes and cells using Pulsed Electron-Electron Double Resonance (PELDOR or DEER) spectroscopy.

TWEET

A new protocol for studying conformational changes and ligand/substrate interactions of bacterial outer membrane proteins in-situ.

COVER TEASER

Studying membrane protein conformations in-situ

ESR Resolves the C Terminus Structure of the Ligand-free Human Glutathione S-Transferase A1-1

Nitroxide- and Cu²⁺-based electron spin resonance (ESR) are combined to provide insight into the conformational states of the functionally important α-helix of the human glutathione S-transferase A1. Distance measurements on various spin-labeled dimeric human glutathione S-transferase A1-1 all result in bimodal distance distributions, indicating that the C-terminus exists in two distinct conformations in solution, one of which closely matches that found in the crystal structure of the ligand-bound enzyme. These measurements permit the generation of a model of the unliganded conformation. Room temperature ESR indicates that the second conformation has high mobility, potentially enabling the enzyme's high degree of substrate promiscuity. This model is then validated using computational modeling and further Cu²⁺-based ESR distance measurements. Cu²⁺-based ESR also provides evidence that the secondary structure of the second conformation is of helical nature. Addition of S-hexyl glutathione results in a shift in relative populations, favoring the state that is similar to the previously known structure of the ligand-bound enzyme.

Structural Characterization of Biomolecules through Atomistic Simulations Guided by DEER Measurements

Double electron-electron resonance (DEER) is a popular technique that exploits attached spin labels to probe the collective dynamics of biomolecules in a native environment. Like most spectroscopic approaches, DEER detects an ensemble of states accounting for biomolecular dynamics as well as the labels' intrinsic flexibility. Hence, the DEER data alone do not provide high-resolution structural information. To disentangle this variability, we introduce a minimally biased simulation method to sample a structural ensemble that reproduces multiple experimental signals within the uncertainty. In contrast to previous approaches, our method targets the raw data themselves, and thereby it brings forth an unbiased molecular interpretation of the experiments. After validation on the T4 lysozyme, we apply this technique to interpret recent DEER experiments on a membrane transporter binding protein (VcSiaP). The results highlight the large-scale conformational movement that occurs upon substrate binding and reveal that the unbound VcSiaP is more open in solution than the X-ray structure.

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.

Electron spin-labelling of the EutC subunit in B12-dependent ethanolamine ammonia-lyase reveals dynamics and a two-state conformational equilibrium in the N-terminal, signal-sequence-associated domain

The B12 (adenosylcobalamin)-dependent ethanolamine ammonia-lyase (EAL) is a product of the ethanolamine utilisation (eut) gene cluster, that is involved in human gut microbiome homeostasis and in disease conditions caused by pathogenic strains of Salmonella and Escherichia coli. Toward elucidation of the molecular basis of EAL catalysis, and its intracellular trafficking and targeting to the Eut biomicrocompartment (BMC), we have applied electron spin-labelling and electron paramagnetic resonance spectroscopy to wild-type (wt) EAL from Salmonella typhimurium, by using the sulphydryl-specific, 4-maleimido-TEMPO (4MT) spin label. One cysteine residue per active site displays exceptional reactivity with 4MT. This site is identified as βC37 on the EutC subunit, by using 4MT-labeling of site-specific cysteine-to-alanine mutants, enzyme kinetics, and accessible surface area calculations. Electron paramagnetic resonance (EPR) spectra of 4MT-labelled wt EAL are collected over 200-265 K in frozen, polycrystalline water-only, and 1% v/v DMSO solvents. EPR simulations reveal two mobility components for each condition. Detectable spin probe reorientational motion of the two components occurs at 215 and 225 K with 1% v/v DMSO, relative to the water-only condition, consistent with formation of an aqueous-DMSO solvent mesodomain around EAL. Parallel trends in fast- and slow-reorientational correlation times and interconversion of the two populations with increasing temperature, indicate 4MT labelling of a single site (βC37). A two-state model is proposed, in which the fast and slow motional populations represent EAL-bound and free conformations of the EutC N-terminal domain. The approximately equal proportion of each state may represent a balance between EutC and EAL protein stability and efficient targeting to the BMC.

Spin labelling for integrative structure modelling: a case study of the polypyrimidine-tract binding protein 1 domains in complexes with short RNAs

A combined method, employing NMR and EPR spectroscopies, has demonstrated its strength in solving structures of protein/RNA and other types of biomolecular complexes. This method works particularly well when the large biomolecular complex consists of a limited number of rigid building blocks, such as RNA-binding protein domains (RBDs). A variety of spin labels is available for such studies, allowing for conventional as well as spectroscopically orthogonal double electron-electron resonance (DEER) measurements in EPR. In this work, we compare different types of nitroxide-based and Gd(iii)-based spin labels attached to isolated RBDs of the polypyrimidine-tract binding protein 1 (PTBP1) and to short RNA fragments. In particular, we demonstrate experiments on spectroscopically orthogonal labelled RBD/RNA complexes. For all experiments we analyse spin labelling, DEER method performance, resulting distance distributions, and their consistency with the predictions from the spin label rotamers analysis. This work provides a set of intra-domain calibration DEER data, which can serve as a basis to start structure determination of the full length PTBP1 complex with an RNA derived from encephalomycarditis virus (EMCV) internal ribosomal entry site (IRES). For a series of tested labelling sites, we discuss their particular advantages and drawbacks in such a structure determination approach.

A Structural Model of a P450-Ferredoxin Complex from Orientation-Selective Double Electron-Electron Resonance Spectroscopy

Cytochrome P450 (CYP) monooxygenases catalyze the oxidation of chemically inert carbon-hydrogen bonds in diverse endogenous and exogenous organic compounds by atmospheric oxygen. This C-H bond oxy-functionalization activity has huge potential in biotechnological applications. Class I CYPs receive the two electrons required for oxygen activation from NAD(P)H via a ferredoxin reductase and ferredoxin. The interaction of Class I CYPs with their cognate ferredoxin is specific. In order to reconstitute the activity of diverse CYPs, structural characterization of CYP-ferredoxin complexes is necessary, but little structural information is available. Here we report a structural model of such a complex (CYP199A2-HaPux) in frozen solution derived from distance and orientation restraints gathered by the EPR technique of orientation-selective double electron-electron resonance (os-DEER). The long-lived oscillations in the os-DEER spectra were well modeled by a single orientation of the CYP199A2-HaPux complex. The structure is different from the two known Class I CYP-Fdx structures: CYP11A1-Adx and CYP101A1-Pdx. At the protein interface, HaPux residues in the [Fe₂S₂] cluster-binding loop and the α3 helix and the C-terminus residue interact with CYP199A2 residues in the proximal loop and the C helix. These residue contacts are consistent with biochemical data on CYP199A2-ferredoxin binding and electron transfer. Electron-tunneling calculations indicate an efficient electron-transfer pathway from the [Fe₂S₂] cluster to the heme. This new structural model of a CYP-Fdx complex provides the basis for tailoring CYP enzymes for which the cognate ferredoxin is not known, to accept electrons from HaPux and display monooxygenase activity.

The contribution of modern EPR to structural biology

Electron paramagnetic resonance (EPR) spectroscopy combined with site-directed spin labelling is applicable to biomolecules and their complexes irrespective of system size and in a broad range of environments. Neither short-range nor long-range order is required to obtain structural restraints on accessibility of sites to water or oxygen, on secondary structure, and on distances between sites. Many of the experiments characterize a static ensemble obtained by shock-freezing. Compared with characterizing the dynamic ensemble at ambient temperature, analysis is simplified and information loss due to overlapping timescales of measurement and system dynamics is avoided. The necessity for labelling leads to sparse restraint sets that require integration with data from other methodologies for building models. The double electron–electron resonance experiment provides distance distributions in the nanometre range that carry information not only on the mean conformation but also on the width of the native ensemble. The distribution widths are often inconsistent with Anfinsen's concept that a sequence encodes a single native conformation defined at atomic resolution under physiological conditions.

Small neutral Gd(iii) tags for distance measurements in proteins by double electron–electron resonance experiments

Small Gd(iii) tags based on DO3A deliver narrow and readily predictable distances by double electron–electron resonance (DEER) measurements.