Molecular dynamics simulation of site-directed spin labeling: experimental validation in muscle fibers

Role of the H-bond between L53 and T56 for Aquaporin-4 epitope in Neuromyelitis Optica

Aquaporin-4 (AQP4) is the CNS water channel organized into well-ordered protein aggregates called Orthogonal Arrays of Particles (OAPs). Neuromyelitis Optica (NMO) is an autoimmune disease caused by anti-OAP autoantibodies (AQP4-IgG). Molecular Dynamics (MD) simulations have identified an H-bond between L53 and T56 as the key for AQP4 epitope and therefore of potential interest for drug design in NMO field. In the present study, we have experimentally tested this MD-prediction using the classic mutagenesis approach. We substituted T56 with V56 and tested this mutant for AQP4 aggregates and AQP4-IgG binding. gSTED super-resolution microscopy showed that the mutation does not affect AQP4 aggregate dimension; immunofluorescence and cytofluorimetric analysis demonstrated its unaltered AQP4-IgG binding, therefore invalidating the MD-prediction. We later investigated whether AQP4, expressed in Sf9 insect and HEK-293F cells, is able to correctly aggregate before and after the purification steps usually applied to obtain AQP4 crystal. The results demonstrated that AQP4-IgG recognizes AQP4 expressed in Sf9 and HEK-293F cells by immunofluorescence even though BN-PAGE analysis showed that AQP4 forms smaller aggregates when expressed in insect cells compared to mammalian cell lines. Notably, after AQP4 purification, from both insect and HEK-293F cells, no aggregates are detectable by BN-PAGE and AQP4-IgG binding is impaired in sandwich ELISA assays. All together these results indicate that 1) the MD prediction under analysis is not supported by experimental data and 2) the procedure to obtain AQP4 crystals might affect its native architecture and, as a consequence, MD simulations. In conclusion, given the complex nature of the AQP4 epitope, MD might not be the suitable for molecular medicine advances in NMO.

A single-stranded junction modulates nanosecond motional ordering of the substrate recognition duplex of a group I ribozyme

Rigid spinning: Site-directed spin-labeling studies using a rigid nitroxide spin label (Ç) reveal that both length and sequence of a single-stranded junction (J1/2) modulate nanosecond motional ordering of the substrate-recognition duplex (P1) of the 120 kD group I ribozyme. The studies demonstrate an approach for experimental measurements of nanosecond dynamics in high-molecular-weight RNA complexes.

P.Arg82Leu von Hippel-Lindau (VHL) gene mutation among three members of a family with familial bilateral pheochromocytoma in India: molecular analysis and in silico characterization

Various missense mutations in the VHL gene have been reported among patients with familial bilateral pheochromocytoma. However, the p.Arg82Leu mutation in the VHL gene described here among patients with familial bilateral pheochromocytoma, has never been reported previously in a germline configuration. Interestingly, long-term follow-up of these patients indicated that the mutation might have had little impact on the normal function of the VHL gene, since all of them have remained asymptomatic. We further attempted to correlate this information with the results obtained by in silico analysis of this mutation using SIFT, PhD-SNP SVM profile, MutPred, PolyPhen2, and SNPs&GO prediction tools. To gain, new mechanistic insight into the structural effect, we mapped the mutation on to 3D structure (PDB ID 1LM8). Further, we analyzed the structural level changes in time scale level with respect to native and mutant protein complexes by using 12 ns molecular dynamics simulation method. Though these methods predict the mutation to have a pathogenic potential, it remains to be seen if these patients will eventually develop symptomatic disease.

How far in-silico computing meets real experiments. A study on the structure and dynamics of spin labeled vinculin tail protein by molecular dynamics simulations and EPR spectroscopy

Background

Investigation of conformational changes in a protein is a prerequisite to understand its biological function. To explore these conformational changes in proteins we developed a strategy with the combination of molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy. The major goal of this work is to investigate how far computer simulations can meet the experiments.

Methods

Vinculin tail protein is chosen as a model system as conformational changes within the vinculin protein are believed to be important for its biological function at the sites of cell adhesion. MD simulations were performed on vinculin tail protein both in water and in vacuo environments. EPR experimental data is compared with those of the simulated data for corresponding spin label positions.

Results

The calculated EPR spectra from MD simulations trajectories of selected spin labelled positions are comparable to experimental EPR spectra. The results show that the information contained in the spin label mobility provides a powerful means of mapping protein folds and their conformational changes.

Conclusions

The results suggest the localization of dynamic and flexible regions of the vinculin tail protein. This study shows MD simulations can be used as a complementary tool to interpret experimental EPR data.

Molecular dynamics simulations and drug discovery

This review discusses the many roles atomistic computer simulations of macromolecular (for example, protein) receptors and their associated small-molecule ligands can play in drug discovery, including the identification of cryptic or allosteric binding sites, the enhancement of traditional virtual-screening methodologies, and the direct prediction of small-molecule binding energies. The limitations of current simulation methodologies, including the high computational costs and approximations of molecular forces required, are also discussed. With constant improvements in both computer power and algorithm design, the future of computer-aided drug design is promising; molecular dynamics simulations are likely to play an increasingly important role.

EPR spectra and molecular dynamics agree that the nucleotide pocket of myosin V is closed and that it opens on binding actin

We have used EPR spectroscopy and computational modeling of nucleotide-analog spin probes to investigate conformational changes at the nucleotide site of myosin V. We find that, in the absence of actin, the mobility of a spin-labeled diphosphate analog [spin-labeled ADP (SLADP)] bound at the active site is strongly hindered, suggesting a closed nucleotide pocket. The mobility of the analog increases when the MV·SLADP complex (MV=myosin V) binds to actin, implying an opening of the active site in the A·MV·SLADP complex (A=actin). The probe mobilities are similar to those seen with myosin II, despite the fact that myosin V has dramatically altered kinetics. Molecular dynamics (MD) simulation was used to understand the EPR spectra in terms of the X-ray database. The X-ray structure of MV·ADP·BeFx shows a closed nucleotide site and has been proposed to be the detached state. The MV·ADP structure shows an open nucleotide site and has been proposed to be the A·MV·ADP state at the end of the working powerstroke. MD simulation of SLADP docked in the closed conformation gave a probe mobility comparable to that seen in the EPR spectrum of the MV·SLADP complex. The simulation of the open conformation gave a probe mobility that was 35-40° greater than that observed experimentally for the A·MV·SLADP state. Thus, EPR, X-ray diffraction, and computational analysis support the closed conformation as a myosin V state that is detached from actin. The MD results indicate that the MV·ADP crystal structure, which may correspond to the strained actin-bound post-powerstroke conformation resulting from head-head interaction in the dimeric processive motor, is superopened.

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.

Multiple conformations of the nucleotide site of Kinesin family motors in the triphosphate state

Identifying conformational changes in kinesin family motors associated with nucleotide and microtubule (MT) binding is essential to determining an atomic-level model for force production and motion by the motors. Using the mobility of nucleotide analog spin probes bound at the active sites of kinesin family motors to monitor conformational changes, we previously demonstrated that, in the ADP state, the open nucleotide site closes upon MT binding [Naber, N., Minehardt, T. J., Rice, S., Chen, X., Grammer, J., Matuska, M., et al. (2003). Closing of the nucleotide pocket of kinesin family motors upon binding to microtubules. Science, 300, 798-801]. We now extend these studies to kinesin-1 (K) and ncd (nonclaret disjunctional protein) motors in ATP and ATP-analog states. Our results reveal structural differences between several triphosphate and transition-state analogs bound to both kinesin and ncd in solution. The spectra of kinesin/ncd in the presence of SLADP•AlFx/BeFx and kinesin, with the mutation E236A (K-E236A; does not hydrolyze ATP) bound to ATP, show an open conformation of the nucleotide pocket similar to that seen in the kinesin/ncd•ADP states. In contrast, the triphosphate analogs K•SLAMPPNP and K-E236A•SLAMPPNP induce a more immobilized component of the electron paramagnetic resonance spectrum, implying closing of the nucleotide site. The MT-bound states of all of the triphosphate analogs reveal two novel spectral components. The equilibrium between these two components is only weakly dependent on temperature. Both components have more restricted mobility than observed in MT-bound diphosphate states. Thus, the closing of the nucleotide pocket when the diphosphate state binds to MTs is amplified in the triphosphate state, perhaps promoting accelerated ATP hydrolysis. Consistent with this idea, molecular dynamics simulations show a good correlation between our spectroscopic data, X-ray crystallography, and the electron microscopy of MT-bound triphosphate-analog states.

Conformational analysis of a nitroxide side chain in an α-helix with density functional theory

In site directed spin labeling, a nitroxide side chain is introduced at a selected site in a protein; the most commonly used is a disulfide-linked side chain designated R1. The electron paramagnetic resonance (EPR) spectra of R1, and the interspin distance between pairs of R1 residues as determined by dipolar EPR spectroscopy, encode a wealth of information on the protein structure and dynamics. However, extracting this information requires structural and dynamical models of the R1 side chain, that is, the favored rotamers, the intraresidue interactions that stabilize them, and the internal modes of motion. X-ray crystal structures of R1 in proteins have revealed a set of preferred rotamers in the crystal lattice. To identify the intraresidue interactions that stabilize particular rotamers of R1 in the absence of interactions with nearby side chains in a helix, and to evaluate models for the internal motion of the side chain, quantum mechanical calculations were performed on a relevant fragment of R1 in a 10-residue α-helix. Relative rotamer energies were determined in the gas phase, and solvation energies were estimated from a continuum solvent model that includes both electrostatic and hydrophobic contributions. The results identified preferred rotamers that are in agreement with the X-ray crystallographic studies. The rotamers are apparently stabilized by intraresidue sulfur-backbone interactions, suggesting that the preferred rotamers may be the same at all solvent-exposed helix sites.

Phosphorylation-induced structural changes in smooth muscle myosin regulatory light chain

We have performed complementary time-resolved fluorescence resonance energy transfer (TR-FRET) experiments and molecular dynamics (MD) simulations to elucidate structural changes in the phosphorylation domain (PD) of smooth muscle regulatory light chain (RLC) bound to myosin. PD is absent in crystal structures, leaving uncertainty about the mechanism of regulation. Donor-acceptor pairs of probes were attached to three site-directed di-Cys mutants of RLC, each having one Cys at position 129 in the C-terminal lobe and the other at position 2, 3, or 7 in the N-terminal PD. Labeled RLC was reconstituted onto myosin subfragment 1 (S1). TR-FRET resolved two simultaneously populated structural states of RLC, closed and open, in both unphosphorylated and phosphorylated biochemical states. All three FRET pairs show that phosphorylation shifts the equilibrium toward the open state, increasing its mol fraction by approximately 20%. MD simulations agree with experiments in remarkable detail, confirming the coexistence of two structural states, with phosphorylation shifting the system toward the more dynamic open structural state. This agreement between experiment and simulation validates the additional structural details provided by MD simulations: In the closed state, PD is bent onto the surface of the C-terminal lobe, stabilized by interdomain salt bridges. In the open state, PD is more helical and straight, resides farther from the C-terminal lobe, and is stabilized by an intradomain salt bridge. The result is a vivid atomic-resolution visualization of the first step in the molecular mechanism by which phosphorylation activates smooth muscle.

Analysis of side chain rotational restrictions of membrane-embedded proteins by spin-label ESR spectroscopy

Site-directed spin-labeling electron spin resonance (SDSL-ESR) is a promising tool for membrane protein structure determination. Here we propose a novel way to translate the local structural constraints gained by SDSL-ESR data into a low-resolution structure of a protein by simulating the restrictions of the local conformational spaces of the spin label attached at different protein sites along the primary structure of the membrane-embedded protein. We test the sensitivity of this approach for membrane-embedded M13 major coat protein decorated with a limited number of strategically placed spin labels employing high-throughput site-directed mutagenesis. We find a reasonably good agreement of the simulated and the experimental data taking a protein conformation close to the one determined by fluorescence resonance energy transfer analysis [P.V. Nazarov, R.B.M. Koehorst, W.L. Vos, V.V. Apanasovich, M.A. Hemminga, FRET study of membrane proteins: determination of the tilt and orientation of the N-terminal domain of M13 major coat protein, Biophys. J. 92 (2007) 1296-1305].

Actin-binding cleft closure in myosin II probed by site-directed spin labeling and pulsed EPR

We present a structurally dynamic model for nucleotide- and actin-induced closure of the actin-binding cleft of myosin, based on site-directed spin labeling and electron paramagnetic resonance (EPR) in Dictyostelium myosin II. The actin-binding cleft is a solvent-filled cavity that extends to the nucleotide-binding pocket and has been predicted to close upon strong actin binding. Single-cysteine labeling sites were engineered to probe mobility and accessibility within the cleft. Addition of ADP and vanadate, which traps the posthydrolysis biochemical state, influenced probe mobility and accessibility slightly, whereas actin binding caused more dramatic changes in accessibility, consistent with cleft closure. We engineered five pairs of cysteine labeling sites to straddle the cleft, each pair having one label on the upper 50-kDa domain and one on the lower 50-kDa domain. Distances between spin-labeled sites were determined from the resulting spin-spin interactions, as measured by continuous wave EPR for distances of 0.7-2 nm or pulsed EPR (double electron-electron resonance) for distances of 1.7-6 nm. Because of the high distance resolution of EPR, at least two distinct structural states of the cleft were resolved. Each of the biochemical states tested (prehydrolysis, posthydrolysis, and rigor), reflects a mixture of these structural states, indicating that the coupling between biochemical and structural states is not rigid. The resulting model is much more dynamic than previously envisioned, with both open and closed conformations of the cleft interconverting, even in the rigor actomyosin complex.

Interactions between branched DNAs and peptide inhibitors of DNA repair

The RecG helicase of Escherichia coli unwinds both Holliday junction (HJ) and replication fork DNA substrates. Our lab previously identified and characterized peptides (WRWYCR and KWWCRW) that block the activity of RecG on these substrates. We determined that the peptides bind HJ DNA and prevent the binding of RecG. Herein, we present further evidence that the peptides are competitive inhibitors of RecG binding to its substrates. We have generated structural models of interactions between WRWYCR and a junction substrate. Using the fluorescent probe 2-aminopurine, we show that inhibitors interact with highest affinity with HJs (K_d = 14 nM) and ∼4- to 9-fold more weakly with replication fork substrates. The fluorescence assay results agree with the structural model, and predict the molecular basis for interactions between HJ-trapping peptides and branched DNA molecules. Specifically, aromatic amino acids in the peptides stack with bases at the center of the DNA substrates. These interactions are stabilized by hydrogen bonds to the DNA and by intrapeptide interactions. These peptides inhibit several proteins involved in DNA repair in addition to RecG, have been useful as tools to dissect recombination, and possess antibiotic activity. Greater understanding of the peptides’ mechanism of action will further increase their utility.

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

A simulated continuous wave electron paramagnetic resonance spectrum of a nitroxide spin label can be obtained from the Fourier transform of a free induction decay. It has been previously shown that the free induction decay can be calculated by solving the time-dependent stochastic Liouville equation for a set of Brownian trajectories defining the rotational dynamics of the label. In this work, a quaternion-based Monte Carlo algorithm has been developed to generate Brownian trajectories describing the global rotational diffusion of a spin-labeled protein. Also, molecular dynamics simulations of two spin-labeled mutants of T4 lysozyme, T4L F153R1, and T4L K65R1 have been used to generate trajectories describing the internal dynamics of the protein and the local dynamics of the spin-label side chain. Trajectories from the molecular dynamics simulations combined with trajectories describing the global rotational diffusion of the protein are used to account for all of the dynamics of a spin-labeled protein. Spectra calculated from these combined trajectories correspond well to the experimental spectra for the buried site T4L F153R1 and the helix surface site T4L K65R1. This work provides a framework to further explore the modeling of the dynamics of the spin-label side chain in the wide variety of labeling environments encountered in site-directed spin labeling studies.

Structure and dynamics of the force-generating domain of myosin probed by multifrequency electron paramagnetic resonance

Spin-labeling and multifrequency EPR spectroscopy were used to probe the dynamic local structure of skeletal myosin in the region of force generation. Subfragment 1 (S1) of rabbit skeletal myosin was labeled with an iodoacetamide spin label at C707 (SH1). X- and W-band EPR spectra were recorded for the apo state and in the presence of ADP and nucleotide analogs. EPR spectra were analyzed in terms of spin-label rotational motion within myosin by fitting them with simulated spectra. Two models were considered: rapid-limit oscillation (spectrum-dependent on the orientational distribution only) and slow restricted motion (spectrum-dependent on the rotational correlation time and the orientational distribution). The global analysis of spectra obtained at two microwave frequencies (9.4 GHz and 94 GHz) produced clear support for the second model and enabled detailed determination of rates and amplitudes of rotational motion and resolution of multiple conformational states. The apo biochemical state is well-described by a single structural state of myosin (M) with very restricted slow motion of the spin label. The ADP-bound biochemical state of myosin also reveals a single structural state (M*, shown previously to be the same as the post-powerstroke ATP-bound state), with less restricted slow motion of the spin label. In contrast, the extra resolution available at 94 GHz reveals that the EPR spectrum of the S1.ADP.V(i)-bound biochemical state of myosin, which presumably mimics the S1.ADP.P(i) state, is resolved clearly into three spectral components (structural states). One state is indistinguishable from that of the ADP-bound state (M*) and is characterized by moderate restriction and slow motion, with a mole fraction of 16%. The remaining 84% (M**) contains two additional components and is characterized by fast rotation about the x axis of the spin label. After analyzing EPR spectra, myosin ATPase activity, and available structural information for myosin II, we conclude that post-powerstroke and pre-powerstroke structural states (M* and M**) coexist in the S1.ADP.V(i) biochemical state. We propose that the pre-powerstroke state M** is characterized by two structural states that could reflect flexibility between the converter and N-terminal domains of myosin.

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

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

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 novel approach to the simulation of nitroxide spin label EPR spectra from a single truncated dynamical trajectory

A simple effective method for calculation of EPR spectra from a single truncated dynamical trajectory of spin probe orientations is reported. It is shown that an accurate simulation can be achieved from the small initial fraction of a dynamical trajectory until the point when the autocorrelation function of re-orientational motion of spin label has relaxed. This substantially reduces the amount of time for spectra simulation compared to previous approaches, which require multiple full length trajectories (normally of several microseconds) to achieve the desired resolution of EPR spectra. Our method is applicable to trajectories generated from both Brownian dynamics and molecular dynamics (MD) calculations. Simulations of EPR spectra from Brownian dynamical trajectories under a variety of motional conditions including bi-modal dynamics with different hopping rates between the modes are compared to those performed by conventional method. Since the relatively short timescales of spin label motions are realistically accessible by modern MD computational methods, our approach, for the first time, opens the prospect of the simulation of EPR spectra entirely from MD trajectories of real proteins structures.

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.

Dynamics of the nitroxide side chain in spin-labeled proteins

The dynamics of the tether linking methanethiosulfonate (MTSSL) spin probes to alpha-helices has been investigated with the purpose of rationalizing its effects on ESR line shapes. Torsional profiles for the chain bonds have been calculated ab initio, and steric interactions with the alpha-helix and the neighboring residues have been introduced at the excluded-volume level. As a consequence of the restrictions deriving from chain geometry and local constraints, a limited number of allowed conformers has been identified that undergo torsional oscillations and conformational jumps. Torsional fluctuations are described as damped oscillations, while transition rates between conformers are calculated according to the Langer multidimensional extension of the Kramers theory. The time scale and amplitude of the different motions are compared; the major role played by rotations of the outermost bonds of the side chain emerges, along with the effects of substituents in the pyrroline ring on the conformer distribution and dynamics. The extent and symmetry of magnetic tensor averaging produced by the side chain motions are estimated, the implications for the ESR spectra of spin-labeled proteins are discussed, and suggestions for the introduction of realistic features of the spin probe dynamics into the line shape simulation are presented.

Modeling a spin-labeled fusion peptide in a membrane: implications for the interpretation of EPR experiments

Site-directed spin-labeling and electron paramagnetic resonance are powerful tools for studying structure and conformational dynamics of proteins, especially in membranes. The position of the spin label is used as an indicator of the position of the site to which it is attached. The interpretation of these experiments is based on the assumptions that the spin label does not affect the peptide configuration and that it has a fixed orientation and distance with respect to the protein backbone. Here, the validity of these assumptions is examined through implicit membrane molecular dynamics simulations of the influenza hemagglutinin fusion peptide that has been labeled with methanethiosulfonate spin label. We find that the methanethiosulfonate spin label can occasionally induce peptide orientations that differ from those adopted by the wild-type peptide. Furthermore, the spin-label resides, on average, several Angstroms deeper in the membrane than the corresponding backbone C(alpha)-atom even at sites pointing toward the solvent. The nitroxide spin label exhibits flexibility and adopts various configurations depending on the surrounding residues.

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.

Motional restrictions of membrane proteins: a site-directed spin labeling study

Site-directed mutagenesis was used to produce 27 single cysteine mutants of bacteriophage M13 major coat protein spanning the whole primary sequence of the protein. Single-cysteine mutants were labeled with nitroxide spin labels and incorporated into phospholipid bilayers with increasing acyl chain length. The SDSL is combined with ESR and CD spectroscopy. CD spectroscopy provided information about the overall protein conformation in different mismatching lipids. The spin label ESR spectra were analyzed in terms of a new spectral simulation approach based on hybrid evolutionary optimization and solution condensation. This method gives the residue-level free rotational space (i.e., the effective space within which the spin label can wobble) and the diffusion constant of the spin label attached to the protein. The results suggest that the coat protein has a large structural flexibility, which facilitates a stable protein-to-membrane association in lipid bilayers with various degrees of hydrophobic mismatch.

Molecular motion of spin labeled side chains in the C-terminal domain of RGL2 protein: A SDSL-EPR and MD study

Five singly spin labeled side chains at surface sites in the C-terminal domain of RGL2 protein have been analyzed to investigate the general relationship between nitroxide side chain mobility and protein structure. At these sites, the structural perturbation produced by replacement of a native residue with a nitroxide side chain appears to be very slight at the level of the backbone fold. The primary determinants of the nitroxide side chain mobility are backbone dynamics and tertiary interactions. On the exposed surfaces of alpha-helices, the side chain mobility is not restricted by tertiary interactions but appears to be determined by backbone dynamics, while in loop sites, the side chain mobility is even higher. For a better understanding of the changes in the EPR spectral line shape, molecular dynamics simulations were performed and found in agreement with EPR spectral data.

Influence of the disulfide bond configuration on the dynamics of the spin label attached to cytochrome c

A series of multi-nanosecond molecular dynamics (MD) simulations of wild-type cytochrome c and its spin-labeled variants with the methanethiosulfonate moiety attached at position C102 were performed (1) to elucidate the effect of the spin probe presence on the protein structure and (2) to describe the structure and dynamics of the spin-label moiety. Comparisons with the reference crystal structure of cytochrome c (PDB entry: 1YCC) indicate that the protein secondary structure is well preserved during simulations of the wild-type cytochrome c but slightly changed in simulations of the cytochrome c labeled at position C102. At the time scale covered in our simulations, the spin label exhibits highly dynamical behavior. The number of observed distinct conformations of the spin label moiety is between 3 and 13. The spin probe was found to form short-lived hydrogen bonds with the protein. Temporary hydrophobic interactions between the probe and the protein were also detected. The MD simulations directly show that the disulfide bond in the tether linking a spin probe with a protein strongly influence the behavior of the nitroxide group. The conformational flexibility and interaction with the protein are different for each of the two low energy conformations of the disulfide bond.

Phospholamban pentamer quaternary conformation determined by in-gel fluorescence anisotropy

We measured in-gel fluorescence anisotropy of phospholamban (PLB) labeled with the biarsenical fluorophore FlAsH at three different sites on the cytoplasmic domain. The 6 kDa monomer bands of FlAsH-tetracysPLB showed high anisotropy (r = 0.29), reflecting null homotransfer and low mobility (S = 0.85) on the nanosecond time scale of the FlAsH fluorescence lifetime. 30 kDa bands (pentameric PLB) within the same lanes exhibited low anisotropy, suggesting intrapentameric fluorescence energy homotransfer between PLB subunits. FlAsH labels positioned at residue -6, 5, or 23 showed a graduated pattern of fluorescence depolarization corresponding to resonance energy transfer radii of 46 +/-2, 38 +/- 4, and <25 A, respectively. Pentamer anisotropy increased with heating or fluorescence photobleaching toward a maximum value similar to that determined for monomeric PLB. Fluorescence resonance energy heterotransfer was also observed in vitro and in vivo within PLB pentamers colabeled with FlAsH and the biarsenical fluorophore ReAsH. In vitro heterotransfer efficiencies were graduated by labeling position, in harmony with homotransfer results. The calculated transfer radii compare favorably to distances predicted by a computer molecular model of the phospholamban pentamer constructed from NMR solution structures. The data support a helical pinwheel model for the PLB pentamer, in which the cytoplasmic domains bend sharply outward from the central bundle of helices.

Force field impact and spin-probe modeling in molecular dynamics simulations of spin-labeled T4 lysozyme

Several attempts have been made to compute electron paramagnetic resonance (EPR) spectra of biomolecules, using motional models or simulated trajectories to describe dynamics. Ideally, the simulated trajectories should capture "fast" (picosecond) snapshots of spin-probe rotations accurately, while being lengthy enough to ensure a proper Fourier integration of the time-domain signal. It is the interplay of the two criteria that poses computational challenges to the method. In this context, an analysis of the spin-probe and protein conformational sampling and equilibration, with different force fields and with explicit solvent, may be a useful attempt. The present work reports a comparative study of the effect of the molecular dynamics (MD) force field on conformational sampling and equilibration in two spin-labeled T4 lysozyme (T4L) variants, N40C and K48C. Ensembles of 10x 3 ns-trajectories per variant and per force field (OPLS/AMBER and AMBER99) are analyzed for a reliable assessment of convergence and sampling. It is found that subtle site-dependent differences in spin-probe rotations and torsions are more readily captured in the AMBER99 trajectories than in the OPLS/AMBER simulations. On the other hand, sampling and equilibration are found to be better with the OPLS/AMBER force field at equal trajectory lengths.

Site-directed spin labeling reveals a conformational switch in the phosphorylation domain of smooth muscle myosin

We have used site-directed spin labeling and EPR spectroscopy to detect structural changes within the regulatory light chain (RLC) of smooth muscle myosin upon phosphorylation. Smooth muscle contraction is activated by phosphorylation of S19 on RLC, but the structural basis of this process is unknown. There is no crystal structure containing a phosphorylated RLC, and there is no crystal structure for the N-terminal region of any RLC. Therefore, we have prepared single-Cys mutations throughout RLC, exchanged each mutant onto smooth muscle heavy meromyosin, verified normal regulatory function, and used EPR to determine dynamics and solvent accessibility at each site. A survey of spin-label sites throughout the RLC revealed that only the N-terminal region (first 24 aa) shows a significant change in dynamics upon phosphorylation, with most of the first 17 residues showing an increase in rotational amplitude. Therefore, we focused on this N-terminal region. Additional structural information was obtained from the pattern of oxygen accessibility along the sequence. In the absence of phosphorylation, little or no periodicity was observed, suggesting a lack of secondary structural order in this region. However, phosphorylation induced a strong helical pattern (3.6-residue periodicity) in the first 17 residues, while increasing accessibility throughout the first 24 residues. We have identified a domain within RLC, the N-terminal phosphorylation domain, in which phosphorylation increases helical order, internal dynamics, and accessibility. These results support a model in which this disorder-to-order transition within the phosphorylation domain results in decreased head-head interactions, activating myosin in smooth muscle.

Solvent Interactions and Protein Dynamics in Spin-labeled T4 Lysozyme

Aspects of T4 lysozyme dynamics and solvent interaction are investigated using atomically detailed Molecular Dynamics (MD) simulations. Two spin-labeled mutants of T4 lysozyme are analyzed (T4L-N40C and T4L-K48C), which have been found from electronic paramagnetic resonance (EPR) experiments to exhibit different mobilities at the site of spin probe attachment (N- and C-terminus of helix B, respectively). Similarities and differences in solvent distribution and diffusion around the spin label, as well as around exposed and buried residues within the protein, are discussed. The purpose is to capture possible strong interactions between the spin label (ring) and solvent molecules, which may affect EPR lineshapes. The effect of backbone motions on the water density profiles is also investigated. The focus is on the domain closure associated with the T4 lysozyme hinge-bending motion, which is analyzed by Essential Dynamics (ED). The N-terminus of helix B is found to be a "hinge" residue, which explains the high degree of flexibility and motional freedom at this site.

Coordination of the two heads of myosin during muscle contraction

We have used luminescence resonance energy transfer between regulatory light chains (RLC) to detect structural changes within the dimeric myosin molecule in contracting muscle fibers. Fully functional scallop muscle fibers were prepared such that each myosin molecule contained a terbium-labeled (luminescent donor) RLC on one head and a rhodamine-labeled (acceptor) RLC on the other. Time-resolved luminescence energy transfer between the two heads increased upon the transition from relaxation (ATP) to contraction (ATP plus Ca) and increased further in rigor (no ATP). Combined with experiments on mutant RLCs labeled specifically at other sites, these results support a model in which the force-generating weak-to-strong transition causes one myosin LC domain to tilt through a 30 degrees angle toward the other, thus acting as a coordinated lever arm.