Single molecule mechanics resolves the earliest events in force generation by cardiac myosin
Abstract
Key steps of cardiac mechanochemistry, including the force-generating working stroke and the release of phosphate (Pi), occur rapidly after myosin-actin attachment. An ultra-high-speed optical trap enabled direct observation of the timing and amplitude of the working stroke, which can occur within <200 μs of actin binding by β-cardiac myosin. The initial actomyosin state can sustain loads of at least 4.5 pN and proceeds directly to the stroke or detaches before releasing ATP hydrolysis products. The rates of these processes depend on the force. The time between binding and stroke is unaffected by 10 mM Pi which, along with other findings, indicates the stroke precedes phosphate release. After Pi release, Pi can rebind enabling reversal of the working stroke. Detecting these rapid events under physiological loads provides definitive indication of the dynamics by which actomyosin converts biochemical energy into mechanical work.
https://doi.org/10.7554/eLife.49266.001Introduction
Myosin is a cytoskeletal motor that uses metabolic energy stored in ATP to do mechanical work. This process underlies muscle contraction, cell motility, and intracellular transport. In muscle, myosin causes displacements of actin filaments connected to work against an external load. The majority of the mechanical work done by myosin is performed in the working stroke, the rotation of the myosin lever arm domain, which occurs soon after actin binding and is closely associated with the release of inorganic phosphate (Pi), a product of ATP hydrolysis. Understanding the connections between the formation of a force-bearing actin-bound state, the biochemical step of Pi release, and the mechanical action of the working stroke is fundamental to understanding energy transduction. Although substantial physiological and biophysical experiments have provided key details associated with these transitions (Bagshaw and Trentham, 1974; Hibberd and Trentham, 1986; Hibberd et al., 1985; Cooke et al., 1988; Dantzig et al., 1992; Araujo and Walker, 1996; White et al., 1997; Takagi et al., 2006; Muretta et al., 2015; Caremani et al., 2015; Rohde et al., 2017; Llinas et al., 2015), controversy remains regarding the order of the force-generating conformational transitions and whether the free energy change from phosphate release precedes or follows the working stroke (Figure 1a). Strengthening of the actomyosin bond due to conformational changes in myosin and actin, tilting of the lever arm, and Pi release all occur within milliseconds of initial actin binding, making the mechanochemical steps difficult to resolve by commonly available biophysical measurements.
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.49266.003
To detect the rapid sequence of events between actin binding and the working stroke, we examined the interactions of recombinant human beta-cardiac myosin (MHY7) with actin using the three-bead optical trapping geometry, where a bead-actin-bead ‘dumbbell’ assembly is brought near to a single myosin molecule attached to a pedestal bead (Finer et al., 1994). The myosin protein was a heavy meromyosin (HMM) construct with a Flag-tag placed at the C-terminus for purification. It was attached non-specifically to nitrocellulose coated silica beads on the coverslip surface. The beta-cardiac isoform not only has great relevance as the dominant myosin in the human ventricles, but its slower kinetic properties compared to fast skeletal myosin facilitated detection of fast rate processes. Although this is a two-headed construct, we did not observe evidence that both heads interacted simultaneously with actin in the presence of ATP (see Discussion). We employed an ultra-fast force clamp (UFFC) which can resolve the onset of actin-myosin binding events to less than 100 microseconds and filament motions at sub-nm resolution (Capitanio et al., 2012). Unlike standard optical trapping techniques which typically have temporal resolutions > 10 ms, UFFC allows observation of the initial myosin working stroke which can occur within 1 ms of actin binding. We applied this system to explore actomyosin attachment durations and displacements in the presence of 1.5–4.5 pN mechanical loads in directions along the long-axis of the actin filament that hinder or assist the stroke. The duration of the events, and the timing and size of filament displacements during interactions were examined in the presence and absence of 10 mM Pi to determine the relationship between actin binding, the working stroke, and Pi release.
Results
Duration of actomyosin interactions observed under load
UFFC uses a feedback loop to maintain a constant force on an actin filament dumbbell by rapid (µs) control of the trapping laser position. The applied force (the geometric sum of the optical trap forces on each bead) causes the filament and two beads to move through the assay solution at a constant velocity set by the bead size and water viscosity and the force is alternated in direction to maintain the dumbbell within 200 nm excursion (Figure 1b–e). When a surface-attached myosin binds to actin and becomes loaded with the applied force, the filament motion stops within ~50 µs (Figure 1d–g). Displacement of the actin filament caused by the myosin working stroke is monitored during the interaction, which ends when myosin dissociates, allowing the actin filament to resume motion under the applied force (Figure 1f). The beginning and end of each actomyosin interaction are detected by examining the velocity of the optical traps as they respond to the feedback system, dropping to zero when myosin binds and increasing again when myosin releases the actin (Figure 1f). When the working stroke displaces the filament toward its pointed end, the beads and the traps move a corresponding amount to maintain a constant force (Figure 1g).
We first analyzed the duration of actomyosin interactions in the absence of Pi under force that hinder or assist the stroke. To distinguish brief, reversible actomyosin attachments from longer attachments when myosins complete the ATPase cycle, most experiments were performed at a low (1 μM), rate-limiting MgATP concentration. Under both assisting and hindering loads, the distributions of event durations span a range from less than 1 ms to greater than 1 s (Figure 2a–e). The actomyosin events with durations less than 10 ms, which represent up to 50% of the observed interactions, have not been observed previously with standard optical trapping techniques because their durations are below the minimum detectable length (deadtime, 10 ms for standard techniques, 0.5–2.7 ms for the current technique) (Capitanio et al., 2012).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.49266.008
The very wide distributions of actomyosin event durations were best described by the linear combination of at least two exponential components, which we quantified using maximum likelihood estimation via the software package MEMLET (Woody et al., 2016). For forces > 2.25 pN, three exponential components were necessary to obtain an adequate fit (Figure 2a and f, Figure 2—figure supplement 1). This result contrasts with the distribution of durations observed with lower time-resolution, non-feedback measurements performed on the same myosin molecules, which are well-described by a single, slow rate of detachment (grey and black, Figure 2a).
The higher detachment rates, kf (Figure 2f, red closed symbols) and kint (green closed symbols), are more than 10-fold faster than the slowest detachment rate, ks (blue closed symbols). kf increases from 1 × 103 s−1 at low loads to up to 3–6 × 103 s−1 with increasing hindering or assisting forces. Events which detach at kf composed 20–40% (Af) of all observed interactions (Supplementary Figure 1b,d). The events which detach at the intermediate rate make up a small fraction of observed interactions (Aint = 10–20%, Figure 2—figure supplement 1) and detach at a rate, kint = 102–103 s−1. This rate has an asymmetric force dependence about zero force, increasing from −3.75 pN to +4.5 pN (Figure 2—figure supplement 2). ks, which is 6–8 s-1 at 1 μM MgATP (consistent with the measured ATP binding rate) (Liu et al., 2015; Deacon et al., 2012; Bloemink et al., 2014), increases to 10–70 s−1 in the presence of 1 mM MgATP (Figure 2—figure supplement 3). This observation is consistent with ADP release limiting the detachment rate at saturating ATP and indicates this population (As) represents motors that have completed the ATPase cycle (Figure 1a) (Liu et al., 2015; Greenberg et al., 2014; Woody et al., 2018a; Sung et al., 2015; Liu et al., 2018).
The actin attachments described by kf are much shorter than expected for molecules that undergo the full, multi-step biochemical reaction while attached to actin that ends with ATP-induced detachment (Figure 1a). These fast detaching interactions do not become faster when MgATP concentration is raised to 1 mM (Figure 2—figure supplement 2; Figure 3), further supporting the suggestion that they detach before completing the ATPase cycle. As shown below, these states do not undergo a working stroke, and are likely the short-lived states described in the physiological and biochemical literature as ‘weak-binding’ states (Brenner, 1991; Stein et al., 1984).
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.49266.011
The fraction of rapid events detaching at kf or kint (Af + Aint) increases with assisting and hindering forces (Figure 2g blue bars, Figure 2h). Although this effect may be partially due to reduced deadtimes at higher forces (x-axis intercepts in Figure 2a-e, see Materials and methods), it also may indicate that higher forces are more likely to pull motor domains off of actin before they can complete their cycle (Capitanio et al., 2012). This hypothesis is further supported by plot of deadtime-corrected fraction Af + Aint, which also increases with force (Figure 2—figure supplement 1).
The myosin working stroke is rapid and more likely for long events
Individual myosin working stroke displacements were observed shortly after the linear motion of the actin stopped when myosin attached (Figure 3a). Upon binding, 5–10 nm displacements were observed in many traces. However, displacements were not readily apparent in all traces, and were rarely seen in interactions with durations below 5 ms (Figure 3a, bottom example). These short duration interactions likely represent myosins that detach before producing a working stroke (see below).
Attempts to objectively quantify these traces using step finding algorithms, hidden Markov chains, or Bayesian non-parametric analysis proved unsatisfactory due to the amplitude and frequency characteristics of the random signal fluctuations. Therefore, ensemble averaging of individual interactions was used to determine the average time and amplitude of the myosin working stroke at forces that hindered the stroke. Interactions were aligned at their detected start times and averaged forward in time (Figure 3b,c). Event displacements were extended in time from their last interacting positions to generate groups with equal durations that could be averaged (Veigel et al., 1999).
The sizes of the average displacements were highly force dependent, decreasing with increasing force (Figure 3b,c). At forces > 1.5 pN interactions with durations < 5 ms showed suppressed displacements (Figure 3d). When short events were excluded from the ensemble averages, the average displacement increased as the minimum event duration increased from 0 to 15 ms (Figure 3e, Figure 3—figure supplement 1). Average displacement of events > 15 ms decreased linearly with increasing hindering load (Figure 3e–g). If the load-dependence of displacement for interactions > 15 ms is due to compliance within the myosin, (Dupuis et al., 1997) then we estimate the stiffness of the motor to be 2.3 pN/nm (68% CI, 2.0–2.7 pN/nm) or greater from the linear fit to the data points in Figure 3g. This value is at the upper end of previous estimates (Kaya and Higuchi, 2010; Lewalle et al., 2008).
Strikingly, ensemble averages show that the myosin working stroke occurs within 5 ms of actin binding, without a detectable time lag at forces > 1.5 pN (Figure 3f, Figure 4). In Figure 4a, the traces are normalized to the amplitude of the initial stroke for visual comparison of their rates. Estimates of the stroke rate from fits of the displacement data to single exponential rise (see Materials and methods) yielded rate constants (kstroke) of 700 to 5,250 s-1, that increase with force from 1.5 pN to 3 pN and then plateau (Figure 4c, blue). These observed rate constants are likely the sum of the rate for proceeding forward (performing the stroke, Figure 1a, k+5) and rate for reversing actomyosin attachment (Figure 1a, k-4, Supplementary file 1 Table 2). In fact, only events in which the stroke is faster than myosin detachment are detected and thus can contribute to the stroke rate in the ensemble average. Therefore, faster strokes are selected as the detachment rate kf increases with force.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.49266.014
A slower 1–1.5 nm further displacement (Figure 3b,e) observed in the ensemble averages was completed in 50–100 ms. The amplitude and time-scale of this slower displacement are consistent with the second phase of the myosin working stroke, previously associated with ADP release (Greenberg et al., 2014; Woody et al., 2018a; Capitanio et al., 2006).
Effects of Pi on the actomyosin attachments and displacement dynamics
The addition of phosphate to contracting demembranated muscle fibers results in a decrease in active force due to the redistribution of myosin states (Hibberd et al., 1985; Cooke et al., 1988; Araujo and Walker, 1996). To probe the effect of phosphate on myosin at the single molecule level, we performed UFFC experiments in the presence of 10 mM free Pi. Adding phosphate did not substantially change the rates of the kf or kint over the range of probed forces (Figure 2f, open red and green symbols); however, their combined amplitudes (Af+Aint) increased in 10 mM Pi at loads < 3 pN, suggesting an increase in the fraction of events that detach from weak-binding states (Figure 2b,c,g,h). The ratio of Af+Aint at 10 mM Pi to Af+Aint at 0 Pi shows there is up to a 50% increase in the number of these rapid detachments (Figure 2h). ks decreased about two-fold at 10 mM Pi and 1 μM MgATP (Figure 2f) in both the UFFC and non-feedback experiments, which is consistent with weak competition between Pi and ATP for binding at the active site, as previously proposed (Amrute-Nayak et al., 2008).
The presence of 10 mM free Pi did not affect the rates of the ensemble displacements for events longer than 25 ms (Figure 4b and Figure 4c (red)). This result strongly implicates a mechanism in which Pi release occurs after the working stroke, as a delay in the stroke kinetics would have been observed with 10 mM Pi if Pi release and rebinding occurred before the working stroke. (Figure 4—figure supplement 1). Comparing the rate of the stroke to the fastest rate of detachment, kf, from the event duration data (Supplementary file 1 Table 2), shows that the average stroke rate is greater than or equal to the fastest phase of detachments. This result implies that molecules in the earliest state resolved in our experiment either perform a stroke or detach, and that Pi release (or any other substantial intermediate step) does not occur between attachment and the stroke.
Many individual displacement records acquired in the presence of 10 mM Pi showed reversals in the working stroke to the pre-stroke level (Figure 5a). However, the individual recordings are variable, and reversals cannot be reliably identified among the Brownian fluctuations. Thus, the dynamic features of the displacements were objectively quantified via ensemble averaging (Figures 4 and 5 b-c).
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.49266.019
The average displacement traces decreased after the peak of the initial working stroke in the presence of 10 mM Pi (Figure 5c, dark colors). We tabulated the total displacement, the size of any declines occurring after the initial positive displacement, and the time to reach the minimum of this dip (Figure 5d–f) (see Materials and methods for definition of the respective quantities). Total displacement was reduced (Figure 5d), and the size of the dip was increased by Pi at all forces above 1.5 pN hindering load (Figure 5e). The changes in total displacement and the amplitude of dip are both greater at higher hindering loads, and the time to reach the minimum of this dip is considerably increased by the presence of 10 mM Pi (Figure 5f). These results are consistent with a model in which Pi is released after the working stroke and then may rebind and cause reversal of the stroke (Figure 4—figure supplement 1). We performed simulations of individual actomyosin interactions using a kinetic model with these features (Figure 5—figure supplement 1). Ensemble averages of these interactions recapitulate the observed effects of Pi (Figure 5g,h).
Discussion
The time resolution afforded by the UFFC system has allowed us to gain novel insights into dynamics of the cardiac myosin initial states following actin attachment and the working stroke. We resolved a short-lived state that is capable of bearing loads up to 4.5 pN and obtained the first estimates of the delay between myosin attachment to actin and the working stroke as a function of applied load from single molecule measurements. We also observed the effects of free phosphate on these events. The results allow us to characterize the order and properties of the rapid processes which follow cardiac myosin association with actin and lead to force generation. Based on the experimental results discussed below, we propose a model (Figure 6) in which an early actomyosin state that is short-lived but nevertheless is capable of bearing load leads promptly to the working stroke and concomitant strong binding, which is then followed by slower Pi release.
Short-lived states are Stereo-Specific, Open-cleft, Pre-stroke, and Pi-bound
Our UFFC experiments reveal short-lived interactions which detach before myosin completes its biochemical cycle via ATP binding. Although we have used a two-headed, HMM construct in our experiments to emulate biological structure, we did not observe changes in the stiffness (measured by the bead covariance signal) during attachments, or lifetime distributions consistent with paired dissociation events that might signal double-headed attachment. Thus, the data are interpreted as resulting from single myosin heads. The short-lived state detected in this study can stop the motion of an actin filament and maintain forces of up to 4.5 pN for hundreds of microseconds. In earlier publications, a weak-binding state of myosin has been described in several ways, and we considered whether the observed short-lived interactions may be in such a weak-binding state.
An electrostatic, non-stereospecific interaction between actin and myosin has been called a weak-binding state in which the myosin head experiences significant orientation fluctuations relative to actin. (Thomas et al., 1995; Thomas et al., 2002; Ostap et al., 1995; Bershitsky et al., 1997; Ferenczi et al., 2005; Wu et al., 2010; Kraft et al., 2005; Lowy et al., 1991; Cooke et al., 1982) The dissociation rate of the short-lived state found here (kf) is substantially more force dependent under hindering loads as compared to assisting loads (Figure 2f, Figure 2—figure supplement 2) and as mentioned above it can briefly support loads. These characteristics suggest that the actomyosin bond in this state is stereo-specific and not a highly disordered interaction.
Weak-binding has also been used to describe a conformation of the myosin head in which a cleft in the actin binding interface is open and is usually associated with ATP or ADP-Pi bound myosin (Sweeney and Houdusse, 2010; Dominguez et al., 1998). This cleft-open state has been proposed to be the initial stereospecific interaction between actin and myosin which can rapidly dissociate from actin unless a relatively slow isomerization leads to ‘strong binding’ (Thomas et al., 2002; Bershitsky et al., 1997; Ferenczi et al., 2005; Kraft et al., 2005; Brenner et al., 2005; Beausang et al., 2013; Brenner et al., 1982). The short-lived states we observed are estimated to dissociate at ~500 s−1 in unloaded conditions, fairly close to the measured ~950–1450 s−1 dissociation rate of actomyosin after actin cleft opening induced by ATP binding (Liu et al., 2015; Deacon et al., 2012; Bloemink et al., 2014). Comparing this to the <0.1 s−1 dissociation rates measured for known cleft-closed states (Woody et al., 2018a), our short-lived state’s lifetime suggests that myosin is in an open-cleft conformation.
Actin-activated Pi release is dependent on a stereo-specific interaction between actin and myosin. If Pi were released prior to the observed short-lived state, there would have to exist an additional stereo-specific state that is too rapid to observe. Based on the estimated rate of dissociation from an open-cleft state of 950 s−1 given above, such a state is highly unlikely. Thus, our short-lived state is likely to be a stereospecific, actin-cleft-open state with Pi bound.
The stroke and cleft closure occur before Pi is released
The model in Figure 6 which places the working stroke before Pi release is supported by the rapid observed rate of the stroke (kstroke) and the similarity between kstroke and kf (Figure 4), as well as the observations of increased and delayed reversals in the presence of free Pi (Figure 5), as discussed further below.
Ensemble averages of actin displacement during actomyosin events show that there is no appreciable working stroke during the short-lived interactions when time resolution is optimal (forces > 1.5 pN, Figure 3d). The working stroke appears to immediately follow this short-lived state, as the observed rate of the working stroke (kstroke, diamonds in Figure 4c) is greater than or equal to the rate of detachment from the short-lived state (kf, circles in Figure 4c). kstroke, surprisingly, shows an increased rate with larger forces in the direction that should hinder the working stroke. If the stroke immediately follows the short-lived state, this result can be readily explained by the slip-bond behavior of kf, as the observed kstroke (≈ k+5 + k-4, Figure 1) should increase when k-4 increases. If a significant intermediate process, such as Pi release, occurred between entering the short-lived state and the stroke, it would have to occur at >10,000 s−1 to obtain the observed value of kstroke (Figure 5—figure supplement 2). In addition, if Pi release and rebinding were to occur before the stroke, kstroke would be slowed by Pi rebinding when it is present in solution (Figure 4—figure supplement 1). The lack of effect from 10 mM Pi on kstroke (Figure 4c, compare open and closed diamonds) is thus further evidence that Pi release and rebinding do not occur before the working stroke.
Biochemical estimates of the Pi release rate for cardiac myosin in solution range from 11.6 to 17 s−1, (Rohde et al., 2017; Liu et al., 2015) which is much slower than the observed stroke rate. The biochemical measurements estimate the microscopic rate of Pi release, while our measurements are observed rates including competing processes (such as detachment) from the same states. If Pi release were to happen before, or simultaneously with, the stroke at a rate of 17 s−1, less than 0.4–4% (high load to zero load) of the molecules occupying the short-lived state would produce a stroke, the rest detaching at ≥500–5000 s−1. The event duration data show that >10–40% (high - low load) of interactions proceed through the entire biochemical cycle (Figure 2g, Figure 2—figure supplement 1). This comparison thus again supports models in which Pi is released after the working stroke and prevents reversal of the stroke and cleft re-opening unless Pi rebinds.
The stroke can be reversed and is reversed more often with Pi present
It has been previously suggested that the myosin actin-binding cleft closure, and a corresponding increase in actin affinity, is tightly coupled to the working stroke (Sweeney and Houdusse, 2010; Málnási-Csizmadia and Kovács, 2010; Xu and Root, 2000), and our data support this hypothesis. If the binding interface cleft remained open after the stroke, a significant fraction of short-lived events would be expected to show a displacement of actin. Under conditions when the time resolution was optimal (forces ≥ 3 pN), negligible displacement was observed in short-lived events.
The data provide evidence that the working stroke and cleft closure are both reversible. The kint detachment rate without added Pi showed an asymmetric force dependence about zero load (Figure 2f, green, Figure 2—figure supplement 2), consistent with a force-dependent rate for the reversal of the working stroke followed immediately by detachment from the short-lived, pre-stroke state. Evidence of these reversals appeared in the ensemble averages of displacement both with and without free Pi (Figure 3c,e and Figure 5b,c) as a dip in the displacement that occurred within a few milliseconds after the initial stroke.
An increase in the fraction of rapidly detaching events was observed with 10 mM added Pi at low forces (Figure 2g,h). This result is consistent with previous single molecule observations (Baker et al., 2002; Debold et al., 2013), and with the scheme of Figure 6 in which Pi rebinding allows more stroke reversals which can lead to premature detachments. At higher forces, it is likely that the increased fraction of direct detachments from the short-lived state (at rate kf) may mask the effects of Pi on attachment durations. Ensemble averages of the long events (>25 ms) show that added Pi leads to a larger and more delayed dip in the displacement (Figure 5c-f), consistent with the scheme of Figure 6 in which Pi rebinding allows reversals to occur both more often and later in the attachment (Figure 4—figure supplement 1). The decreased displacement in the presence of 10 mM Pi compared to control (Figure 5d) is consistent with myosin detaching from a pre-stroke conformation after a delayed reversal of the working stroke.
It should be noted that for this model to describe the ensemble average data (Figure 5g,h), the lifetime of the pre-stroke state following a reversal was approximately 1 ms under load (Figure 5—figure supplement 1). This is longer than the estimated 300 μs lifetime of the initial pre-stroke state under load. This might suggest that after Pi is released, the Pi rebinding does not cause the actin cleft to reopen fully, even if it does allow power stroke reversal. Additional studies and modeling will be required to test this speculation, but it has been suggested previously that Pi rebinding may induce transitions off of the canonical pathway (Debold et al., 2013).
The model in Figure 6 can explain all our observations. In this model, electrostatically mediated interactions between actin and myosin transition to the short-lived, stereo-specific state detected in our study. In this state, myosin can briefly bear force with an open actin binding cleft and Pi bound. The working stroke of myosin occurs directly from this state and the closure of the cleft is tightly coupled to the stroke. Both the stroke and cleft-closure can be reversed under zero load, with a rate (kint) that increases with hindering load. Phosphate release occurs after the stroke and rebinding of phosphate allows reversals of the stroke.
Relation to previous work (Walker et al., 1992)
The scheme of Figure 6 is fully consistent with a previously proposed model developed from muscle fiber experiments for skeletal myosin II (Dantzig et al., 1992; Walker et al., 1992). That model included a force bearing state followed by reversible force generation (the working stroke) and then phosphate release (Dantzig et al., 1992). Similar studies on cardiac fibers showed this model also is applicable to cardiac myosin (Araujo and Walker, 1996). More recent studies of the rate of the working stroke in skeletal (Muretta et al., 2015) and cardiac (Rohde et al., 2017) myosin II, and myosin V (Trivedi et al., 2015) using FRET probes and time-resolved fluorescence have found evidence of a conformational change in the motor domain, likely a rotation of the lever arm, that precedes phosphate release. The data from previous experiments on fast-skeletal myosin with the UFFC system support that the model of Figure 6 also applies to this myosin isoform (Capitanio et al., 2012). The present data provide the most direct evidence for this ordering of events after actin attachment in cardiac muscle myosin.
High resolution crystal structures of myosin have been used to suggest that phosphate must be released before the stroke occurs (Llinas et al., 2015; Sweeney and Houdusse, 2010). In addition, some models of contraction developed to explain a wide variety of muscle fiber data have suggested that phosphate release and the working stroke are not tightly coupled (Caremani et al., 2015; Caremani et al., 2013). While our data cannot exclude that Pi may occasionally be released before the stroke, if this occurred more than 20% of the time, there would be a discernable delay in the stroke rate in the presence of Pi which was not observed (Figure 4). The possibility remains that other myosin isoforms which exhibit different kinetic properties than muscle myosin, such as non-muscle myosins V and VI, may proceed through actin binding, the working stroke, and phosphate release with different dynamics or even different sequences (Llinas et al., 2015; Sweeney and Houdusse, 2010; Rosenfeld and Lee Sweeney, 2004). Specific experiments on those systems are needed to clarify the mechanisms.
Materials and methods
Protein purification and source
Request a detailed protocolA heavy meromyosin (HMM) construct of human β-cardiac myosin (MHY7) was expressed in C2C12 myoblasts and purified as previously described (Winkelmann et al., 2015). The HMM protein has 1146 residues that include residues 1–1138 of the human MYH7 gene and a FLAG tag on the C-terminus (res. 1139–1146). The sequence of the HMM preparations used in this study were confirmed by LC/MS/MS of protein digests. Bound light chains are those that are constitutively expressed in the C2C12 cells (MLC1/MLC3 and rLC2).
Experimental protocol
Request a detailed protocolExperiments were conducted in flow cells constructed and prepared as previously described (Woody et al., 2018a). The assay solutions contained 60 mM MOPs, 10 mM DTT, 1 mM EGTA, 1 mM MgCl2, and 1 mg/mL BSA. The calculated ionic strength of the experimental solutions was held constant across ATP and Pi conditions. The values of the components and total ionic strength of the solutions for the various Pi and MgATP conditions are given in Supplementary file 1 Table 3.
Myosin was added to nitrocellulose-coated chambers and allowed to attach nonspecifically to 2.5 μm silica pedestal beads on the surface and then the chambers were blocked with 1 mg/mL BSA as previously described (Woody et al., 2018a). The concentration of myosin (~0.2 ng/mL)) led to about 1 out of 10 bead locations showing interactions between myosin and actin. The final assay solution contained pM concentrations of 10–15% biotinylated actin (Cytoskeleton Inc,) from rabbit skeletal muscle stabilized by rhodamine phalloidin, followed by ~3 microliters of the final assay solution containing 5 ng mL−1 500 nm diameter beads (PolySciences Inc,), coated nonspecifically with neutravidin (Woody et al., 2018a) (ThermoFisher Inc).
Experiments were performed with trap stiffness 0.06–0.08 pN/nm and 4–5.5 pN pretension applied to the actin dumbbells. When a myosin molecule was found on a pedestal bead, a stage stabilization system was engaged (Woody et al., 2018b; Capitanio et al., 2005) and actomyosin events were recorded without force-clamp feedback to determine the size of the myosin working stroke and directionality of the filament (Woody et al., 2018a). The ultra-fast force clamp was then engaged, and 5–10 min of data were recorded for each applied force. Occasional filament breakage prevented collection at all forces (1.5–4.5 pN) for some dumbbells.
Optical trapping implementation
Request a detailed protocolThe optical trapping setup used in these experiments was described previously in detail, (Woody et al., 2017) including optimizations to increase the update response time of the feedback loop to less than 8 μs. Briefly, polarization-split 1064 nm beams are steered by two 1D electro-optical deflectors (EODs, Conoptics, Inc) into a 60x water immersion objective (Nikon). The laser light collected from the chamber of an oil immersion condenser is projected onto two quadrant photodiodes conjugate to the back focal plane of the condenser for direct force detection (Gittes and Schmidt, 1998). The feedback loop and data acquisition used a National Instruments Field Programmable Gate Array controlled by custom LabVIEW virtual instruments.
Experiments were conducted otherwise as previously described (Capitanio et al., 2012; Woody et al., 2018b; Gardini et al., 2018) with applied forces of 1.5–4.5 pN in directions alternating toward each end of the actin filament. The trap moved under the applied force for 200 nm before the applied force was switched to the opposite direction and the cycle of motions back and forth repeated. To eliminate effects that tension on the actin filament might have on myosin interaction, the forces were distributed between the two beads so that the pretension on the filament was constant even as the applied force changed magnitude (Figure 5—figure supplement 3). The force signals and feedback output signals to the EODs were digitized at 250 kHz through a 125 kHz anti-aliasing filter in the sum and difference amplifiers for the x- and y-direction difference signals from the photodiode currents. The positions of the traps were determined from the EOD driving signals and optical calibration measurements of corresponding bead displacements.
We added a novel real-time drift and slope correction (DSC) system to our UFFC instrument to improve stability and accuracy of the force signal as a result of thermal changes and other variations. This system is described in the below and Figure 5—figure supplement 3.
Real-time drift and slope correction (DSC)
Request a detailed protocolRelative to acousto-optic deflectors (AODs) more commonly used in dynamic optical trap instruments, the electro-optic deflectors (EODs) used here provide a much smoother angular deflection of the IR beams (constant slope of deflection vs. input voltage) at the expense of maximum deflection angle (Woody et al., 2018b; Valentine et al., 2008). The EOD crystals are more narrow (2 x 2 mm) and much longer (100 mm) than typical AODs (~4 x 4 x 10 mm) and it is essential to keep the laser beams wholly within the crystals, so alignment is more critical and was checked and/or realigned daily.
Even with regular alignment, thermal drifts and or air currents produced instrumental 0.25–2 pN fluctuations over an hour period of recording. To obtain accurate and reliable force signals, we utilized an real time drift and slope correction (DSC) system which periodically pauses the feedback system during an experiment to probe the true ‘zero load’ (F0) signal values and correct for the drift. While the feedback system is briefly paused, the trap positions are not updated and (when myosin is not interacting with actin) the only force on the beads is from the pretension applied to the filament. The force signals observed during each pause can be used to reset the zero-force voltage and calculate the setpoint in the feedback loop to the recorded F0 level plus the applied force.
We implemented this system on the FPGA feedback controller such that every 10 ms the feedback loop pauses twice, first when the trap position reaches the top of the set excursion distance, and again when the trap position is at the bottom of the excursion (Figure 5—figure supplement 3c). Each pause lasts 2 ms, and the force signals from the last 1.5 ms of the pause are averaged for each bead to provide an estimate of the F0 signals. These individual values are averaged over 20 pauses to calculate the setpoint during the experiment. Either the top or bottom pause allows correction for the slow thermal drifts. Additionally, by using both the top and bottom F0 values and the known excursion distance, the value of the EOD slope is be calculated (Figure 5—figure supplement 3).
This slope value can then be used in the FPGA to update the setpoint value as the trap positions change in order to account for the slope in the force signals.
Examples of the drift offset and the calculated slope values over a 30 s trace are shown in Figure 5—figure supplement 3d.
When the data are analyzed, the recorded force signals can be transformed into the actual force signals by using a similar algorithm in the analysis software that was used in real time on the FPGA (Figure 5—figure supplement 3e). When detecting and characterizing binding events, the forces and positions of the traps when the feedback was paused are ignored. This DSC system markedly improved reliability and consistency within and across experiments, reduced the width of the velocity distributions, lowered effective deadtimes, and improved data reliability at lower forces.
Data analysis
Event detection
Request a detailed protocolEvents were detected using a modified version of software written in Matlab (Mathworks, Natick, MA) used previously for UFFC experiments (Capitanio et al., 2012; Gardini et al., 2018). Briefly, the velocity of the traps is calculated and smoothed using a 160–800 μs width Gaussian filter, and a threshold is determined for distinguishing periods when the trap motion has stopped that balances false binding and false unbinding events. The amount of smoothing is set so that the number of false events is estimated to be <1% of all detected events based on statistical estimates (Colquhoun and Sigworth, 1995). The smoothing varies with the applied force but is constant across experimental conditions. Because higher loads cause faster motion of the actin, less smoothing is required and shorter events can be detected at higher loads. Small corrections are applied to the start- and end-times of the events to account for delay based on the Gaussian filter, as previously described (Capitanio et al., 2012). The estimated uncertainty of the localization of binding time (σt) can be calculated and ranged from 500 μs at 1.5 pN to <60 μs at 4.5 pN applied load. From each data trace, a theoretical minimum detectable event length (deadtime) was calculated based on the velocity distributions and the threshold value as previously described (Capitanio et al., 2012). Events were designated as occurring either during hindering or assisting load based on the direction of actin motion when binding occurred and the polarity of the actin as determined by initial non-feedback experiments.
Duration analysis
Request a detailed protocolFor a given force, all events which were shorter than the largest calculated deadtime for traces recorded at that force were excluded from the analysis. The deadtimes, set to be constant across different biochemical conditions, were 2.7, 1.25, 0.73, 0.65, and 0.5 ms for forces 1.5, 2.25, 3, 3.75, and 4.5 pN respectively. Events from molecules at the same magnitudes and directions of force were pooled for each given condition (e.g. forces, concentrations of MgATP and free Pi). Parameters for deadtime corrected single, double, and triple exponential distributions were estimated in MEMLET (Woody et al., 2016) for each data set and the log-likelihood ratio test was used with a p-value cutoff of 0.05 to determine if the triple and/or double exponential distributions were statistically justified. Each molecule contributed equally to the fitted values as described previously (Woody et al., 2018a), although this did not lead to considerable differences compared to the fits weighted by the number of interactions. The observed amplitudes of each phase are reported as well as the deadtime corrected amplitudes (Figure 2—figure supplement 1), which consider how many events are likely missed due to the experimental deadtime and assuming exponential duration components (Woody et al., 2016).
Ensemble averaging
Request a detailed protocolEnsemble averages were carried out as have been described previously (Woody et al., 2018a). The beginnings of the events were aligned based on their detected time of binding from the crossing of the smoothed velocity trace through the detection threshold as previously described (Capitanio et al., 2012). The trap position for the leading bead (the bead with the higher magnitude of force applied) was taken as the position of the actin filament, as there is less influence of non-linear end compliance on this more highly loaded bead. The extension point was taken to be the displacement at time σt (uncertainty of binding time) before the detected end of the event, meaning this value was used in the averages for all time points longer than the actual event duration. The value of σt depends on the smoothing of the velocity for event detection, and thus was smaller (more accurate) for larger forces. The average value of σt varied from 0.5 ms for 1.5 pN to ~72 μs for 4.5 pN. Ensembles were weighted such that each molecule contributed equally to the final average.
Events longer than 15 ms were analyzed together for the 0 Pi conditions in Figure 3 to isolate events which went through the entire biochemical cycle, including ATP binding to terminate the actomyosin interaction, as the total observed step did not change >10% by excluding any additional events longer than 15 ms (Figure 3—figure supplement 1a). For the 10 mM Pi data and for analysis of the 0 Pi data which was directly compared to that with 10 mM Pi (Figures 4 and 5), events longer than 25 ms were included because the total working stroke size continued to increase as events between 15 ms and 25 ms were excluded (Figure 3—figure supplement 1b), which could be due to reversals and detachments occurring later in an interaction when Pi is present.
Quantifications of ensemble averages
Request a detailed protocolQuantification of the ensemble average was done programmatically, using the methods and parameters reported below. The minimum position of the ensemble average within the first millisecond of detected actin binding was set to be zero displacement. The size of the initial displacement was calculated from the raw (unfiltered) ensemble averages by finding the time (tinit) when the ensemble average reached its maximum within the first 1.2 ms following the minimum position, with the exception of the 1.5 pN data, which were searched over the first 5 ms. The displacement from a window of 200 µs centered on tinit was averaged to determine the maximum initial position. The minimum and maximum positions were used to normalize the displacements to range from 0 to 1, as plotted in Figure 4a,b. The normalized displacement data between the minimum position and tinit was used to fit a single exponential rise with an amplitude of 1 to determine the rate of the working stroke as reported in Figure 4c.
For quantification of the dip size and time, the initial stroke size was determined by averaging 100 μs on either side of tinit. The location and position of the dip were determined by taking a 1 ms moving average of the ensemble and finding the minimum point between tinit and 15 ms after the detected binding time. The size and timing of the dip are reported as the difference between the time and displacement of initial stroke and the minimum dip respectively. The total displacement is quantified as an average over the last 200 μs of the ensemble average.
Modeling
Request a detailed protocolSimulations of individual traces were performed using a Monte-Carlo based method (Gillespie, 1977) using the rates and transitions given in Figure 5—figure supplement 1. For each simulated state (i), and for each transition out of that state (), times, , were randomly selected from an exponential distribution based on the transition rate (). The shortest was used as the lifetime of the state i and the transition associated with that time was used to assign the next state (j). This method leads to the expected results that the flux through a transition () will be the rate of that transition divided by the sum of all rates exiting the current state . The effective rate of the transition will be . Once all states were assigned for a given simulated interaction (the last state is always a detached state), a position value was assigned to each state based on the mechanical properties of the states. The simulated position values were drawn for each interaction from a Gaussian distribution with a 0.1 nm standard deviation from their set mean value. The position trace before the initial interaction was a constant slope with a velocity similar to what is observed in the data for the simulated applied force. Upon initial binding, an exponential function was used to simulated stretching of the myosin under the applied load using a linear stiffness of 2 pN/nm. 3000 interactions were simulated, and the interactions lasting longer than a 25 ms cutoff were included for averaging with the same routine used for the data analysis.
Data availability
The primary optical trapping data files are > 1TB, and are available upon request. Source data files have been provided for Figures 1-5.
References
-
Inorganic phosphate binds to the empty nucleotide binding pocket of conventional myosin IIJournal of Biological Chemistry 283:3773–3781.https://doi.org/10.1074/jbc.M706779200
-
The hypertrophic cardiomyopathy myosin mutation R453C alters ATP binding and hydrolysis of human cardiac β-MyosinJournal of Biological Chemistry 289:5158–5167.https://doi.org/10.1074/jbc.M113.511204
-
BookDriving Filament Sliding: Weak binding cross-bridge states, strong binding cross-bridge states, and the power strokeIn: Sugi H, editors. Sliding Filament Mechanism in Muscle Contraction, 565. Springer. pp. 75–91.https://doi.org/10.1007/0-387-24990-7_7
-
Position control and optical manipulation for nanotechnology applicationsThe European Physical Journal B 46:1–8.https://doi.org/10.1140/epjb/e2005-00232-2
-
Force and number of myosin motors during muscle shortening and the coupling with the release of the ATP hydrolysis productsThe Journal of Physiology 593:3313–3332.https://doi.org/10.1113/JP270265
-
BookFitting and Statistical Analysis of Single-Channel RecordsIn: Sakmann B, Neher E, editors. Single-Channel Recording. Springer. pp. 483–587.https://doi.org/10.1007/978-1-4419-1229-9_19
-
The inhibition of rabbit skeletal muscle contraction by hydrogen ions and phosphateThe Journal of Physiology 395:77–97.https://doi.org/10.1113/jphysiol.1988.sp016909
-
Identification of functional differences between recombinant human α and β cardiac myosin motorsCellular and Molecular Life Sciences 69:2261–2277.https://doi.org/10.1007/s00018-012-0927-3
-
Actin filament mechanics in the laser trapJournal of Muscle Research and Cell Motility 18:17–30.https://doi.org/10.1023/A:1018672631256
-
High-Speed optical tweezers for the study of single molecular motorsMethods Mol. Biol. Clifton NJ 1805:151–184.https://doi.org/10.1007/978-1-4939-8556-2_9
-
Exact stochastic simulation of coupled chemical reactionsThe Journal of Physical Chemistry 81:2340–2361.https://doi.org/10.1021/j100540a008
-
BookMyosin Motors: Kinetics of MyosinIn: Egelman E. H, editors. Comprehensive Biophysics, 4. Elsevier. pp. 151–169.https://doi.org/10.1016/B978-0-12-374920-8.00411-2
-
Relationships between chemical and mechanical events during muscular contractionAnnual Review of Biophysics and Biophysical Chemistry 15:119–161.https://doi.org/10.1146/annurev.bb.15.060186.001003
-
Single-Molecule measurement of the stiffness of the rigor myosin headBiophysical Journal 94:2160–2169.https://doi.org/10.1529/biophysj.107.119396
-
Controlling load-dependent kinetics of β-cardiac myosin at the single-molecule levelNature Structural & Molecular Biology 25:505–514.https://doi.org/10.1038/s41594-018-0069-x
-
How actin initiates the motor activity of myosinDevelopmental Cell 33:401–412.https://doi.org/10.1016/j.devcel.2015.03.025
-
Emerging complex pathways of the actomyosin powerstrokeTrends in Biochemical Sciences 35:684–690.https://doi.org/10.1016/j.tibs.2010.07.012
-
A model of myosin V processivityJournal of Biological Chemistry 279:40100–40111.https://doi.org/10.1074/jbc.M402583200
-
Rate-limiting step in the actomyosin adenosine triphosphatase cycleBiochemistry 23:1555–1563.https://doi.org/10.1021/bi00302a033
-
Structural and functional insights into the myosin motor mechanismAnnual Review of Biophysics 39:539–557.https://doi.org/10.1146/annurev.biophys.050708.133751
-
Force generation in single conventional actomyosin complexes under high dynamic loadBiophysical Journal 90:1295–1307.https://doi.org/10.1529/biophysj.105.068429
-
The mechanism of force generation in myosin: a disorder-to-order transition, coupled to internal structural changesBiophysical Journal 68:135–141.
-
Changes in actin and myosin structural dynamics due to their weak and strong interactionsResults and Problems in Cell Differentiation 36:7–19.https://doi.org/10.1007/978-3-540-46558-4_2
-
Precision steering of an optical trap by electro-optic deflectionOptics Letters 33:599–601.https://doi.org/10.1364/OL.33.000599
-
Effects of Ca2+ on the kinetics of phosphate release in skeletal muscleThe Journal of Biological Chemistry 267:2459–2466.
-
ConferenceAn ultra-fast EOD-based force-clamp detects rapid biomechanical transitions. in 10347SPIE Nanoscience + Engineering.https://doi.org/10.1117/12.2277130
-
Conformational selection during weak binding at the actin and myosin interfaceBiophysical Journal 79:1498–1510.https://doi.org/10.1016/S0006-3495(00)76401-0
Article and author information
Author details
Funding
National Institutes of Health (P01-GM087253)
- E Michael Ostap
- Yale E Goldman
National Institutes of Health (R01-HL133863)
- Donald Winkelmann
- E Michael Ostap
National Science Foundation (CMMI: 15-48571)
- E Michael Ostap
- Yale E Goldman
National Science Foundation (NSF GRFP)
- Michael Woody
National Institutes of Health (R35-GM118139)
- Yale E Goldman
Horizon 2020 Framework Programme (654148)
- Marco Capitanio
Ente Cassa di Risparmio di Firenze
- Marco Capitanio
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2019, Woody et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,937
- views
-
- 341
- downloads
-
- 73
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Lipopolysaccharides (LPS) confer resistance against harsh conditions, including antibiotics, in Gram-negative bacteria. The lipopolysaccharide transport (Lpt) complex, consisting of seven proteins (A-G), exports LPS across the cellular envelope. LptB2FG forms an ATP-binding cassette transporter that transfers LPS to LptC. How LptB2FG couples ATP binding and hydrolysis with LPS transport to LptC remains unclear. We observed the conformational heterogeneity of LptB2FG and LptB2FGC in micelles and/or proteoliposomes using pulsed dipolar electron spin resonance spectroscopy. Additionally, we monitored LPS binding and release using laser-induced liquid bead ion desorption mass spectrometry. The β-jellyroll domain of LptF stably interacts with the LptG and LptC β-jellyrolls in both the apo and vanadate-trapped states. ATP binding at the cytoplasmic side is allosterically coupled to the selective opening of the periplasmic LptF β-jellyroll domain. In LptB2FG, ATP binding closes the nucleotide binding domains, causing a collapse of the first lateral gate as observed in structures. However, the second lateral gate, which forms the putative entry site for LPS, exhibits a heterogeneous conformation. LptC binding limits the flexibility of this gate to two conformations, likely representing the helix of LptC as either released from or inserted into the transmembrane domains. Our results reveal the regulation of the LPS entry gate through the dynamic behavior of the LptC transmembrane helix, while its β-jellyroll domain is anchored in the periplasm. This, combined with long-range ATP-dependent allosteric gating of the LptF β-jellyroll domain, may ensure efficient and unidirectional transport of LPS across the periplasm.
-
- Biochemistry and Chemical Biology
- Cancer Biology
Many cells in high glucose repress mitochondrial respiration, as observed in the Crabtree and Warburg effects. Our understanding of biochemical constraints for mitochondrial activation is limited. Using a Saccharomyces cerevisiae screen, we identified the conserved deubiquitinase Ubp3 (Usp10), as necessary for mitochondrial repression. Ubp3 mutants have increased mitochondrial activity despite abundant glucose, along with decreased glycolytic enzymes, and a rewired glucose metabolic network with increased trehalose production. Utilizing ∆ubp3 cells, along with orthogonal approaches, we establish that the high glycolytic flux in glucose continuously consumes free Pi. This restricts mitochondrial access to inorganic phosphate (Pi), and prevents mitochondrial activation. Contrastingly, rewired glucose metabolism with enhanced trehalose production and reduced GAPDH (as in ∆ubp3 cells) restores Pi. This collectively results in increased mitochondrial Pi and derepression, while restricting mitochondrial Pi transport prevents activation. We therefore suggest that glycolytic flux-dependent intracellular Pi budgeting is a key constraint for mitochondrial repression.