A structure-based simulation approach for electron paramagnetic resonance spectra using molecular and stochastic dynamics simulations

Electron paramagnetic resonance (EPR) spectroscopy using site-directed spin-labeling is an appropriate technique to analyze the structure and dynamics of flexible protein regions as well as protein-protein interactions under native conditions. The analysis of a set of protein mutants with consecutive spin-label positions leads to the identification of secondary and tertiary structure elements. In the first place, continuous-wave EPR spectra reflect the motional freedom of the spin-label specifically linked to a desired site within the protein. EPR spectra calculations based on molecular dynamics (MD) and stochastic dynamics simulations facilitate verification or refinement of predicted computer-aided models of local protein conformations. The presented spectra simulation algorithm implies a specialized in vacuo MD simulation at 600 K with additional restrictions to sample the entire accessible space of the bound spin-label without large temporal effort. It is shown that the distribution of spin-label orientations obtained from such MD simulations at 600 K agrees well with the extrapolated motion behavior during a long timescale MD at 300 K with explicit water. The following potential-dependent stochastic dynamics simulation combines the MD data about the site-specific orientation probabilities of the spin-label with a realistic rotational diffusion coefficient yielding a set of trajectories, each more than 700 ns long, essential to calculate the EPR spectrum. Analyses of a structural model of the loop between helices E and F of bacteriorhodopsin are illustrated to demonstrate the applicability and potentials of the reported simulation approach. Furthermore, effects on the motional freedom of bound spin-labels induced by solubilization of bacteriorhodopsin with Triton X-100 are examined.     

High-field EPR studies of the structure and conformational changes of site-directed spin labeled bacteriorhodopsin

Cw and pulsed high-field EPR (95 GHz, 3.4 T) are performed on site-directed spin labeled bacteriorhodopsin (BR) mutants. The enhanced Zeeman splitting leads to spectra with resolved g-tensor components of the nitroxide spin label. The g(xx) component shift determined for 10 spin labels located in the cytoplasmic loop region and in the protein interior along the BR proton channel reveals a maximum close to position 46 between the proton donor D96 and the retinal. A plot of g(xx) versus A(zz) of the nitrogen discloses grouping of 12 spin labeled sites in protic and aprotic sites. Spin labels at positions 46, 167 and 171 show the aprotic character of the cytoplasmic moiety of the proton channel whereas nitroxides at positions 53, 194 and 129 reveal the protic environment in the extracellular channel. The enhanced sensitivity of high-field EPR with respect to anisotropic reorientational motion of nitroxides allows the characterization of different motional modes for spin labels bound to positions 167 and 170. The motional restriction of the nitroxide at position 167 of the double mutant V167C/D96N is decreased in the M(N) photo-intermediate. An outward shift of the cytoplasmic moiety of helix F in the M(N) intermediate would account for the high-field EPR results and is in agreement with diffraction and recent X-band EPR data.     

Librational motion of an "immobilized" spin label: hemoglobin spin labeled by a maleimide derivative.

The spin label Tempo-maleimide, when "immobilized" in hemoglobin, is shown to exhibit motional fluctuation whose amplitude and/or frequency depend on temperature and solution conditions. These motional fluctuations are observable by several electron spin resonance techniques. For desalted hemoglobin the fluctuations are detectable at approximately -15 degrees C using saturation transfer techniques and at approximately +25 degrees C using line-width measurements of normal absorption spectra. In ammonium sulfate precipitated hemoglobin, however, motional fluctuations are not detectable by either technique up to at least 40 degrees C. The most probable mechanism for spin-label motion appears to be either fluctuations in protein conformation which affect the label binding site or conformational transitions of the nitroxide ring itself. These motional fluctuations are shown to introduce a librational character to the overall label motion during hemoglobin rotational diffusion, with the librational motion significantly affecting the use of spin-label spectral shapes to calculate hemoglobin rotational correlation times.     

Saturation-transfer electron paramagnetic resonance detection of anisotropic motion by sickle hemoglobin molecules in the polymer state

The motional behavior of spin-labeled deoxygenated sickle hemoglobin has been studied by using both 9- and 35-GHz saturation-transfer electron paramagnetic resonance (EPR). Using spectral subtraction techniques and saturation-transfer EPR parameter correlation plots, we find that the saturation-transfer EPR spectra for the sickle hemoglobin gel state at high temperature and high hemoglobin concentration cannot be described as a simple superposition of spectra from immobilized hemoglobin plus solution-state hemoglobin but instead suggest that the individual sickle hemoglobin molecules exhibit limited, anisotropic, rotational oscillation within the polymer fiber. The spectra also imply that the symmetry axis for sickle hemoglobin rotational oscillation is approximately coincident with the nitroxide z axis of the covalently attached spin-label. We suggest that this anisotropic rotational motion may be produced by one or two of the known intermolecular contact sites within the sickle hemoglobin fiber acting as strong intermolecular binding sites, and producing "motional alignment" within the fiber; determining the location of the strong binding site should be important in focusing the future development of antisickling agents.     

Apolipoprotein A-I assumes a "looped belt" conformation on reconstituted high density lipoprotein

Apolipoprotein A-I (apoA-I) plays a central role in the reverse cholesterol transport pathway; however, the structural basis for its antiatherogenic effects remains poorly understood. Here we employ EPR spectroscopy and fluorescence resonance energy transfer to elucidate the conformation and relative alignment of apoA-I monomers on discoidal (9.4 nm) reconstituted high density lipoprotein (rHDL). EPR spectroscopy provided evidence for an extended helical secondary structure. Position 139 since it was the only residue examined to display a dynamic motional character consistent with a flexible loop structure. The EPR spectra of nitroxide probes at positions 133 and 146 exhibit spin coupling, indicating that these positions are proximal to an apoA-I paired counterpart on the perimeter of rHDL. fluorescence resonance energy transfer studies employing engineered apoA-I variants possessing a single tryptophan (energy donor) and/or a single cysteine (whose thiol moiety was covalently labeled with an extrinsic energy acceptor) provided evidence that paired apoA-I molecules around the perimeter of rHDL align in an extended antiparallel conformation. Taken together with the observation that the EPR spectra of nitroxide probes positioned at intervening sequence positions (134-145) do not exhibit spin coupling, this has led us to propose a "looped belt" model, wherein residues 133-146 comprise a flexible loop segment that confers to apoA-I an intrinsic ability to adapt its structure to accommodate changing particle lipid content. Specifically, in the looped belt model, with the exception of amino acids 134-145, apoA-I aligns with its counterpart in a helix 5-helix 5 registry, centered at position 139.     

Different conformations of the 3' termini of initiator and elongator transfer ribonucleic acids. An EPR study

The 3' ends of transfer ribonucleic acids were covalently labeled with a nitroxide spin label. The 3' end of initiator tRNA (tRNAMetf) from Escherichia coli shows different motional behavior than the 3' terminus of elongator tRNAs as monitored by EPR. The line shapes of the EPR spectra are quite sensitive to the buffer conditions, as shown by measurements in 4 different buffers. The data are consistent with a constrained or folded back 3' terminus in the initiator tRNA as opposed to the freely rotating elongator 3' terminus. The EPR spectra are also sensitive to aggregation of the tRNA.     

Progress towards a DNA-based model system for the study of electron spin-spin interactions

A DNA-based model system is described for studying electron spin-spin interactions between a paramagnetic metal ion and a nitroxide spin label. The modified base deoxythymidine-EDTA (dT-EDTA) chelates the divalent or trivalent metal ion and produces a new feature in the circular dichroism (CD) spectra that makes it possible to monitor local DNA melting. Based on the results of optical and electron paramagnetic resonance (EPR) experiments, we find that the terminus of the DNA duplex that incorporates dT-EDTA and the spin-label melts at a higher temperature than the rest of the DNA duplex. EPR microwave progressive power saturation experiments performed at 77 K are consistent with the specific binding of Dy(III) at the EDTA site and an intramolecular dipole-dipole interaction between the nitroxide spin-label and the chelated Dy(III). This model system should be suitable for studying the relaxation properties of metal ions by saturation-recovery EPR.     

Two-dimensional 1H NMR of three spin-labeled derivatives of bovine pancreatic trypsin inhibitor

Three nitroxide spin-labeled monoderivatives of bovine pancreatic trypsin inhibitor were prepared with the amino-specific reagent succinimidyl 1-oxy-2,2,5,5-tetramethyl-3-pyrroline-3-carboxylate. The monoderivatives were purified by ion-exchange and affinity chromatography. Thin-layer maps of tryptic peptides of the monoderivatives showed that the spin-label was incorporated at either the alpha-amino group, Lys-15, or Lys-26. Two-dimensional J-correlated 1H NMR spectra of the monoderivatives were recorded. Spectra were also recorded after reduction by ascorbic acid of the nitroxide label to hydroxylamine. With the nitroxide label present, significant line-broadening effects on many of the cross peaks in the spectra were observed. The extent of line broadening for the C alpha H-NH cross peaks was qualitatively correlated with the distance between the labeled amino group and the average C alpha H-NH position in the crystal structure. The spin-label affects cross peaks of protons within approximately 15 A. This study suggests that it is feasible to accumulate sufficient intramolecular distances in order to determine protein solution structures with the aid of distance geometry algorithms.     

Dynamics and interactions of the anion channel in intact human erythrocytes: an electron paramagnetic resonance spectroscopic study employing a new membrane-impermeant bifunctional spin-label

We have developed a new membrane-impermeant, bifunctional spin-labeling reagent, bis-(sulfo-N-succinimidyl) doxyl-2-spiro-4'-pimelate (BSSDP), and employed it in an electron paramagnetic resonance (EPR) study of the rotational diffusion of the anion-exchange channel (band 3) in intact human erythrocytes. BSSDP reacts in a covalent manner and with high specificity with the extracytoplasmic domain of band 3, forming a complex in which the spin-label is immobilized on the protein. The linear EPR spectrum of BSSDP-labeled intact erythrocytes is characteristic of a highly immobilized, spatially isolated nitroxide probe. The saturation-transfer EPR spectrum of the same sample indicates that the anion channel in intact erythrocytes exhibits rotational dynamics in the 0.1-1 ms correlation time range at 20 degrees C. Rotational dynamics in this motional domain are consistent with a strong interaction of the anion-exchange channel with the erythrocyte cytoskeleton. The saturation-transfer EPR spectrum of ghosts prepared from BSSDP-labeled erythrocytes indicates a significant increase in rotational mobility of the anion channel, suggesting a significant disruption on lysis of interactions between the anion channel and the cytoskeleton.     

Electron spin resonance study of human alpha 1-acid glycoprotein interaction with a spin labelled steroid.

The interaction of human alpha 1-acid glycoprotein (AAG) with a corticosteroid was studied using nitroxide labeled deoxycorticosterone and electron spin resonance (ESR) spectroscopy. The ESR spectra of the spin labeled steroid in the presence of AAG could be used to characterize the ligand-protein interaction at equilibrium without the need of a separation between bound and free species. An association constant Ka of 6.10(5) M-1 at 20 degrees C and a binding capacity of one site per mole protein were found. ESR spectra recorded at equilibrium at various temperatures allowed the calculation of enthalpy and entropy variations for the steroid-protein interaction; these thermodynamic parameters exhibited a rapid change above 45 degrees C which may be related to a protein conformational modification above this temperature, as detected by circular dichroism study. The ESR spectra width could be used to define a polar character for the spin label environment in the steroid binding site of AAG and to calculate an apparent rotational correlation time of 2.8 x 10(-8) sec for the steroid-protein complex in aqueous solution at 20 degrees C. It can be concluded that spin labeling and ESR methodology is of value in the study of steroid-protein interactions of biological significance above all because it can provide direct physico-chemical information concerning the local environment of the ligand in its binding site at equilibrium.     

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.     

Effect of phospholipid substitution on the mobility of spin labels bound to the ATPase of sarcoplasmic reticulum

The mobility of spin labels covalently bound to the Ca2+-transport ATPase (ATP phosphohydrolase [EC 3.y.1.3]) was studied by electron spin-resonance spectroscopy in purified ATPase and reconstituted vesicles. The purified ATPase of sarcoplasmic reticulum of rabbit skeletal muscle was covalently labeled with maleimide spin-labels of different chain length and the phospholipids were exchanged for dipalmitoylphosphatidylcholine. The spectrum of the short-chain maleimide spin-label, bound to purified ATPase indicates reduced mobility after substitution of endogenous phospholipids with dipalmitoylphosphatidylcholine. With the long-chain maleimide derivative no difference was detected in the spectra, measured at 20-35 degrees C temperature before and after substitution with dipalmitoylphosphatidylcholine. Below 10 degrees C temperature the substitution with dipalmitoylphosphatidylcholine decreased the mobility of the prove, indicating that the microviscosity of environment in the vicinity of nitroxide groups was influenced by changes in the fatty acid composition. With both short and long chain spin-labels bound to purified ATPase adn sarcoplasmic reticulum vesicles the amplitude of weakly immobilized component sharply decreased in media containing 20-50% glycerol. Therefore, the mobility of covalently bound nitroxide group in short or long chain maleimide derivatives is also sensitive to the viscosity of the water phase.     

Preparation and characterization of spin-labeled derivatives of epidermal growth factor (EGF) for investigations of the interactions of EGF with its receptor by electron paramagnetic resonance spectroscopy.

We prepared, purified, and characterized derivatives of epidermal growth factor (EGF) having a nitroxide spin-label attached covalently at the amino terminus. Characterization of these derivatives with regard to the positions of attachment of the spin-label was accomplished by a combination of peptide mapping, protein sequencing, and fast atom bombardment-mass spectrometry. One derivative was chosen for use in initial investigations by electron paramagnetic resonance (EPR) spectroscopy of receptor-bound EGF and its dissociation kinetics. This derivative was found to be equipotent with the native hormone in competitive binding assays, in activating the EGF receptor kinase, and in stimulating the formation of EGF receptor dimers in solubilized cell extracts. Upon binding to solubilized EGF receptor, the spin-labeled EGF derivative became immobilized, giving rise to a visually distinct slow-motion EPR spectrum. The resulting spectrum showed no detectable dipolar interaction between nitroxides, indicating that the nitroxide moieties of spin-labels reacted at the amino termini of receptor-bound spin-labeled EGF molecules are separated by a distance of at least 16 A. An EPR study of the kinetics of dissociation of spin-labeled EGF in the presence of excess unlabeled EGF revealed a rapid component with a k off approximately 2 x 10(-2) s-1 and a less well resolved slow component.     

Integrated analysis of the conformation of a protein-linked spin label by crystallography,EPR and NMR spectroscopy

Long-range structural information derived from paramagnetic relaxation enhancement observed in the presence of a paramagnetic nitroxide radical is highly useful for structural characterization of globular, modular and intrinsically disordered proteins, as well as protein–protein and protein-DNA complexes. Here we characterized the conformation of a spin-label attached to the homodimeric protein CylR2 using a combination of X-ray crystallography, electron paramagnetic resonance (EPR) and NMR spectroscopy. Close agreement was found between the conformation of the spin label observed in the crystal structure with interspin distances measured by EPR and signal broadening in NMR spectra, suggesting that the conformation seen in the crystal structure is also preferred in solution. In contrast, conformations of the spin label observed in crystal structures of T4 lysozyme are not in agreement with the paramagnetic relaxation enhancement observed for spin-labeled CylR2 in solution. Our data demonstrate that accurate positioning of the paramagnetic center is essential for high-resolution structure determination.     

Electron paramagnetic resonance probing of macromolecules: a comparison of structure-function relationships in chymotrypsinogen, -chymotrypsin and anhydrochymotrypsin

Chymotrypsinogen, chymotrypsin and anhydrochymotrypsin have been covalently spin-labeled by an analog of bromoacetamide, and the latter two proteins have been labeled by an analog of 1-chloro-3-tosylamido-4-phenyl butanone. The electron paramagnetic resonance spectra of the labeled proteins indicate protein conformational changes accompanying (1) activation of the zymogen and (2) the binding of protons and substrates by the native and anhydro enzymes, and tertiary structural differences between these protein forms which are at once informative and predictable. A spin-label linked to the thioether side-chain of methionine 192 in Chymotrypsinogen may be in contact with a hydrophobic surface. This interaction is lost upon zymogen activation with little change in the isotropic rotational freedom of the nitroxide group. The rotational freedom of the group increases sigmoidally with pH; a spectral dependence upon an ionizing group (pK_a = 8.9) is demonstrated. The binding of indole to the labeled enzyme raises the pK_a of the ionizing group to 10.2. A spin-label linked to histidine 57 in chymotrypsin senses both indole binding and pH changes directly; the same label in anhydrochymotrypsin responds directly only to changes in pH. Neither histidine-labeled derivative exhibits enzymic activity. The electron paramagnetic resonance spectra of these two labeled proteins at high pH indicate a decrease in the motional freedom of the spin label. The spectral data show that the conformational state of the labeled zymogen is not similar to the high-pH conformational state of the labeled enzyme. Furthermore, the pH-dependent conformational transition of labeled chymotrypsin requires neither the serine 195 hydroxyl nor the histidine 57 imidazole, since the transition occurs normally in derivatized and chemically modified protein forms. The chemical reactivity of histidine 57 in anhydrochymotrypsin is evaluated and the catalytic activities of two histidine alkylated enzymes are compared.     

Mapping local protein electrostatics by EPR of pH-sensitive thiol-specific nitroxide

A first thiol-specific pH-sensitive nitroxide spin-label of the imidazolidine series, methanethiosulfonic acid S-(1-oxyl-2,2,3,5,5-pentamethylimidazolidin-4-ylmethyl) ester (IMTSL), has been synthesized and characterized. X-Band (9 GHz) and W-band (94 GHz) EPR spectral parameters of the new spin-label in its free form and covalently attached to an amino acid cysteine and a tripeptide glutathione were studied as a function of pH and solvent polarity. The pKa value of the protonatable tertiary amino group of the spin-label was found to be unaffected by other ionizable groups present in side chains of unstructured small peptides. The W-band EPR spectra were shown to allow for pKa determination from precise g-factor measurements. Is has been demonstrated that the high accuracy of pKa determination for pH-sensitive nitroxides could be achieved regardless of the frequency of measurements or the regime of spin exchange: fast at X-band and slow at W-band. IMTSL was found to react specifically with a model protein, iso-1-cytochrome c from the yeast Saccharomyces cerevisiae, giving EPR spectra very similar to those of the most commonly employed cysteine-specific label MTSL. CD data indicated no perturbations to the overall protein structure upon IMTSL labeling. It was found that for IMTSL, g iso correlates linearly with A iso, but the slopes are different for the neutral and charged forms of the nitroxide. This finding was attributed to the solvent effects on the spin density at the oxygen atom of the NO group and on the excitation energy of the oxygen lone-pair orbital.     

Selective detection of the rotational dynamics of the protein-associated lipid hydrocarbon chains in sarcoplasmic reticulum membranes.

We have developed a saturation transfer EPR (ST-EPR) method to measure selectively the rotational dynamics of those lipids that are motionally restricted by integral membrane proteins and have applied this methodology to measure lipid-protein interactions in native sarcoplasmic reticulum (SR) membranes. This analysis involves the measurement of spectral saturation using a series of six stearic acid spin labels that are labeled with a nitroxide at different carbon atom positions. A large amount of spectral saturation is observed for spin labels in native SR membranes, but not for spin labels in dispersions of extracted SR lipids, implying that the motional properties of those lipids interacting with the Ca-ATPase, i.e., the boundary or annular lipid, can be directly measured without the need for spectral subtraction procedures. A comparison of the motional properties of the restricted lipid, measured by ST-EPR, with those measured by digital subtraction of conventional EPR spectra qualitatively agree, for in both cases the Ca-ATPase restricts the rotational mobility of a population of lipids, whose rotational mobility increases as the nitroxide is positioned toward the center of the bilayer. However, the ability of ST-EPR to directly measure the motionally restricted lipid in a model-independent means provides the greater precision necessary to measure small changes in the rotational dynamics of the lipid at the protein-lipid interface, providing a valuable tool in clarifying the relationship between the physical nature of the protein-lipid interface and membrane function.