Simulation of nitroxide electron paramagnetic resonance spectra from brownian trajectories and molecular dynamics simulations

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.

Site directed spin labeling to elucidating the mechanism of the cyanobacterial circadian clock

Electron Paramagnetic Resonance (EPR) is a spectroscopic technique that provides structural and dynamic information on unpaired spins and their surrounding environments. Introduction of exogenous spin labels via site directed spin labeling (SDSL) enables characterization of systems of interests lacking intrinsic unpaired spins. This chapter describes the use of SDSL in quantifying KaiB-KaiC binding in the cyanobacterial circadian clock (Kai Clock), exploiting the changes in mobility of the local environment around the spin label on KaiB-KaiC interactions. While the Kai system serves as our model system to demonstrate SDSL-EPR utility in quantifying protein-protein interactions, this technique is readily amenable to other systems of interest whenever specific protein-protein interactions need to be isolated. We first present a protocol for spin labeling KaiB. Then, we detail the sample preparation and acquisition processes to maximize signal-to-noise for downstream analysis. We close this chapter by highlighting recent advances in SDSL technology to incorporate spin labels into proteins of interest and in EPR technology to improve detection sensitivity that may allow greater flexibilities to the types of experiments possible.

Monitoring Protein-Protein Interactions in the Cyanobacterial Circadian Clock in Real Time via Electron Paramagnetic Resonance Spectroscopy

The cyanobacterial circadian clock in Synechococcus elongatus consists of three proteins, KaiA, KaiB, and KaiC. KaiA and KaiB rhythmically interact with KaiC to generate stable oscillations of KaiC phosphorylation with a period of 24 h. The observation of stable circadian oscillations when the three clock proteins are reconstituted and combined in vitro makes it an ideal system for understanding its underlying molecular mechanisms and circadian clocks in general. These oscillations were historically monitored in vitro by gel electrophoresis of reaction mixtures based on the differing electrophoretic mobilities between various phosphostates of KaiC. As the KaiC phospho-distribution represents only one facet of the oscillations, orthogonal tools are necessary to explore other interactions to generate a full description of the system. However, previous biochemical assays are discontinuous or qualitative. To circumvent these limitations, we developed a spin-labeled KaiB mutant that can differentiate KaiC-bound KaiB from free KaiB using continuous-wave electron paramagnetic resonance spectroscopy that is minimally sensitive to KaiA. Similar to wild-type (WT-KaiB), this labeled mutant, in combination with KaiA, sustains robust circadian rhythms of KaiC phosphorylation. This labeled mutant is hence a functional surrogate of WT-KaiB and thus participates in and reports on autonomous macroscopic circadian rhythms generated by mixtures that include KaiA, KaiC, and ATP. Quantitative kinetics could be extracted with improved precision and time resolution. We describe design principles, data analysis, and limitations of this quantitative binding assay and discuss future research necessary to overcome these challenges.

Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories

Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

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.

Trajectory-Based Simulation of EPR Spectra: Models of Rotational Motion for Spin Labels on Proteins

Direct time-domain simulation of continuous-wave (CW) electron paramagnetic resonance (EPR) spectra from molecular dynamics (MD) trajectories has become increasingly popular, especially for proteins labeled with nitroxide spin labels. Due to the time-consuming nature of simulating adequately long MD trajectories, two approximate methods have been developed to reduce the MD-trajectory length required for modeling EPR spectra: hindered Brownian diffusion (HBD) and hidden Markov models (HMMs). Here, we assess the accuracy of these two approximate methods relative to direct simulations from MD trajectories for three spin-labeled protein systems (a simple helical peptide, a soluble protein, and a membrane protein) and two nitroxide spin labels with differing mobilities (R1 and 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC)). We find that the HMMs generally outperform HBD. Although R1 dynamics partially resembles hindered Brownian diffusion, HMMs accommodate the multiple dynamic time scales for the transitions between rotameric states of R1 that cannot be captured accurately by a HBD model. The MD trajectories of the TOAC-labeled proteins show that its dynamics closely resembles slow multisite exchange between twist-boat and chair ring puckering states. This motion is modeled well by HMM but not by HBD. All MD-trajectory data processing, stochastic trajectory simulations, and CW EPR spectral simulations are implemented in EasySpin, a free software package for MATLAB.

Spin Probe Determination of Molecular Orientation Distribution and Rotational Mobility in Liquid Crystals: Model-Free Approach

A model-free approach for simulation of EPR spectra of nitroxide spin probes in liquid-crystalline materials was suggested and used to obtain parameters of molecular orientation and rotational mobility. The developed method is based on experimental recording and numerical simulation of the angular dependence of EPR spectra, which is shown to be much more informative in comparison with a single EPR spectrum. Quantitative spectral simulations considering both local orientational ordering and distribution of local directors in the sample were used for discrimination of models of rotational mobility and orientational alignment. The method was applied for detailed quantitative characterization of axial, orthorhombic, and low-symmetry non-orthorhombic molecular orientation distributions. It is shown that the ordinarily used model of rotational diffusion in a mean-field potential is suitable for the description of molecular mobility and orientational ordering only for relatively low sample temperatures and low-mobility probe molecules with large sizes. In cases of high molecular mobility, the more realistic jump mechanism of molecular moves can be approximately described as quasi-librations. For ordered liquid crystals it was found that mostly the order parameters up to rank 12-14 are essential and easily determined. When well-aligned materials are described, the order parameters up to 18th rank or even higher become meaningful. Both molecular and sample biaxiality is analyzed and quantitatively characterized. The local molecular ordering and sample orientational alignment are quantitatively characterized separately.

Direct Prediction of EPR Spectra from Lipid Bilayers: Understanding Structure and Dynamics in Biological Membranes

Of the many biophysical techniques now being brought to bear on studies of membranes, electron paramagnetic resonance (EPR) of nitroxide spin probes was the first to provide information about both mobility and ordering in lipid membranes. Here, we report the first prediction of variable temperature EPR spectra of model lipid bilayers in the presence and absence of cholesterol from the results of large scale fully atomistic molecular dynamics (MD) simulations. Three types of structurally different spin probes were employed in order to study different parts of the bilayer. Our results demonstrate very good agreement with experiment and thus confirm the accuracy of the latest lipid force fields. The atomic resolution of the simulations allows the interpretation of the molecular motions and interactions in terms of their impact on the sensitive EPR line shapes. Direct versus indirect effects of cholesterol on the dynamics of spin probes are analysed. Given the complexity of structural organisation in lipid bilayers, the advantage of using a combined MD-EPR simulation approach is two-fold. Firstly, prediction of EPR line shapes directly from MD trajectories of actual phospholipid structures allows unambiguous interpretation of EPR spectra of biological membranes in terms of complex motions. Secondly, such an approach provides an ultimate test bed for the up-to-date MD simulation models employed in the studies of biological membranes, an area that currently attracts great attention.

All-atom molecular dynamics simulations of spin labelled double and single-strand DNA for EPR studies

Prediction of motional EPR spectra of spin labelled DNA structures from fully atomistic MD simulations.

Prediction of EPR Spectra of Lyotropic Liquid Crystals using a Combination of Molecular Dynamics Simulations and the Model-Free Approach

We report the first application of fully atomistic molecular dynamics (MD) simulations to the prediction of the motional electron paramagnetic resonance (EPR) spectra of lyotropic liquid crystals in different aggregation states doped with a paramagnetic spin probe. The purpose of this study is twofold. First, given that EPR spectra are highly sensitive to the motions and order of the spin probes doped within lyotropic aggregates, simulation of EPR line shapes from the results of MD modelling provides an ultimate test bed for the force fields currently employed to model such systems. Second, the EPR line shapes are simulated using the motional parameters extracted from MD trajectories using the Model-Free (MF) approach. Thus a combined MD-EPR methodology allowed us to test directly the validity of the application of the MF approach to systems with multi-component molecular motions. All-atom MD simulations using the General AMBER Force Field (GAFF) have been performed on sodium dodecyl sulfate (SDS) and dodecyltrimethylammonium chloride (DTAC) liquid crystals. The resulting MD trajectories were used to predict and interpret the EPR spectra of pre-micellar, micellar, rod and lamellar aggregates. The predicted EPR spectra demonstrate good agreement with most of experimental line shapes thus confirming the validity of both the force fields employed and the MF approach for the studied systems. At the same time simulation results confirm that GAFF tends to overestimate the packing and the order of the carbonyl chains of the surfactant molecules.

The Synergetic Effects of Combining Structural Biology and EPR Spectroscopy on Membrane Proteins

Protein structures as provided by structural biology such as X-ray crystallography, cryo-electron microscopy and NMR spectroscopy are key elements to understand the function of a protein on the molecular level. Nonetheless, they might be error-prone due to crystallization artifacts or, in particular in case of membrane-imbedded proteins, a mostly artificial environment. In this review, we will introduce different EPR spectroscopy methods as powerful tools to complement and validate structural data gaining insights in the dynamics of proteins and protein complexes such that functional cycles can be derived. We will highlight the use of EPR spectroscopy on membrane-embedded proteins and protein complexes ranging from receptors to secondary active transporters as structural information is still limited in this field and the lipid environment is a particular challenge.

A combined EPR and MD simulation study of a nitroxyl spin label with restricted internal mobility sensitive to protein dynamics

EPR studies combined with fully atomistic Molecular Dynamics (MD) simulations and an MD-EPR simulation method provide evidence for intrinsic low rotameric mobility of a nitroxyl spin label, Rn, compared to the more widely employed label MTSL (R1). Both experimental and modelling results using two structurally different sites of attachment to Myoglobin show that the EPR spectra of Rn are more sensitive to the local protein environment than that of MTSL. This study reveals the potential of using the Rn spin label as a reporter of protein motions.

Determination of nitroxide spin label conformations via PELDOR and X-ray crystallography

PELDOR is used to unravel the position and orientation of MTSSL in six singly-labelled azurin mutants. A comparison with X-ray structures of the mutants shows good agreement with respect to the position and orientation of the nitroxide group.

Characterizing the conformational dynamics of metal-free PsaA using molecular dynamics simulations and electron paramagnetic resonance spectroscopy

Prokaryotic metal-ion receptor proteins, or solute-binding proteins, facilitate the acquisition of metal ions from the extracellular environment. Pneumococcal surface antigen A (PsaA) is the primary Mn(2+)-recruiting protein of the human pathogen Streptococcus pneumoniae and is essential for its in vivo colonization and virulence. The recently reported high-resolution structures of metal-free and metal-bound PsaA have provided the first insights into the mechanism of PsaA-facilitated metal binding. However, the conformational dynamics of metal-free PsaA in solution remain unknown. Here, we use continuous wave electron paramagnetic resonance (EPR) spectroscopy and molecular dynamics (MD) simulations to study the relative flexibility of the structural domains in metal-free PsaA and its distribution of conformations in solution. The results show that the crystal structure of metal-free PsaA is a good representation of the dominant conformation in solution, but the protein also samples structurally distinct conformations that are not captured by the crystal structure. Further, these results suggest that the metal binding site is both larger and more solvent exposed than indicated by the metal-free crystal structure. Collectively, this study provides atomic-resolution insight into the conformational dynamics of PsaA prior to metal binding and lays the groundwork for future EPR and MD based studies of PsaA in solution.

Full Atom Simulations of Spin Label Conformations

Interpretation of EPR measurables from spin labels in terms of structure and dynamics requires knowledge of label behavior. General strategies were developed for simulation of labels used in EPR of proteins. The criteria for those simulations are (a) exhaustive sampling of rotamer space, (b) consensus of results independent of starting points, and (c) inclusion of entropy. These criteria are satisfied only when the number of transitions in any dihedral angle exceeds 100 and the simulation maintains thermodynamic equilibrium. Methods such as conventional MD do not efficiently cross energetic barriers, simulated anealing, Monte Carlo or popular Rotamer Library methodologies are potential energy based and ignore entropy (in addition to their specific shortcomings: environment fluctuations, fixed environment, or electrostatics). The Simulated scaling method avoids the above flaws by modulating the force fields between a reduced (allowing crossing energy barriers) and full potential (sampling minima). Spin label diffuses on this surface while remaining in thermodynamic equilibrium. Simulations show that (a) adopting a single conformation is rare, often there are two to four populated rotamers and (b) position of the NO varies up to 16 Å. These results illustrate necessity for caution when interpreting EPR signals in terms of molecular structure. For example, the 10-16 Å distance change in DEER should not be interpreted as a large conformational change, it can well be a flip about Cα-Cβ bond. Rigorous exploration of possible rotamer structures of a spin label is paramount in signal interpretation. We advocate use of bifunctional labels, motion of which is restricted 10,000-fold and the NO position is restricted to 2-5 Å.

Automated structure refinement for a protein heterodimer complex using limited EPR spectroscopic data and a rigid-body docking algorithm: a three-dimensional model for an ankyrin-CDB3 complex

We report here specialized functions incorporated recently in the rigid-body docking software toolkit TagDock to utilize electron paramagnetic resonance derived (EPR-derived) interresidue distance measurements and spin-label accessibility data. The TagDock package extensions include a custom methanethiosulfonate spin label rotamer library to enable explicit, all-atom spin-label side-chain modeling and scripts to evaluate spin-label surface accessibility. These software enhancements enable us to better utilize the biophysical data routinely available from various spin-labeling experiments. To illustrate the power and utility of these tools, we report the refinement of an ankyrin:CDB3 complex model that exhibits much improved agreement with the EPR distance measurements, compared to model structures published previously.

Simulation of electron paramagnetic resonance spectra of spin-labeled molecules from replica-exchange molecular dynamics

We report a general approach for the simulation of the electron paramagnetic resonance (EPR) spectra of spin labels attached to peptides and proteins directly from replica-exchange molecular dynamics (REMD) trajectories. Conventional MD trajectories are generally inadequate for the prediction of EPR line shapes since the label can become trapped in one or more of a set of rotameric states, thus preventing both conformational sampling and accurate estimates of the exchange rates between different rotamers. The advantage of using REMD is its ability to provide both efficient conformational sampling and kinetic information for spin-label dynamics. Our approach is illustrated with spin-labeled peptide. This approach to REMD-EPR simulation paves the way for the wider application of MD modeling to the simulation and interpretation of EPR spectra of spin-labeled molecules.

Conformational dynamics and distribution of nitroxide spin labels

Long-range distance measurements based on paramagnetic relaxation enhancement (PRE) in NMR, quantification of surface water dynamics near biomacromolecules by Overhauser dynamic nuclear polarization (DNP) and sensitivity enhancement by solid-state DNP all depend on introducing paramagnetic species into an otherwise diamagnetic NMR sample. The species can be introduced by site-directed spin labeling, which offers precise control for positioning the label in the sequence of a biopolymer. However, internal flexibility of the spin label gives rise to dynamic processes that potentially influence PRE and DNP behavior and leads to a spatial distribution of the electron spin even in solid samples. Internal dynamics of spin labels and their static conformational distributions have been studied mainly by electron paramagnetic resonance spectroscopy and molecular dynamics simulations, with a large body of results for the most widely applied methanethiosulfonate spin label MTSL. These results are critically discussed in a unifying picture based on rotameric states of the group that carries the spin label. Deficiencies in our current understanding of dynamics and conformations of spin labeled groups and of their influence on NMR observables are highlighted and directions for further research suggested.

Tunable, mixed-resolution modeling using library-based Monte Carlo and graphics processing units

Building on our recently introduced library-based Monte Carlo (LBMC) approach, we describe a flexible protocol for mixed coarse-grained (CG)/all-atom (AA) simulation of proteins and ligands. In the present implementation of LBMC, protein side chain configurations are pre-calculated and stored in libraries, while bonded interactions along the backbone are treated explicitly. Because the AA side chain coordinates are maintained at minimal run-time cost, arbitrary sites and interaction terms can be turned on to create mixed-resolution models. For example, an AA region of interest such as a binding site can be coupled to a CG model for the rest of the protein. We have additionally developed a hybrid implementation of the generalized Born/surface area (GBSA) implicit solvent model suitable for mixed-resolution models, which in turn was ported to a graphics processing unit (GPU) for faster calculation. The new software was applied to study two systems: (i) the behavior of spin labels on the B1 domain of protein G (GB1) and (ii) docking of randomly initialized estradiol configurations to the ligand binding domain of the estrogen receptor (ERα). The performance of the GPU version of the code was also benchmarked in a number of additional systems.

Direct simulation of magnetic resonance relaxation rates and line shapes from molecular trajectories

We simulate spin relaxation processes, which may be measured by either continuous wave or pulsed magnetic resonance techniques, using trajectory-based simulation methodologies. The spin-lattice relaxation rates are extracted numerically from the relaxation simulations. The rates obtained from the numerical fitting of the relaxation curves are compared to those obtained by direct simulation from the relaxation Bloch-Wangsness-Abragam-Redfield theory (BWART). We have restricted our study to anisotropic rigid-body rotational processes, and to the chemical shift anisotropy (CSA) and a single spin-spin dipolar (END) coupling mechanisms. Examples using electron paramagnetic resonance (EPR) nitroxide and nuclear magnetic resonance (NMR) deuterium quadrupolar systems are provided. The objective is to compare those rates obtained by numerical simulations with the rates obtained by BWART. There is excellent agreement between the simulated and BWART rates for a Hamiltonian describing a single spin (an electron) interacting with the bath through the chemical shift anisotropy (CSA) mechanism undergoing anisotropic rotational diffusion. In contrast, when the Hamiltonian contains both the chemical shift anisotropy (CSA) and the spin-spin dipolar (END) mechanisms, the decay rate of a single exponential fit of the simulated spin-lattice relaxation rate is up to a factor of 0.2 smaller than that predicted by BWART. When the relaxation curves are fit to a double exponential, the slow and fast rates extracted from the decay curves bound the BWART prediction. An extended BWART theory, in the literature, includes the need for multiple relaxation rates and indicates that the multiexponential decay is due to the combined effects of direct and cross-relaxation mechanisms.

Nitroxide sensing of a DNA microenvironment: mechanistic insights from EPR spectroscopy and molecular dynamics simulations

The behavior of the nitroxide spin labels 1-oxyl-4-bromo-2,2,5,5-tetramethylpyrroline (R5a) and 1-oxyl-2,2,5,5-tetramethylpyrroline (R5) attached at a phosphorothioate-substituted site in a DNA duplex is modulated by the DNA in a site- and stereospecific manner. A better understanding of the mechanisms of R5a/R5 sensing of the DNA microenvironment will enhance our capability to relate information from nitroxide spectra to sequence-dependent properties of DNA. Toward this goal, electron paramagnetic resonance (EPR) spectroscopy and molecular dynamics (MD) simulations were used to investigate R5 and R5a attached as R(p) and S(p) diastereomers at phosphorothioate (pS)C(7) of d(CTACTG(pS)C(7)Y(8)TTAG). d(CTAAAGCAGTAG) (Y = T or U). X-band continuous-wave EPR spectra revealed that the dT(8) to dU(8) change alters nanosecond rotational motions of R(p)-R5a but produces no detectable differences for S(p)-R5a, R(p)-R5, and S(p)-R5. MD simulations were able to qualitatively account for these spectral variations and provide a plausible physical basis for the R5/R5a behavior. The simulations also revealed a correlation between DNA backbone B(I)/B(II) conformations and R5/R5a rotational diffusion, thus suggesting a direct connection between DNA local backbone dynamics and EPR-detectable R5/R5a motion. These results advance our understanding of how a DNA microenvironment influences nitroxide motion and the observed EPR spectra. This may enable use of R5/R5a for a quantitative description of the sequence-dependent properties of large biologically relevant DNA molecules.

Simulating the dynamics and orientations of spin-labeled side chains in a protein-DNA complex

Site-directed spin labeling, wherein a nitroxide side chain is introduced into a protein at a selected mutant site, is increasingly employed to investigate biological systems by electron spin resonance (ESR) spectroscopy. An understanding of the packing and dynamics of the spin label is needed to extract the biologically relevant information about the macromolecule from ESR measurements. In this work, molecular dynamics (MD) simulations were performed on the spin-labeled restriction endonuclease, EcoRI in complex with DNA. Mutants of this homodimeric enzyme were previously constructed, and distance measurements were performed using the double electron electron resonance experiment. These correlated distance constraints have been leveraged with MD simulations to learn about side chain packing and preferred conformers of the spin label on sites in an α-helix and a β-strand. We found three dihedral angles of the spin label side chain to be most sensitive to the secondary structure where the spin label was located. Conformers sampled by the spin label differed between secondary structures as well. C(α)-C(α) distance distributions were constructed and used to extract details about the protein backbone mobility at the two spin labeled sites. These simulation studies enhance our understanding of the behavior of spin labels in proteins and thus expand the ability of ESR spectroscopy to contribute to knowledge of protein structure and dynamics.

MtsslWizard: In Silico Spin-Labeling and Generation of Distance Distributions in PyMOL

MtsslWizard is a computer program, which operates as a plugin for the PyMOL molecular graphics system. MtsslWizard estimates distances between spin labels on proteins quickly with user-configurable options through a simple graphical interface. In default mode, the program searches for ensembles of possible MTSSL conformations that do not clash with a static model of the protein. Once conformations are assigned, distance distributions between two or more ensembles are calculated, displayed, and can be exported to other software. The program's use is evaluated in a number of challenging test cases and its strengths and weaknesses evaluated. The benefits of the program are its accuracy and simplicity. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00723-012-0314-0) contains supplementary material, which is available to authorized users.

Distance measurements on a dual-labeled TOAC AChR M2δ peptide in mechanically aligned DMPC bilayers via dipolar broadening CW-EPR spectroscopy

A membrane alignment technique has been used to measure the distance between two TOAC nitroxide spin labels on the membrane-spanning M2δ, peptide of the nicotinic acetylcholine receptor (AChR), via CW-EPR spectroscopy. The TOAC-labeled M2δ peptides were mechanically aligned using DMPC lipids on a planar quartz support, and CW-EPR spectra were recorded at specific orientations. Global analysis in combination with rigorous spectral simulation was used to simultaneously analyze data from two different sample orientations for both single- and double-labeled peptides. We measured an internitroxide distance of 14.6 Å from a dual TOAC-labeled AChR M2δ peptide at positions 7 and 13 that closely matches with the 14.5 Å distance obtained from a model of the labeled AChR M2δ peptide. In addition, the angles determining the relative orientation of the two nitroxides have been determined, and the results compare favorably with molecular modeling. The global analysis of the data from the aligned samples gives much more precise estimates of the parameters defining the geometry of the two labels than can be obtained from a randomly dispersed sample.

Rotamer libraries of spin labelled cysteines for protein studies

Studies of structure and dynamics of proteins using site-directed spin labelling rely on explicit modelling of spin label conformations. The large computational effort associated with such modelling with molecular dynamics (MD) simulations can be avoided by a rotamer library approach based on a coarse-grained representation of the conformational space of the spin label. We show here that libraries of about 200 rotamers, obtained by iterative projection of a long MD trajectory of the free spin label onto a set of canonical dihedral angles, provide a representation of the underlying trajectory adequate for EPR distance measurements. Rotamer analysis was performed on selected X-ray structures of spin labelled T4 lysozyme mutants to characterize the spin label rotamer ensemble on a single protein site. Furthermore, predictions based on the rotamer library approach are shown to be in nearly quantitative agreement with electron paramagnetic resonance (EPR) distance data on the Na(+)/H(+) antiporter NhaA and on the light-harvesting complex LHCII whose structures are known from independent cryo electron microscopy and X-ray studies, respectively. Suggestions for the selection of labelling sites in proteins are given, limitations of the approach discussed, and requirements for further development are outlined.

Conformational changes in BAK, a pore-forming proapoptotic Bcl-2 family member, upon membrane insertion and direct evidence for the existence of BH3-BH3 contact interface in BAK homo-oligomers

During apoptosis, the pro-apoptotic Bcl-2 family proteins BAK and BAX form large oligomeric pores in the mitochondrial outer membrane. Apoptotic factors, including cytochrome c, are released through these pores from the mitochondrial intermembrane space into the cytoplasm where they initiate the cascade of events leading to cell death. To better understand this pivotal step toward apoptosis, a method was developed to induce membrane permeabilization by BAK in the membrane without using the full-length protein. Using a soluble form of BAK with a hexahistidine tag at the C terminus and a liposomal system containing the Ni(2+)-nitrilotriacetic acid lipid analog that can bind hexahistidine-tagged proteins, BAK oligomers were formed in the presence of the activator protein p7/p15Bid. In this system, we determined the conformational changes in BAK upon membrane insertion by applying the site-directed spin labeling method of EPR to 13 different amino acid locations. Upon membrane insertion, the BH3 domains were reorganized, and the alpha5-alpha6 helical hairpin structure was partially exposed to the membrane environment. The monomer-monomer interface in the oligomeric structure was also mapped by measuring the distance-dependent spin-spin interactions for each residue location. Spin labels attached in the BH3 domain were juxtaposed within 5-10 A distance in the oligomeric form in the membrane. These results are consistent with the current hypothesis that BAK or BAX forms homodimers, and these homodimers assemble into a higher order oligomeric pore. Detailed analyses of the data provide new insights into the structure of the BAX or BAK homodimer.

Probing the effect of transport inhibitors on the conformation of the mitochondrial citrate transport protein via a site-directed spin labeling approach

The present investigation utilized the site-directed spin labeling method of electron paramagnetic resonance (EPR) spectroscopy to identify the effect of citrate, the natural ligand, and transport inhibitors on the conformation of the yeast mitochondrial citrate transport protein (CTP) reconstituted in liposomal vesicles. Spin label was placed at six different locations within the CTP in order to monitor conformational changes that occurred near each of the transporter's two substrate binding sites, as well as at more distant domains within the CTP architecture. We observed that citrate caused little change in the EPR spectra. In contrast the transport inhibitors 1,2,3-benzenetricarboxylate (BTC), pyridoxal 5'-phosphate (PLP), and compound 792949 resulted in spectral changes that indicated a decrease in the flexibility of the attached spin label at each of the six locations tested. The rank order of the immobilizing effect was compound 792949 > PLP > BTC. The four spin-label locations that report on the CTP substrate binding sites displayed the greatest changes in the EPR spectra upon addition of inhibitor. Furthermore, we found that when compound 792949 was added vectorially (i.e., extra- and/or intra-liposomally), the immobilizing effect was mediated nearly exclusively by external reagent. In contrast, upon addition of PLP vectorially, the effect was mediated to a similar extent from both the external and the internal compartments. In combination our data indicate that: i) citrate binding to the CTP substrate binding sites does not alter side-chain and/or backbone mobility in a global manner and is consistent with our expectation that both in the absence and presence of substrate the CTP displays the flexibility required of a membrane transporter; and ii) binding of each of the transport inhibitors tested locked multiple CTP domains into more rigid conformations, thereby exhibiting long-range inter-domain conformational communication. The differential vectorial effects of compound 792949 and PLP are discussed in the context of the CTP homology-modeled structure and potential mechanistic molecular explanations are given.

Spin labeling EPR

Site-directed spin labeling in combination with electron paramagnetic resonance spectroscopy has emerged as an efficient tool to elucidate the structure and conformational dynamics of biomolecules under native-like conditions. This article summarizes the basics as well as recent progress of site-directed spin labeling. Continuous wave EPR spectra analyses and pulse EPR techniques are reviewed with special emphasis on applications to the sensory rhodopsin-transducer complex mediating the photophobic response of the halophilic archaeum Natronomonas pharaonis and the photosynthetic reaction center from Rhodobacter sphaeroides R26.

Molecular dynamics simulations reveal that AEDANS is an inert fluorescent probe for the study of membrane proteins

Computer simulations were carried out of a number of AEDANS-labeled single cysteine mutants of a small reference membrane protein, M13 major coat protein, covering 60% of its primary sequence. M13 major coat protein is a single membrane-spanning, α-helical membrane protein with a relatively large water-exposed region in the N-terminus. In 10-ns molecular dynamics simulations, we analyze the behavior of the AEDANS label and the native tryptophan, which were used as acceptor and donor in previous FRET experiments. The results indicate that AEDANS is a relatively inert environmental probe that can move unhindered through the lipid membrane when attached to a membrane protein.

Multifrequency electron spin resonance spectra of a spin-labeled protein calculated from molecular dynamics simulations

Multifrequency electron spin resonance (ESR) spectra provide a wealth of structural and dynamic information about the local environment of the spin label and, indirectly, about the spin-labeled protein. Relating the features of the observed spectra to the underlying molecular motions and interactions is, however, challenging. To make progress toward a rigorous interpretation of ESR spectra, we perform extensive molecular dynamics (MD) simulations of fully solvated T4 Lysozyme, labeled with the spin label MTSSL at positions 72 and 131. These two sites have been the object of numerous experimental studies and are generally considered as prototypical solvent-exposed sites on the surfaces of alpha-helices. To extend the time window afforded by the MD simulations, stochastic Markov models reflecting the dynamics of the spin label side chains in terms of their rotameric states are constructed from the trajectories. The calculated multifrequency ESR spectra are in very good agreement with experiment for three different magnetic field strengths without adjusting any parameters. During the trajectories, the spin labels interconvert among a fairly large number of conformations and display a propensity to form interactions with protein residues other than their nearest neighbors along the helix. The detailed picture of the spin label emerging from the MD simulations provides useful insight into the molecular origins of the available spectroscopic and crystallographic data.

Parametrization, molecular dynamics simulation, and calculation of electron spin resonance spectra of a nitroxide spin label on a polyalanine alpha-helix

The nitroxide spin label 1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulfonate (MTSSL), commonly used in site-directed spin labeling of proteins, is studied with molecular dynamics (MD) simulations. After developing force field parameters for the nitroxide moiety and the spin label linker, we simulate MTSSL attached to a polyalanine alpha-helix in explicit solvent to elucidate the factors affecting its conformational dynamics. Electron spin resonance spectra at 9 and 250 GHz are simulated in the time domain using the MD trajectories and including global rotational diffusion appropriate for the tumbling of T4 Lysozyme in solution. Analysis of the MD simulations reveals the presence of significant hydrophobic interactions of the spin label with the alanine side chains.

Site-directed spin labeling studies on nucleic acid structure and dynamics

Site-directed spin labeling (SDSL) uses electron paramagnetic resonance (EPR) spectroscopy to monitor the behavior of a stable nitroxide radical attached at specific locations within a macromolecule such as protein, DNA, or RNA. Parameters obtained from EPR measurements, such as internitroxide distances and descriptions of the rotational motion of a nitroxide, provide unique information on features near the labeling site. With recent advances in solid-phase synthesis of nucleic acids and developments in EPR methodologies, particularly pulsed EPR technologies, SDSL has been increasingly used to study the structure and dynamics of DNA and RNA at the level of the individual nucleotides. This chapter summarizes the current SDSL studies on nucleic acids, with discussions focusing on literature from the last decade.