chiLife: An open-source Python package for in silico spin labeling and integrative protein modeling

Non-canonical amino acids for site-directed spin labeling of membrane proteins

Membrane proteins remain challenging targets for conventional structural biology techniques because they need to reside within complex hydrophobic lipid environments to maintain proper structure and function. Magnetic resonance combined with site-directed spin labeling is an alternative method that provides atomic-level structural and dynamical information from effects introduced by an electron- or nuclear-based spin label. With the advent of bioorthogonal click chemistries and genetically engineered non-canonical amino acids (ncAAs), options for linking spin probes to biomolecules have substantially broadened outside the conventional cysteine-based labeling scheme. Here, we highlight current strategies to spin-label membrane proteins through ncAAs for nuclear and electron paramagnetic resonance applications. Such advances are critical for developing bioorthogonal spin labeling schemes to achieve in-cell labeling and in-cell measurements of membrane protein conformational dynamics.

Spin labels for 19F ENDOR distance determination: resolution, sensitivity and distance predictability

19F electron-nuclear double resonance (ENDOR) has emerged as an attractive method for determining distance distributions in biomolecules in the range of 0.7-2 nm, which is not easily accessible by pulsed electron dipolar spectroscopy. The 19F ENDOR approach relies on spin labeling, and in this work, we compare various labels' performance. Four protein variants of GB1 and ubiquitin bearing fluorinated residues were labeled at the same site with nitroxide and trityl radicals and a Gd(III) chelate. Additionally, a double-histidine variant of GB1 was labeled with a Cu(II) nitrilotriacetic acid chelate. ENDOR measurements were carried out at W-band (95 GHz) where 19F signals are well separated from 1H signals. Differences in sensitivity were observed, with Gd(III) chelates providing the highest signal-to-noise ratio. The new trityl label, OXMA, devoid of methyl groups, exhibited a sufficiently long phase memory time to provide an acceptable sensitivity. However, the longer tether of this label effectively reduces the maximum accessible distance between the 19F and the Cα of the spin-labeling site. The nitroxide and Cu(II) labels provide valuable additional geometric insights via orientation selection. Prediction of electron-nuclear distances based on the known structures of the proteins were the closest to the experimental values for Gd(III) labels, and distances obtained for Cu(II) labeled GB1 are in good agreement with previously published NMR results. Overall, our results offer valuable guidance for selecting optimal spin labels for 19F ENDOR distance measurement in proteins.

Ligand-Coupled Conformational Changes in a Cyclic Nucleotide-Gated Ion Channel Revealed by Time-Resolved Transition Metal Ion FRET

Cyclic nucleotide-binding domain (CNBD) ion channels play crucial roles in cellular-signaling and excitability and are regulated by the direct binding of cyclic adenosine- or guanosine-monophosphate (cAMP, cGMP). However, the precise allosteric mechanism governing channel activation upon ligand binding, particularly the energetic changes within domains, remains poorly understood. The prokaryotic CNBD channel SthK offers a valuable model for investigating this allosteric mechanism. In this study, we investigated the conformational dynamics and energetics of the SthK C-terminal region using a combination of steady-state and time-resolved transition metal ion Förster resonance energy transfer (tmFRET) experiments. We engineered donor-acceptor pairs at specific sites within a SthK C-terminal fragment by incorporating a fluorescent noncanonical amino acid donor and metal ion acceptors. Measuring tmFRET with fluorescence lifetimes, we determined intramolecular distance distributions in the absence and presence of cAMP or cGMP. The probability distributions between conformational states without and with ligand were used to calculate the changes in free energy ( ΔG ) and differences in free energy change ( ΔΔG ) in the context of a simple four-state model. Our findings reveal that cAMP binding produces large structural changes, with a very favorable ΔΔG . In contrast to cAMP, cGMP behaved as a partial agonist and only weakly promoted the active state. Furthermore, we assessed the impact of protein oligomerization and ionic strength on the structure and energetics of the conformational states. This study demonstrates the effectiveness of time-resolved tmFRET in determining the conformational states and the ligand-dependent energetics of the SthK C-terminal region.

Alternating access of a bacterial homolog of neurotransmitter: sodium symporters determined from AlphaFold2 ensembles and DEER spectroscopy

Neurotransmitter:sodium symporters (NSSs) play critical roles in neural signaling by regulating neurotransmitter uptake into cells powered by sodium electrochemical gradients. Bacterial NSSs orthologs, including MhsT from Bacillus halodurans, have emerged as model systems to understand the structural motifs of alternating access in NSSs and the extent of conservation of these motifs across the family. Here, we apply a computational/experimental methodology to illuminate the conformational landscape of MhsT alternating access. Capitalizing on our recently developed method, Sampling Protein Ensembles and Conformational Heterogeneity with AlphaFold2 (SPEACH_AF), we derived clusters of MhsT models spanning the transition from inward-facing to outward-facing conformations. Systematic application of double electron-electron resonance (DEER) spectroscopy revealed ligand-dependent movements of multiple structural motifs that underpin MhsT's conformational cycle. Remarkably, comparative DEER analysis in detergent micelles and lipid nanodiscs highlights the profound effect of the environment on the energetics of conformational changes. Through experimentally derived selection of collective variables, we present a model of ion and substrate-powered transport by MhsT consistent with the conformational cycle derived from DEER. Our findings not only advance the understanding of MhsT's function but also uncover motifs of conformational dynamics conserved within the broader context of the NSS family and within the LeuT-fold class of transporters. Importantly, our methodological blueprint introduces an approach that can be applied across a diverse spectrum of transporters to describe their conformational landscapes.

Exploring protein structural ensembles: Integration of sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling methods

Under physiological conditions, proteins continuously undergo structural fluctuations on different timescales. Some conformations are only sparsely populated, but still play a key role in protein function. Thus, meaningful structure-function frameworks must include structural ensembles rather than only the most populated protein conformations. To detail protein plasticity, modern structural biology combines complementary experimental and computational approaches. In this review, we survey available computational approaches that integrate sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling techniques to derive all-atom structural models of rare protein conformations. We also propose strategies to increase the reliability and improve efficiency using deep learning approaches, thus advancing the field of integrative structural biology.

Field-domain rapid-scan EPR at 240GHz for studies of protein functional dynamics at room temperature

We present field-domain rapid-scan (RS) electron paramagnetic resonance (EPR) at 8.6T and 240GHz. To enable this technique, we upgraded a home-built EPR spectrometer with an FPGA-enabled digitizer and real-time processing software. The software leverages the Hilbert transform to recover the in-phase (I) and quadrature (Q) channels, and therefore the raw absorptive and dispersive signals, χ' and χ'', from their combined magnitude (I2+Q2). Averaging a magnitude is simpler than real-time coherent averaging and has the added benefit of permitting long-timescale signal averaging (up to at least 2.5×106 scans) because it eliminates the effects of source-receiver phase drift. Our rapid-scan (RS) EPR provides a signal-to-noise ratio that is approximately twice that of continuous wave (CW) EPR under the same experimental conditions, after scaling by the square root of acquisition time. We apply our RS EPR as an extension of the recently reported time-resolved Gd-Gd EPR (TiGGER) [Maity et al., 2023], which is able to monitor inter-residue distance changes during the photocycle of a photoresponsive protein through changes in the Gd-Gd dipolar couplings. RS, opposed to CW, returns field-swept spectra as a function of time with 10ms time resolution, and thus, adds a second dimension to the static field transients recorded by TiGGER. We were able to use RS TiGGER to track time-dependent and temperature-dependent kinetics of AsLOV2, a light-activated phototropin domain found in oats. The results presented here combine the benefits of RS EPR with the improved spectral resolution and sensitivity of Gd chelates at high magnetic fields. In the future, field-domain RS EPR at high magnetic fields may enable studies of other real-time kinetic processes with time resolutions that are otherwise difficult to access in the solution state.

Measuring conformational equilibria in allosteric proteins with time-resolved tmFRET

Proteins are the workhorses of biology, orchestrating a myriad of cellular functions through intricate conformational changes. Protein allostery, the phenomenon where binding of ligands or environmental changes induce conformational rearrangements in the protein, is fundamental to these processes. We have previously shown that transition metal Förster resonance energy transfer (tmFRET) can be used to interrogate the conformational rearrangements associated with protein allostery and have recently introduced novel FRET acceptors utilizing metal-bipyridyl derivatives to measure long (>20 Å) intramolecular distances in proteins. Here, we combine our tmFRET system with fluorescence lifetime measurements to measure the distances, conformational heterogeneity, and energetics of maltose-binding protein, a model allosteric protein. Time-resolved tmFRET captures near-instantaneous snapshots of distance distributions, offering insights into protein dynamics. We show that time-resolved tmFRET can accurately determine distance distributions and conformational heterogeneity of proteins. Our results demonstrate the sensitivity of time-resolved tmFRET in detecting subtle conformational or energetic changes in protein conformations, which are crucial for understanding allostery. In addition, we extend the use of metal-bipyridyl compounds, showing that Cu(phen)2+ can serve as a spin label for pulse dipolar electron paramagnetic resonance (EPR) spectroscopy, a method that also reveals distance distributions and conformational heterogeneity. The EPR studies both establish Cu(phen)2+ as a useful spin label for pulse dipolar EPR and validate our time-resolved tmFRET measurements. Our approach offers a versatile tool for deciphering conformational landscapes and understanding the regulatory mechanisms governing biological processes.

Long-distance tmFRET using bipyridyl- and phenanthroline-based ligands

With the great progress on determining protein structures over the last decade comes a renewed appreciation that structures must be combined with dynamics and energetics to understand function. Fluorescence spectroscopy, specifically Förster resonance energy transfer (FRET), provides a great window into dynamics and energetics due to its application at physiological temperatures and ability to measure dynamics on the ångström scale. We have recently advanced transition metal FRET (tmFRET) to study allosteric regulation of maltose binding protein and have reported measurements of maltose-dependent distance changes with an accuracy of ∼1.5 Å. When paired with the noncanonical amino acid Acd as a donor, our previous tmFRET acceptors were useful over a working distance of 10 to 20 Å. Here, we use cysteine-reactive bipyridyl and phenanthroline compounds as chelators for Fe2+ and Ru2+ to produce novel tmFRET acceptors to expand the working distance to as long as 50 Å, while preserving our ability to resolve even small maltose-dependent changes in distance. We compare our measured FRET efficiencies to predictions based on models using rotameric ensembles of the donors and acceptors to demonstrate that steady-state measurements of tmFRET with our new probes have unprecedented ability to measure conformational rearrangements under physiological conditions.

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.

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

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

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.

Measuring conformational equilibria in allosteric proteins with time-resolved tmFRET

Proteins are the workhorses of biology, orchestrating a myriad of cellular functions through intricate conformational changes. Protein allostery, the phenomenon where binding of ligands or environmental changes induce conformational rearrangements in the protein, is fundamental to these processes. We have previously shown that transition metal Förster resonance energy transfer (tmFRET) can be used to interrogate the conformational rearrangements associated with protein allostery and have recently introduced novel FRET acceptors utilizing metal-bipyridyl derivatives to measure long (>20 Å) intramolecular distances in proteins. Here, we combine our tmFRET system with fluorescence lifetime measurements to measure the distances, conformational heterogeneity, and energetics of maltose binding protein (MBP), a model allosteric protein. Time-resolved tmFRET captures near-instantaneous snapshots of distance distributions, offering insights into protein dynamics. We show that time-resolved tmFRET can accurately determine distance distributions and conformational heterogeneity of proteins. Our results demonstrate the sensitivity of time-resolved tmFRET in detecting subtle conformational or energetic changes in protein conformations, which are crucial for understanding allostery. In addition, we extend the use of metal-bipyridyl compounds, showing Cu(phen)2+ can serve as a spin label for pulse dipolar electron paramagnetic resonance (EPR) spectroscopy, a method which also reveals distance distributions and conformational heterogeneity. The EPR studies both establish Cu(phen)2+ as a useful spin label for pulse dipolar EPR and validate our time-resolved tmFRET measurements. Our approach offers a versatile tool for deciphering conformational landscapes and understanding the regulatory mechanisms governing biological processes.

Long-distance tmFRET using bipyridyl- and phenanthroline-based ligands

With the great progress on determining protein structures over the last decade comes a renewed appreciation that structures must be combined with dynamics and energetics to understand function. Fluorescence spectroscopy, specifically Förster resonance energy transfer (FRET), provides a great window into dynamics and energetics due to its application at physiological temperatures and ability to measure dynamics on the ångström scale. We have recently advanced transition metal FRET (tmFRET) to study allosteric regulation of maltose binding protein and have reported measurements of maltose-dependent distance changes with an accuracy of ~1.5 Å. When paired with the noncanonical amino acid Acd as a donor, our previous tmFRET acceptors were useful over a working distance of 10 Å to 20 Å. Here, we use cysteine-reactive bipyridyl and phenanthroline compounds as chelators for Fe2+ and Ru2+ to produce novel tmFRET acceptors to expand the working distance to as long as 50 Å, while preserving our ability to resolve even small maltose-dependent changes in distance. We compare our measured FRET efficiencies to predictions based on models using rotameric ensembles of the donors and acceptors to demonstrate that steady-state measurements of tmFRET with our new probes have unprecedented ability to measure conformational rearrangements under physiological conditions.

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.

Revealing the DNA Binding Modes of CsoR by EPR Spectroscopy

In pathogens, a unique class of metalloregulator proteins, called gene regulatory proteins, sense specific metal ions that initiate gene transcription of proteins that export metal ions from the cell, thereby preventing toxicity and cell death. CsoR is a metalloregulator protein found in various bacterial systems that "sense" Cu(I) ions with high affinity. Upon copper binding, CsoR dissociates from the DNA promoter region, resulting in initiation of gene transcription. Crystal structures of CsoR in the presence and absence of Cu(I) from various bacterial systems have been reported, suggesting either a dimeric or tetrameric structure of these helical proteins. However, structural information about the CsoR-DNA complex is missing. Here, we applied electron paramagnetic resonance (EPR) spectroscopy to follow the conformational and dynamical changes that Mycobacterium tuberculosis CsoR undergoes upon DNA binding in solution. We showed that the quaternary structure is predominantly dimeric in solution, and only minor conformational and dynamical changes occur in the DNA bound state. Also, labeling of the unresolved C- terminus revealed no significant change in dynamics upon DNA binding. These observations are unique, since for other bacterial copper metalloregulators, such as the MerR and CopY families, major conformational changes were observed upon DNA binding, indicating a different mode of action for this protein family.

Spectroscopically Orthogonal Labelling to Disentangle Site-Specific Nitroxide Label Distributions

Biomolecular applications of pulse dipolar electron paramagnetic resonance spectroscopy (PDS) are becoming increasingly valuable in structural biology. Site-directed spin labelling of proteins is routinely performed using nitroxides, with paramagnetic metal ions and other organic radicals gaining popularity as alternative spin centres. Spectroscopically orthogonal spin labelling using different types of labels potentially increases the information content available from a single sample. When analysing experimental distance distributions between two nitroxide spin labels, the site-specific rotamer information has been projected into the distance and is not readily available, and the contributions of individual labelling sites to the width of the distance distribution are not obvious from the PDS data. Here, we exploit the exquisite precision of labelling double-histidine (dHis) motifs with CuII chelate complexes. The contribution of this label to the distance distribution widths in model protein GB1 has been shown to be negligible. By combining a dHis CuII labelling site with cysteine-specific nitroxide labelling, we gather insights on the label rotamers at two distinct sites, comparing their contributions to distance distributions based on different in silico modelling approaches and structural models. From this study, it seems advisable to consider discrepancies between different in silico modelling approaches when selecting labelling sites for PDS studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00723-023-01611-1.

DEERefiner-assisted structural refinement using pulsed dipolar spectroscopy: a study on multidrug transporter LmrP

Pulsed dipolar spectroscopy, such as double electron-electron resonance (DEER), has been underutilized in protein structure determination, despite its ability to provide valuable spatial information. In this study, we present DEERefiner, a user-friendly MATLAB-based GUI program that enables the modeling of protein structures by combining an initial structure and DEER distance restraints. We illustrate the effectiveness of DEERefiner by successfully modeling the ligand-dependent conformational changes of the proton-drug antiporter LmrP to an extracellular-open-like conformation with an impressive precision of 0.76 Å. Additionally, DEERefiner was able to uncover a previously hypothesized but experimentally unresolved proton-dependent conformation of LmrP, characterized as an extracellular-closed/partially intracellular-open conformation, with a precision of 1.16 Å. Our work not only highlights the ability of DEER spectroscopy to model protein structures but also reveals the potential of DEERefiner to advance the field by providing an accessible and applicable tool for precise protein structure modeling, thereby paving the way for deeper insights into protein function.

Revealing light-induced structural shifts in G-quadruplex-porphyrin complexes: a pulsed dipolar EPR study

The binding of G-quadruplex structures (G4s) with photosensitizers is of considerable importance in medicinal chemistry and drug discovery due to their promising potential in photodynamic therapy applications. G4s can experience structural changes as a result of ligand interactions and light exposure. Understanding these modifications is essential to uncover the fundamental biological roles of the complexes and optimize their therapeutic potential. The structural diversity of G4s makes it challenging to study their complexes with ligands, necessitating the use of various complementary methods to fully understand these interactions. In this study, we introduce, for the first time, the application of laser-induced dipolar EPR as a method to characterize G-quadruplex DNA complexes containing photosensitizers and to investigate light-induced structural modifications in these systems. To demonstrate the feasibility of this approach, we studied complexes of the human telomeric G-quadruplex (HTel-22) with cationic 5,10,15,20-tetrakis(1-methyl-4-pyridinio) porphyrin tetra(p-toluenesulfonate) (TMPyP4). In addition to showcasing a new methodology, we also aimed to provide insights into the mechanisms underlying photoinduced HTel-22/TMPyP4 structural changes, thereby aiding in the advancement of approaches targeting G4s in photodynamic therapy. EPR revealed G-quadruplex unfolding and dimer formation upon light exposure. Our findings demonstrate the potential of EPR spectroscopy for examining G4 complexes with photosensitizers and contribute to a better understanding of G4s' interactions with ligands under light.

Ultrafast Bioorthogonal Spin-Labeling and Distance Measurements in Mammalian Cells Using Small, Genetically Encoded Tetrazine Amino Acids

Site-directed spin-labeling (SDSL)─in combination with double electron-electron resonance (DEER) spectroscopy─has emerged as a powerful technique for determining both the structural states and the conformational equilibria of biomacromolecules. DEER combined with in situ SDSL in live cells is challenging since current bioorthogonal labeling approaches are too slow to allow for complete labeling with low concentrations of spin label prior to loss of signal from cellular reduction. Here, we overcome this limitation by genetically encoding a novel family of small, tetrazine-bearing noncanonical amino acids (Tet-v4.0) at multiple sites in proteins expressed in Escherichia coli and in human HEK293T cells. We achieved specific and quantitative spin-labeling of Tet-v4.0-containing proteins by developing a series of strained trans-cyclooctene (sTCO)-functionalized nitroxides─including a gem-diethyl-substituted nitroxide with enhanced stability in cells─with rate constants that can exceed 106 M-1 s-1. The remarkable speed of the Tet-v4.0/sTCO reaction allowed efficient spin-labeling of proteins in live cells within minutes, requiring only sub-micromolar concentrations of sTCO-nitroxide. DEER recorded from intact cells revealed distance distributions in good agreement with those measured from proteins purified and labeled in vitro. Furthermore, DEER was able to resolve the maltose-dependent conformational change of Tet-v4.0-incorporated and spin-labeled MBP in vitro and support assignment of the conformational state of an MBP mutant within HEK293T cells. We anticipate the exceptional reaction rates of this system, combined with the relatively short and rigid side chains of the resulting spin labels, will enable structure/function studies of proteins directly in cells, without any requirements for protein purification.

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.