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.     

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.     

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.     

Calculating slow-motional electron paramagnetic resonance spectra from molecular dynamics using a diffusion operator approach

A number of groups have utilized molecular dynamics (MD) to calculate slow-motional electron paramagnetic resonance (EPR) spectra of spin labels attached to biomolecules. Nearly all such calculations have been based on some variant of the trajectory method introduced by Robinson, Slutsky and Auteri (J. Chem. Phys. 1992,96, 2609-2616). Here we present an alternative approach that is specifically adapted to the diffusion operator-based stochastic Liouville equation (SLE) formalism that is also widely used to calculate slow-motional EPR line shapes. Specifically, the method utilizes MD trajectories to derive diffusion parameters such as the rotational diffusion tensor, diffusion tilt angles, and expansion coefficients of the orienting potential, which are then used as direct inputs to the SLE line shape program. This approach leads to a considerable improvement in computational efficiency over trajectory-based methods, particularly for high frequency, high field EPR. It also provides a basis for deconvoluting the effects of local spin label motion and overall motion of the labeled molecule or domain: once the local motion has been characterized by this approach, the label diffusion parameters may be used in conjunction with line shape analysis at lower EPR frequencies to characterize global motions. The method is validated by comparison of the MD predicted line shapes to experimental high frequency (250 GHz) EPR spectra.     

On the interaction of bovine serum albumin with ionic surfactants: temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations, HPS presented a smaller effect at concentrations above 1.5mM. At 10mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue.     

Site-directed electrostatic measurements with a thiol-specific ph-sensitive nitroxide: differentiating local pK and polarity effects by high-field EPR

This communication describes the use of a methanethiosulfonate derivative of an imidazolidine nitroxide, methanethiosulfonic acid S-(1-oxyl-2,2,3,5,5-pentamethyl-imidazolidin-4-ylmethyl) ester, IMTSL, for site-directed pKa determination of peptides by electron paramagnetic resonance. This spin label is covalently attached to the thiol group of unique cysteines incorporated into peptide structures. The tertiary amine nitrogen N3 of the label readily participates in proton exchange reactions, which are monitored through changes in EPR spectra of nitroxide moiety. Using EPR at 95 GHz (W-band) isotropic magnetic parameters of this nitroxide, both Aiso and giso, were calibrated in solvents of different polarity and pH. Two different linear correlations between Aiso and giso for acidic and basic forms of IMTSL were observed, making it possible to differentiate effects of local polarity from N3 protonation on nitroxide EPR spectra. Titration of a synthetic P11 peptide fragment of the laminin B1 chain illustrates the utility of this method.     

Motional heterogeneity and phase separation of functionalized polyester polyurethanes

The polyester polyurethanes, PU based on isophoronediisocyanate, polycaprolactone, and 1,4-butanediol with different amounts of functional groups introduced into the hard segments via second chain extender, 2,2′-bis-(hydroxymethyl) propionic acid, were investigated by electron spin resonance, ESR, spin label method, wide-angle X-ray diffraction, WAXD, optical microscopy and differential scanning calorimetry, DSC. The objective of this study is to clarify the effect of functional groups on the motional heterogeneity, microphase separation and crystallisation of the polyurethanes. The concentration of carboxylic groups varied from 0 to 0.45 mmol g⁻¹. The temperature-dependent ESR spectra of spin labelled PU hard segments chain ends with stable nitroxide radical 2,2,6,6-tetramethyl-4-aminopiperidin-1-yloxyl are sensitive to the amount of functional groups attached to the hard segments. Composite ESR spectra of functionalized PU, with fast and slow component, suggest that PU hard segments are partitioned in two motionally different environments. According to the ratio of fast and slow component motional heterogeneity increases with an increase of functional groups up to 0.35 mmol g⁻¹ and above this concentration slow component decreases indicating higher degree of phase mixing and stronger effect of soft segments. Polarized micrographs and the extent of ordering from WAXD measurements reveal the changes of phase morphology with the carboxylic groups content in a similar way as shown from the motional behaviour of spin label on the segmental level. The degree of crystallinity and the separation of spherulitic rings are decreasing above a certain concentration of functional groups. The effect of functional groups in PU on the hard and soft segment mixing is discussed in terms of additional noncovalent interactions and chain structure which at critical level of interactions lead to a formation of more open hard segment structure accessible to interaction with the soft segment.     

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.     

Conformationally Restricted Isoindoline‐Derived Spin Labels in Duplex DNA: Distances and Rotational Flexibility by Pulsed Electron–Electron Double Resonance Spectroscopy

Three structurally related isoindoline‐derived spin labels that have different mobilities were incorporated into duplex DNA to systematically study the effect of motion on orientation‐dependent pulsed electron–electron double resonance (PELDOR) measurements. To that end, a new nitroxide spin label, ^ExIm U , was synthesized and incorporated into DNA oligonucleotides. ^ExIm U is the first example of a conformationally unambiguous spin label for nucleic acids, in which the nitroxide N?O bond lies on the same axis as the three single bonds used to attach the otherwise rigid isoindoline‐based spin label to a uridine base. Continuous‐wave (CW) EPR measurements of ^ExIm U confirm a very high rotational mobility of the spin label in duplex DNA relative to the structurally related spin label ^Im U , which has restricted mobility due to an intramolecular hydrogen bond. The X‐band CW‐EPR spectra of ^ExIm U can be used to identify mismatches in duplex DNA. PELDOR distance measurements between pairs of the spin labels ^Im U , ^Ox U , and ^ExIm U in duplex DNA showed a strong angular dependence for ^Im U , a medium dependence for ^Ox U , and no orientation effect for ^ExIm U . Thus, precise distances can be extracted from ^ExIm U without having to take orientational effects into account.     

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.     

Triplet Fullerenes as Prospective Spin Labels for Nanoscale Distance Measurements by Pulsed Dipolar EPR Spectroscopy

Precise nanoscale distance measurements by pulsed electron paramagnetic resonance (EPR) spectroscopy play a crucial role in structural studies of biomolecules. The properties of the spin labels used in this approach determine the sensitivity limits, attainable distances, and proximity to biological conditions. Herein, we propose and validate the use of photoexcited fullerenes as spin labels for pulsed dipolar (PD) EPR distance measurements. Hyperpolarization and the narrower spectrum of fullerenes compared to other triplets (e.g., porphyrins) boost the sensitivity, and superior relaxation properties allow PD EPR measurements up to a near‐room temperature. This approach is demonstrated using fullerene–nitroxide and fullerene–triarylmethyl pairs, as well as a supramolecular complex of fullerene with nitroxide‐labeled protein. Photoexcited triplet fullerenes can be considered as new spin labels with outstanding spectroscopic properties for future structural studies of biomolecules.     

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.     

A Spirocyclohexyl Nitroxide Amino Acid Spin Label for Pulsed EPR Spectroscopy Distance Measurements

Site‐directed spin labeling and EPR spectroscopy offer accurate, sensitive tools for the characterization of structure and function of macromolecules and their assemblies. A new rigid spin label, spirocyclohexyl nitroxide α‐amino acid and its N‐(9‐fluorenylmethoxycarbonyl) derivative, have been synthesized, which exhibit slow enough spin‐echo dephasing to permit accurate distance measurements by pulsed EPR spectroscopy at temperatures up to 125 K in 1:1 water/glycerol and at higher temperatures in matrices with higher glass transition temperatures. Distance measurements in the liquid nitrogen temperature range are less expensive than those that require liquid helium, which will greatly facilitate applications of pulsed EPR spectroscopy to the study of structure and conformation of peptides and proteins.     

Low operational stability of enzymes in dry organic solvents: changes in the active site might affect catalysis

The potential of enzyme catalysis in organic solvents for synthetic applications has been overshadowed by the fact that their catalytic properties are affected by organic solvents. In addition, it has recently been shown that an enzyme's initial activity diminishes considerably after prolonged exposure to organic media. Studies geared towards understanding this last drawback have yielded unclear results. In the present work we decided to use electron paramagnetic resonance spectroscopy (EPR) to study the motion of an active site spin label (a nitroxide free radical) during 96 h of exposure of the serine protease subtilisin Carlsberg to four different organic solvents. Our EPR data shows a typical two component spectra that was quantified by the ratio of the anisotropic and isotropic signals. The isotropic component, associated with a mobile nitroxide free radical, increases during prolonged exposure to all solvents used in the study. The maximum increase (of 43%) was observed in 1,4-dioxane. Based on these and previous studies we suggest that prolonged exposure of the enzyme to these solvents provokes a cascade of events that could induce substrates to adopt different binding conformations. This is the first EPR study of the motion of an active-site spin label during prolonged exposure of an enzyme to organic solvents ever reported.     

EPR characterization of cellulose triacetate fibers used for enzyme immobilization

EPR studies of a nitroxide spin label and of the nitroxide spin-labeled albumin entrapped in cellulose triacetate fibers were carried out. The EPR spectra have shown that within the fiber only two phases are present: a liquid one of medium viscosity trapped inside microcavities, and a polymeric one surrounding them. After entrapment, spin-labeled albumin is distributed mainly in the liquid phase, though a not negligible amount of it remains within the polymeric matrix. The EPR studies have shown that, after the standard procedure of drying, the albumin is almost completely precipitated, but about 85% of it returns to solution when the fiber is again placed in the solution. The behavior of the albumin dissolved inside the microcavities toward denaturating agents and pH change, and that of the free albumin in solution is similar; the minor differences noticed indicate a second-order interaction between the fiber and the protein.     

Molecular latent space simulators

Small integration time steps limit molecular dynamics (MD) simulations to millisecond time scales. Markov state models (MSMs) and equation-free approaches learn low-dimensional kinetic models from MD simulation data by performing configurational or dynamical coarse-graining of the state space. The learned kinetic models enable the efficient generation of dynamical trajectories over vastly longer time scales than are accessible by MD, but the discretization of configurational space and/or absence of a means to reconstruct molecular configurations precludes the generation of continuous atomistic molecular trajectories. We propose latent space simulators (LSS) to learn kinetic models for continuous atomistic simulation trajectories by training three deep learning networks to (i) learn the slow collective variables of the molecular system, (ii) propagate the system dynamics within this slow latent space, and (iii) generatively reconstruct molecular configurations. We demonstrate the approach in an application to Trp-cage miniprotein to produce novel ultra-long synthetic folding trajectories that accurately reproduce atomistic molecular structure, thermodynamics, and kinetics at six orders of magnitude lower cost than MD. The dramatically lower cost of trajectory generation enables greatly improved sampling and greatly reduced statistical uncertainties in estimated thermodynamic averages and kinetic rates.

Latent space simulators learn kinetic models for atomistic simulations and generate novel trajectories at six orders of magnitude lower cost.     

Structural dynamics of supercooled water from quasielastic neutron scattering and molecular simulations

One of the outstanding challenges presented by liquid water is to understand how molecules can move on a picosecond time scale despite being incorporated in a three-dimensional network of relatively strong H-bonds. This challenge is exacerbated in the supercooled state, where the dramatic slowing down of structural dynamics is reminiscent of the, equally poorly understood, generic behavior of liquids near the glass transition temperature. By probing single-molecule dynamics on a wide range of time and length scales, quasielastic neutron scattering (QENS) can potentially reveal the mechanistic details of water's structural dynamics, but because of interpretational ambiguities this potential has not been fully realized. To resolve these issues, we present here an extensive set of high-quality QENS data from water in the range 253-293 K and a corresponding set of molecular dynamics (MD) simulations to facilitate and validate the interpretation. Using a model-free approach, we analyze the QENS data in terms of two motional components. Based on the dynamical clustering observed in MD trajectories, we identify these components with two distinct types of structural dynamics: picosecond local (L) structural fluctuations within dynamical basins and slower interbasin jumps (J). The Q-dependence of the dominant QENS component, associated with J dynamics, can be quantitatively rationalized with a continuous-time random walk (CTRW) model with an apparent jump length that depends on low-order moments of the jump length and waiting time distributions. Using a simple coarse-graining algorithm to quantitatively identify dynamical basins, we map the newtonian MD trajectory on a CTRW trajectory, from which the jump length and waiting time distributions are computed. The jump length distribution is gaussian and the rms jump length increases from 1.5 to 1.9 A? as the temperature increases from 253 to 293 K. The rms basin radius increases from 0.71 to 0.75 A? over the same range. The waiting time distribution is exponential at all investigated temperatures, ruling out significant dynamical heterogeneity. However, a simulation at 238 K reveals a small but significant dynamical heterogeneity. The macroscopic diffusion coefficient deduced from the QENS data agrees quantitatively with NMR and tracer results. We compare our QENS analysis with existing approaches, arguing that the apparent dynamical heterogeneity implied by stretched exponential fitting functions results from the failure to distinguish intrabasin (L) from interbasin (J) structural dynamics. We propose that the apparent dynamical singularity at ～220 K corresponds to freezing out of J dynamics, while the calorimetric glass transition corresponds to freezing out of L dynamics.     

EPR analysis and DFT computations of a series of polynitroxides

Polynitroxides with varying numbers of nitroxide groups (one to four) derived from different aromatic core structures show intramolecular electron spin-spin coupling. The scope of this study is to establish an easy methodology for extracting structural, dynamical, and thermodynamical information from the EPR spectra of these polynitroxides which might find use as spin probes in complex systems, such as biological and host/guest systems, and as polarizing agents in dynamic nuclear polarization (DNP) applications. Density functional theory (DFT) calculations at the B3LYP/6-31G(d) level provided information on the structural details such as bond lengths and angles in the gas phase, which were compared with the single crystal X-ray diffraction data in the solid state. Polarizable continuum model (PCM) calculations were performed to account for solvent influences. The electron paramagnetic resonance (EPR) spectra of the polynitroxides in chloroform were analyzed in detail to extract information such as the percentages of different conformers, hyperfine coupling constants a, and rotational correlation times τ(c). The temperature dependence on the line shape of the EPR spectra gave thermodynamic parameters ΔH and ΔS for the conformational transitions. These parameters were found to depend on the number and relative positions of the nitroxide and other polar groups.     

自旋标记法研究STVPh/PEO共混体系的链运动

利用自旋标记法将氮氧自由基连接在聚氧化乙烯分子末端,测得其与不同羟基含量的苯乙烯-4-乙烯基苯酚共聚物组成高分子共混体系的ESR波谱,研究了各组分的分子链运动.自旋标记聚氧化乙烯的ESR谱图在整个温度范围内只显示快运动或慢运动的一种运动成分,表明氮氧自由基处于单一的环境中.高分子共混体系的ESR谱在一定温度范围内同时存在快运动和慢运动的两种运动成分,且两成分相对比例随温度的变化而变化,表明体系中氮氧自由基处于不同的微相环境中.由ESR谱得到的T_5mT,τ_c,E_a的值都随着体系中羟基含量的增加而变大,显示标记分子的链运动在共混体系中的活动逐步受阻.与此同时,随着该共混物中酚羟基和醚氧基间氢键相互作用强度的增大,共混物的相容性得到逐步改善.     

Restricted orientational motion of nitroxides in molecular glasses: direct estimation of the motional time scale basing on the comparative study of primary and stimulated electron spin echo decays

A comparative study of anisotropic relaxation in two-pulse primary and three-pulse stimulated electron spin echo decays provides a direct way to distinguish fast (correlation time tau(c)<10(-6) s) and slow (tau(c)>10(-6) s) motions. Anisotropic relaxation is detected as a difference of the decay rates for different resonance field positions in anisotropic electron paramagnetic resonance spectra. For fast motion anisotropic relaxation influences the primary echo decay and does not influence the stimulated echo decay. For slow motion it is seen in both two-pulse echo and three-pulse stimulated echo decays. For nitroxide spin probes dissolved in glassy glycerol only fast motion was found below 200 K. Increase of temperature above 200 K results in the appearance of slow motion. Its amplitude increases rapidly with temperature increase. While in glycerol glass slow motion appears above glass transition temperature T(g), in ethanol glass it is observable below T(g). The scenario of motional dynamics in glasses is proposed which involves the broadening of the correlation time distribution with increasing temperature.