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.     

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.     

A Reactive,Rigid GdIII Labeling Tag for In‐Cell EPR Distance Measurements in Proteins

The cellular environment of proteins differs considerably from in vitro conditions under which most studies of protein structures are carried out. Therefore, there is a growing interest in determining dynamics and structures of proteins in the cell. A key factor for in‐cell distance measurements by the double electron–electron resonance (DEER) method in proteins is the nature of the used spin label. Here we present a newly designed Gd^III spin label, a thiol‐specific DOTA‐derivative (DO3MA‐3BrPy), which features chemical stability and kinetic inertness, high efficiency in protein labelling, a short rigid tether, as well as favorable spectroscopic properties, all are particularly suitable for in‐cell distance measurements by the DEER method carried out at W‐band frequencies. The high performance of DO3MA‐3BrPy‐Gd^III is demonstrated on doubly labelled ubiquitin D39C/E64C, both in vitro and in HeLa cells. High‐quality DEER data could be obtained in HeLa cells up to 12 h after protein delivery at in‐cell protein concentrations as low as 5–10 μm .     

A Reactive,Rigid GdIII Labeling Tag for In‐Cell EPR Distance Measurements in Proteins

The cellular environment of proteins differs considerably from in vitro conditions under which most studies of protein structures are carried out. Therefore, there is a growing interest in determining dynamics and structures of proteins in the cell. A key factor for in‐cell distance measurements by the double electron–electron resonance (DEER) method in proteins is the nature of the used spin label. Here we present a newly designed Gd^III spin label, a thiol‐specific DOTA‐derivative (DO3MA‐3BrPy), which features chemical stability and kinetic inertness, high efficiency in protein labelling, a short rigid tether, as well as favorable spectroscopic properties, all are particularly suitable for in‐cell distance measurements by the DEER method carried out at W‐band frequencies. The high performance of DO3MA‐3BrPy‐Gd^III is demonstrated on doubly labelled ubiquitin D39C/E64C, both in vitro and in HeLa cells. High‐quality DEER data could be obtained in HeLa cells up to 12 h after protein delivery at in‐cell protein concentrations as low as 5–10 μm .     

In‐Situ Spin Labeling of His‐Tagged Proteins: Distance Measurements under In‐Cell Conditions

New spin labeling strategies have immense potential in studying protein structure and dynamics under physiological conditions with electron paramagnetic resonance (EPR) spectroscopy. Here, a new spin‐labeled chemical recognition unit for switchable and concomitantly high affinity binding to His‐tagged proteins was synthesized. In combination with an orthogonal site‐directed spin label, this novel spin probe, Proxyl‐trisNTA (P‐trisNTA) allows the extraction of structural constraints within proteins and macromolecular complexes by EPR. By using the multisubunit maltose import system of E. coli: 1) the topology of the substrate‐binding protein, 2) its substrate‐dependent conformational change, and 3) the formation of the membrane multiprotein complex can be extracted. Notably, the same distance information was retrieved both in vitro and in situ allowing for site‐specific spin labeling in cell lysates under in‐cell conditions. This approach will open new avenues towards in‐cell EPR.     

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.     

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.     

A GRID-derived water network stabilizes molecular dynamics computer simulations of a protease

Structural water molecules are crucial for the stability and function of proteins. Recently, we presented a molecular dynamics (MD) study on blood coagulation factor Xa (fXa) to investigate the effect of water molecules on the flexibility of the protein structure. We showed that neglecting important water positions at the outset of the simulation leads to severe structural distortions during the MD simulations: A stable trajectory was obtained with a water set that was derived from all 73 X-ray structures of the protein. However, for many proteins of interest, only limited structural data is available, which precludes the merging of information from many X-ray structures. Here, we show that an in silico assembled water network, derived from molecular interaction fields generated with the GRID program, is a viable alternative to X-ray data. MD simulations with the GRID water set show a significantly improved stability over alternative setups without water or the X-ray resolved water molecules in the starting structure. The performance is comparable to a water setup derived from a recently presented clustering approach.     

Ox-SLIM: Synthesis of and Site-Specific Labelling with a Highly Hydrophilic Trityl Spin Label

The combination of pulsed dipolar electron paramagnetic resonance spectroscopy (PDS) with site-directed spin labelling is a powerful tool in structural biology. Rational design of trityl-based spin labels has enabled studying biomolecular structures at room temperature and within cells. However, most current trityl spin labels suffer either from aggregation with proteins due to their hydrophobicity, or from bioconjugation groups not suitable for in-cell measurements. Therefore, we introduce here the highly hydrophilic trityl spin label Ox-SLIM. Engineered as a short-linked maleimide, it combines the most recent developments in one single molecule, as it does not aggregate with proteins, exhibits high resistance under in-cell conditions, provides a short linker, and allows for selective and efficient spin labelling via cysteines. Beyond establishing synthetic access to Ox-SLIM, its suitability as a spin label is illustrated and ultimately, highly sensitive PDS measurements are presented down to protein concentrations as low as 45 nm resolving interspin distances of up to 5.5 nm.     

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

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

A Bioresistant Nitroxide Spin Label for In‐Cell EPR Spectroscopy: In Vitro and In Oocytes Protein Structural Dynamics Studies

Approaching protein structural dynamics and protein–protein interactions in the cellular environment is a fundamental challenge. Owing to its absolute sensitivity and to its selectivity to paramagnetic species, site‐directed spin labeling (SDSL) combined with electron paramagnetic resonance (EPR) has the potential to evolve into an efficient method to follow conformational changes in proteins directly inside cells. Until now, the use of nitroxide‐based spin labels for in‐cell studies has represented a major hurdle because of their short persistence in the cellular context. The design and synthesis of the first maleimido‐proxyl‐based spin label (M‐TETPO) resistant towards reduction and being efficient to probe protein dynamics by continuous wave and pulsed EPR is presented. In particular, the extended lifetime of M‐TETPO enabled the study of structural features of a chaperone in the absence and presence of its binding partner at endogenous concentration directly inside cells.     

Long Distance Measurements up to 160 Å in the GroEL Tetradecamer Using Q‐Band DEER EPR Spectroscopy

Current distance measurements between spin‐labels on multimeric protonated proteins using double electron–electron resonance (DEER) EPR spectroscopy are generally limited to the 15–60 Å range. Here we show how DEER experiments can be extended to dipolar evolution times of ca. 80 μs, permitting distances up to 170 Å to be accessed in multimeric proteins. The method relies on sparse spin‐labeling, supplemented by deuteration of protein and solvent, to minimize the deleterious impact of multispin effects and substantially increase the apparent spin‐label phase memory relaxation time, complemented by high sensitivity afforded by measurements at Q‐band. We demonstrate the approach using the tetradecameric molecular machine GroEL as an example. Two engineered surface‐exposed mutants, R268C and E315C, are used to measure pairwise distance distributions with mean values ranging from 20 to 100 Å and from 30 to 160 Å, respectively, both within and between the two heptameric rings of GroEL. The measured distance distributions are consistent with the known crystal structure of apo GroEL. The methodology presented here should significantly expand the use of DEER for the structural characterization of conformational changes in higher order oligomers.     

A generalized Langevin dynamics approach to model solvent dynamics effects on proteins via a solvent-accessible surface. The carboxypeptidase A inhibitor protein as a model

A generalized Langevin dynamics (GLD) scheme is derived for (bio)macromolecules having internal structure, arbitrary shapes and a size larger than solvent molecules (i.e. proteins). The concept of solvent-accessible surface area (SASA) is used to incorporate solvent effects via external forces thereby avoiding its explicit molecular representation. A simulation algorithm is implemented in the GROMOS molecular dynamics (MD) program including random forces and memory effects, while solvation effects enter via derivatives of the surface area. The potato carboxypeptidase inhibitor (PCI), a small protein, is used to numerically test the approach. This molecule has N- and C-terminal tails whose structure and fluctuations are solvent dependent. A 1-ns MD trajectory was analyzed in depth. X-ray and NMR structures are used in conjunction with MD simulations with and without explicit solvent to gauge the quality of the results. All the analyses showed that the GLD simulation approached the results obtained for the MD simulation with explicit simple-point-charge-model water molecules. The SASAs of the polar atoms show a natural exposure towards the solvent direction. A FLS solvent simulation was completed in order to sense memory effects. The approach and results presented here could be of great value for developing alternatives to the use of explicit solvent molecules in the MD simulation of proteins, expanding its use and the time-scale explored. Received: 2 February 2000 / Revised: 12 March 2000 / Accepted: 26 May 2000 / Published online: 2 November 2000     

Synthesis of μ2-Oxo-Bridged Iron(III) Tetraphenylporphyrin–Spacer–Nitroxide Dimers and their Structural and Dynamics Characterization by using EPR and MD Simulations

Iron(III) porphyrins have the propensity to form μ₂-oxo-dimers, the structures of which resemble two wheels on an axle. Whereas their crystal structure is known, their solution structure and internal dynamics is not. In the present work, the structure and dynamics of such dimers were studied by means of electron paramagnetic resonance (EPR) spectroscopy and quantum chemistry based molecular dynamics (MD) simulations by using the semiempirical tight-binding method (GFN-xTB). To enable EPR investigation of the dimers, a nitroxide was attached to each of the tetraphenylporphyrin cores through a linear and a bent linker. The inter-nitroxide distance distributions within the dimers were determined by continuous-wave (cw)-EPR and pulsed electron–electron double resonance (PELDOR or DEER) experiments and, with the help of MD, interpreted in terms of the rotation of the porphyrin planes with respect to each other around the Fe–O–Fe axis. It was found that such rotation is restricted to the four registers defined by the phenyl substituents. Within the registers, the rotation angle swings between 30° and 60° in the proximal and between 125° and 145° in the distal register. With EPR, all four angles were found to be equally populated, whereas the 30° and 145° angles are strongly favored to the expense of the 60° and 125° angles in the MD simulation. In either case, the internal dynamics of these dimers thus resemble the motion of a step motor.     

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.     

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.     

Post‐synthetic Spin‐Labeling of RNA through Click Chemistry for PELDOR Measurements

Site‐directed spin labeling of RNA based on click chemistry is used in combination with pulsed electron‐electron double resonance (PELDOR) to benchmark a nitroxide spin label, called here d? . We compare this approach with another established method that employs the rigid spin label Çm for RNA labeling. By using CD spectroscopy, thermal denaturation measurements, CW‐EPR as well as PELDOR we analyzed and compared the influence of d? and Çm on a self‐complementary RNA duplex. Our results demonstrate that the conformational diversity of d? is significantly reduced near the freezing temperature of a phosphate buffer, resulting in strongly orientation‐selective PELDOR time traces of the d? ‐labeled RNA duplex.     

Light‐Induced Porphyrin‐Based Spectroscopic Ruler for Nanometer Distance Measurements

We present a novel pulsed electron paramagnetic resonance (EPR) spectroscopic ruler to test the performance of a recently developed spin‐labeling method based on the photoexcited triplet state (S=1). Four‐pulse electron double resonance (PELDOR) experiments are carried out on a series of helical peptides, labeled at the N‐terminal end with the porphyrin moiety, which can be excited to the triplet state, and with the nitroxide at various sequence positions, spanning distances in the range 1.8–8 nm. The PELDOR traces provide accurate distance measurements for all the ruler series, showing deep envelope modulations at frequencies varying in a progressive way according to the increasing distance between the spin labels. The upper limit is evaluated and found to be around 8 nm. The PELDOR‐derived distances are in excellent agreement with theoretical predictions. We demonstrate that high sensitivity is acquired using the triplet state as a spin label by comparison with Cu(II)–porphyrin analogues. The new labeling approach has a high potential for measuring nanometer distances in more complex biological systems due to the properties of the porphyrin triplet state.     

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.     

Measurement of Angstrom to Nanometer Molecular Distances with 19F Nuclear Spins by EPR/ENDOR Spectroscopy

Spectroscopic and biophysical methods for structural determination at atomic resolution are fundamental in studies of biological function. Here we introduce an approach to measure molecular distances in bio‐macromolecules using ¹⁹F nuclear spins and nitroxide radicals in combination with high‐frequency (94 GHz/3.4 T) electron–nuclear double resonance (ENDOR). The small size and large gyromagnetic ratio of the ¹⁹F label enables to access distances up to about 1.5 nm with an accuracy of 0.1–1 Å. The experiment is not limited by the size of the bio‐macromolecule. Performance is illustrated on synthesized fluorinated model compounds as well as spin‐labelled RNA duplexes. The results demonstrate that our simple but strategic spin‐labelling procedure combined with state‐of‐the‐art spectroscopy accesses a distance range crucial to elucidate active sites of nucleic acids or proteins in the solution state.