Structural dynamics of the actomyosin complex probed by a bifunctional spin label that cross-links SH1 and SH2

Structural Dynamics of Protein Interactions Using Site-Directed Spin Labeling of Cysteines to Measure Distances and Rotational Dynamics with EPR Spectroscopy

Here we review applications of site-directed spin labeling (SDSL) with engineered cysteines in proteins, to study the structural dynamics of muscle and non-muscle proteins, using and developing the electron paramagnetic resonance (EPR) spectroscopic techniques of dipolar EPR, double electron electron resonance (DEER), saturation transfer EPR (STEPR), and orientation measured by EPR. The SDSL technology pioneered by Wayne Hubbell and collaborators has greatly expanded the use of EPR, including the measurement of distances between spin labels covalently attached to proteins and peptides. The Thomas lab and collaborators have applied these techniques to elucidate dynamic interactions in the myosin-actin complex, myosin-binding protein C, calmodulin, ryanodine receptor, phospholamban, utrophin, dystrophin, β-III-spectrin, and Aurora kinase. The ability to design and engineer cysteines in proteins for site-directed covalent labeling has enabled the use of these powerful EPR techniques to measure distances, while showing that they are complementary with optical spectroscopy measurements.

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins

Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is an established tool for exploring protein structure and dynamics. Although nitroxide side chains attached to a single cysteine via a disulfide linkage are commonly employed in SDSL-EPR, their internal flexibility complicates applications to monitor slow internal motions in proteins and to structure determination by distance mapping. Moreover, the labile disulfide linkage prohibits the use of reducing agents often needed for protein stability. To enable the application of SDSL-EPR to the measurement of slow internal dynamics, new spin labels with hindered internal motion are desired. Here, we introduce a highly ordered nitroxide side chain, designated R9, attached at a single cysteine residue via a non-reducible thioether linkage. The reaction to introduce R9 is highly selective for solvent-exposed cysteine residues. Structures of R9 at two helical sites in T4 Lysozyme were determined by X-ray crystallography and the mobility in helical sequences was characterized by EPR spectral lineshape analysis, Saturation Transfer EPR, and Saturation Recovery EPR. In addition, interspin distance measurements between pairs of R9 residues are reported. Collectively, all data indicate that R9 will be useful for monitoring slow internal structural fluctuations, and applications to distance mapping via dipolar spectroscopy and relaxation enhancement methods are anticipated.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00723-023-01618-8.

Characterization of Bifunctional Spin Labels for Investigating the Structural and Dynamic Properties of Membrane Proteins Using EPR Spectroscopy

Site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is a very powerful technique to study structural and dynamic properties of membrane proteins. The most widely used spin label is methanthiosulfonate (MTSL). However, the flexibility of this spin label introduces greater uncertainties in EPR measurements obtained for determining structures, side-chain dynamics, and backbone motion of membrane protein systems. Recently, a newer bifunctional spin label (BSL), 3,4-bis(methanethiosulfonylmethyl)-2,2,5,5-tetramethyl-2,5-dihydro-1H-pyrrol-1-yloxy, has been introduced to overcome the dynamic limitations associated with the MTSL spin label and has been invaluable in determining protein backbone dynamics and inter-residue distances due to its restricted internal motion and fewer size restrictions. While BSL has been successful in providing more accurate information about the structure and dynamics of several proteins, a detailed characterization of the spin label is still lacking. In this study, we characterized BSLs by performing CW-EPR spectral line shape analysis as a function of temperature on spin-labeled sites inside and outside of the membrane for the integral membrane protein KCNE1 in POPC/POPG lipid bilayers and POPC/POPG lipodisq nanoparticles. The experimental data revealed a powder pattern spectral line shape for all of the KCNE1-BSL samples at 296 K, suggesting the motion of BSLs approaches the rigid limit regime for these series of samples. BSLs were further utilized to report for the first time the distance measurement between two BSLs attached on an integral membrane protein KCNE1 in POPC/POPG lipid bilayers at room temperature using dipolar line broadening CW-EPR spectroscopy. The CW dipolar line broadening EPR data revealed a 15 ± 2 Å distance between doubly attached BSLs on KCNE1 (53/57-63/67) which is consistent with molecular dynamics modeling and the solution NMR structure of KCNE1 which yielded a distance of 17 Å. This study demonstrates the utility of investigating the structural and dynamic properties of membrane proteins in physiologically relevant membrane mimetics using BSLs.

A Drosophila model of dominant inclusion body myopathy type 3 shows diminished myosin kinetics that reduce muscle power and yield myofibrillar defects

Individuals with inclusion body myopathy type 3 (IBM3) display congenital joint contractures with early-onset muscle weakness that becomes more severe in adulthood. The disease arises from an autosomal dominant point mutation causing an E706K substitution in myosin heavy chain type IIa. We have previously expressed the corresponding myosin mutation (E701K) in homozygous Drosophila indirect flight muscles and recapitulated the myofibrillar degeneration and inclusion bodies observed in the human disease. We have also found that purified E701K myosin has dramatically reduced actin-sliding velocity and ATPase levels. Since IBM3 is a dominant condition, we now examine the disease state in heterozygote Drosophila in order to gain a mechanistic understanding of E701K pathogenicity. Myosin ATPase activities in heterozygotes suggest that approximately equimolar levels of myosin accumulate from each allele. In vitro actin sliding velocity rates for myosin isolated from the heterozygotes were lower than the control, but higher than for the pure mutant isoform. Although sarcomeric ultrastructure was nearly wild type in young adults, mechanical analysis of skinned indirect flight muscle fibers revealed a 59% decrease in maximum oscillatory power generation and an approximately 20% reduction in the frequency at which maximum power was produced. Rate constant analyses suggest a decrease in the rate of myosin attachment to actin, with myosin spending decreased time in the strongly bound state. These mechanical alterations result in a one-third decrease in wing beat frequency and marginal flight ability. With aging, muscle ultrastructure and function progressively declined. Aged myofibrils showed Z-line streaming, consistent with the human heterozygote phenotype. Based upon the mechanical studies, we hypothesize that the mutation decreases the probability of the power stroke occurring and/or alters the degree of movement of the myosin lever arm, resulting in decreased in vitro motility, reduced muscle power output and focal myofibrillar disorganization similar to that seen in individuals with IBM3.

Summary: Reduced muscle power output and progressive myofibrillar defects in a Drosophila model of inclusion body myopathy 3 arise from the decreased rate of weak to strong actin-binding transition of myosin.

A combined EPR and MD simulation study of a nitroxyl spin label with restricted internal mobility sensitive to protein dynamics

EPR studies combined with fully atomistic Molecular Dynamics (MD) simulations and an MD-EPR simulation method provide evidence for intrinsic low rotameric mobility of a nitroxyl spin label, Rn, compared to the more widely employed label MTSL (R1). Both experimental and modelling results using two structurally different sites of attachment to Myoglobin show that the EPR spectra of Rn are more sensitive to the local protein environment than that of MTSL. This study reveals the potential of using the Rn spin label as a reporter of protein motions.

The use of the Rx spin label in orientation measurement on proteins, by EPR

The bipedal spin label Rx is more restricted in its conformation and dynamics than its monopodal counterpart R1. To systematically investigate the utility of the Rx label, we have attempted to comprehensively survey the attachment of Rx to protein secondary structures. We have examined the formation, structure and dynamics of the spin label in relation to the underlying protein in order to determine feasibility and optimum conditions for distance and orientation measurement by pulsed EPR. The labeled proteins have been studied using molecular dynamics, CW EPR, pulsed EPR distance measurement at X-band and orientation measurement at W-band. The utility of different modes and positions of attachment have been compared and contrasted.

A bifunctional spin label reports the structural topology of phospholamban in magnetically-aligned bicelles

We have applied a bifunctional spin label and EPR spectroscopy to determine membrane protein structural topology in magnetically-aligned bicelles, using monomeric phospholamban (PLB) as a model system. Bicelles are a powerful tool for studying membrane proteins by NMR and EPR spectroscopies, where magnetic alignment yields topological constraints by resolving the anisotropic spectral properties of nuclear and electron spins. However, EPR bicelle studies are often hindered by the rotational mobility of monofunctional Cys-linked spin labels, which obscures their orientation relative to the protein backbone. The rigid and stereospecific TOAC label provides high orientational sensitivity but must be introduced via solid-phase peptide synthesis, precluding its use in large proteins. Here we show that a bifunctional methanethiosulfonate spin label attaches rigidly and stereospecifically to Cys residues at i and i+4 positions along PLB's transmembrane helix, thus providing orientational resolution similar to that of TOAC, while being applicable to larger membrane proteins for which synthesis is impractical. Computational modeling and comparison with NMR data shows that these EPR experiments provide accurate information about helix tilt relative to the membrane normal, thus establishing a robust method for determining structural topology in large membrane proteins with a substantial advantage in sensitivity over NMR.

Bifunctional Spin Labeling of Muscle Proteins: Accurate Rotational Dynamics, Orientation, and Distance by EPR

While EPR allows for the characterization of protein structure and function due to its exquisite sensitivity to spin label dynamics, orientation, and distance, these measurements are often limited in sensitivity due to the use of labels that are attached via flexible monofunctional bonds, incurring additional disorder and nanosecond dynamics. In this chapter, we present methods for using a bifunctional spin label (BSL) to measure muscle protein structure and dynamics. We demonstrate that bifunctional attachment eliminates nanosecond internal rotation of the spin label, thereby allowing the accurate measurement of protein backbone rotational dynamics, including microsecond-to-millisecond motions by saturation transfer EPR. BSL also allows for accurate determination of helix orientation and disorder in mechanically and magnetically aligned systems, due to the label's stereospecific attachment. Similarly, labeling with a pair of BSL greatly enhances the resolution and accuracy of distance measurements measured by double electron-electron resonance (DEER). Finally, when BSL is applied to a protein with high helical content in an assembly with high orientational order (e.g., muscle fiber or membrane), two-probe DEER experiments can be combined with single-probe EPR experiments on an oriented sample in a process we call BEER, which has the potential for ab initio high-resolution structure determination.

High-resolution helix orientation in actin-bound myosin determined with a bifunctional spin label

Using electron paramagnetic resonance (EPR) of a bifunctional spin label (BSL) bound stereospecifically to Dictyostelium myosin II, we determined with high resolution the orientation of individual structural elements in the catalytic domain while myosin is in complex with actin. BSL was attached to a pair of engineered cysteine side chains four residues apart on known α-helical segments, within a construct of the myosin catalytic domain that lacks other reactive cysteines. EPR spectra of BSL-myosin bound to actin in oriented muscle fibers showed sharp three-line spectra, indicating a well-defined orientation relative to the actin filament axis. Spectral analysis indicated that orientation of the spin label can be determined within <2.1° accuracy, and comparison with existing structural data in the absence of nucleotide indicates that helix orientation can also be determined with <4.2° accuracy. We used this approach to examine the crucial ADP release step in myosin's catalytic cycle and detected reversible rotations of two helices in actin-bound myosin in response to ADP binding and dissociation. One of these rotations has not been observed in myosin-only crystal structures.

A nucleotide-independent cyclic nitroxide label for monitoring segmental motions in nucleic acids

Background

Spin labels, which are chemically stable radicals attached at specific sites of a bio-molecule, enable investigations on structure and dynamics of proteins and nucleic acids using techniques such as site-directed spin labeling and paramagnetic NMR. Among spin labels developed, the class of rigid labels have limited or no independent motions between the radical bearing moiety and the target, and afford a number of advantages in measuring distances and monitoring local dynamics within the parent bio-molecule. However, a general method for attaching a rigid label to nucleic acids in a nucleotide-independent manner has not been reported.

Results

We developed an approach for installing a nearly rigid nitroxide spin label, designated as R5c, at a specific site of the nucleic acid backbone in a nucleotide-independent manner. The method uses a post-synthesis approach to covalently attach the nitroxide moiety in a cyclic fashion to phosphorothioate groups introduced at two consecutive nucleotides of the target strand. R5c-labeled nucleic acids are capable of pairing with their respective complementary strands, and the cyclic nature of R5c attachment significantly reduced independence motions of the label with respect to the parent duplex, although it may cause distortion of the local environment at the site of labeling. R5c yields enhanced sensitivity to the collective motions of the duplex, as demonstrated by its capability to reveal changes in collective motions of the substrate recognition duplex of the 120-kDa Tetrahymena group I ribozyme, which elude detection by a flexible label.

Conclusions

The cyclic R5c nitroxide can be efficiently attached to a target nucleic acid site using a post-synthetic coupling approach conducted under mild biochemical conditions, and serves as a viable label for experimental investigation of segmental motions in nucleic acids, including large folded RNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-015-0019-5) contains supplementary material, which is available to authorized users.

Electron paramagnetic resonance resolves effects of oxidative stress on muscle proteins

We have used site-directed spin labeling and electron paramagnetic resonance (EPR) to explore the effects of oxidation on muscle function, with particular focus on the actin-myosin interaction. EPR measurements show that aging or oxidative modification causes a decrease in the fraction of myosins in the strong-binding state, which can be traced to the actin-binding cleft of the myosin catalytic domain.

Coupling of lever arm swing and biased Brownian motion in actomyosin

An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement. In this study, we address the long-standing controversy between the lever-arm hypothesis and the Brownian motion hypothesis through in silico observations of an actomyosin system. We study a system composed of myosin II and actin filament by calculating free-energy landscapes of actin-myosin interactions using the molecular dynamics method and by simulating transitions among dynamically changing free-energy landscapes using the Monte Carlo method. The results obtained by this combined multi-scale calculation show that myosin with inorganic phosphate (P_i) and ADP weakly binds to actin and that after releasing P_i and ADP, myosin moves along the actin filament toward the strong-binding site by exhibiting the biased Brownian motion, a behavior consistent with the observed single-molecular behavior of myosin. Conformational flexibility of loops at the actin-interface of myosin and the N-terminus of actin subunit is necessary for the distinct bias in the Brownian motion. Both the 5.5–11 nm displacement due to the biased Brownian motion and the 3–5 nm displacement due to lever-arm swing contribute to the net displacement of myosin. The calculated results further suggest that the recovery stroke of the lever arm plays an important role in enhancing the displacement of myosin through multiple cycles of ATP hydrolysis, suggesting a unified movement mechanism for various members of the myosin family.

Myosin II is a molecular motor that is fueled by ATP hydrolysis and generates mechanical force by interacting with actin filament. Comparison among various myosin structures obtained by X-ray and electron microscope analyses has led to the hypothesis that structural change of myosin in ATP hydrolysis cycle is the driving mechanism of force generation. However, single-molecule experiments have suggested an alternative mechanism in which myosin moves stochastically in a biased direction along actin filament. Computer simulation serves as a platform for assessing these hypotheses by revealing the prominent features of the dynamically changing landscape of actin-myosin interaction. The calculated results show that myosin binds to actin at different locations of actin filament in the weak- and strong-binding states and that the free energy has a global gradient from the weak-binding site to the strong-binding site. Myosin relaxing into the strong-binding state therefore necessarily shows the biased Brownian motion toward the strong-binding site. Lever-arm swing is induced during this relaxation process; therefore, lever-arm swing and the biased Brownian motion are coupled to contribute to the net displacement of myosin. This coupling should affect the dynamical behaviors of muscle and cardiac systems.

Output-input ratio in thermally fluctuating biomolecular machines

Biological molecular machines are proteins that operate under isothermal conditions and hence are referred to as free energy transducers. They can be formally considered as enzymes that simultaneously catalyze two chemical reactions: the free energy-donating (input) reaction and the free energy-accepting (output) one. Most if not all biologically active proteins display a slow stochastic dynamics of transitions between a variety of conformational substates composing their native state. This makes the description of the enzymatic reaction kinetics in terms of conventional rate constants insufficient. In the steady state, upon taking advantage of the assumption that each reaction proceeds through a single pair (the gate) of transition conformational substates of the enzyme-substrates complex, the degree of coupling between the output and the input reaction fluxes has been expressed in terms of the mean first-passage times on a conformational transition network between the distinguished substates. The theory is confronted with the results of random-walk simulations on the five-dimensional hypercube. The formal proof is given that, for single input and output gates, the output-input degree of coupling cannot exceed unity. As some experiments suggest such exceeding, looking for the conditions for increasing the degree of coupling value over unity challenges the theory. Performed simulations of random walks on several model networks involving more extended gates indicate that the case of the degree of coupling value higher than 1 is realized in a natural way on critical branching trees extended by long-range shortcuts. Such networks are scale-free and display the property of the small world. For short-range shortcuts, the networks are scale-free and fractal, representing a reasonable model for biomolecular machines displaying tight coupling, i.e., the degree of coupling equal exactly to unity. A hypothesis is stated that the protein conformational transition networks, as just as higher-level biological networks, the protein interaction network, and the metabolic network, have evolved in the process of self-organized criticality.

Novel configuration of a myosin II transient intermediate analogue revealed by quick-freeze deep-etch replica electron microscopy

In the present paper, we described our attempt to characterize the rough three-dimensional features of the structural analogue of the key intermediate of myosin's cross-bridge cycle. Using quick-freeze deep-etch replica electron microscopy, we observed that actin-attached myosin during in vitro sliding was bent superficially as postulated by the conventional hypothesis, but in the opposite direction of the putative pre-power-stroke configuration, as for ADP·Vi (inorganic vanadate)-bound myosin. We searched for the conformational species with a similar appearance and found that SH1-SH2 (thiols 1 and 2)-cross-linked myosin is a good candidate. To characterize such small asymmetric structures, we employed a new pattern-recognition procedure that accommodates the metal-replicated samples. In this method, the best-matched views of the target microscopic images were selected from a comprehensive set of images simulated from known atomic co-ordinates of relevant proteins. Together with effective morphological filtering, we could define the conformational species and the view angles of the catalytic domain and the lever arm cropped from averaged images of disulfide-cross-linked myosin. Whereas the catalytic domain of the new conformer closely resembled the pPDM (N,N'-p-phenylenedimaleimide)-treated, but SH2 Lys705-cross-linked, structure (PDB code 1L2O), a minor product of the same cross-linking reaction, the lever arm projected differently. Using separately determined view angles of the catalytic domain and the lever arm, we built a model of disulfide-cross-linked myosin. Further combination with the 'displacement-mapping' procedure enabled us to reconstruct the global three-dimensional envelope of the unusual structure whose lever arm orientation is compatible with our reports on the actin-sliding cross-bridge structure. Assuming this conformer as the structural analogue of the transient intermediate during actin sliding, the power stroke of the lever arm might accompany the reversal of the disorganized SH1 helix.

Conformationally trapping the actin-binding cleft of myosin with a bifunctional spin label

We have trapped the catalytic domain of Dictyostelium (Dicty) myosin II in a weak actin-binding conformation by chemically crosslinking two engineered cysteines across the actin-binding cleft, using a bifunctional spin label (BSL). By connecting the lower and upper 50 kDa domains of myosin, the crosslink restricts the conformation of the actin-binding cleft. Crosslinking has no effect on the basal ATPase activity of isolated myosin, but it impairs rigor actin binding and actin-activation of myosin ATPase. EPR spectra of BSL provide insight into actomyosin structural dynamics. BSL is highly immobilized within the actin-binding cleft and is thus exquisitely sensitive to the global orientation and rotational motions of the myosin head. Conventional EPR shows that myosin heads bound to oriented actin filaments are highly disordered with respect to the actin filament axis, in contrast to the nearly crystalline order of myosin heads in rigor. This disorder is similar to that of weakly bound heads induced by ATP, but saturation transfer EPR shows that the disorder of crosslinked myosin is at least 100 times slower. Thus this cleft-crosslinked myosin is remarkably similar, in both actin affinity and rotational dynamics, to SH1-SH2 crosslinked BSL-myosin S1. We conclude that, whether myosin is trapped at the actin-myosin interface or in the force-generating region between the active site and lever arm, the structural state of myosin is intermediate between the weak-binding state preceding phosphate release and the strong-binding state that succeeds it. We propose that it represents the threshold of force generation.

Temperature-enhanced association of proteins due to electrostatic interaction: a coarse-grained simulation of actin-myosin binding

Association of protein molecules constitutes the basis for the interaction network in a cell. Despite its fundamental importance, the thermodynamic aspect of protein-protein binding, particularly the issues relating to the entropy change upon binding, remains elusive. The binding of actin and myosin, which are vital proteins in motility, is a typical example, in which two different binding mechanisms have been argued: the binding affinity increases with increasing temperature and with decreasing salt-concentration, indicating the entropy-driven binding and the enthalpy-driven binding, respectively. How can these thermodynamically different binding mechanisms coexist? To address this question, which is of general importance in understanding protein-protein bindings, we conducted an in silico titration of the actin-myosin system by molecular dynamics simulation using a residue-level coarse-grained model, with particular focus on the role of the electrostatic interaction. We found a good agreement between in silico and in vitro experiments on the salt-concentration dependence and the temperature dependence of the binding affinity. We then figured out how the two binding mechanisms can coexist: the enthalpy (due to electrostatic interaction between actin and myosin) provides the basal binding affinity, and the entropy (due to the orientational disorder of water molecules) enhances it at higher temperatures. In addition, we analyzed the actin-myosin complex structures observed during the simulation and obtained a variety of weak-binding complex structures, among which were found an unusual binding mode suggested by an earlier experiment and precursor structures of the strong-binding complex proposed by electron microscopy. These results collectively indicate the potential capability of a residue-level coarse-grained model to simulate the association-dissociation dynamics (particularly for transient weak-bindings) exhibited by larger and more complicated systems, as in a cell.

Expression of the inclusion body myopathy 3 mutation in Drosophila depresses myosin function and stability and recapitulates muscle inclusions and weakness

A Drosophila model of myosin-based inclusion body myopathy type 3 is presented. Muscle function, ATPase activity, and actin sliding velocity were dramatically reduced. The mutant myosin is prone to aggregate, likely accounting for the observed cytoplasmic inclusions and disorganized muscle filaments reminiscent of the human disease.

Hereditary myosin myopathies are characterized by variable clinical features. Inclusion body myopathy 3 (IBM-3) is an autosomal dominant disease associated with a missense mutation (E706K) in the myosin heavy chain IIa gene. Adult patients experience progressive muscle weakness. Biopsies reveal dystrophic changes, rimmed vacuoles with cytoplasmic inclusions, and focal disorganization of myofilaments. We constructed a transgene encoding E706K myosin and expressed it in Drosophila (E701K) indirect flight and jump muscles to establish a novel homozygous organism with homogeneous populations of fast IBM-3 myosin and muscle fibers. Flight and jump abilities were severely reduced in homozygotes. ATPase and actin sliding velocity of the mutant myosin were depressed >80% compared with wild-type myosin. Light scattering experiments and electron microscopy revealed that mutant myosin heads bear a dramatic propensity to collapse and aggregate. Thus E706K (E701K) myosin appears far more labile than wild-type myosin. Furthermore, mutant fly fibers exhibit ultrastructural hallmarks seen in patients, including cytoplasmic inclusions containing aberrant proteinaceous structures and disorganized muscle filaments. Our Drosophila model reveals the unambiguous consequences of the IBM-3 lesion on fast muscle myosin and fibers. The abnormalities observed in myosin function and muscle ultrastructure likely contribute to muscle weakness observed in our flies and patients.

Three distinct actin-attached structural states of myosin in muscle fibers

We have used thiol cross-linking and electron paramagnetic resonance (EPR) to resolve structural transitions of myosin's light chain domain (LCD) and catalytic domain (CD) that are associated with force generation. Spin labels were incorporated into the LCD of muscle fibers by exchanging spin-labeled regulatory light chain for endogenous regulatory light chain, with full retention of function. To trap myosin in a structural state analogous to the elusive posthydrolysis ternary complex A.M'.D.P, we used pPDM to cross-link SH1 (Cys(707)) to SH2 (Cys(697)) on the CD. LCD orientation and dynamics were measured in three biochemical states: relaxation (A.M.T), SH1-SH2 cross-linked (A.M'.D.P analog), and rigor (A.M.D). EPR showed that the LCD of cross-linked fibers has an orientational distribution intermediate between relaxation and rigor, and saturation transfer EPR revealed slow rotational dynamics indistinguishable from that of rigor. Similar results were obtained for the CD using a bifunctional spin label to cross-link SH1-SH2, but the CD was more disordered than the LCD. We conclude that SH1-SH2 cross-linking traps a state in which both the CD and LCD are intermediate between relaxation (highly disordered and microsecond dynamics) and rigor (highly ordered and rigid), supporting the hypothesis that the cross-linked state is an A.M'D.P analog on the force generation pathway.

Structure and dynamics of a conformationally constrained nitroxide side chain and applications in EPR spectroscopy

A disulfide-linked nitroxide side chain (R1) is the most widely used spin label for determining protein topology, mapping structural changes, and characterizing nanosecond backbone motions by site-directed spin labeling. Although the internal motion of R1 and the number of preferred rotamers are limited, translating interspin distance measurements and spatial orientation information into structural constraints is challenging. Here, we introduce a highly constrained nitroxide side chain designated RX as an alternative to R1 for these applications. RX is formed by a facile cross-linking reaction of a bifunctional methanethiosulfonate reagent with pairs of cysteine residues at i and i + 3 or i and i + 4 in an α-helix, at i and i + 2 in a β-strand, or with cysteine residues in adjacent strands in a β-sheet. Analysis of EPR spectra, a crystal structure of RX in T4 lysozyme, and pulsed electron-electron double resonance (ELDOR) spectroscopy on an immobilized protein containing RX all reveal a highly constrained internal motion of the side chain. Consistent with the constrained geometry, interspin distance distributions between pairs of RX side chains are narrower than those from analogous R1 pairs. As an important consequence of the constrained internal motion of RX, spectral diffusion detected with ELDOR reveals microsecond internal motions of the protein. Collectively, the data suggest that the RX side chain will be useful for distance mapping by EPR spectroscopy, determining spatial orientation of helical segments in oriented specimens, and measuring structural fluctuations on the microsecond time scale.

Structural kinetics of myosin by transient time-resolved FRET

For many proteins, especially for molecular motors and other enzymes, the functional mechanisms remain unsolved due to a gap between static structural data and kinetics. We have filled this gap by detecting structure and kinetics simultaneously. This structural kinetics experiment is made possible by a new technique, (TR)(2)FRET (transient time-resolved FRET), which resolves protein structural states on the submillisecond timescale during the transient phase of a biochemical reaction. (TR)(2)FRET is accomplished with a fluorescence instrument that uses a pulsed laser and direct waveform recording to acquire an accurate subnanosecond time-resolved fluorescence decay every 0.1 ms after stopped flow. To apply this method to myosin, we labeled the force-generating region site specifically with two probes, mixed rapidly with ATP to initiate the recovery stroke, and measured the interprobe distance by (TR)(2)FRET with high resolution in both space and time. We found that the relay helix bends during the recovery stroke, most of which occurs before ATP is hydrolyzed, and two structural states (relay helix straight and bent) are resolved in each nucleotide-bound biochemical state. Thus the structural transition of the force-generating region of myosin is only loosely coupled to the ATPase reaction, with conformational selection driving the motor mechanism.

Multiscale modeling of structural dynamics underlying force generation and product release in actomyosin complex

To decrypt the mechanistic basis of myosin motor function, it is essential to probe the conformational changes in actomyosin with high spatial and temporal resolutions. In a computational effort to meet this challenge, we have performed a multiscale modeling of the allosteric couplings and transition pathway of actomyosin complex by combining coarse-grained modeling of the entire complex with all-atom molecular dynamics simulations of the active site. Our modeling of allosteric couplings at the pre-powerstroke state has pinpointed key actin-activated couplings to distant myosin parts which are critical to force generation and the sequential release of phosphate and ADP. At the post-powerstroke state, we have identified isoform-dependent couplings which underlie the reciprocal coupling between actin binding and nucleotide binding in fast Myosin II, and load-dependent ADP release in Myosin V. Our modeling of transition pathway during powerstroke has outlined a clear sequence of structural events triggered by actin binding, which lead to subsequent force generation, twisting of central beta-sheet, and the sequential release of phosphate and ADP. Finally we have performed atomistic simulations of active-site dynamics based on an on-path "transition-state" myosin conformation, which has revealed significantly weakened coordination of phosphate by Switch II, and a disrupted key salt bridge between Switch I and II. Meanwhile, the coordination of MgADP by Switch I and P loop is less perturbed. As a result, the phosphate can be released prior to MgADP. This study has shed new lights on the controversy over the structural mechanism of actin-activated phosphate release and force generation in myosin motor.

Site-directed spectroscopic probes of actomyosin structural dynamics

Spectroscopy of myosin and actin has entered a golden age. High-resolution crystal structures of isolated actin and myosin have been used to construct detailed models for the dynamic actomyosin interactions that move muscle. Improved protein mutagenesis and expression technologies have facilitated site-directed labeling with fluorescent and spin probes. Spectroscopic instrumentation has achieved impressive advances in sensitivity and resolution. Here we highlight the contributions of site-directed spectroscopic probes to understanding the structural dynamics of myosin II and its actin complexes in solution and muscle fibers. We emphasize studies that probe directly the movements of structural elements within the myosin catalytic and light-chain domains, and changes in the dynamics of both actin and myosin due to their alternating strong and weak interactions in the ATPase cycle. A moving picture emerges in which single biochemical states produce multiple structural states, and transitions between states of order and dynamic disorder power the actomyosin engine.

Mechanistic role of movement and strain sensitivity in muscle contraction

Tension generation can be studied by applying step perturbations to contracting muscle fibers and subdividing the mechanical response into exponential phases. The de novo tension-generating isomerization is associated with one of these phases. Earlier work has shown that a temperature jump perturbs the equilibrium constant directly to increase tension. Here, we show that a length jump functions quite differently. A step release (relative movement of thick and thin filaments) appears to release a steric constraint on an ensemble of noncompetent postphosphate release actomyosin cross-bridges, enabling them to generate tension, a concentration jump in effect. Structural studies [Taylor KA, et al. (1999) Tomographic 3D reconstruction of quick-frozen, Ca(2+)-activated contracting insect flight muscle. Cell 99:421-431] that map to these kinetics indicate that both catalytic and lever arm domains of noncompetent myosin heads change angle on actin, whereas lever arm movement alone mediates the power stroke. Together, these kinetic and structural observations show a 13-nm overall interaction distance of myosin with actin, including a final 4- to 6-nm power stroke when the catalytic domain is fixed on actin. Raising fiber temperature with both perturbation techniques accelerates the forward, but slows the reverse rate constant of tension generation, kinetics akin to the unfolding/folding of small proteins. Decreasing strain, however, causes both forward and reverse rate constants to increase. Despite these changes in rate, the equilibrium constant is strain-insensitive. Activation enthalpy and entropy data show this invariance to be the result of enthalpy-entropy compensation. Reaction amplitudes confirm a strain-invariant equilibrium constant and thus a strain-insensitive ratio of pretension- to tension-generating states as work is done.