Simulating electron spin resonance spectra of nitroxide spin labels from molecular dynamics and stochastic trajectories

Simulating electron spin resonance spectra of nitroxide spin labels from motional models is necessary for the quantitative analysis of experimental spectra. We present a framework for modeling the spin label dynamics by using trajectories such as those from molecular dynamics (MD) simulations combined with stochastic treatment of the global protein tumbling. This is achieved in the time domain after two efficient numerical integrators are developed: One for the quantal dynamics of the spins and the other for the classical rotational diffusion. For the quantal dynamics, we propagate the relevant part of the spin density matrix in Hilbert space. For the diffusional tumbling, we work with quaternions, which enables the treatment of anisotropic diffusion in a potential expanded as a sum of spherical harmonics. Time-averaging arguments are invoked to bridge the gap between the smaller time step of the MD trajectories and the larger time steps appropriate for the rotational diffusion and/or quantal spin dynamics.     

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.     

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.     

Stearic acid spin labels in lipid bilayers: insight through atomistic simulations

Spin-labeled stearic acid species are commonly used for electron paramagnetic resonance (EPR) studies of cell membranes to investigate phase transitions, fluidity, and other physical properties. In this paper, we use large-scale molecular dynamics simulations to investigate the position and behavior of nitroxide spin labels attached to stearic acid molecules in dipalmitoylphosphatidylcholine (DPPC) bilayers. The results of these studies are potentially very important for the interpretation of EPR spectra, which rely on assumptions about the position of the label in the membrane. Additionally, we investigate the effect of chirality and ionization of the carboxyl group of the label. For a non-ionized species, we observe that spin-label molecules are even able to make flip-flop transitions between the leaflets of the bilayer. Such transitions have been previously observed only in very rare cases in molecular simulations.     

Spectroscopic selection of distance measurements in a protein dimer with mixed nitroxide and Gd3+ spin labels

The pulse DEER (Double Electron-Electron Resonance) technique is frequently applied for measuring nanometer distances between specific sites in biological macromolecules. In this work we extend the applicability of this method to high field distance measurements in a protein assembly with mixed spin labels, i.e. a nitroxide spin label and a Gd(3+) tag. We demonstrate the possibility of spectroscopic selection of distance distributions between two nitroxide spin labels, a nitroxide spin label and a Gd(3+) ion, and two Gd(3+) ions. Gd(3+)-nitroxide DEER measurements possess high potential for W-band long range distance measurements (6 nm) by combining high sensitivity with ease of data analysis, subject to some instrumental improvements.     

Versatile Trityl Spin Labels for Nanometer Distance Measurements on Biomolecules In Vitro and within Cells

Structure determination of biomacromolecules under in‐cell conditions is a relevant yet challenging task. Electron paramagnetic resonance (EPR) distance measurements in combination with site‐directed spin labeling (SDSL) are a valuable tool in this endeavor but the usually used nitroxide spin labels are not well‐suited for in‐cell measurements. In contrast, triarylmethyl (trityl) radicals are highly persistent, exhibit a long relaxation time and a narrow spectral width. Here, the synthesis of a versatile collection of trityl spin labels and their application in in vitro and in‐cell trityl–iron distance measurements on a cytochrome P450 protein are described. The trityl labels show similar labeling efficiencies and better signal‐to‐noise ratios (SNR) as compared to the popular methanethiosulfonate spin label (MTSSL) and enabled a successful in‐cell measurement.     

High-field EPR and ESEEM investigation of the nitrogen quadrupole interaction of nitroxide spin labels in disordered solids: toward differentiation between polarity and proticity matrix effects on protein function

The combination of high-field electron paramagnetic resonance (EPR) with site-directed spin labeling (SDSL) techniques employing nitroxide radicals has turned out to be particularly powerful in revealing subtle changes of the polarity and proticity profiles in proteins enbedded in membranes. This information can be obtained by orientation-selective high-field EPR resolving principal components of the nitroxide Zeeman (g) and hyperfine ( A) tensors of the spin labels attached to specific molecular sites. In contrast to the g- and A-tensors, the (14)N ( I = 1) quadrupole interaction tensor of the nitroxide spin label has not been exploited in EPR for probing effects of the microenvironment of functional protein sites. In this work it is shown that the W-band (95 GHz) high-field electron spin echo envelope modulation (ESEEM) method is well suited for determining with high accuracy the (14)N quadrupole tensor principal components of a nitroxide spin label in disordered frozen solution. By W-band ESEEM the quadrupole components of a five-ring pyrroline-type nitroxide radical in glassy ortho-terphenyl and glycerol solutions have been determined. This radical is the headgroup of the MTS spin label widely used in SDSL protein studies. By DFT calulations and W-band ESEEM experiments it is demonstrated that the Q(yy) value is especially sensitive to the proticity and polarity of the nitroxide environment in H-bonding and nonbonding situations. The quadrupole tensor is shown to be rather insensitive to structural variations of the nitroxide label itself. When using Q(yy) as a testing probe of the environment, its ruggedness toward temperature changes represents an important advantage over the g xx and A(zz) parameters which are usually employed for probing matrix effects on the spin labeled molecular site. Thus, beyond measurenments of g xx and A(zz) of spin labeled protein sites in disordered solids, W-band high-field ESEEM studies of (14)N quadrupole interactions open a new avenue to reliably probe subtle environmental effects on the electronic structure. This is a significant step forward on the way to differentiate between effects from matrix polarity and hydrogen-bond formation.     

Using Markov models to simulate electron spin resonance spectra from molecular dynamics trajectories

Simulating electron spin resonance (ESR) spectra directly from molecular dynamics simulations of a spin-labeled protein necessitates a large number (hundreds or thousands) of relatively long (hundreds of nanoseconds) trajectories. To meet this challenge, we explore the possibility of constructing accurate stochastic models of the spin label dynamics from atomistic trajectories. A systematic, two-step procedure, based on the probabilistic framework of hidden Markov models, is developed to build a discrete-time Markov chain process that faithfully captures the internal spin label dynamics on time scales longer than about 150 ps. The constructed Markov model is used both to gain insight into the long-lived conformations of the spin label and to generate the stochastic trajectories required for the simulation of ESR spectra. The methodology is illustrated with an application to the case of a spin-labeled poly alanine alpha helix in explicit solvent.     

Exploring the local conformational space of a membrane protein by site-directed spin labeling

Molecular modeling based on a hybrid evolutionary optimization and an information condensation algorithm, called GHOST, of spin label ESR spectra was applied to study the structure and dynamics of membrane proteins. The new method is capable of providing detailed molecular information about the conformational space of the spin-labeled segment of the protein in a membrane system. The method is applied to spin-labeled bacteriophage M13 major coat protein, which is used as a model membrane protein. Single cysteine mutants of the coat protein were labeled with nitroxide spin labels and incorporated in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayers. The new computational method allows us to monitor distributions of local spatial constraints and molecular mobility, in addition to information about the location of the protein in a membrane. Furthermore, the results suggest that different local conformations may coexist in the membrane protein. The knowledge of different local conformations may help us to better understand the function-structure relationship of membrane proteins.     

Inclusion complexes of PBN-type nitrone spin traps and their superoxide spin adducts with cyclodextrin derivatives: parallel determination of the association constants by NMR titrations and 2D-EPR simulations

(1)H NMR and electron paramagnetic resonance (EPR) titrations were used to determine the association constants of the complexes of alpha-phenyl-N-tert-butylnitrone (PBN) analogues and their superoxide spin adducts, respectively, with methylated beta-cyclodextrins. A 1:1 stoichiometry for the nitrones with randomly methylated beta-cyclodextrin and 2,6-di-O-methyl-beta-cyclodextrin and 1:1 and 1:2 stoichiometries for the corresponding cyclodextrin-nitroxide complexes were observed. After the superoxide radical spin trapping reaction, EPR titrations afforded the association constants of the corresponding cyclodextrin-nitroxide complexes. Two-dimensional EPR simulations indicated a bimodal inclusion of the nitroxide free radical spin adducts into the cyclodextrins. For all the nitrone-cyclodextrin and nitroxide-cyclodextrin complexes, the association constants were always higher for the nitroxide complexes than for the nitrone complexes. A cooperative system concerning the complexation of the nitroxide spin adduct with a cyclodextrin was evidenced by EPR titrations. The efficiency of the cyclodextrin inclusion technique to trap superoxide and to resist bioreduction by sodium l-ascorbate was also investigated.     

Differential pulse polarography and electron spin resonance spectroscopy of nitroxyl free radicals used as ESR-spin probes

Differential pulse polarography (DPP) and electron spin resonance (ESR) were used to study the influence of substituents and of the pH of the medium on DPP peak potentials (electrochemical reduction) resp. kreduction (chemical reduction) of nitroxyl free radicals. The DPP peak potentials can be used to select the appropriate nitroxide spin label for relevant biochemical and biophysical applications.     

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

We detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton (¹H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels‐Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double‐quantum‐filtered ¹H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD‐proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 62–78     

Site-specific hydration dynamics in the nonpolar core of a molten globule by dynamic nuclear polarization of water

Water-protein interactions play a direct role in protein folding. The chain collapse that accompanies protein folding involves extrusion of water from the nonpolar core. For many proteins, including apomyoglobin (apoMb), hydrophobic interactions drive an initial collapse to an intermediate state before folding to the final structure. However, the debate continues as to whether the core of the collapsed intermediate state is hydrated and, if so, what the dynamic nature of this water is. A key challenge is that protein hydration dynamics is significantly heterogeneous, yet suitable experimental techniques for measuring hydration dynamics with site-specificity are lacking. Here, we introduce Overhauser dynamic nuclear polarization at 0.35 T via site-specific nitroxide spin labels as a unique tool to probe internal and surface protein hydration dynamics with site-specific resolution in the molten globular, native, and unfolded protein states. The (1)H NMR signal enhancement of water carries information about the local dynamics of the solvent within ～10 ? of a spin label. EPR is used synergistically to gain insights on local polarity and mobility of the spin-labeled protein. Several buried and solvent-exposed sites of apoMb are examined, each bearing a covalently bound nitroxide spin label. We find that the nonpoloar core of the apoMb molten globule is hydrated with water bearing significant translational dynamics, only 4-6-fold slower than that of bulk water. The hydration dynamics of the native state is heterogeneous, while the acid-unfolded state bears fast-diffusing hydration water. This study provides a high-resolution glimpse at the folding-dependent nature of protein hydration dynamics.     

Pulsed eldor measurement of nitrogen T1 in spin labels

A 180° pulse is delivered to one hyperfine line of a nitroxide spin label, and the arrival and disappearance of saturation at another hyperfine line is monitored with a second microwave field. Electron and nitrogen nuclear relaxation times are found to be in poor agreement ,vith the electron-nuclear dipolar (END) mechanism.     

Synthesis of new, BODIPY-based sensors and labels

New, 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY) dye based thiol-reactive fluorescent label, fluorescent amino acid, and fluoroionophore compounds with 540-560 nm emission are described. Combination of a BODIPY dye with a nitronyl nitroxide or an imino nitroxide or a bifunctional pyrroline nitroxide furnished a nitric oxide, a redox sensitive molecule and a double (spin and fluorescence) label, respectively.     

Mapping electron paramagnetic resonance spin label conformations by the simulated scaling method

In order to efficiently simulate spin label behavior when attached to the protein backbone we developed a novel approach that enhances local conformational sampling. The simulated scaling (SS) approach (Li, H., et al. J. Chem. Phys. 2007, 126, 24106) couples the random walk of a potential scaling parameter and molecular dynamics in the framework of hybrid Monte Carlo. This approach allows efficient barrier crossings between conformations. The method retains the thermodynamic detailed balance allowing for determination of relative free energies between various conformations. The accuracy of our method was validated by comparison with the recently resolved X-ray crystal structure of a spin labeled T4 lysozyme in which the spin label was in the interior of the protein. Consistent potentials of mean force (PMF) are obtained for the spin label torsion angles to illustrate their behavior in various protein environments: surface, semiburied, and buried. These PMFs reflect the experimentally observed trends and provide the rationale for the spin label dynamics. We have used this method to compare an implicit and explicit solvent model in spin label modeling. The implicit model, which is computationally faster, was found to be in excellent agreement with the explicit solvent treatment. Based on this collection of results, we believe that the presented approach has great potential in the general strategy of describing the behavior of the spin label using molecular modeling and using this information in the interpretation of EPR measurements in terms of protein conformation and dynamics.     

Recent progress in EPR study of spin labeled polymers and spin probed polymer systems

The EPR technique is commonly used for the detection and characterization of paramagnetic centers in chemical science. This method can provides a lot of information, such as identity, structure, dynamics, interaction, orientation, glass transition temperature, adsorption behavior, functionality, phase behavior, nano-inhomogeneities, and conformation of the free-radical portion of the polymer chain. Most polymers intrinsically possess diamagnetic properties, so in order to study polymers with EPR, paramagnetic centers need to be incorporated into the polymer systems. Spin labeling and spin probing are main methods of covalently attaching paramagnetic centers to polymer chains or embedding them in polymer matrices through non-covalent interactions, respectively. Spin labeling and spin probing techniques for polymers and polymer systems (especially with nitroxide radicals) have also been studied, which have a profound impact on polymer science. This review focuses on the continuous wave EPR technique and introduces the recent advances in spin labeled polymers and spin probed polymer systems in EPR research.     

Dynamic nuclear polarization enhanced nuclear magnetic resonance and electron spin resonance studies of hydration and local water dynamics in micelle and vesicle assemblies

We present a unique analysis tool for the selective detection of local water inside soft molecular assemblies (hydrophobic cores, vesicular bilayers, and micellar structures) suspended in bulk water. Through the use of dynamic nuclear polarization (DNP), the (1)H NMR signal of water is amplified, as it interacts with stable radicals that possess approximately 658 times higher spin polarization. We utilized stable nitroxide radicals covalently attached along the hydrophobic tail of stearic acid molecules that incorporate themselves into surfactant-based micelle or vesicle structures. Here, we present a study of local water content and fluid viscosity inside oleate micelles and vesicles and Triton X-100 micelles to serve as model systems for soft molecular assemblies. This approach is unique because the amplification of the NMR signal is performed in bulk solution and under ambient conditions with site-specific spin labels that only detect the water that is directly interacting with the localized spin labels. Continuous wave (cw) electron spin resonance (ESR) analysis provides rotational dynamics of the spin-labeled molecular chain segments and local polarity parameters that can be related to hydration properties, whereas we show that DNP-enhanced (1)H NMR analysis of fluid samples directly provides translational water dynamics and permeability of the local environment probed by the spin label. Our technique therefore has the potential to become a powerful analysis tool, complementary to cw ESR, to study hydration characteristics of surfactant assemblies, lipid bilayers, or protein aggregates, where water dynamics is a key parameter of their structure and function. In this study, we find that there is significant penetration of water inside the oleate micelles with a higher average local water viscosity ( approximately 1.8 cP) than in bulk water, and Triton X-100 micelles and oleate vesicle bilayers mostly exclude water while allowing for considerable surfactant chain motion and measurable water permeation through the soft structure.     

Triradical cation of p-phenylenediamine having two nitroxide radical groups: spin alignment mediated by delocalized spin

The cationic state of a p-phenylenediamine (PDA) molecule having two nitroxide radical groups was prepared and characterized using electrochemical, electron spin resonance (ESR) spectroscopic, and absorption spectroscopic methods. The delocalized intervalence state of the p-phenylendiamine (PDA) moiety was detected in the cationic state. From the pulsed ESR measurements, it was confirmed that the delocalized spin induces parallel spin alignment between the localized two nitroxide groups which are magnetically weakly coupled in the neutral state. It was found that the resulting high-spin alignment does not seriously affect the delocalized intervalence state of the PDA radical cation.     

A distance ruler for RNA using EPR and site-directed spin labeling

As a basic model study for measuring distances in RNA molecules using continuous wave (CW) EPR spectroscopy, site-directed spin-labeled 10-mer RNA duplexes and HIV-1 TAR RNA motifs with various interspin distances were examined. The spin labels were attached to the 2'-NH2 positions of appropriately placed uridines in the duplexes, and interspin distances were measured from both molecular dynamics simulations (MD) and Fourier deconvolution methods (FD). The 10-mer duplexes have interspin distances ranging from 10 A to 30 A based on MD; however, dipolar line broadening of the CW EPR spectrum is only observed for the RNAs for predicted interspin distances of 10-21 A and not for distances over 25 A. The conformational changes in TAR (transactivating responsive region) RNA in the presence and in the absence of different divalent metal ions were monitored by measuring distances between two nucleotides in the bulge region. The predicted interspin distances obtained from the FD method and those from MD calculations match well for both the model RNA duplexes and the structural changes predicted for TAR RNA. These results demonstrate that distance measurement using EPR spectroscopy is a potentially powerful method to help predict the structures of RNA molecules.