Structural determination of spin label immobilization and orientation: a Monte Carlo minimization approach

Changes in the fraction of strongly attached cross bridges in mouse atrophic and hypertrophic muscles as revealed by continuous wave electron paramagnetic resonance

Electron paramagnetic resonance (EPR), coupled with site-directed spin labeling, has been proven to be a particularly suitable technique to extract information on the fraction of myosin heads strongly bound to actin upon muscle contraction. The approach can be used to investigate possible structural changes occurring in myosin of fiber s altered by diseases and aging. In this work, we labeled myosin at position Cys707, located in the SH1-SH2 helix in the myosin head cleft, with iodoacetamide spin label, a spin label that is sensitive to the reorientational motion of this protein during the ATPase cycle and characterized the biochemical states of the labeled myosin head by means of continuous wave EPR. After checking the sensitivity and the power of the technique on different muscles and species, we investigated whether changes in the fraction of strongly bound myosin heads might explain the contractile alterations observed in atrophic and hypertrophic murine muscles. In both conditions, the difference in contractile force could not be justified simply by the difference in muscle mass. Our results showed that in atrophic muscles the decrease in force generation was attributable to a lower fraction of strongly bound cross bridges during maximal activation. In contrast in hypertrophic muscles, the increase in force generation was likely due to several factors, as pointed out by the comparison of the EPR experiments with the tension measurements on single skinned fibers.

Characterization of the Domain Orientations of E. coli 5'-Nucleotidase by Fitting an Ensemble of Conformers to DEER Distance Distributions

Escherichia coli 5'-nucleotidase is a two-domain enzyme exhibiting a unique 96° domain motion that is required for catalysis. Here we present an integrated structural biology study that combines DEER distance distributions with structural information from X-ray crystallography and computational biology to describe the population of presumably almost isoenergetic open and closed states in solution. Ensembles of models that best represent the experimental distance distributions are determined by a Monte Carlo search algorithm. As a result, predominantly open conformations are observed in the unliganded state indicating that the majority of enzyme molecules await substrate binding for the catalytic cycle. The addition of a substrate analog yields ensembles with an almost equal mixture of open and closed states. Thus, in the presence of substrate, efficient catalysis is provided by the simultaneous appearance of open conformers (binding substrate or releasing product) and closed conformers (enabling the turnover of the substrate).

Exploring Structure, Dynamics, and Topology of Nitroxide Spin-Labeled Proteins Using Continuous-Wave Electron Paramagnetic Resonance Spectroscopy

Structural and dynamical characterization of proteins is of central importance in understanding the mechanisms underlying their biological functions. Site-directed spin labeling (SDSL) combined with continuous-wave electron paramagnetic resonance (CW EPR) spectroscopy has shown the capability of providing this information with site-specific resolution under physiological conditions for proteins of any degree of complexity, including those associated with membranes. This chapter introduces methods commonly employed for SDSL and describes selected CW EPR-based methods that can be applied to (1) map secondary and tertiary protein structure, (2) determine membrane protein topology, (3) measure protein backbone flexibility, and (4) reveal the existence of conformational exchange at equilibrium.

Multiple Drug Transport Pathways through Human P-Glycoprotein

P-Glycoprotein (P-gp) is a plasma membrane efflux pump that is commonly associated with therapy resistances in cancers and infectious diseases. P-gp can lower the intracellular concentrations of many drugs to subtherapeutic levels by translocating them out of the cell. Because of the broad range of substrates transported by P-gp, overexpression of P-gp causes multidrug resistance. We reported previously on dynamic transitions of P-gp as it moved through conformations based on crystal structures of homologous ABCB1 proteins using in silico targeted molecular dynamics techniques. We expanded these studies here by docking transport substrates to drug binding sites of P-gp in conformations open to the cytoplasm, followed by cycling the pump through conformations that opened to the extracellular space. We observed reproducible transport of two substrates, daunorubicin and verapamil, by an average of 11-12 Å through the plane of the membrane as P-gp progressed through a catalytic cycle. Methylpyrophosphate, a ligand that should not be transported by P-gp, did not show this movement through P-gp. Drug binding to either of two subsites on P-gp appeared to determine the initial pathway used for drug movement through the membrane. The specific side-chain interactions with drugs within each pathway seemed to be, at least in part, stochastic. The docking and transport properties of a P-gp inhibitor, tariquidar, were also studied. A mechanism of inhibition by tariquidar that involves stabilization of an outward open conformation with tariquidar bound in intracellular loops or at the drug binding domain of P-gp is presented.

Full Atom Simulations of Spin Label Conformations

Interpretation of EPR measurables from spin labels in terms of structure and dynamics requires knowledge of label behavior. General strategies were developed for simulation of labels used in EPR of proteins. The criteria for those simulations are (a) exhaustive sampling of rotamer space, (b) consensus of results independent of starting points, and (c) inclusion of entropy. These criteria are satisfied only when the number of transitions in any dihedral angle exceeds 100 and the simulation maintains thermodynamic equilibrium. Methods such as conventional MD do not efficiently cross energetic barriers, simulated anealing, Monte Carlo or popular Rotamer Library methodologies are potential energy based and ignore entropy (in addition to their specific shortcomings: environment fluctuations, fixed environment, or electrostatics). The Simulated scaling method avoids the above flaws by modulating the force fields between a reduced (allowing crossing energy barriers) and full potential (sampling minima). Spin label diffuses on this surface while remaining in thermodynamic equilibrium. Simulations show that (a) adopting a single conformation is rare, often there are two to four populated rotamers and (b) position of the NO varies up to 16 Å. These results illustrate necessity for caution when interpreting EPR signals in terms of molecular structure. For example, the 10-16 Å distance change in DEER should not be interpreted as a large conformational change, it can well be a flip about Cα-Cβ bond. Rigorous exploration of possible rotamer structures of a spin label is paramount in signal interpretation. We advocate use of bifunctional labels, motion of which is restricted 10,000-fold and the NO position is restricted to 2-5 Å.

RosettaEPR: rotamer library for spin label structure and dynamics

An increasingly used parameter in structural biology is the measurement of distances between spin labels bound to a protein. One limitation to these measurements is the unknown position of the spin label relative to the protein backbone. To overcome this drawback, we introduce a rotamer library of the methanethiosulfonate spin label (MTSSL) into the protein modeling program Rosetta. Spin label rotamers were derived from conformations observed in crystal structures of spin labeled T4 lysozyme and previously published molecular dynamics simulations. Rosetta’s ability to accurately recover spin label conformations and EPR measured distance distributions was evaluated against 19 experimentally determined MTSSL labeled structures of T4 lysozyme and the membrane protein LeuT and 73 distance distributions from T4 lysozyme and the membrane protein MsbA. For a site in the core of T4 lysozyme, the correct spin label conformation (Χ₁ and Χ₂) is recovered in 99.8% of trials. In surface positions 53% of the trajectories agree with crystallized conformations in Χ₁ and Χ₂. This level of recovery is on par with Rosetta performance for the 20 natural amino acids. In addition, Rosetta predicts the distance between two spin labels with a mean error of 4.4 Å. The width of the experimental distance distribution, which reflects the flexibility of the two spin labels, is predicted with a mean error of 1.3 Å. RosettaEPR makes full-atom spin label modeling available to a wide scientific community in conjunction with the powerful suite of modeling methods within Rosetta.

Tunable, mixed-resolution modeling using library-based Monte Carlo and graphics processing units

Building on our recently introduced library-based Monte Carlo (LBMC) approach, we describe a flexible protocol for mixed coarse-grained (CG)/all-atom (AA) simulation of proteins and ligands. In the present implementation of LBMC, protein side chain configurations are pre-calculated and stored in libraries, while bonded interactions along the backbone are treated explicitly. Because the AA side chain coordinates are maintained at minimal run-time cost, arbitrary sites and interaction terms can be turned on to create mixed-resolution models. For example, an AA region of interest such as a binding site can be coupled to a CG model for the rest of the protein. We have additionally developed a hybrid implementation of the generalized Born/surface area (GBSA) implicit solvent model suitable for mixed-resolution models, which in turn was ported to a graphics processing unit (GPU) for faster calculation. The new software was applied to study two systems: (i) the behavior of spin labels on the B1 domain of protein G (GB1) and (ii) docking of randomly initialized estradiol configurations to the ligand binding domain of the estrogen receptor (ERα). The performance of the GPU version of the code was also benchmarked in a number of additional systems.

Orientation and rotational motions of single molecules by polarized total internal reflection fluorescence microscopy (polTIRFM)

In this article, we describe methods to detect the spatial orientation and rotational dynamics of single molecules using polarized total internal reflection fluorescence microscopy (polTIRFM). polTIRFM determines the three-dimensional angular orientation and the extent of wobble of a fluorescent probe bound to the macromolecule of interest. We discuss single-molecule versus ensemble measurements, as well as single-molecule techniques for orientation and rotation, and fluorescent probes for orientation studies. Using calmodulin (CaM) as an example of a target protein, we describe a method for labeling CaM with bifunctional rhodamine (BR). We also describe the physical principles and experimental setup of polTIRFM. We conclude with a brief introduction to assays using polTIRFM to assess the interaction of actin and myosin.

Simulating the dynamics and orientations of spin-labeled side chains in a protein-DNA complex

Site-directed spin labeling, wherein a nitroxide side chain is introduced into a protein at a selected mutant site, is increasingly employed to investigate biological systems by electron spin resonance (ESR) spectroscopy. An understanding of the packing and dynamics of the spin label is needed to extract the biologically relevant information about the macromolecule from ESR measurements. In this work, molecular dynamics (MD) simulations were performed on the spin-labeled restriction endonuclease, EcoRI in complex with DNA. Mutants of this homodimeric enzyme were previously constructed, and distance measurements were performed using the double electron electron resonance experiment. These correlated distance constraints have been leveraged with MD simulations to learn about side chain packing and preferred conformers of the spin label on sites in an α-helix and a β-strand. We found three dihedral angles of the spin label side chain to be most sensitive to the secondary structure where the spin label was located. Conformers sampled by the spin label differed between secondary structures as well. C(α)-C(α) distance distributions were constructed and used to extract details about the protein backbone mobility at the two spin labeled sites. These simulation studies enhance our understanding of the behavior of spin labels in proteins and thus expand the ability of ESR spectroscopy to contribute to knowledge of protein structure and dynamics.

Computer modeling of nitroxide spin labels on proteins

Electron paramagnetic resonance using site-directed spin labeling can be used as an approach for determination of protein structures that are difficult to solve by other methods. One important aspect of this approach is the measurement of interlabel distances using the double electron-electron resonance (DEER) method. Interpretation of experimental data could be facilitated by a computational approach to calculation of interlabel distances. We describe an algorithm, PRONOX, for rapid computation of interlabel distances based on calculation of spin label conformer distributions at any site of a protein. The program incorporates features of the label distribution established experimentally, including weighting of favorable conformers of the label. Distances calculated by PRONOX were compared with new DEER distances for amphiphysin and annexin B12 and with published data for FCHo2 (F-BAR), endophilin, and α-synuclein, a total of 44 interlabel distances. The program reproduced these distances accurately (r(2) = 0.94, slope = 0.98). For 9 of the 11 distances for amphiphysin, PRONOX reproduced the experimental data to within 2.5 Å. The speed and accuracy of PRONOX suggest that the algorithm can be used for fitting to DEER data for determination of protein tertiary structure.

Broad disorder and the allosteric mechanism of myosin II regulation by phosphorylation

Double electron electron resonance EPR methods was used to measure the effects of the allosteric modulators, phosphorylation, and ATP, on the distances and distance distributions between the two regulatory light chain of myosin (RLC). Three different states of smooth muscle myosin (SMM) were studied: monomers, the short-tailed subfragment heavy meromyosin, and SMM filaments. We reconstituted myosin with nine single cysteine spin-labeled RLC. For all mutants we found a broad distribution of distances that could not be explained by spin-label rotamer diversity. For SMM and heavy meromyosin, several sites showed two heterogeneous populations in the unphosphorylated samples, whereas only one was observed after phosphorylation. The data were consistent with the presence of two coexisting heterogeneous populations of structures in the unphosphorylated samples. The two populations were attributed to an on and off state by comparing data from unphosphorylated and phosphorylated samples. Models of these two states were generated using a rigid body docking approach derived from EM [Wendt T, Taylor D, Trybus KM, Taylor K (2001) Proc Natl Acad Sci USA 98:4361-4366] (PNAS, 2001, 98:4361-4366), but our data revealed a new feature of the off-state, which is heterogeneity in the orientation of the two RLC. Our average off-state structure was very similar to the Wendt model reveal a new feature of the off state, which is heterogeneity in the orientations of the two RLC. As found previously in the EM study, our on-state structure was completely different from the off-state structure. The heads are splayed out and there is even more heterogeneity in the orientations of the two RLC.

EPR in protein science : intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) form a unique protein category characterized by the absence of a well-defined structure and by remarkable conformational flexibility. Electron Paramagnetic Resonance (EPR) spectroscopy combined with site-directed spin labeling (SDSL) is amongst the most suitable methods to unravel their structure and dynamics. This review summarizes the tremendous methodological developments in the area of SDSL EPR and its applications in protein research. Recent results on the intrinsically disordered Parkinson's disease protein α-synuclein illustrate that the method has gained increasing attention in IDP research. SDSL EPR has now reached a level where broad application in this rapidly advancing field is feasible.

The structure of p85ni in class IA phosphoinositide 3-kinase exhibits interdomain disorder

Regulation of the class IA PI 3-kinase involves inhibition and stabilization of the catalytic subunit (p110) by the regulatory subunit (p85). Regulation is achieved by two major contacts: a stable interface involving the adapter-binding domain (ABD) of p110 and the inter-SH2 (iSH2) domain of p85 and a regulatory interaction between the N-terminal SH2 (nSH2) domain of p85 and the helical domain of p110. In the present study, we have examined the relative orientation of the nSH2 and iSH2 of p85alpha using site-directed spin labeling and pulsed EPR. Surprisingly, both distance measurements and distance distributions suggest that the nSH2 domain is highly disordered relative to the iSH2 domain. Molecular modeling based on EPR distance restraints suggests that the nSH2 domain moves in a hinge-like manner, sampling a torus space around the proximal end of the iSH2 domain. These data have important implications for the mechanism by which p85/p110 dimers are regulated by phosphopeptides.

Osmolyte perturbation reveals conformational equilibria in spin-labeled proteins

Recent evidence suggests that proteins at equilibrium can exist in a manifold of conformational substates, and that these substates play important roles in protein function. Therefore, there is great interest in identifying regions in proteins that are in conformational exchange. Electron paramagnetic resonance spectra of spin-labeled proteins containing the nitroxide side chain (R1) often consist of two (or more) components that may arise from slow exchange between conformational substates (lifetimes > 100 ns). However, crystal structures of proteins containing R1 have shown that multicomponent spectra can also arise from equilibria between rotamers of the side chain itself. In this report, it is shown that these scenarios can be distinguished by the response of the system to solvent perturbation with stabilizing osmolytes such as sucrose. Thus, site-directed spin labeling (SDSL) emerges as a new tool to explore slow conformational exchange in proteins of arbitrary size, including membrane proteins in a native-like environment. Moreover, equilibrium between substates with even modest differences in conformation is revealed, and the simplicity of the method makes it suitable for facile screening of multiple proteins. Together with previously developed strategies for monitoring picosecond to millisecond backbone dynamics, the results presented here expand the timescale over which SDSL can be used to explore protein flexibility.

Accommodating discontinuities in dimeric left-handed coiled coils in ATP synthase external stalks

ATP synthases from coupling membranes are complex rotary motors that convert the energy of proton gradients across coupling membranes into the chemical potential of the beta-gamma anhydride bond of ATP. Proton movement within the ring of c subunits localized in the F(0)-sector drives gamma and epsilon rotation within the F(1)alpha(3)beta(3) catalytic core where substrates are bound and products are released. An external stalk composed of homodimeric subunits b(2) in Escherichia coli or heterodimeric bb' in photosynthetic synthases connects F(0) subunit a with F(1) subunits delta and most likely alpha. The external stalk resists rotation, and is of interest both functionally and structurally. Hypotheses that the external stalk contributes to the overall efficiency of the reaction through elastic coupling of rotational substeps, and that stalks form staggered, right-handed coiled coils, are investigated here. We report on different structures that accommodate heptad discontinuities with either local or global underwinding. Analyses of the knob-and-hole packing of the E. coli b(2) and Synechocystis bb' stalks strongly support the possibility that these proteins can adopt conventional left-handed coiled coils.

Nucleotide-induced flexibility change in neck linkers of dimeric kinesin as detected by distance Measurements using spin-labeling EPR

Using dipolar continuous-wave and pulsed electron paramagnetic resonance methods, we have determined the distribution of the distances between two spin labels placed on the middle of each of the neck linkers of dimeric kinesin. In the absence of microtubules, the distance was centered at 3.3 nm, but displayed a broad distribution with a width of 2.7 nm. This broad distribution implies that the linkers are random coils and extend well beyond the 2.5-nm distance expected of crystal structures. In the presence of microtubules, two linker populations were found: one similar to that observed in the absence of microtubules (a broad distribution centered at 3.3 nm), and the second population with a narrower distribution centered at 1.3-2.5 nm. In the absence of nucleotide but in the presence of microtubules, approximately 40% of the linkers were at a distance centered at 1.9 nm with a 1.2-nm width; the remaining fraction was at 3.3 nm, as before. This suggests that neck linkers exhibit dynamics covering a wide distance range between 1.0 and 5.0 nm. In the presence of ATP analogs adenosine 5'-(beta,gamma-imido)triphosphate and adenosine 5'-(gamma-thio)triphosphate, 40-50% of the spins showed a very narrow distribution centered at 1.6 nm, with a width of 0.4-0.5 nm. The remaining population displayed the broad 3.3-nm distribution. Under these conditions, a large fraction of linkers are docked firmly onto a motor core or microtubule, while the remainder is disordered. We propose that large nucleotide-dependent flexibility changes in the linkers contribute to the directional bias of the kinesin molecule stepping 8 nm along the microtubule.

Accurate long-range distance measurements in a doubly spin-labeled protein by a four-pulse, double electron-electron resonance method

Distance determination in disordered systems by a four-pulse double electron-electron resonance method (DEER or PELDOR) is becoming increasingly popular because long distances (several nanometers) and their distributions can be measured. From the distance distributions eventual heterogeneities and dynamics can be deduced. To make full use of the method, typical distance distributions for structurally well-defined systems are needed. Here, the structurally well-characterized protein azurin is investigated by attaching two (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl) methanethiosulfonate spin labels (MTSL) by site-directed mutagenesis. Mutations at the surface sites of the protein Q12, K27, and N42 are combined in the double mutants Q12C/K27C and K27C/N42C. A distance of 4.3 nm is found for Q12C/K27C and 4.6 nm for K27C/N42C. For Q12C/K27C the width of the distribution (0.24 nm) is smaller than for the K27C/N42C mutant (0.36 nm). The shapes of the distributions are close to Gaussian. These distance distributions agree well with those derived from a model to determine the maximally accessible conformational space of the spin-label linker. Additionally, the expected distribution for the shorter distance variant Q12C/N42C was modeled. The width is larger than the calculated one for Q12C/K27C by 21%, revealing the effect of the different orientation and shorter distance. The widths and the shapes of the distributions are suited as a reference for two unperturbed MTSL labels at structurally well-defined sites.

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.

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.

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.

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.

Structure of the cytosolic part of the subunit b-dimer of Escherichia coli F0F1-ATP synthase

The structure of the external stalk and its function in the catalytic mechanism of the F(0)F(1)-ATP synthase remains one of the important questions in bioenergetics. The external stalk has been proposed to be either a rigid stator that binds F(1) or an elastic structural element that transmits energy from the small rotational steps of subunits c to the F(1) sector during catalysis. We employed proteomics, sequence-based structure prediction, molecular modeling, and electron spin resonance spectroscopy using site-directed spin labeling to understand the structure and interfacial packing of the Escherichia coli b-subunit homodimer external stalk. Comparisons of bacterial, cyanobacterial, and plant b-subunits demonstrated little sequence similarity. Supersecondary structure predictions, however, show that all compared b-sequences have extensive heptad repeats, suggesting that the proteins all are capable of packing as left-handed coiled-coils. Molecular modeling subsequently indicated that b(2) from the E. coli ATP synthase could pack into stable left-handed coiled-coils. Thirty-eight substitutions to cysteine in soluble b-constructs allowed the introduction of spin labels and the determination of intersubunit distances by ESR. These distances correlated well with molecular modeling results and strongly suggest that the E. coli subunit b-dimer can stably exist as a left-handed coiled-coil.

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.

Computation of nitroxide-nitroxide distances in spin-labeled DNA duplexes

Nanometer distances in nucleic acids can be measured by EPR using two 1-oxyl-2,2,5,5-tetramethylpyrroline radicals, with each label attached via a methylene group to a phosphorothioate-substituted backbone position as one of two phosphorothioate diastereomers (R(P) and S(P)). Correlating the internitroxide distance to the geometry of the parent molecule requires computational analysis of the label conformers. Here, we report sixteen 4-ns MD simulations on a DNA duplex d(CTACTGCTTTAG) .d(CTAAAGCAGTAG) with label pairs at C7/C19, T5/A17, and T2/T14, respectively. For each labeled duplex, four simulations were performed with S(P)/S(P), R(P)/R(P), S(P)/R(P), and R(P)/S(P) labels, with initial all trans label conformations. Another set of four simulations was performed for the 7/19-labeled duplex using a different label starting conformation. The average internitroxide distance r(MD) was within 0.2 A for the two sets of simulations for the 7/19-labeled duplex, indicating sufficient sampling of conformational space. For all three labeled duplexes studied, r(MD) agreed with experimental values, as well as with average distances obtained from an efficient conformer search algorithm (NASNOX). The simulations also showed that the labels have conformational preferences determined by the linker chemistry and label-DNA interactions. These results establish computational algorithms that allow use of the 1-oxyl-2,2,5,5-tetramethylpyrroline label for mapping global structures of nucleic acids.

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.

Distance measurements on spin-labelled biomacromolecules by pulsed electron paramagnetic resonance

The biological function of protein, DNA, and RNA molecules often depends on relative movements of domains with dimensions of a few nanometers. This length scale can be accessed by distance measurements between spin labels if pulsed electron paramagnetic resonance (EPR) techniques such as electron-electron double resonance (ELDOR) and double-quantum EPR are used. The approach does not require crystalline samples and is well suited to biomacromolecules with an intrinsic flexibility as distributions of distances can be measured. Furthermore, oligomerization or complexation of biomacromolecules can also be studied, even if it is incomplete. The sensitivity of the technique and the reliability of the measured distance distribution depend on careful optimization of the experimental conditions and procedures for data analysis. Interpretation of spin-to-spin distance distributions in terms of the structure of the biomacromolecules furthermore requires a model for the conformational distribution of the spin labels.

Measuring distances by pulsed dipolar ESR spectroscopy: spin-labeled histidine kinases

Applications of dipolar ESR spectroscopy to structural biology are rapidly expanding, and it has become a useful method that is aimed at resolving protein structure and functional mechanisms. The method of pulsed dipolar ESR spectroscopy (PDS) is outlined in the first half of the chapter, and it illustrates the simplicity and potential of this developing technology with applications to various biological systems. A more detailed description is presented of the implementation of PDS to reconstruct the ternary structure of a large dimeric protein complex from Thermotoga maritima, formed by the histidine kinase CheA and the coupling protein CheW. This protein complex is a building block of an extensive array composed of coupled supramolecular structures assembled from CheA/CheW proteins and transmembrane signaling chemoreceptors, which make up a sensor that is key to controlling the motility in bacterial chemotaxis. The reconstruction of the CheA/CheW complex has employed several techniques, including X-ray crystallography and pulsed ESR. Emphasis is on the role of PDS, which is part of a larger effort to reconstruct the entire signaling complex, including chemoreceptor, by means of PDS structural mapping. In order to precisely establish the mode of coupling of CheW to CheA and to globally map the complex, approximately 70 distances have already been determined and processed into molecular coordinates by readily available methods of distance geometry constraints.

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.

Solution structure of the cytoplasmic domain of erythrocyte membrane band 3 determined by site-directed spin labeling

The cytoplasmic domain of the anion exchange protein (cdb3) serves as a critical organizing center for protein-protein interactions that stabilize the erythrocyte membrane. The structure of the central core of cdb3, determined by X-ray crystallography from crystals grown at pH 4.8, revealed a compact dimer for residues 55-356 and unresolved N- and C-termini on each monomer [Zhang et al. (2000) Blood 96, 2925-2933]. Given that previous studies had suggested a highly asymmetric structure for cdb3 and that pH dependent structural transitions of cdb3 have been reported, the structure of cdb3 in solution at neutral pH was investigated via site-directed spin labeling in combination with conventional electron paramagnetic resonance (EPR) and double electron electron resonance (DEER) spectroscopies. These studies show that the structure of the central compact dimer (residues 55-356) is indistinguishable from the crystal structure determined at pH 4.8. N-Terminal residues 1-54 and C-terminal residues 357-379 are dynamically disordered and show no indications of stable secondary structure. These results establish a structural model for cdb3 in solution at neutral pH which represents an important next step in characterizing structural details of the protein-protein interactions that stabilize the erythrocyte membrane.

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.

Rotational motions of macro-molecules by single-molecule fluorescence microscopy

Several complementary techniques have been developed to determine average orientation, dynamics on multiple time scales, and concerted rotational motions of individual fluorescent probes bound to biological macromolecules. In both protein domains and nucleic acids, tilting and wobble are relevant to their functional mechanisms. Here we briefly review methods to detect angles and rotational motions of single fluorophores and give an example of three-dimensional, total internal reflection, single-molecule fluorescence polarization applied to actin as it is translocated by conventional muscle myosin.

Explicit Treatment of Spin Labels in Modeling of Distance Constraints from Dipolar EPR and DEER

Current SDSL-EPR methods allow measurement of dipolar distances in the 8-70 A range; however, the use of extrinsic probes complicates the interpretation of these distances in modeling macromolecular structure and conformational changes. The data presented here show that interprobe distances correlate only weakly with Cbeta-Cbeta distances, especially for distances that are on the order of the spin label tether lengths. Explicitly incorporating the spin label into the modeling process increases the experiment/model correlation 4-fold and reduces the distance error from 6 A to 3 A.

Conformational switching in muscle

Muscles can be studied as complex systems of many interacting proteins and investigated at many different levels of organization. This talk will describe how we modeled the mechanism of Ca activation, the structure of the muscle proteins, and protein complexes (including actin monomers, tropomyosin and troponin complexes, and myosin) to examine two different scientific problems: the mechano-chemical energy transduction mechanism, and the control system of that mechanism. The methods we used--saturation transfer electron paramagnetic resonance, phosphorescence anisotropy, and fluorescence resonance energy transfer--reveal two specific structures: a hinge between the motor and regulatory domains, and a stiff regulatory domain. This indicates that the structure of the myosin head is capable of generating translating conformational changes within the motor domain to the swing of the regulatory domain, and that the regulatory domain is rigid enough to act as a lever arm.

Myosin regulatory domain orientation in skeletal muscle fibers: application of novel electron paramagnetic resonance spectral decomposition and molecular modeling methods

Reorientation of the regulatory domain of the myosin head is a feature of all current models of force generation in muscle. We have determined the orientation of the myosin regulatory light chain (RLC) using a spin-label bound rigidly and stereospecifically to the single Cys-154 of a mutant skeletal isoform. Labeled RLC was reconstituted into skeletal muscle fibers using a modified method that results in near-stoichiometric levels of RLC and fully functional muscle. Complex electron paramagnetic resonance spectra obtained in rigor necessitated the development of a novel decomposition technique. The strength of this method is that no specific model for a complex orientational distribution was presumed. The global analysis of a series of spectra, from fibers tilted with respect to the magnetic field, revealed two populations: one well-ordered (+/-15 degrees ) with the spin-label z axis parallel to actin, and a second population with a large distribution (+/-60 degrees ). A lack of order in relaxed or nonoverlap fibers demonstrated that regulatory domain ordering was defined by interaction with actin rather than the thick filament surface. No order was observed in the regulatory domain during isometric contraction, consistent with the substantial reorientation that occurs during force generation. For the first time, spin-label orientation has been interpreted in terms of the orientation of a labeled domain. A Monte Carlo conformational search technique was used to determine the orientation of the spin-label with respect to the protein. This in turn allows determination of the absolute orientation of the regulatory domain with respect to the actin axis. The comparison with the electron microscopy reconstructions verified the accuracy of the method; the electron paramagnetic resonance determined that axial orientation was within 10 degrees of the electron microscopy model.

Binding of Calcium Ions to an Avian Flight Muscle Troponin T†

Numerous troponin T (TnT) isoforms are generated by alternative RNA splicing primarily in its N-terminal hypervariable region, but the functions of these isoforms are not completely understood. Here for the first time, we discovered that a chicken fast TnT isoform with a unique Tx motif (HEEAH)(n) binds calcium. The metal binding behavior of this TnT isoform was first investigated using terbium as a calcium analogue due to its more readily detectable fluorescence variation upon TnT binding. Both intact TnT and TnT N-terminal fragment (TnT N47) bound terbium with high affinity indicating that the N-terminal sequence was the site of binding. Since terbium often substitutes at calcium-binding sites, radioactive calcium was tested and found to bind both intact TnT and TnT N47. Fluorescence measurements using the calcium-sensitive fluorescent dye, calcium green 5N, confirmed that calcium bound to the tertiary complex of TnT and the tropomyosin dimer with a fast on-rate (10(6)-10(7) M(-1) s(-1)) as detected in stopped-flow analysis. Consistent with these observations, computational predictions suggest that TnT N47 might fold into an elongated structure with at least one high-affinity metal ion binding pocket comprised primarily of the Tx motif sequence and several lower affinity binding sites. These results suggest that TnT may play a role in modulating the calcium-mediated regulatory process of striated muscle contraction.

Exploring the conformational space of membrane protein folds matching distance constraints

Herein we present a computational technique for generating helix-membrane protein folds matching a predefined set of distance constraints, such as those obtained from NMR NOE, chemical cross-linking, dipolar EPR, and FRET experiments. The purpose of the technique is to provide initial structures for local conformational searches based on either energetic considerations or ad-hoc scoring criteria. In order to properly screen the conformational space, the technique generates an exhaustive list of conformations within a specified root-mean-square deviation (RMSD) where the helices are positioned in order to match the provided distances. Our results indicate that the number of structures decreases exponentially as the number of distances increases, and increases exponentially as the errors associated with the distances increases. We also found the number of solutions to be smaller when all the distances share one helix in common, compared to the case where the distances connect helices in a daisy-chain manner. We found that for 7 helices, at least 15 distances with errors up to 8 A are needed to produce a number of solutions that is not too large to be processed by local search refinement procedures. Finally, without energetic considerations, our enumeration technique retrieved the transmembrane domains of Bacteriorhodopsin (PDB entry1c3w), Halorhodopsin (1e12), Rhodopsin (1f88), Aquaporin-1 (1fqy), Glycerol uptake facilitator protein (1fx8), Sensory Rhodopsin (1jgj), and a subunit of Fumarate reductase flavoprotein (1qlaC) with Calpha level RMSDs of 3.0 A, 2.3 A, 3.2 A, 4.6 A, 6.0 A, 3.7 A, and 4.4 A, respectively.

Structure of the inhibitory region of troponin by site directed spin labeling electron paramagnetic resonance

Site-directed spin labeling EPR (SDSL-EPR) was used to determine the structure of the inhibitory region of TnI in the intact cardiac troponin ternary complex. Maeda and collaborators have modeled the inhibitory region of TnI (skeletal 96-112: the structural motif that communicates the Ca(2+) signal to actin) as a kinked alpha-helix [Vassylyev, D., Takeda, S., Wakatsuki, S., Maeda, K. & Maeda, Y. (1998) Proc. Natl. Acad. Sci. USA 95, 4847-4852), whereas Trewhella and collaborators have proposed the same region to be a flexible beta-hairpin [Tung, C. S., Wall, M. E., Gallagher, S. C. & Trewhella, J. (2000) Protein Sci. 9, 1312-1326]. To distinguish between the two models, residues 129-145 of cardiac TnI were mutated sequentially to cysteines and labeled with the extrinsic spin probe, MTSSL. Sequence-dependent solvent accessibility was measured as a change in power saturation of the spin probe in the presence of the relaxation agent. In the ternary complex, the 129-137 region followed a pattern characteristic of a regular 3.6 residues/turn alpha-helix. The following region, residues 138-145, showed no regular pattern in solvent accessibility. Measurements of 4 intradomain distances within the inhibitory sequence, using dipolar EPR, were consistent with an alpha-helical structure. The difference in side-chain mobility between the ternary (C.I.T) and binary (C.I) complexes revealed a region of interaction of TnT located at the N-terminal end of the inhibitory sequence, residues 130-135. The above findings for the troponin complex in solution do not support either of the computational models of the binary complex; however, they are in very good agreement with a preliminary report of the x-ray structure of the cardiac ternary complex [Takeda, S. Yamashita, A., Maeda, K. & Maeda, Y. (2002) Biophys. J. 82, 832].