Actomyosin-ADP states, interhead cooperativity, and the force-velocity relation of skeletal muscle

Insights into Muscle Contraction Derived from the Effects of Small-Molecular Actomyosin-Modulating Compounds

Bottom-up mechanokinetic models predict ensemble function of actin and myosin based on parameter values derived from studies using isolated proteins. To be generally useful, e.g., to analyze disease effects, such models must also be able to predict ensemble function when actomyosin interaction kinetics are modified differently from normal. Here, we test this capability for a model recently shown to predict several physiological phenomena along with the effects of the small molecular compound blebbistatin. We demonstrate that this model also qualitatively predicts effects of other well-characterized drugs as well as varied concentrations of MgATP. However, the effects of one compound, amrinone, are not well accounted for quantitatively. We therefore systematically varied key model parameters to address this issue, leading to the increased amplitude of the second sub-stroke of the power stroke from 1 nm to 2.2 nm, an unchanged first sub-stroke (5.3−5.5 nm), and an effective cross-bridge attachment rate that more than doubled. In addition to better accounting for the effects of amrinone, the modified model also accounts well for normal physiological ensemble function. Moreover, a Monte Carlo simulation-based version of the model was used to evaluate force−velocity data from small myosin ensembles. We discuss our findings in relation to key aspects of actin−myosin operation mechanisms causing a non-hyperbolic shape of the force−velocity relationship at high loads. We also discuss remaining limitations of the model, including uncertainty of whether the cross-bridge elasticity is linear or not, the capability to account for contractile properties of very small actomyosin ensembles (<20 myosin heads), and the mechanism for requirements of a higher cross-bridge attachment rate during shortening compared to during isometric contraction.

Sensitivity of muscle force response of a two-state cross-bridge model to variations in model parameters

Muscle models based on the cross-bridge theory (Huxley-type models) are frequently used to calculate muscle forces for different contractile conditions. Dynamic and nonlinear characteristics of muscle forces produced during isometric, concentric, and eccentric contractions can be represented to a limited extent by using cross-bridge models. Cross-bridge models use various parameters to simulate force responses. However, there remains uncertainty as to the effect of changes in model parameters on force responses in Huxley-type models. In this study, we aimed to analyze the sensitivity of force response to changes in model parameters in Huxley-type models. A two-state Huxley model was used to determine the cross-bridge attachment distributions and forces for shortening and lengthening contractions. Sensitivity of muscle force to changes in attachment rate, detachment rate, and cross-bridge binding distance was examined within a range of ±20% of the nominal value using Monte Carlo simulations. Changes in the detachment rate influenced the predicted muscle forces the most for lengthening contractions, while changes in attachment rate and binding distance affected forces the most for shortening contractions. These results show once more the asymmetry between shortening and lengthening contractions and the difficulty in using a single cross-bridge model to predict forces during shortening and elongation accurately.

Multistep orthophosphate release tunes actomyosin energy transduction

Muscle contraction and a range of critical cellular functions rely on force-producing interactions between myosin motors and actin filaments, powered by turnover of adenosine triphosphate (ATP). The relationship between release of the ATP hydrolysis product ortophosphate (Pi) from the myosin active site and the force-generating structural change, the power-stroke, remains enigmatic despite its central role in energy transduction. Here, we present a model with multistep Pi-release that unifies current conflicting views while also revealing additional complexities of potential functional importance. The model is based on our evidence from kinetics, molecular modelling and single molecule fluorescence studies of Pi binding outside the active site. It is also consistent with high-speed atomic force microscopy movies of single myosin II molecules without Pi at the active site, showing consecutive snapshots of pre- and post-power stroke conformations. In addition to revealing critical features of energy transduction by actomyosin, the results suggest enzymatic mechanisms of potentially general relevance.

Release of the ATP hydrolysis product orthophosphate (Pi) from the myosin active site is central in force generation but is poorly understood. Here, Moretto et al. present evidence for multistep Pi-release reconciling apparently contradictory results.

Critical Evaluation of Current Hypotheses for the Pathogenesis of Hypertrophic Cardiomyopathy

Hereditary hypertrophic cardiomyopathy (HCM), due to mutations in sarcomere proteins, occurs in more than 1/500 individuals and is the leading cause of sudden cardiac death in young people. The clinical course exhibits appreciable variability. However, typically, heart morphology and function are normal at birth, with pathological remodeling developing over years to decades, leading to a phenotype characterized by asymmetric ventricular hypertrophy, scattered fibrosis and myofibrillar/cellular disarray with ultimate mechanical heart failure and/or severe arrhythmias. The identity of the primary mutation-induced changes in sarcomere function and how they trigger debilitating remodeling are poorly understood. Support for the importance of mutation-induced hypercontractility, e.g., increased calcium sensitivity and/or increased power output, has been strengthened in recent years. However, other ideas that mutation-induced hypocontractility or non-uniformities with contractile instabilities, instead, constitute primary triggers cannot yet be discarded. Here, we review evidence for and criticism against the mentioned hypotheses. In this process, we find support for previous ideas that inefficient energy usage and a blunted Frank-Starling mechanism have central roles in pathogenesis, although presumably representing effects secondary to the primary mutation-induced changes. While first trying to reconcile apparently diverging evidence for the different hypotheses in one unified model, we also identify key remaining questions and suggest how experimental systems that are built around isolated primarily expressed proteins could be useful.

The Synergic Role of Actomyosin Architecture and Biased Detachment in Muscle Energetics: Insights in Cross Bridge Mechanism Beyond the Lever-Arm Swing

Muscle energetics reflects the ability of myosin motors to convert chemical energy into mechanical energy. How this process takes place remains one of the most elusive questions in the field. Here, we combined experimental measurements of in vitro sliding velocity based on DNA-origami built filaments carrying myosins with different lever arm length and Monte Carlo simulations based on a model which accounts for three basic components: (i) the geometrical hindrance, (ii) the mechano-sensing mechanism, and (iii) the biased kinetics for stretched or compressed motors. The model simulations showed that the geometrical hindrance due to acto-myosin spatial mismatching and the preferential detachment of compressed motors are synergic in generating the rapid increase in the ATP-ase rate from isometric to moderate velocities of contraction, thus acting as an energy-conservation strategy in muscle contraction. The velocity measurements on a DNA-origami filament that preserves the motors’ distribution showed that geometrical hindrance and biased detachment generate a non-zero sliding velocity even without rotation of the myosin lever-arm, which is widely recognized as the basic event in muscle contraction. Because biased detachment is a mechanism for the rectification of thermal fluctuations, in the Brownian-ratchet framework, we predict that it requires a non-negligible amount of energy to preserve the second law of thermodynamics. Taken together, our theoretical and experimental results elucidate less considered components in the chemo-mechanical energy transduction in muscle.

Force-velocity and tension transient measurements from Drosophila jump muscle reveal the necessity of both weakly-bound cross-bridges and series elasticity in models of muscle contraction

Muscle contraction is a fundamental biological process where molecular interactions between the myosin molecular motor and actin filaments result in contraction of a whole muscle, a process spanning size scales differing in eight orders of magnitude. Since unique behavior is observed at every scale in between these two extremes, to fully understand muscle function it is vital to develop multi-scale models. Based on simulations of classic measurements of muscle heat generation as a function of work, and shortening rate as a function of applied force, we hypothesize that a model based on molecular measurements must be modified to include a weakly-bound interaction between myosin and actin in order to fit measurements at the muscle fiber or whole muscle scales. This hypothesis is further supported by the model's need for a weakly-bound state in order to qualitatively reproduce the force response that occurs when a muscle fiber is rapidly stretched a small distance. We tested this hypothesis by measuring steady-state force as a function of shortening velocity, and the force transient caused by a rapid length step in Drosophila jump muscle fibers. Then, by performing global parameter optimization, we quantitatively compared the predictions of two mathematical models, one lacking a weakly-bound state and one with a weakly-bound state, to these measurements. Both models could reproduce our force-velocity measurements, but only the model with a weakly-bound state could reproduce our force transient measurements. However, neither model could concurrently fit both measurements. We find that only a model that includes weakly-bound cross-bridges with force-dependent detachment and an elastic element in series with the cross-bridges is able to fit both of our measurements. This result suggests that the force response after stretch is not a reflection of distinct steps in the cross-bridge cycle, but rather arises from the interaction of cross-bridges with a series elastic element. Additionally, the model suggests that the curvature of the force-velocity relationship arises from a combination of the force-dependence of weakly- and strongly-bound cross-bridges. Overall, this work presents a minimal cross-bridge model that has predictive power at the fiber level.

Switch-1 instability at the active site decouples ATP hydrolysis from force generation in myosin II

Myosin active site elements (i.e., switch‐1) bind both ATP and a divalent metal to coordinate ATP hydrolysis. ATP hydrolysis at the active site is linked via allosteric communication to the actin polymer binding site and lever arm movement, thus coupling the free energy of ATP hydrolysis to force generation. How active site motifs are functionally linked to actin binding and the power stroke is still poorly understood. We hypothesize that destabilizing switch‐1 movement at the active site will negatively affect the tight coupling of the ATPase catalytic cycle to force production. Using a metal‐switch system, we tested the effect of interfering with switch‐1 coordination of the divalent metal cofactor on force generation. We found that while ATPase activity increased, motility was inhibited. Our results demonstrate that a single atom change that affects the switch‐1 interaction with the divalent metal directly affects actin binding and productive force generation. Even slight modification of the switch‐1 divalent metal coordination can decouple ATP hydrolysis from motility. Switch‐1 movement is therefore critical for both structural communication with the actin binding site, as well as coupling the energy of ATP hydrolysis to force generation.

Hypothesis: Single Actomyosin Properties Account for Ensemble Behavior in Active Muscle Shortening and Isometric Contraction

Muscle contraction results from cyclic interactions between myosin II motors and actin with two sets of proteins organized in overlapping thick and thin filaments, respectively, in a nearly crystalline lattice in a muscle sarcomere. However, a sarcomere contains a huge number of other proteins, some with important roles in muscle contraction. In particular, these include thin filament proteins, troponin and tropomyosin; thick filament proteins, myosin binding protein C; and the elastic protein, titin, that connects the thin and thick filaments. Furthermore, the order and 3D organization of the myofilament lattice may be important per se for contractile function. It is possible to model muscle contraction based on actin and myosin alone with properties derived in studies using single molecules and biochemical solution kinetics. It is also possible to reproduce several features of muscle contraction in experiments using only isolated actin and myosin, arguing against the importance of order and accessory proteins. Therefore, in this paper, it is hypothesized that “single molecule actomyosin properties account for the contractile properties of a half sarcomere during shortening and isometric contraction at almost saturating Ca concentrations”. In this paper, existing evidence for and against this hypothesis is reviewed and new modeling results to support the arguments are presented. Finally, further experimental tests are proposed, which if they corroborate, at least approximately, the hypothesis, should significantly benefit future effective analysis of a range of experimental studies, as well as drug discovery efforts.

Active Force Generation in Cardiac Muscle Cells: Mathematical Modeling and Numerical Simulation of the Actin-Myosin Interaction

Cardiac in silico numerical simulations are based on mathematical models describing the physical processes involved in the heart function. In this review paper, we critically survey biophysically-detailed mathematical models describing the subcellular mechanisms behind the generation of active force, that is the process by which the chemical energy of ATP (adenosine triphosphate) is transformed into mechanical work, thus making the muscle tissue contract. While presenting these models, that feature different levels of biophysical detail, we analyze the trade-off between the accuracy in the description of the subcellular mechanisms and the number of parameters that need to be estimated from experiments. Then, we focus on a generalized version of the classic Huxley model, built on the basis of models available in the literature, that is able to reproduce the main experimental characterizations associated to the time scales typical of a heartbeat-such as the force-velocity relationship and the tissue stiffness in response to small steps-featuring only four independent parameters. Finally, we show how those parameters can be calibrated starting from macroscopic measurements available from experiments.

Hierarchical modeling of force generation in cardiac muscle

Performing physiologically relevant simulations of the beating heart in clinical context requires to develop detailed models of the microscale force generation process. These models, however, may reveal difficult to implement in practice due to their high computational costs and complex calibration. We propose a hierarchy of three interconnected muscle contraction models-from the more refined to the more simplified-that are rigorously and systematically related to each other, offering a way to select, for a specific application, the model that yields a good trade-off between physiological fidelity, computational cost and calibration complexity. The three model families are compared to the same set of experimental data to systematically assess what physiological indicators can be reproduced or not and how these indicators constrain the model parameters. Finally, we discuss the applicability of these models for heart simulation.

Cardiac muscle regulatory units are predicted to interact stronger than neighboring cross-bridges

Strong interactions between cross-bridges (XB) and regulatory units (RU) lead to a steep response of cardiac muscle to an increase in intracellular calcium. We developed a model to quantitatively assess the influence of different types of interactions within the sarcomere on the properties of cardiac muscle. In the model, the ensembles consisting of cross-bridge groups connected by elastic tropomyosin are introduced, and their dynamics is described by a set of partial differential equations. Through large scans in the free energy landscape, we demonstrate the different influence of RU-RU, XB-XB, and XB-RU interactions on the cooperativity coefficient of calcium binding, developed maximal force, and calcium sensitivity. The model solution was fitted to reproduce experimental data on force development during isometric contraction, shortening in physiological contraction, and ATP consumption by acto-myosin. On the basis of the fits, we quantified the free energy change introduced through RU-RU and XB-XB interactions and showed that RU-RU interaction leads to ~ 5 times larger change in the free energy profile of the reaction than XB-XB interaction. Due to the deterministic description of muscle contraction and its thermodynamic consistency, we envision that the developed model can be used to study heart muscle biophysics on tissue and organ levels.

The effects of inorganic phosphate on muscle force development and energetics: challenges in modelling related to experimental uncertainties

Muscle force and power are developed by myosin cross-bridges, which cyclically attach to actin, undergo a force-generating transition and detach under turnover of ATP. The force-generating transition is intimately associated with release of inorganic phosphate (Pi) but the exact sequence of events in relation to the actual Pi release step is controversial. Details of this process are reflected in the relationships between [Pi] and the developed force and shortening velocity. In order to account for these relationships, models have proposed branched kinetic pathways or loose coupling between biochemical and force-generating transitions. A key hypothesis underlying the present study is that such complexities are not required to explain changes in the force–velocity relationship and ATP turnover rate with altered [Pi]. We therefore set out to test if models without branched kinetic paths and Pi-release occurring before the main force-generating transition can account for effects of varied [Pi] (0.1–25 mM). The models tested, one assuming either linear or non-linear cross-bridge elasticity, account well for critical aspects of muscle contraction at 0.5 mM Pi but their capacity to account for the maximum power output vary. We find that the models, within experimental uncertainties, account for the relationship between [Pi] and isometric force as well as between [Pi] and the velocity of shortening at low loads. However, in apparent contradiction with available experimental findings, the tested models produce an anomalous force–velocity relationship at elevated [Pi] and high loads with more than one possible velocity for a given load. Nevertheless, considering experimental uncertainties and effects of sarcomere non-uniformities, these discrepancies are insufficient to refute the tested models in favour of more complex alternatives.

Blebbistatin Effects Expose Hidden Secrets in the Force-Generating Cycle of Actin and Myosin

Cyclic interactions between myosin II motors and actin filaments driven by ATP turnover underlie muscle contraction and have key roles in the motility of nonmuscle cells. A remaining enigma in the understanding of this interaction is the relationship between the force-generating structural change and the release of the ATP-hydrolysis product, inorganic phosphate (Pi), from the active site of myosin. Here, we use the small molecular compound blebbistatin to probe otherwise hidden states and transitions in this process. Different hypotheses for the Pi release mechanism are tested by interpreting experimental results from in vitro motility assays and isolated muscle fibers in terms of mechanokinetic actomyosin models. The data fit with ideas that actomyosin force generation is preceded by Pi release, which in turn is preceded by two serial transitions after/coincident with cross-bridge attachment. Blebbistatin changes the rate limitation of the cycle from the first to the second of these transitions, uncovering functional roles of an otherwise short-lived pre-power stroke state that has been implicated by structural data.

A chemical kinetic model for Ca2+ induced spontaneous oscillatory contraction of myocardium

The Ca²⁺ induced Spontaneous Oscillatory Contraction (Ca-SPOC) of cardiac myofibrils oscillate with a period similar to resting heartbeat of several animal species, and its auto-oscillatory properties set the basic rhythm of cardiac contraction. To explain the dynamics of Ca-SPOC, the present paper constructs a novel chemical kinetical model based upon the cooperative behavior between the two heads of myosin II dimer, also considering the reaction-diffusion effect of ATP inside myocardial fibers. The simulation results show that the concentration of ATP inside myocardial fibers oscillates over time under some special conditions, together with the proportions of myosin II dimers in different states periodically changing with time, which contributes to produce the sustained oscillations of contractive tension. These results indicate that the SPOC of muscles may be partly due to chemical oscillation involved in the actomyosin ATPase cycle, which has been ignored by the previous theoretical studies.

On the Shape of the Force-Velocity Relationship in Skeletal Muscles: The Linear, the Hyperbolic, and the Double-Hyperbolic

The shape of the force-velocity (F-V) relationship has important implications for different aspects of muscle physiology, such as muscle efficiency and fatigue, the understanding of the pathophysiology of several myopathies or the mechanisms of muscle contraction per se, and may be of relevance for other fields, such as the development of robotics and prosthetic applications featuring natural muscle-like properties. However, different opinions regarding the shape of the F-V relationship and the underlying mechanisms exist in the literature. In this review, we summarize relevant evidence on the shape of the F-V relationship obtained over the last century. Studies performed at multiple scales ranging from the sarcomere to the organism level have described the concentric F-V relationship as linear, hyperbolic or double-hyperbolic. While the F-V relationship has most frequently been described as a rectangular hyperbola, a large number of studies have found deviations from the hyperbolic function at both ends of the F-V relation. Indeed, current evidence suggests that the F-V relation in skeletal muscles follows a double-hyperbolic pattern, with a breakpoint located at very high forces/low velocities, which may be a direct consequence of the kinetic properties of myofilament cross-bridge formation. Deviations at low forces/high velocities, by contrast, may be related to a recently discovered, calcium-independent regulatory mechanism of muscle contraction, which may also explain the low metabolic cost of very fast muscle shortening contractions. Controversial results have also been reported regarding the eccentric F-V relationship, with studies in prepared muscle specimens suggesting that maximum eccentric force is substantially greater than isometric force, whereas in vivo studies in humans show only a modest increase, no change, or even a decrease in force in lengthening contractions. This review discusses possible reasons reported in the literature for these discrepant findings, including the testing procedures (familiarization, pre-load condition, and temperature) and a potential neural inhibition at higher lengthening velocities. Finally, some unresolved questions and recommendations for F-V testing in humans are reported at the end of this document.

Comparing models with one versus multiple myosin-binding sites per actin target zone: The power of simplicity

Mechanokinetic statistical models describe the mechanisms of muscle contraction on the basis of the average behavior of a large ensemble of actin-myosin motors. Such models often assume that myosin II motor domains bind to regularly spaced, discrete target zones along the actin-based thin filaments and develop force in a series of strain-dependent transitions under the turnover of ATP. The simplest models assume that there is just one myosin-binding site per target zone and a uniform spatial distribution of the myosin motor domains in relation to each site. However, most of the recently developed models assume three myosin-binding sites per target zone, and some models include a spatially explicit 3-D treatment of the myofilament lattice and thereby of the geometry of the actin-myosin contact points. Here, I show that the predictions for steady-state contractile behavior of muscle are very similar whether one or three myosin-binding sites per target zone is assumed, provided that the model responses are appropriately scaled to the number of sites. Comparison of the model predictions for isometrically contracting mammalian muscle cells suggests that each target zone contains three or more myosin-binding sites. Finally, I discuss the strengths and weaknesses of one-site spatially inexplicit models in relation to three-site models, including those that take into account the detailed 3-D geometry of the myofilament lattice. The results of this study suggest that single-site models, with reduced computational cost compared with multisite models, are useful for several purposes, e.g., facilitated molecular mechanistic insights.

Nonlinear Actomyosin Elasticity in Muscle?

Cyclic interactions between myosin II motor domains and actin filaments that are powered by turnover of ATP underlie muscle contraction and have key roles in motility of nonmuscle cells. The elastic characteristics of actin-myosin cross-bridges are central in the force-generating process, and disturbances in these properties may lead to disease. Although the prevailing paradigm is that the cross-bridge elasticity is linear (Hookean), recent single-molecule studies suggest otherwise. Despite convincing evidence for substantial nonlinearity of the cross-bridge elasticity in the single-molecule work, this finding has had limited influence on muscle physiology and physiology of other ordered cellular actin-myosin ensembles. Here, we use a biophysical modeling approach to close the gap between single molecules and physiology. The model is used for analysis of available experimental results in the light of possible nonlinearity of the cross-bridge elasticity. We consider results obtained both under rigor conditions (in the absence of ATP) and during active muscle contraction. Our results suggest that a wide range of experimental findings from mechanical experiments on muscle cells are consistent with nonlinear actin-myosin elasticity similar to that previously found in single molecules. Indeed, the introduction of nonlinear cross-bridge elasticity into the model improves the reproduction of key experimental results and eliminates the need for force dependence of the ATP-induced detachment rate, consistent with observations in other single-molecule studies. The findings have significant implications for the understanding of key features of actin-myosin-based production of force and motion in living cells, particularly in muscle, and for the interpretation of experimental results that rely on stiffness measurements on cells or myofibrils.

A new myofilament contraction model with ATP consumption for ventricular cell model

A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.

Do Actomyosin Single-Molecule Mechanics Data Predict Mechanics of Contracting Muscle?

In muscle, but not in single-molecule mechanics studies, actin, myosin and accessory proteins are incorporated into a highly ordered myofilament lattice. In view of this difference we compare results from single-molecule studies and muscle mechanics and analyze to what degree data from the two types of studies agree with each other. There is reasonable correspondence in estimates of the cross-bridge power-stroke distance (7–13 nm), cross-bridge stiffness (~2 pN/nm) and average isometric force per cross-bridge (6–9 pN). Furthermore, models defined on the basis of single-molecule mechanics and solution biochemistry give good fits to experimental data from muscle. This suggests that the ordered myofilament lattice, accessory proteins and emergent effects of the sarcomere organization have only minor modulatory roles. However, such factors may be of greater importance under e.g., disease conditions. We also identify areas where single-molecule and muscle data are conflicting: (1) whether force generation is an Eyring or Kramers process with just one major power-stroke or several sub-strokes; (2) whether the myofilaments and the cross-bridges have Hookean or non-linear elasticity; (3) if individual myosin heads slip between actin sites under certain conditions, e.g., in lengthening; or (4) if the two heads of myosin cooperate.

Temperature Effects on Force and Actin⁻Myosin Interaction in Muscle: A Look Back on Some Experimental Findings

Observations made in temperature studies on mammalian muscle during force development, shortening, and lengthening, are re-examined. The isometric force in active muscle goes up substantially on warming from less than 10 °C to temperatures closer to physiological (>30 °C), and the sigmoidal temperature dependence of this force has a half-maximum at ~10 °C. During steady shortening, when force is decreased to a steady level, the sigmoidal curve is more pronounced and shifted to higher temperatures, whereas, in lengthening muscle, the curve is shifted to lower temperatures, and there is a less marked increase with temperature. Even with a small rapid temperature-jump (T-jump), force in active muscle rises in a definitive way. The rate of tension rise is slower with adenosine diphosphate (ADP) and faster with increased phosphate. Analysis showed that a T-jump enhances an early, pre-phosphate release step in the acto-myosin (crossbridge) ATPase cycle, thus inducing a force-rise. The sigmoidal dependence of steady force on temperature is due to this endothermic nature of crossbridge force generation. During shortening, the force-generating step and the ATPase cycle are accelerated, whereas during lengthening, they are inhibited. The endothermic force generation is seen in different muscle types (fast, slow, and cardiac). The underlying mechanism may involve a structural change in attached myosin heads and/or their attachments on heat absorption.

The force-generation process in active muscle is strain sensitive and endothermic: a temperature-perturbation study

In experiments on active muscle, we examined the tension decline and its temperature sensitivity at the onset of ramp shortening and at a range of velocities. A segment (∼1.5 mm long) of a skinned muscle fibre isolated from rabbit psoas muscle was held isometrically (sarcomere length ∼2.5 µm) at 8-9°C, maximally Ca²⁺-activated and a ramp shortening applied. The tension decline with a ramp shortening showed an early decrease of slope (the P₁ transition) followed by a slower decrease in slope (the P₂ transition) to the steady (isotonic) force. The tension level at the initial P₁ transition and the time to that transition decreased as the velocity was increased; the length change to this transition increased with shortening velocity to a steady value of ∼8 nm half-sarcomere^-1 A small, rapid, temperature jump (T-jump) (3-4°C, <0.2 ms) applied coincident with the onset of ramp shortening showed force enhancement by T-jump and changed the tension decline markedly. Analyses showed that the rate of T-jump-induced force rise increased linearly with increase of shortening velocity. These results provide crucial evidence that the strain-sensitive cross-bridge force generation, or a step closely coupled to it, is endothermic.

Robust mechanobiological behavior emerges in heterogeneous myosin systems

Biological complexity presents challenges for understanding natural phenomenon and engineering new technologies, particularly in systems with molecular heterogeneity. Such complexity is present in myosin motor protein systems, and computational modeling is essential for determining how collective myosin interactions produce emergent system behavior. We develop a computational approach for altering myosin isoform parameters and their collective organization, and support predictions with in vitro experiments of motility assays with α-actinins as molecular force sensors. The computational approach models variations in single myosin molecular structure, system organization, and force stimuli to predict system behavior for filament velocity, energy consumption, and robustness. Robustness is the range of forces where a filament is expected to have continuous velocity and depends on used myosin system energy. Myosin systems are shown to have highly nonlinear behavior across force conditions that may be exploited at a systems level by combining slow and fast myosin isoforms heterogeneously. Results suggest some heterogeneous systems have lower energy use near stall conditions and greater energy consumption when unloaded, therefore promoting robustness. These heterogeneous system capabilities are unique in comparison with homogenous systems and potentially advantageous for high performance bionanotechnologies. Findings open doors at the intersections of mechanics and biology, particularly for understanding and treating myosin-related diseases and developing approaches for motor molecule-based technologies.

Modulators of actin-myosin dissociation: basis for muscle type functional differences during fatigue

The muscle types present with variable fatigue tolerance, in part due to the myosin isoform expressed. However, the critical steps that define “fatigability” in vivo of fast vs. slow myosin isoforms, at the molecular level, are not yet fully understood. We examined the modulation of the ATP-induced myosin subfragment 1 (S1) dissociation from pyrene-actin by inorganic phosphate (P_i), pH, and temperature using a specially modified stopped-flow system that allowed fast kinetics measurements at physiological temperature. We contrasted the properties of rabbit psoas (fast) and bovine masseter (slow) myosins (obtained from samples collected from New Zealand rabbits and from a licensed abattoir, respectively, according to institutional and national ethics permits). To identify ATP cycling biochemical intermediates, we assessed ATP binding to a preequilibrated mixture of actomyosin and variable [ADP], pH (pH 7 vs. pH 6.2), and P_i (zero, 15, or 30 added mM P_i) in a range of temperatures (5 to 45°C). Temperature and pH variations had little, if any, effect on the ADP dissociation constant (K_ADP) for fast S1, but for slow S1, K_ADP was weakened with increasing temperature or low pH. In the absence of ADP, the dissociation constant for phosphate (K_Pi) was weakened with increasing temperature for fast S1. In the presence of ADP, myosin type differences were revealed at the apparent phosphate affinity, depending on pH and temperature. Overall, the newly revealed kinetic differences between myosin types could help explain the in vivo observed muscle type functional differences at rest and during fatigue.

Reinterpretation of the Tension Response of Muscle to Stretches and Releases

We have reexamined the experimental time courses of tension in frog muscle after rapid length steps. The early tension recoveries are biexponential. After 3 nm/hs stretches and releases, the rates of the immediate rapid tension changes are similar but the subsequent tension fall after a stretch is much slower than the rise after a release. After 1.5 nm/hs length steps, the entire tension responses are more nearly mirror images. To identify the underlying processes, we used a model of the muscle cross-bridge cycle with two tension-generating (tensing) steps. Analysis of the time course of the tension, the rates of the steps in the cycle, and their contributions to tension provided insights into previously puzzling features of the experimental response. After a stretch, the initial rapid tension fall in the model is caused principally by the reversal of the first tensing step, but after a few milliseconds the tensing step resumes its forward direction. We conclude that the remaining response should not be included in phase 2, the period of early tension recovery. With this exclusion, T₂, the tension at the end of this period, rises with an increase of stretch. The rate of early tension recovery also increases with stretch size, showing that the reversal of the first tensing step is strain sensitive. After small length steps, the fast and slow components of the early tension recovery are both caused mainly by the first tensing step. The fast component is triggered by the initial sliding of the filaments, and the slow component is due to further sliding that occurs as the tension recovers. With small length steps (<0.5 nm/hs), the time course of the response to a stretch is the reverse of that to a release.

Actomyosin based contraction: one mechanokinetic model from single molecules to muscle?

Bridging the gaps between experimental systems on different hierarchical scales is needed to overcome remaining challenges in the understanding of muscle contraction. Here, a mathematical model with well-characterized structural and biochemical actomyosin states is developed to that end. We hypothesize that this model accounts for generation of force and motion from single motor molecules to the large ensembles of muscle. In partial support of this idea, a wide range of contractile phenomena are reproduced without the need to invoke cooperative interactions or ad hoc states/transitions. However, remaining limitations exist, associated with ambiguities in available data for model definition e.g.: (1) the affinity of weakly bound cross-bridges, (2) the characteristics of the cross-bridge elasticity and (3) the exact mechanistic relationship between the force-generating transition and phosphate release in the actomyosin ATPase. Further, the simulated number of attached myosin heads in the in vitro motility assay differs several-fold from duty ratios, (fraction of strongly attached ATPase cycle times) derived in standard analysis. After addressing the mentioned issues the model should be useful in fundamental studies, for engineering of myosin motors as well as for studies of muscle disease and drug development.

Mechanosensing in Myosin Filament Solves a 60 Years Old Conflict in Skeletal Muscle Modeling between High Power Output and Slow Rise in Tension

Almost 60 years ago Andrew Huxley with his seminal paper (Huxley, 1957) laid the foundation of modern muscle modeling, linking chemical to mechanical events. He described mechanics and energetics of muscle contraction through the cyclical attachment and detachment of myosin motors to the actin filament with ad-hoc assumptions on the dependence of the rate constants on the strain of the myosin motors. That relatively simple hypothesis is still present in recent models, even though with several modifications to adapt the model to the different experimental constraints which became subsequently available. However, already in that paper, one controversial aspect of the model became clear. Relatively high attachment and detachment rates of myosin to the actin filament were needed to simulate the high power output at intermediate velocity of shortening. However, these rates were incompatible with the relatively slow rise in tension upon activation, despite the rise should be generated by the same rate functions. This discrepancy has not been fully solved till today, despite several hypotheses have been forwarded to reconcile the two aspects. Here, using a conventional muscle model, we show that the recently revealed mechanosensing mechanism of recruitment of myosin motors (Linari et al., 2015) can solve this long standing problem without any further ad-hoc hypotheses.

Including Thermal Fluctuations in Actomyosin Stable States Increases the Predicted Force per Motor and Macroscopic Efficiency in Muscle Modelling

Muscle contractions are generated by cyclical interactions of myosin heads with actin filaments to form the actomyosin complex. To simulate actomyosin complex stable states, mathematical models usually define an energy landscape with a corresponding number of wells. The jumps between these wells are defined through rate constants. Almost all previous models assign these wells an infinite sharpness by imposing a relatively simple expression for the detailed balance, i.e., the ratio of the rate constants depends exponentially on the sole myosin elastic energy. Physically, this assumption corresponds to neglecting thermal fluctuations in the actomyosin complex stable states. By comparing three mathematical models, we examine the extent to which this hypothesis affects muscle model predictions at the single cross-bridge, single fiber, and organ levels in a ceteris paribus analysis. We show that including fluctuations in stable states allows the lever arm of the myosin to easily and dynamically explore all possible minima in the energy landscape, generating several backward and forward jumps between states during the lifetime of the actomyosin complex, whereas the infinitely sharp minima case is characterized by fewer jumps between states. Moreover, the analysis predicts that thermal fluctuations enable a more efficient contraction mechanism, in which a higher force is sustained by fewer attached cross-bridges.

Mathematical models are of fundamental importance in the quantitative verification of biological hypotheses. Muscle contraction models assume the existence of several stable states for the myosin head, whereas the transition rates between states are defined to fit experimental data. The ratio of the forward and backward rates is linked to the ratio of the probabilities of being in one or other stable state at equilibrium through a detailed balance condition. A commonly used assumption leads to a relatively simple expression for this balance condition that depends only on the values of the energy at the minima and not on the minima shape. Mathematically, this hypothesis corresponds to infinite sharpness at these minima; physically, it neglects the small thermal fluctuations within actomyosin stable states. In this work, we compare this classical approach with a model that includes thermal fluctuations within wide minima, and quantitatively assess how much this hypothesis affects the model outcomes at the single molecule, single fiber, and whole heart levels. It is shown that, using parameters compatible with known behavior in muscle mechanics, relaxing the infinitely sharp minima hypothesis improves the predicted force generation and efficiency at the macroscopic level.

Weightlifting and the actomyosin cycle

How does a human lift a weight? Can we relate the dynamics of the lift to the molecular actin–myosin interactions responsible for muscle contraction? We address these questions with bench press experiments that we analyse with a theoretical model, based on the sliding filament theory. The agreement is fair, and we discuss its possible extension to medical diagnostics.

Direct Measurements of Local Coupling between Myosin Molecules Are Consistent with a Model of Muscle Activation

Muscle contracts due to ATP-dependent interactions of myosin motors with thin filaments composed of the proteins actin, troponin, and tropomyosin. Contraction is initiated when calcium binds to troponin, which changes conformation and displaces tropomyosin, a filamentous protein that wraps around the actin filament, thereby exposing myosin binding sites on actin. Myosin motors interact with each other indirectly via tropomyosin, since myosin binding to actin locally displaces tropomyosin and thereby facilitates binding of nearby myosin. Defining and modeling this local coupling between myosin motors is an open problem in muscle modeling and, more broadly, a requirement to understanding the connection between muscle contraction at the molecular and macro scale. It is challenging to directly observe this coupling, and such measurements have only recently been made. Analysis of these data suggests that two myosin heads are required to activate the thin filament. This result contrasts with a theoretical model, which reproduces several indirect measurements of coupling between myosin, that assumes a single myosin head can activate the thin filament. To understand this apparent discrepancy, we incorporated the model into stochastic simulations of the experiments, which generated simulated data that were then analyzed identically to the experimental measurements. By varying a single parameter, good agreement between simulation and experiment was established. The conclusion that two myosin molecules are required to activate the thin filament arises from an assumption, made during data analysis, that the intensity of the fluorescent tags attached to myosin varies depending on experimental condition. We provide an alternative explanation that reconciles theory and experiment without assuming that the intensity of the fluorescent tags varies.

Despite decades of study, there is no clear connection between muscle contraction at the molecular and the macroscopic scale. For example, we cannot yet predict how a genetic defect in a muscle protein will result in a physiological change in the heart. This multi-scale understanding is difficult, in part, because molecules cooperate during muscle contraction; that is, one molecule’s behavior is influenced by the behavior of its neighbors. It is difficult to make direct measurements from such coupled molecular systems and also difficult to describe them quantitatively. Despite these obstacles, we recently published experimental measurements and theoretical models of this coupling, but there were apparent discrepancies between the two. Here, we use detailed computer simulations of these experiments to show that, in fact, the measurements agree with the model to a remarkable extent. This agreement suggests that the model captures the essential molecular events that underlie the coupling between muscle molecules. This removes a major obstacle to a multi-scale understanding of muscle contraction and, while more work is necessary, suggests that a connection between the molecular and macroscopic scale is within reach.

Cross-Bridge Group Ensembles Describing Cooperativity in Thermodynamically Consistent Way

The aim of this work is to incorporate cooperativity into Huxley-type cross-bridge model in thermodynamically consistent way. While the Huxley-type models assume that cross-bridges act independently from each other, we take into account that each cross-bridge is influenced by its neighbors and cooperativity is induced by tropomyosin movement. For that, we introduce ensembles of cross-bridge groups connected by elastic tropomyosin. By taking into account that the mechanical displacement of tropomyosin induces free energy change of the cross-bridge group ensemble, we develop the formalism for thermodynamically consistent description of the cooperativity in muscle contraction. An example model was composed to test the approach. The model parameters were found by optimization from the linear relation between oxygen consumption and stress-strain area as well as experimentally measured stress dynamics of rat trabecula. We have found a good agreement between the optimized model solution and experimental data. Simulations also showed that it is possible to study cooperativity with the approach developed in this work.

Isoforms Confer Characteristic Force Generation and Mechanosensation by Myosin II Filaments

Myosin II isoforms with varying mechanochemistry and filament size interact with filamentous actin (F-actin) arrays to generate contractile forces in muscle and nonmuscle cells. How myosin II force production is shaped by isoform-specific motor properties and environmental stiffness remains poorly understood. Here, we used computer simulations to analyze force production by an ensemble of myosin motors against an elastically tethered actin filament. We found that force output depends on two timescales: the duration of F-actin attachment, which varies sharply with the ensemble size, motor duty ratio, and external load; and the time to build force, which scales with the ensemble stall force, gliding speed, and environmental stiffness. Although force-dependent kinetics were not required to sense changes in stiffness, the myosin catch bond produced positive feedback between the attachment time and force to trigger switch-like transitions from transient attachments, generating small forces, to high-force-generating runs. Using parameters representative of skeletal muscle myosin, nonmuscle myosin IIB, and nonmuscle myosin IIA revealed three distinct regimes of behavior, respectively: 1) large assemblies of fast, low-duty ratio motors rapidly build stable forces over a large range of environmental stiffness; 2) ensembles of slow, high-duty ratio motors serve as high-affinity cross-links with force buildup times that exceed physiological timescales; and 3) small assemblies of low-duty ratio motors operating at intermediate speeds are poised to respond sharply to changes in mechanical context-at low force or stiffness, they serve as low-affinity cross-links, but they can transition to force production via the positive-feedback mechanism described above. Together, these results reveal how myosin isoform properties may be tuned to produce force and respond to mechanical cues in their environment.

Poorly understood aspects of striated muscle contraction

Muscle contraction results from cyclic interactions between the contractile proteins myosin and actin, driven by the turnover of adenosine triphosphate (ATP). Despite intense studies, several molecular events in the contraction process are poorly understood, including the relationship between force-generation and phosphate-release in the ATP-turnover. Different aspects of the force-generating transition are reflected in the changes in tension development by muscle cells, myofibrils and single molecules upon changes in temperature, altered phosphate concentration, or length perturbations. It has been notoriously difficult to explain all these events within a given theoretical framework and to unequivocally correlate observed events with the atomic structures of the myosin motor. Other incompletely understood issues include the role of the two heads of myosin II and structural changes in the actin filaments as well as the importance of the three-dimensional order. We here review these issues in relation to controversies regarding basic physiological properties of striated muscle. We also briefly consider actomyosin mutation effects in cardiac and skeletal muscle function and the possibility to treat these defects by drugs.

Nanowire-imposed geometrical control in studies of actomyosin motor function

Recently, molecular motor gliding assays with actin and myosin from muscle have been realized on semiconductor nanowires coated with Al2O3. This opens for unique nanotechnological applications and novel fundamental studies of actomyosin motor function. Here, we provide a comparison of myosin-driven actin filament motility on Al2O3 to both nitrocellulose and trimethylchlorosilane derivatized surfaces. We also show that actomyosin motility on the less than 200 nm wide tips of arrays of Al2O3-coated nanowires can be used to control the number, and density, of myosin-actin attachment points. Results obtained using nanowire arrays with different inter-wire spacing are consistent with the idea that the actin filament sliding velocity is determined both by the total number and the average density of attached myosin heads along the actin filament. Further, the results are consistent with buckling of long myosin-free segments of the filaments as a factor underlying reduced velocity. On the other hand, the findings do not support a mechanistic role in decreasing velocity, of increased nearest neighbor distance between available myosin heads. Our results open up for more advanced studies that may use nanowire-based structures for fundamental investigations of molecular motors, including the possibility to create a nanowire-templated bottom-up assembly of 3D, muscle-like structures.

Muscle activation described with a differential equation model for large ensembles of locally coupled molecular motors

Molecular motors, by turning chemical energy into mechanical work, are responsible for active cellular processes. Often groups of these motors work together to perform their biological role. Motors in an ensemble are coupled and exhibit complex emergent behavior. Although large motor ensembles can be modeled with partial differential equations (PDEs) by assuming that molecules function independently of their neighbors, this assumption is violated when motors are coupled locally. It is therefore unclear how to describe the ensemble behavior of the locally coupled motors responsible for biological processes such as calcium-dependent skeletal muscle activation. Here we develop a theory to describe locally coupled motor ensembles and apply the theory to skeletal muscle activation. The central idea is that a muscle filament can be divided into two phases: an active and an inactive phase. Dynamic changes in the relative size of these phases are described by a set of linear ordinary differential equations (ODEs). As the dynamics of the active phase are described by PDEs, muscle activation is governed by a set of coupled ODEs and PDEs, building on previous PDE models. With comparison to Monte Carlo simulations, we demonstrate that the theory captures the behavior of locally coupled ensembles. The theory also plausibly describes and predicts muscle experiments from molecular to whole muscle scales, suggesting that a micro- to macroscale muscle model is within reach.

Hypothesis and theory: mechanical instabilities and non-uniformities in hereditary sarcomere myopathies

Familial hypertrophic cardiomyopathy (HCM), due to point mutations in genes for sarcomere proteins such as myosin, occurs in 1/500 people and is the most common cause of sudden death in young individuals. Similar mutations in skeletal muscle, e.g., in the MYH7 gene for slow myosin found in both the cardiac ventricle and slow skeletal muscle, may also cause severe disease but the severity and the morphological changes are often different. In HCM, the modified protein function leads, over years to decades, to secondary remodeling with substantial morphological changes, such as hypertrophy, myofibrillar disarray, and extensive fibrosis associated with severe functional deterioration. Despite intense studies, it is unclear how the moderate mutation-induced changes in protein function cause the long-term effects. In hypertrophy of the heart due to pressure overload (e.g., hypertension), mechanical stress in the myocyte is believed to be major initiating stimulus for activation of relevant cell signaling cascades. Here it is considered how expression of mutated proteins, such as myosin or regulatory proteins, could have similar consequences through one or both of the following mechanisms: (1) contractile instabilities within each sarcomere (with more than one stable velocity for a given load), (2) different tension generating capacities of cells in series. These mechanisms would have the potential to cause increased tension and/or stretch of certain cells during parts of the cardiac cycle. Modeling studies are used to illustrate these ideas and experimental tests are proposed. The applicability of similar ideas to skeletal muscle is also postulated, and differences between heart and skeletal muscle are discussed.

Dynamic coupling of regulated binding sites and cycling myosin heads in striated muscle

Myosim, a computer model of muscle contraction, includes molecular-level effects and incorporates dynamic coupling of myosin heads and binding sites.

In an activated muscle, binding sites on the thin filament and myosin heads switch frequently between different states. Because the status of the binding sites influences the status of the heads, and vice versa, the binding sites and myosin heads are dynamically coupled. The functional consequences of this coupling were investigated using MyoSim, a new computer model of muscle. MyoSim extends existing models based on Huxley-type distribution techniques by incorporating Ca²⁺ activation and cooperative effects. It can also simulate arbitrary cross-bridge schemes set by the researcher. Initial calculations investigated the effects of altering the relative speeds of binding-site and cross-bridge kinetics, and of manipulating cooperative processes. Subsequent tests fitted simulated force records to experimental data recorded using permeabilized myocardial preparations. These calculations suggest that the rate of force development at maximum activation is limited by myosin cycling kinetics, whereas the rate at lower levels of activation is limited by how quickly binding sites become available. Additional tests investigated the behavior of transiently activated cells by driving simulations with experimentally recorded Ca²⁺ signals. The unloaded shortening profile of a twitching myocyte could be reproduced using a model with two myosin states, cooperative activation, and strain-dependent kinetics. Collectively, these results demonstrate that dynamic coupling of binding sites and myosin heads is important for contractile function.

Power-stroke-driven actomyosin contractility

In ratchet-based models describing actomyosin contraction the activity is usually associated with actin binding potential while the power-stroke mechanism, residing inside myosin heads, is viewed as passive. To show that contraction can be propelled directly through a conformational change, we propose an alternative model where the power stroke is the only active mechanism. The asymmetry, ensuring directional motion, resides in steric interaction between the externally driven power-stroke element and the passive nonpolar actin filament. The proposed model can reproduce all four discrete states of the minimal actomyosin catalytic cycle even though it is formulated in terms of continuous Langevin dynamics. We build a conceptual bridge between processive and nonprocessive molecular motors by demonstrating that not only the former but also the latter can use structural transformation as the main driving force.

Hill's equation of muscle performance and its hidden insight on molecular mechanisms

Muscles shorten faster against light loads than they do against heavy loads. The hyperbolic equation first used by A.V. Hill over seven decades ago to illustrate the relationship between shortening velocity and load is still the predominant method used to characterize muscle performance, even though it has been regarded as purely empirical and lacking precision in predicting velocities at high and low loads. Popularity of the Hill equation has been sustained perhaps because of historical reasons, but its simplicity is certainly attractive. The descriptive nature of the equation does not diminish its role as a useful tool in our quest to understand animal locomotion and optimal design of muscle-powered devices like bicycles. In this Review, an analysis is presented to illustrate the connection between the historic Hill equation and the kinetics of myosin cross-bridge cycle based on the latest findings on myosin motor interaction with actin filaments within the structural confines of a sarcomere. In light of the new data and perspective, some previous studies of force-velocity relations of muscle are revisited to further our understanding of muscle mechanics and the underlying biochemical events, specifically how extracellular and intracellular environment, protein isoform expression, and posttranslational modification of contractile and regulatory proteins change the interaction between myosin and actin that in turn alter muscle force, shortening velocity, and the relationship between them.

Nonlinear cross-bridge elasticity and post-power-stroke events in fast skeletal muscle actomyosin

Generation of force and movement by actomyosin cross-bridges is the molecular basis of muscle contraction, but generally accepted ideas about cross-bridge properties have recently been questioned. Of the utmost significance, evidence for nonlinear cross-bridge elasticity has been presented. We here investigate how this and other newly discovered or postulated phenomena would modify cross-bridge operation, with focus on post-power-stroke events. First, as an experimental basis, we present evidence for a hyperbolic [MgATP]-velocity relationship of heavy-meromyosin-propelled actin filaments in the in vitro motility assay using fast rabbit skeletal muscle myosin (28-29°C). As the hyperbolic [MgATP]-velocity relationship was not consistent with interhead cooperativity, we developed a cross-bridge model with independent myosin heads and strain-dependent interstate transition rates. The model, implemented with inclusion of MgATP-independent detachment from the rigor state, as suggested by previous single-molecule mechanics experiments, accounts well for the [MgATP]-velocity relationship if nonlinear cross-bridge elasticity is assumed, but not if linear cross-bridge elasticity is assumed. In addition, a better fit is obtained with load-independent than with load-dependent MgATP-induced detachment rate. We discuss our results in relation to previous data showing a nonhyperbolic [MgATP]-velocity relationship when actin filaments are propelled by myosin subfragment 1 or full-length myosin. We also consider the implications of our results for characterization of the cross-bridge elasticity in the filament lattice of muscle.

The working stroke of the myosin II motor in muscle is not tightly coupled to release of orthophosphate from its active site

Skeletal muscle shortens faster against a lower load. This force-velocity relationship is the fundamental determinant of muscle performance in vivo and is due to ATP-driven working strokes of myosin II motors, during their cyclic interactions with the actin filament in each half-sarcomere. Crystallographic studies suggest that the working stroke is associated with the release of phosphate (Pi) and consists of 70 deg tilting of a light-chain domain that connects the catalytic domain of the myosin motor to the myosin tail and filament. However, the coupling of the working stroke with Pi release is still an unsolved question. Using nanometre-microsecond mechanics on skinned muscle fibres, we impose stepwise drops in force on an otherwise isometric contraction and record the isotonic velocity transient, to measure the mechanical manifestation of the working stroke of myosin motors and the rate of its regeneration in relation to the half-sarcomere load and [Pi]. We show that the rate constant of the working stroke is unaffected by [Pi], while the subsequent transition to steady velocity shortening is accelerated. We propose a new chemo-mechanical model that reproduces the transient and steady state responses by assuming that: (i) the release of Pi from the catalytic site of a myosin motor can occur at any stage of the working stroke, and (ii) a myosin motor, in an intermediate state of the working stroke, can slip to the next actin monomer during filament sliding. This model explains the efficient action of muscle molecular motors working as an ensemble in the half-sarcomere.

A cross-bridge cycle with two tension-generating steps simulates skeletal muscle mechanics

We examined whether cross-bridge cycle models with one or two tension-generating steps can account for the force-velocity relation of and tension response to length steps of frog skeletal muscle. Transition-state theory defined the strain dependence of the rate constants. The filament stiffness was non-Hookean. Models were refined against experimental data by simulated annealing and downhill simplex runs. Models with one tension-generating step were rejected, as they had a low efficiency and fitted the experimental data relatively poorly. The best model with two tension-generating steps (stroke distances 5.6 and 4.6 nm) and a cross-bridge stiffness of 1.7 pN/nm gave a good account of the experimental data. The two tensing steps allow an efficiency of up to 38% during shortening. In an isometric contraction, 54.7% of the attached heads were in a pre-tension-generating state, 44.5% of the attached heads had undergone the first tension-generating step, and only 0.8% had undergone both tension-generating steps; they bore 34%, 64%, and 2%, respectively, of the isometric tension. During slow shortening, the second tensing step made a greater contribution. During lengthening, up to 93% of the attached heads were in a pre-tension-generating state yet bore elevated tension by being dragged to high strains before detaching.

Global fit analysis of myosin-5b motility reveals thermodynamics of Mg2+-sensitive acto-myosin-ADP states

Kinetic and thermodynamic studies of the mechanochemical cycle of myosin motors are essential for understanding the mechanism of energy conversion. Here, we report our investigation of temperature and free Mg²⁺-ion dependencies of sliding velocities of a high duty ratio class-5 myosin motor, myosin-5b from D. discoideum using in vitro motility assays. Previous studies have shown that the sliding velocity of class-5 myosins obeys modulation by free Mg²⁺-ions. Free Mg²⁺-ions affect ADP release kinetics and the dwell time of actin-attached states. The latter determines the maximal velocity of actin translocation in the sliding filament assay. We measured the temperature dependence of sliding velocity in the range from 5 to 55°C at two limiting free Mg²⁺-ion concentrations. Arrhenius plots demonstrated non-linear behavior. Based on this observation we propose a kinetic model, which explains both sensitivity towards free Mg²⁺-ions and non-linearity of the temperature dependence of sliding velocity. According to this model, velocity is represented as a simple analytical function of temperature and free Mg²⁺-ion concentrations. This function has been applied to global non-linear fit analysis of three data sets including temperature and magnesium (at 20°C) dependence of sliding velocity. As a result we obtain thermodynamic parameters (ΔH_Mg and ΔS_Mg) of a fast equilibrium between magnesium free (AM·D) and magnesium bound acto-myosin-ADP (AM· Mg²⁺D) states and the corresponding enthalpic barriers associated with ADP release (ΔH₁^‡ and ΔH₂^‡). The herein presented integrative approach of data analysis based on global fitting can be applied to the remaining steps of the acto-myosin ATPase cycle facilitating the determination of energetic parameters and thermodynamics of acto-myosin interactions.

Muscle contraction and the elasticity-mediated crosstalk effect

Cooperative action of molecular motors is essential for many cellular processes. One possible regulator of motor coordination is the elasticity-mediated crosstalk (EMC) coupling between myosin II motors whose origin is the tensile stress that they collectively generate in actin filaments. Here, we use a statistical mechanical analysis to investigate the influence of the EMC effect on the sarcomere -- the basic contractile unit of skeletal muscles. We demonstrate that the EMC effect leads to an increase in the attachment probability of motors located near the end of the sarcomere while simultaneously decreasing the attachment probability of the motors in the central part. Such a polarized attachment probability would impair the motors' ability to cooperate efficiently. Interestingly, this undesired phenomenon becomes significant only when the system size exceeds that of the sarcomere in skeletal muscles, which provides an explanation for the remarkable lack of sarcomere variability in vertebrates. Another phenomenon that we investigate is the recently observed increase in the duty ratio of the motors with the tension in muscle. We reveal that the celebrated Hill's equation for muscle contraction is very closely related to this observation.

Mechanical coupling between myosin molecules causes differences between ensemble and single-molecule measurements

In contracting muscle, individual myosin molecules function as part of a large ensemble, hydrolyzing ATP to power the relative sliding of actin filaments. The technological advances that have enabled direct observation and manipulation of single molecules, including recent experiments that have explored myosin's force-dependent properties, provide detailed insight into the kinetics of myosin's mechanochemical interaction with actin. However, it has been difficult to reconcile these single-molecule observations with the behavior of myosin in an ensemble. Here, using a combination of simulations and theory, we show that the kinetic mechanism derived from single-molecule experiments describes ensemble behavior; but the connection between single molecule and ensemble is complex. In particular, even in the absence of external force, internal forces generated between myosin molecules in a large ensemble accelerate ADP release and increase how far actin moves during a single myosin attachment. These myosin-induced changes in strong binding lifetime and attachment distance cause measurable properties, such as actin speed in the motility assay, to vary depending on the number of myosin molecules interacting with an actin filament. This ensemble-size effect challenges the simple detachment limited model of motility, because even when motility speed is limited by ADP release, increasing attachment rate can increase motility speed.

From single molecule fluctuations to muscle contraction: a Brownian model of A.F. Huxley's hypotheses

Muscular force generation in response to external stimuli is the result of thermally fluctuating, cyclical interactions between myosin and actin, which together form the actomyosin complex. Normally, these fluctuations are modelled using transition rate functions that are based on muscle fiber behaviour, in a phenomenological fashion. However, such a basis reduces the predictive power of these models. As an alternative, we propose a model which uses direct single molecule observations of actomyosin fluctuations reported in the literature. We precisely estimate the actomyosin potential bias and use diffusion theory to obtain a brownian ratchet model that reproduces the complete cross-bridge cycle. The model is validated by simulating several macroscopic experimental conditions, while its interpretation is compatible with two different force-generating scenarios.

Actomyosin kinetics of pure fast and slow rat myosin isoforms studied by in vitro motility assay approach

An in vitro motility assay approach was used to investigate the mechanisms of the functional differences between myosin isoforms, by studying the effect of MgATP and MgADP on actin sliding velocity (V(f)) of pure slow and fast rat skeletal myosin at different temperatures. The value of V(f) depended on [MgATP] according to Michaelis-Menten kinetics, with an apparent constant (K(m)) of 54.2, 64.4 and 200 μm for the fast isoform and 18.6, 36.5 and 45.5 μM for the slow isoform at 20, 25 and 35°C, respectively. The presence of 2 mM MgADP decreased V(f) and yielded an inhibition constant (K(i)) of 377, 463 and 533 μM for the fast isoform at 20, 25 and 35°C, respectively, and 120 and 355 μM for the slow isoform at 25 and 35°C, respectively. The analysis of K(m) and K(i) suggested that slow and fast isoforms differ in the kinetics limiting V(f). Moreover, the higher sensitivity of the fast myosin isoform to a drop in [MgATP] is consistent with the higher fatigability of fast fibres than slow fibres. From the Michaelis-Menten relation in the absence of MgADP, we calculated the rate of actomyosin dissociation by MgATP (k(+ATP)) and the rate of MgADP release (k(-ADP)). We found values of k(+ATP) of 4.8 × 10(6), 6.5 × 10(6) and 6.6 × 10(6) M(-1) s(-1) for the fast isoform and 3.3 × 10(6), 2.9 × 10(6) and 6.7 × 10(6) M(-1) s(-1) for the slow isoform and values of k(-ADP) of 263, 420 and 1320 s(-1) for the fast isoform and 62, 107 and 306 s(-1) for the slow isoform at 20, 25 and 35°C, respectively. The results suggest that k(-ADP) could be the major determinant of functional differences between the fast and slow myosin isoforms at physiological temperatures.

Mechanism of force enhancement during and after lengthening of active muscle: a temperature dependence study

The aim of the present study was to examine the temperature dependence of active force in lengthening and shortening muscle. Experiments were done, in vitro, on bundles of intact fibres (fibre length L(0) ~2 mm; sarcomere length ~2.5 μm) isolated from a rat fast muscle (flexor hallucis brevis) and a ramp length change of 5-7% L(0) was applied on the plateau of an isometric tetanic contraction. Ramp lengthening increased and ramp shortening decreased the muscle tension to new approximately steady levels in a velocity-dependent way. The isometric tension and the lower steady tension reached at a given shortening velocity, increased with warming from 10 to 35 °C and the relation between tension and reciprocal absolute temperature was sigmoidal. However, the tension-temperature curve of shortening muscle was sharper and shifted to higher temperature with increased velocity. In contrast, the enhanced steady tension during lengthening at a given velocity was largely temperature-insensitive within the same temperature range; we hypothesize that the tension-temperature curve may be shifted to lower temperatures in lengthening muscle. Consequently, when normalised to the isometric tension at each temperature, the tension during lengthening at a given velocity decreased exponentially with increase of temperature. The residual force enhancement that remains after ramp lengthening showed a similar behaviour and was markedly reduced in warming from 10 to 35 °C. The findings are consistent with the thesis that active force generation in muscle is endothermic and strain-sensitive; during shortening with a faster crossbridge cycle it becomes more pronounced, but during lengthening it becomes depressed as the cycle slows in a velocity-dependent way. The residual force enhancement may be caused by the same process in addition to non-crossbridge mechanism(s).