Cross-bridge model of muscle contraction. Quantitative analysis

Modeling thick filament activation suggests a molecular basis for force depression

Multiscale models aiming to connect muscle's molecular and cellular function have been difficult to develop, in part due to a lack of self-consistent multiscale data. To address this gap, we measured the force response from single, skinned rabbit psoas muscle fibers to ramp shortenings and step stretches performed on the plateau region of the force-length relationship. We isolated myosin from the same muscles and, under similar conditions, performed single-molecule and ensemble measurements of myosin's ATP-dependent interaction with actin using laser trapping and in vitro motility assays. We fit the fiber data by developing a partial differential equation model that includes thick filament activation, whereby an increase in force on the thick filament pulls myosin out of an inhibited state. The model also includes a series elastic element and a parallel elastic element. This parallel elastic element models a titin-actin interaction proposed to account for the increase in isometric force after stretch (residual force enhancement). By optimizing the model fit to a subset of our fiber measurements, we specified seven unknown parameters. The model then successfully predicted the remainder of our fiber measurements and also our molecular measurements from the laser trap and in vitro motility. The success of the model suggests that our multiscale data are self-consistent and can serve as a testbed for other multiscale models. Moreover, the model captures the decrease in isometric force observed in our muscle fibers after active shortening (force depression), suggesting a molecular mechanism for force depression, whereby a parallel elastic element combines with thick filament activation to decrease the number of cycling cross-bridges.

Design Principles and Benefits of Spatially Explicit Models of Myofilament Function

Spatially explicit models of muscle contraction include fine-scale details about the spatial, kinetic, and/or mechanical properties of the biological processes being represented within the model network. Over the past 25 years, this has primarily consisted of a set of mathematical and computational algorithms representing myosin cross-bridge activity, Ca2+-activation of contraction, and ensemble force production within a half-sarcomere representation of the myofilament network. Herein we discuss basic design principles associated with creating spatially explicit models of myofilament function, as well as model assumptions underlying model development. A brief overview of computational approaches is introduced. Opportunities for new model directions that could investigate coupled regulatory pathways between the thick-filament and thin-filaments are also presented. Given the modular design and flexibility associated with spatially explicit models, we highlight some advantages of this approach compared to other model formulations.

Modeling Thick Filament Activation Suggests a Molecular Basis for Force Depression

Multiscale models aiming to connect muscle's molecular and cellular function have been difficult to develop, in part, due to a lack of self-consistent multiscale data. To address this gap, we measured the force response from single skinned rabbit psoas muscle fibers to ramp shortenings and step stretches performed on the plateau region of the force-length relationship. We isolated myosin from the same muscles and, under similar conditions, performed single molecule and ensemble measurements of myosin's ATP-dependent interaction with actin using laser trapping and in vitro motility assays. We fit the fiber data by developing a partial differential equation model that includes thick filament activation, whereby an increase in force on the thick filament pulls myosin out of an inhibited state. The model also includes a series elastic element and a parallel elastic element. This parallel elastic element models a titin-actin interaction proposed to account for the increase in isometric force following stretch (residual force enhancement). By optimizing the model fit to a subset of our fiber measurements, we specified seven unknown parameters. The model then successfully predicted the remainder of our fiber measurements and also our molecular measurements from the laser trap and in vitro motility. The success of the model suggests that our multiscale data are self-consistent and can serve as a testbed for other multiscale models. Moreover, the model captures the decrease in isometric force observed in our muscle fibers after active shortening (force depression), suggesting a molecular mechanism for force depression, whereby a parallel elastic element combines with thick filament activation to decrease the number of cycling cross-bridges.

New paradigms in actomyosin energy transduction: Critical evaluation of non-traditional models for orthophosphate release

Release of the ATP hydrolysis product ortophosphate (Pi) from the active site of myosin is central in chemo-mechanical energy transduction and closely associated with the main force-generating structural change, the power-stroke. Despite intense investigations, the relative timing between Pi-release and the power-stroke remains poorly understood. This hampers in depth understanding of force production by myosin in health and disease and our understanding of myosin-active drugs. Since the 1990s and up to today, models that incorporate the Pi-release either distinctly before or after the power-stroke, in unbranched kinetic schemes, have dominated the literature. However, in recent years, alternative models have emerged to explain apparently contradictory findings. Here, we first compare and critically analyze three influential alternative models proposed previously. These are either characterized by a branched kinetic scheme or by partial uncoupling of Pi-release and the power-stroke. Finally, we suggest critical tests of the models aiming for a unified picture.

Varying thin filament activation in the framework of the Huxley'57 model

Muscle contraction is triggered by the activation of the actin sites of the thin filament by calcium ions. It results that the thin filament activation level varies over time. Moreover, this activation process is also used as a regulation mechanism of the developed force. Our objective is to build a model of varying actin site activation level within the classical Huxley'57 two-state framework. This new model is obtained as an enhancement of a previously proposed formulation of the varying thick filament activation within the same framework. We assume that the state of an actin site depends on whether it is activated and whether it forms a cross-bridge with the associated myosin head, which results in four possible states. The transitions between the actin site states are controlled by the global actin sites activation level and the dynamics of these transitions is coupled with the attachment-detachment process. A preliminary calibration of the model with experimental twitch contraction data obtained at varying sarcomere lengths is performed.

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.

Passive viscoelastic response of striated muscles

Muscle cells with sarcomeric structure exhibit highly non trivial passive mechanical response. The difficulty of its continuum modeling is due to the presence of long-range interactions transmitted by extended protein skeleton. To build a rheological model for muscle 'material', we use a stochastic micromodel, and derive a linear response theory for a half-sarcomere, which can be extended to the whole fibre. Instead of the first order rheological equation, anticipated by Hill on the phenomenological grounds, we obtain a novel second order equation which shows that tension depends not only on its current length and the velocity of stretching, but also on its acceleration. Expressing the model in terms of elementary rheological elements, we show that one contribution to the visco-elastic properties of the fibre originates in cross-bridges, while the other can be linked to inert elements which move in the sarcoplasm. We apply this model to explain the striking qualitative difference between the relaxation in experiments involving perturbation of length vs. those involving perturbation of force, and we use the values of the microscopic parameters for frog muscles to show that the model is in excellent quantitative agreement with physiological experiments.

A multiscale sliding filament model of lymphatic muscle pumping

The lymphatics maintain fluid balance by returning interstitial fluid to veins via contraction/compression of vessel segments with check valves. Disruption of lymphatic pumping can result in a condition called lymphedema with interstitial fluid accumulation. Lymphedema treatments are often ineffective, which is partially attributable to insufficient understanding of specialized lymphatic muscle lining the vessels. This muscle exhibits cardiac-like phasic contractions and smooth muscle-like tonic contractions to generate and regulate flow. To understand the relationship between this sub-cellular contractile machinery and organ-level pumping, we have developed a multiscale computational model of phasic and tonic contractions in lymphatic muscle and coupled it to a lymphangion pumping model. Our model uses the sliding filament model (Huxley in Prog Biophys Biophys Chem 7:255-318, 1957) and its adaptation for smooth muscle (Mijailovich in Biophys J 79(5):2667-2681, 2000). Multiple structural arrangements of contractile components and viscoelastic elements were trialed but only one provided physiologic results. We then coupled this model with our previous lumped parameter model of the lymphangion to relate results to experiments. We show that the model produces similar pressure, diameter, and flow tracings to experiments on rat mesenteric lymphatics. This model provides the first estimates of lymphatic muscle contraction energetics and the ability to assess the potential effects of sub-cellular level phenomena such as calcium oscillations on lymphangion outflow. The maximum efficiency value predicted (40%) is at the upper end of estimates for other muscle types. Spontaneous calcium oscillations during diastole were found to increase outflow up to approximately 50% in the range of frequencies and amplitudes tested.

A Physiology-Guided Classification of Active-Stress and Active-Strain Approaches for Continuum-Mechanical Modeling of Skeletal Muscle Tissue

The well-established sliding filament and cross-bridge theory explain the major biophysical mechanism responsible for a skeletal muscle's active behavior on a cellular level. However, the biomechanical function of skeletal muscles on the tissue scale, which is caused by the complex interplay of muscle fibers and extracellular connective tissue, is much less understood. Mathematical models provide one possibility to investigate physiological hypotheses. Continuum-mechanical models have hereby proven themselves to be very suitable to study the biomechanical behavior of whole muscles or entire limbs. Existing continuum-mechanical skeletal muscle models use either an active-stress or an active-strain approach to phenomenologically describe the mechanical behavior of active contractions. While any macroscopic constitutive model can be judged by it's ability to accurately replicate experimental data, the evaluation of muscle-specific material descriptions is difficult as suitable data is, unfortunately, currently not available. Thus, the discussions become more philosophical rather than following rigid methodological criteria. Within this work, we provide a extensive discussion on the underlying modeling assumptions of both the active-stress and the active-strain approach in the context of existing hypotheses of skeletal muscle physiology. We conclude that the active-stress approach resolves an idealized tissue transmitting active stresses through an independent pathway. In contrast, the active-strain approach reflects an idealized tissue employing an indirect, coupled pathway for active stress transmission. Finally the physiological hypothesis that skeletal muscles exhibit redundant pathways of intramuscular stress transmission represents the basis for considering a mixed-active-stress-active-strain constitutive framework.

Chronic Red Bull Consumption during Adolescence: Effect on Mesocortical and Mesolimbic Dopamine Transmission and Cardiovascular System in Adult Rats

Energy drinks are very popular nonalcoholic beverages among adolescents and young adults for their stimulant effects. Our study aimed to investigate the effect of repeated intraoral Red Bull (RB) infusion on dopamine transmission in the nucleus accumbens shell and core and in the medial prefrontal cortex and on cardiac contractility in adult rats exposed to chronic RB consumption. Rats were subjected to 4 weeks of RB voluntary consumption from adolescence to adulthood. Monitoring of in vivo dopamine was carried out by brain microdialysis. In vitro cardiac contractility was studied on biomechanical properties of isolated left-ventricular papillary muscle. The main finding of the study was that, in treated animals, RB increased shell dopamine via a nonadaptive mechanism, a pattern similar to that of drugs of abuse. No changes in isometric and isotonic mechanical parameters were associated with chronic RB consumption. However, a prolonged time to peak tension and half-time of relaxation and a slower peak rate of tension fall were observed in RB-treated rats. It is likely that RB treatment affects left-ventricular papillary muscle contraction. The neurochemical results here obtained can explain the addictive properties of RB, while the cardiovascular investigation findings suggest a hidden papillary contractility impairment.

STZ-diabetic rat heart maintains developed tension amplitude by increasing sarcomere length and crossbridge density

New Findings

What is the central question of this study?

In the papillary muscle from type I diabetic rats, does diabetes‐associated altered ventricular function result from changes of acto‐myosin interactions and are these modifications attributable to a possible sarcomere rearrangement?

What is the main finding and its importance?

For the first time, we showed that type‐I diabetes altered sarcomeric ultrastructure, as seen by transmission electron microscopy, consistent with physiological parameters. The diabetic condition induced slower timing parameters, which is compatible with a diastolic dysfunction. At the sarcomeric level, augmented β‐myosin heavy chain content and increased sarcomere length and crossbridges' number preserve myocardial stroke and could concur to maintain the ejection fraction.

Abstract

We investigated whether diabetes‐associated altered ventricular function, in a type I diabetes animal model, results from a modification of acto‐myosin interactions, through the in vitro recording of left papillary muscle mechanical parameters and examination of sarcomere morphology by transmission electron microscopy (TEM). Experiments were performed on streptozotocin‐induced diabetic and age‐matched control female Wistar rats. Mechanical isometric and isotonic indexes and timing parameters were determined. Using Huxley's equations, we calculated mechanics, kinetics and energetics of myosin crossbridges. Sarcomere length and A‐band length were measured on TEM images. Type I and III collagen and β‐myosin heavy chain (MHC) expression were determined by immunoblotting. No variation in resting and developed tension or maximum extent of shortening was evident between groups, but diabetic rats showed lower maximum shortening velocity and prolonged timing parameters. Compared to controls, diabetics also displayed a higher number of crossbridges with lower unitary force. Moreover, no change in type I and III collagen was associated to diabetes, but pathological rats showed a two‐fold enhancement of β‐MHC content and longer sarcomeres and A‐band, detected by ultrastructural morphometry. Overall, these data address whether a preserved systolic function accompanied by an altered diastolic phase results from a recruitment of super‐relaxed myosin heads or the phosphorylation of the regulatory light chain site in myosin. Although the early signs of diabetic cardiomyopathy were well expressed, the striking finding of our study was that, in diabetics, sarcomere modification may be a possible compensatory mechanism that preserves systolic function.

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.

The spatial distribution of thin filament activation influences force development and myosin activity in computational models of muscle contraction

Striated muscle contraction is initiated by Ca²⁺ binding to, and activating, thin filament regulatory units (RU) within the sarcomere, which then allows myosin cross-bridges from the opposing thick filament to bind actin and generate force. The amount of overlap between the filaments dictates how many potential cross-bridges are capable of binding, and thus how force is generated by the sarcomere. Myopathies and atrophy can impair muscle function by limiting cross-bridge interactions between the filaments, which can occur when the length of the thin filament is reduced or when RU function is disrupted. To investigate how variations in thin filament length and RU density affect ensemble cross-bridge behavior and force production, we simulated muscle contraction using a spatially explicit computational model of the half-sarcomere. Thin filament RUs were disabled either uniformly from the pointed end of the filament (to model shorter thin filament length) or randomly throughout the length of the half-sarcomere. Both uniform and random RU 'knockout' schemes decreased overall force generation during maximal and submaximal activation. The random knockout scheme also led to decreased calcium sensitivity and cooperativity of the force-pCa relationship. We also found that the rate of force development slowed with the random RU knockout, compared to the uniform RU knockout or conditions of normal RU activation. These findings imply that the relationship between RU density and force production within the sarcomere involves more complex coordination than simply the raw number of RUs available for myosin cross-bridge binding, and that the spatial pattern in which activatable RU are distributed throughout the sarcomere influences the dynamics of force production.

Control for multifunctionality: bioinspired control based on feeding in Aplysia californica

Animals exhibit remarkable feats of behavioral flexibility and multifunctional control that remain challenging for robotic systems. The neural and morphological basis of multifunctionality in animals can provide a source of bioinspiration for robotic controllers. However, many existing approaches to modeling biological neural networks rely on computationally expensive models and tend to focus solely on the nervous system, often neglecting the biomechanics of the periphery. As a consequence, while these models are excellent tools for neuroscience, they fail to predict functional behavior in real time, which is a critical capability for robotic control. To meet the need for real-time multifunctional control, we have developed a hybrid Boolean model framework capable of modeling neural bursting activity and simple biomechanics at speeds faster than real time. Using this approach, we present a multifunctional model of Aplysia californica feeding that qualitatively reproduces three key feeding behaviors (biting, swallowing, and rejection), demonstrates behavioral switching in response to external sensory cues, and incorporates both known neural connectivity and a simple bioinspired mechanical model of the feeding apparatus. We demonstrate that the model can be used for formulating testable hypotheses and discuss the implications of this approach for robotic control and neuroscience.

Mechanical contribution to muscle thin filament activation

Vertebrate striated muscle thin filaments are thought to be thermodynamically activated in response to an increase in Ca²⁺ concentration. We tested this hypothesis by measuring time intervals for gliding runs and pauses of individual skeletal muscle thin filaments in cycling myosin motility assays. A classic thermodynamic mechanism predicts that if chemical potential is constant, transitions between runs and pauses of gliding thin filaments will occur at constant rate as given by a Poisson distribution. In this scenario, rate is given by the odds of a pause, and hence, run times between pauses fit an exponential distribution that slopes negatively for all observable run times. However, we determined that relative density of observed run times fits an exponential only at low Ca²⁺ levels that activate filament gliding. Further titration with Ca²⁺, or adding excess regulatory proteins tropomyosin and troponin, shifted the relative density of short run times to fit the positive slope of a gamma distribution, which derives from waiting times between Poisson events. Events that arise during a run and prevent the chance of ending a run for a random interval of time account for the observed run time distributions, suggesting that the events originate with cycling myosin. We propose that regulatory proteins of the thin filament require the mechanical force of cycling myosin to achieve the transition state for activation. During activation, combinations of cycling myosin that contribute insufficient activation energy delay deactivation.

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.

An Equivocal Final Link - Quantitative Determination of the Thermodynamic Efficiency of ATP Hydrolysis - Sullies the Chain of Electric, Ionic, Mechanical and Metabolic Steps Underlying Cardiac Contraction

Each beat of the heart completes the final step in a sequence of events commencing with electrical excitation-triggered release of Ca2+ from the sarcoplasmic reticulum which, in turn, triggers ATP-hydrolysis-dependent mechanical contraction. Given that Thermodynamics is inherently detail-independent, the heart can be thus be viewed as a mechanical pump - the generator of pressure that drives blood through the systemic and pulmonary circulations. The beat-to-beat pressure-volume work (W) of the heart is relatively straightforward to measure experimentally. Given an ability to measure, simultaneously, the accompanying heat production or oxygen consumption, it is trivial to calculate the mechanical efficiency: ε = W/ΔH where ΔH is the change of enthalpy: (W + Q), Q representing the accompanying production of heat. But it is much less straightforward to measure the thermodynamic efficiency: η = W/ΔG ATP , where ΔG ATP signifies the Gibbs Free Energy of ATP hydrolysis. The difficulty arises because of uncertain quantification of the substrate-dependent yield of ATP - conveniently expressed as the P/O2 ratio. P/O2 ratios, originally ("classically") inferred from thermal studies, have been considerably reduced over the past several decades by re-analysis of the stoichiometric coefficients separating sequential steps in the electron transport system - in particular, dropping the requirement that the coefficients have integer values. Since the early classical values are incompatible with the more recent estimates, we aim to probe this discrepancy with a view to its reconciliation. Our probe consists of a simple, thermodynamically constrained, algebraic model of cardiac mechano-energetics. Our analysis fails to reconcile recent and classical estimates of PO2 ratios; hence, we are left with a conundrum.

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.

Myosin and Other Energy-Transducing ATPases: Structural Dynamics Studied by Electron Paramagnetic Resonance

The objective of this article was to document the energy-transducing and regulatory interactions in supramolecular complexes such as motor, pump, and clock ATPases. The dynamics and structural features were characterized by motion and distance measurements using spin-labeling electron paramagnetic resonance (EPR) spectroscopy. In particular, we focused on myosin ATPase with actin-troponin-tropomyosin, neural kinesin ATPase with microtubule, P-type ion-motive ATPase, and cyanobacterial clock ATPase. Finally, we have described the relationships or common principles among the molecular mechanisms of various energy-transducing systems and how the large-scale thermal structural transition of flexible elements from one state to the other precedes the subsequent irreversible chemical reactions.

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.

From Rough to Precise: Human-Inspired Phased Target Learning Framework for Redundant Musculoskeletal Systems

Redundant muscles in human-like musculoskeletal robots provide additional dimensions to the solution space. Consequently, the computation of muscle excitations remains an open question. Conventional methods like dynamic optimization and reinforcement learning usually have high computational costs or unstable learning processes when applied to a complex musculoskeletal system. Inspired by human learning, we propose a phased target learning framework that provides different targets to learners at varying levels, to guide their training process and to avoid local optima. By introducing an extra layer of neurons reflecting a preference, we improve the Q-network method to generate continuous excitations. In addition, based on information transmission in the human nervous system, two kinds of biological noise sources are introduced into our framework to enhance exploration over the solution space. Tracking experiments based on a simplified musculoskeletal arm model indicate that under guidance of phased targets, the proposed framework prevents divergence of excitations, thus stabilizing training. Moreover, the enhanced exploration of solutions results in smaller motion errors. The phased target learning framework can be expanded for general-purpose reinforcement learning, and it provides a preliminary interpretation for modeling the mechanisms of human motion learning.

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

Models of muscle contraction often assume that a myosin motor interacts with areas of one to three distinct binding sites along an actin filament. Månsson shows that computational models employing different numbers of sites return similar predictions, suggesting one-site models are often sufficient.

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.

A short history of the development of mathematical models of cardiac mechanics

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

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.

Thermodynamic analysis questions claims of improved cardiac efficiency by dietary fish oil

Studies in the literature describe the ability of dietary supplementation by omega-3 fish oil to increase the pumping efficiency of the left ventricle. Here we attempt to reconcile such studies with our own null results. We undertake a quantitative analysis of the improvement that could be expected theoretically, subject to physiological constraints, by posing the following question: By how much could efficiency be expected to increase if inefficiencies could be eliminated? Our approach utilizes thermodynamic analyses to investigate the contributions, both singly and collectively, of the major components of cardiac energetics to total cardiac efficiency. We conclude that it is unlikely that fish oils could achieve the required diminution of inefficiencies without greatly compromising cardiac performance.

Skeletal muscle mechanics, energetics and plasticity

The following papers by Richard Lieber (Skeletal Muscle as an Actuator), Thomas Roberts (Elastic Mechanisms and Muscle Function), Silvia Blemker (Skeletal Muscle has a Mind of its Own: a Computational Framework to Model the Complex Process of Muscle Adaptation) and Sabrina Lee (Muscle Properties of Spastic Muscle (Stroke and CP) are summaries of their representative contributions for the session on skeletal muscle mechanics, energetics and plasticity at the 2016 Biomechanics and Neural Control of Movement Conference (BANCOM 2016). Dr. Lieber revisits the topic of sarcomere length as a fundamental property of skeletal muscle contraction. Specifically, problems associated with sarcomere length non-uniformity and the role of sarcomerogenesis in diseases such as cerebral palsy are critically discussed. Dr. Roberts then makes us aware of the (often neglected) role of the passive tissues in muscles and discusses the properties of parallel elasticity and series elasticity, and their role in muscle function. Specifically, he identifies the merits of analyzing muscle deformations in three dimensions (rather than just two), because of the potential decoupling of the parallel elastic element length from the contractile element length, and reviews the associated implications for the architectural gear ratio of skeletal muscle contraction. Dr. Blemker then tackles muscle adaptation using a novel way of looking at adaptive processes and what might drive adaptation. She argues that cells do not have pre-programmed behaviors that are controlled by the nervous system. Rather, the adaptive responses of muscle fibers are determined by sub-cellular signaling pathways that are affected by mechanical and biochemical stimuli; an exciting framework with lots of potential. Finally, Dr. Lee takes on the challenging task of determining human muscle properties in vivo. She identifies the dilemma of how we can demonstrate the effectiveness of a treatment, specifically in cases of muscle spasticity following stroke or in children with cerebral palsy. She then discusses the merits of ultrasound based elastography, and the clinical possibilities this technique might hold. Overall, we are treated to a vast array of basic and clinical problems in skeletal muscle mechanics and physiology, with some solutions, and many suggestions for future research.

Force Responses and Sarcomere Dynamics of Cardiac Myofibrils Induced by Rapid Changes in [Pi]

The second phase of the biphasic force decay upon release of phosphate from caged phosphate was previously interpreted as a signature of kinetics of the force-generating step in the cross-bridge cycle. To test this hypothesis without using caged compounds, force responses and individual sarcomere dynamics upon rapid increases or decreases in concentration of inorganic phosphate [P_i] were investigated in calcium-activated cardiac myofibrils. Rapid increases in [P_i] induced a biphasic force decay with an initial slow decline (phase 1) and a subsequent 3-5-fold faster major decay (phase 2). Phase 2 started with the distinct elongation of a single sarcomere, the so-called sarcomere "give". "Give" then propagated from sarcomere to sarcomere along the myofibril. Propagation speed and rate constant of phase 2 (k_+Pi(2)) had a similar [P_i]-dependence, indicating that the kinetics of the major force decay (phase 2) upon rapid increase in [P_i] is determined by sarcomere dynamics. In contrast, no "give" was observed during phase 1 after rapid [P_i]-increase (rate constant k_+Pi(1)) and during the single-exponential force rise (rate constant k_-Pi) after rapid [P_i]-decrease. The values of k_+Pi(1) and k_-Pi were similar to the rate constant of mechanically induced force redevelopment (k_TR) and Ca²⁺-induced force development (k_ACT) measured at same [P_i]. These results indicate that the major phase 2 of force decay upon a P_i-jump does not reflect kinetics of the force-generating step but results from sarcomere "give". The other phases of P_i-induced force kinetics that occur in the absence of "give" yield the same information as mechanically and Ca²⁺-induced force kinetics (k_+Pi(1) ∼ k_-Pi ∼ k_TR ∼ k_ACT). Model simulations indicate that P_i-induced force kinetics neither enable the separation of P_i-release from the rate-limiting transition f into force states nor differentiate whether the "force-generating step" occurs before, along, or after the P_i-release.

Chemomechanical regulation of myosin Ic cross-bridges: Deducing the elastic properties of an ensemble from single-molecule mechanisms

Myosin Ic is thought to be the principal constituent of the motor that adjusts mechanical responsiveness during adaptation to prolonged stimuli by hair cells, the sensory receptors of the inner ear. In this context myosin molecules operate neither as filaments, as occurs in muscles, nor as single or few molecules, as characterizes intracellular transport. Instead, myosin Ic molecules occur in a complex cluster in which they may exhibit cooperative properties. To better understand the motor’s remarkable function, we introduce a theoretical description of myosin Ic’s chemomechanical cycle based on experimental data from recent single-molecule studies. The cycle consists of distinct chemical states that the myosin molecule stochastically occupies. We explicitly calculate the probabilities of the occupancy of these states and show their dependence on the external force, the availability of actin, and the nucleotide concentrations as required by thermodynamic constraints. This analysis highlights that the strong binding of myosin Ic to actin is dominated by the ADP state for small external forces and by the ATP state for large forces. Our approach shows how specific parameter values of the chemomechanical cycle for myosin Ic result in behaviors distinct from those of other members of the myosin family. Integrating this single-molecule cycle into a simplified ensemble description, we predict that the average number of bound myosin heads is regulated by the external force and nucleotide concentrations. The elastic properties of such an ensemble are determined by the average number of myosin cross-bridges. Changing the binding probabilities and myosin’s stiffness under a constant force results in a mechanical relaxation which is large enough to account for fast adaptation in hair cells.

Myosin molecules are biological nanomachines that transduce chemical energy into mechanical work and thus produce directed motion in living cells. These molecules proceed through cyclic reactions in which they change their conformational states upon the binding and release of nucleotides while attaching to and detaching from filaments. The myosin family consists of many distinct members with diverse functions such as muscle contraction, cargo transport, cell migration, and sensory adaptation. How these functions emerge from the biophysical properties of the individual molecules is an open question. We present an approach that integrates recent findings from single-molecule experiments into a thermodynamically consistent description of myosin Ic and demonstrate how the specific parameter values of the cycle result in a distinct function. The free variables of our description are the chemical input and external force, both of which are experimentally accessible and define the cellular environment in which these proteins function. We use this description to predict the elastic properties of an ensemble of molecules and discuss the implications for myosin Ic’s function in the inner ear as a tension regulator mediating adaptation, a hallmark of biological sensory systems. In this situation myosin molecules cooperate in an intermediate regime, neither as a large ensemble as in muscle nor as a single or a few molecules as in intracellular transport.

Minimum number of myosin motors accounting for shortening velocity under zero load in skeletal muscle

KEY POINTS

Myosin filament mechanosensing determines the efficiency of the contraction by adapting the number of switched ON motors to the load. Accordingly, the unloaded shortening velocity (V₀ ) is already set at the end of latency relaxation (LR), ∼10 ms after the start of stimulation, when the myosin filament is still in the OFF state. Here the number of actin-attached motors per half-myosin filament (n) during V₀ shortening imposed either at the end of LR or at the plateau of the isometric contraction is estimated from the relation between half-sarcomere compliance and force during the force redevelopment after shortening. The value of n decreases progressively with shortening and, during V₀ shortening starting at the end of LR, is 1-4. Reduction of n is accounted for by a constant duty ratio of 0.05 and a parallel switching OFF of motors, explaining the very low rate of ATP utilization found during unloaded shortening.

ABSTRACT

The maximum velocity at which a skeletal muscle can shorten (i.e. the velocity of sliding between the myosin filament and the actin filament under zero load, V₀ ) is already set at the end of the latency relaxation (LR) preceding isometric force generation, ∼10 ms after the start of electrical stimulation in frog muscle fibres at 4°C. At this time, Ca²⁺ -induced activation of the actin filament is maximal, while the myosin filament is in the OFF state characterized by most of the myosin motors lying on helical tracks on the filament surface, making them unavailable for actin binding and ATP hydrolysis. Here, the number of actin-attached motors per half-thick filament during V₀ shortening (n) is estimated by imposing, on tetanized single fibres from Rana esculenta (at 4°C and sarcomere length 2.15 μm), small 4 kHz oscillations and determining the relation between half-sarcomere (hs) compliance and force during the force development following V₀ shortening. When V₀ shortening is superimposed on the maximum isometric force T₀ , n decreases progressively with the increase of shortening (range 30-80 nm per hs) and, when V₀ shortening is imposed at the end of LR, n can be as low as 1-4. Reduction of n is accounted for by a constant duty ratio of the myosin motor of ∼0.05 and a parallel switching OFF of the thick filament, providing an explanation for the very low rate of ATP utilization during extended V₀ shortening.

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.

Measuring Rotations of a Few Cross-Bridges in Skeletal Muscle

The ability to measure properties of a single cross-bridge in working muscle is important because it avoids averaging the signal from a large number of molecules and because it probes cross-bridges in their native crowded environment. Because the concentration of myosin in muscle is large, observing the kinetics of a single myosin molecule requires that the signal be collected from small volumes. The introduction of small observational volumes defined by diffraction-limited laser beams and confocal detection has made it possible to limit the observational volume to a femtoliter (10 15 liter). By restraining labeling to 1 fluorophore per 100 myosin molecules, we were able to follow the kinetics of approximately 400 cross-bridges. To reduce this number further, we used two-photon (2P) microscopy. The focal plane in which the laser power density was high enough to produce 2P absorption was thinner than in confocal microscopy. Using 2P microscopy, we were able to observe approximately 200 cross-bridges during contraction. The novel method of confocal total internal reflection (CTIR) provides a method to reduce the observational volume even further, to approximately 1 attoliter (10 18 liter), and to measure fluorescence with a high signal-to-noise (S/N) ratio. In this method, the observational volume is made shallow by illuminating the sample with an evanescent field produced by total internal reflection (TIR) of the incident laser beam. To guarantee the small lateral dimensions of the observational volume, a confocal aperture is inserted in the conjugate-image plane of the objective. With a 3.5-μm confocal aperture, we achieved a volume of 1.5 attoliter. Association-dissociation of the myosin head was probed with rhodamine attached at cys707 of the heavy chain of myosin. Signal was contributed by one to five fluorescent myosin molecules. Fluorescence decayed in a series of discrete steps, corresponding to bleaching of individual molecules of rhodamine. The S/N ratio was sufficiently large to make statistically significant comparisons from rigor and contracting myofibrils.

Computational Modeling of Muscle Regeneration and Adaptation to Advance Muscle Tissue Regeneration Strategies

Skeletal muscle has an exceptional ability to regenerate and adapt following injury. Tissue engineering approaches (e.g. cell therapy, scaffolds, and pharmaceutics) aimed at enhancing or promoting muscle regeneration from severe injuries are a promising and active field of research. Computational models are beginning to advance the field by providing insight into regeneration mechanisms and therapies. In this paper, we summarize the contributions computational models have made to understanding muscle remodeling and the functional implications thereof. Next, we describe a new agent-based computational model of skeletal muscle inflammation and regeneration following acute muscle injury. Our computational model simulates the recruitment and cellular behaviors of key inflammatory cells (e.g. neutrophils and M1 and M2 macrophages) and their interactions with native muscle cells (muscle fibers, satellite stem cells, and fibroblasts) that result in the clearance of necrotic tissue and muscle fiber regeneration. We demonstrate the ability of the model to track key regeneration metrics during both unencumbered regeneration and in the case of impaired macrophage function. We also use the model to simulate regeneration enhancement when muscle is primed with inflammatory cells prior to injury, which is a putative therapeutic intervention that has not yet been investigated experimentally. Computational modeling of muscle regeneration, pursued in combination with experimental analyses, provides a quantitative framework for evaluating and predicting muscle regeneration and enables the rational design of therapeutic strategies for muscle recovery.

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.

Mechanics of Locomotion

Energetic and mechanical principles of walking and running are reviewed, using information available from force-plate studies. A mathematical model of walking is described that conserves the sum of the kinetic and gravitational potential energies of the body. In running, energy is stored transiently in the elastic deformations of stretched muscles and tendons. Theory and experiments are described using these principles and others to find the range of stiffness values for a running track that both lowers the potential for injuries and increases running speed.

Temperature effect on the chemomechanical regulation of substeps within the power stroke of a single Myosin II

Myosin IIs in the skeletal muscle are highly efficient nanoscale machines evolved in nature. Understanding how they function can not only bring insights into various biological processes but also provide guidelines to engineer synthetic nanoscale motors working in the vicinity of thermal noise. Though it was clearly demonstrated that the behavior of a skeletal muscle fiber, or that of a single myosin was strongly affected by the temperature, how exactly the temperature affects the kinetics of a single myosin is not fully understood. By adapting the newly developed transitional state model, which successfully explained the intriguing motor force regulation during skeletal muscle contraction, here we systematically explain how exactly the power stroke of a single myosin proceeds, with the consideration of the chemomechanical regulation of sub-steps within the stroke. The adapted theory is then utilized to investigate the temperature effect on various aspects of the power stroke. Our analysis suggests that, though swing rates, the isometric force, and the maximal stroke size all strongly vary with the temperature, the temperature can have a very small effect on the releasable elastic energy within the power stroke.

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.

Myosin MgADP release rate decreases at longer sarcomere length to prolong myosin attachment time in skinned rat myocardium

Cardiac contractility increases as sarcomere length increases, suggesting that intrinsic molecular mechanisms underlie the Frank-Starling relationship to confer increased cardiac output with greater ventricular filling. The capacity of myosin to bind with actin and generate force in a muscle cell is Ca(2+) regulated by thin-filament proteins and spatially regulated by sarcomere length as thick-to-thin filament overlap varies. One mechanism underlying greater cardiac contractility as sarcomere length increases could involve longer myosin attachment time (ton) due to slowed myosin kinetics at longer sarcomere length. To test this idea, we used stochastic length-perturbation analysis in skinned rat papillary muscle strips to measure ton as [MgATP] varied (0.05-5 mM) at 1.9 and 2.2 μm sarcomere lengths. From this ton-MgATP relationship, we calculated cross-bridge MgADP release rate and MgATP binding rates. As MgATP increased, ton decreased for both sarcomere lengths, but ton was roughly 70% longer for 2.2 vs. 1.9 μm sarcomere length at maximally activated conditions. These ton differences were driven by a slower MgADP release rate at 2.2 μm sarcomere length (41 ± 3 vs. 74 ± 7 s(-1)), since MgATP binding rate was not different between the two sarcomere lengths. At submaximal activation levels near the pCa50 value of the tension-pCa relationship for each sarcomere length, length-dependent increases in ton were roughly 15% longer for 2.2 vs. 1.9 μm sarcomere length. These changes in cross-bridge kinetics could amplify cooperative cross-bridge contributions to force production and thin-filament activation at longer sarcomere length and suggest that length-dependent changes in myosin MgADP release rate may contribute to the Frank-Starling relationship in the heart.

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.

Constriction model of actomyosin ring for cytokinesis by fission yeast using a two-state sliding filament mechanism

We developed a model describing the structure and contractile mechanism of the actomyosin ring in fission yeast, Schizosaccharomyces pombe. The proposed ring includes actin, myosin, and α-actinin, and is organized into a structure similar to that of muscle sarcomeres. This structure justifies the use of the sliding-filament mechanism developed by Huxley and Hill, but it is probably less organized relative to that of muscle sarcomeres. Ring contraction tension was generated via the same fundamental mechanism used to generate muscle tension, but some physicochemical parameters were adjusted to be consistent with the proposed ring structure. Simulations allowed an estimate of ring constriction tension that reproduced the observed ring constriction velocity using a physiologically possible, self-consistent set of parameters. Proposed molecular-level properties responsible for the thousand-fold slower constriction velocity of the ring relative to that of muscle sarcomeres include fewer myosin molecules involved, a less organized contractile configuration, a low α-actinin concentration, and a high resistance membrane tension. Ring constriction velocity is demonstrated as an exponential function of time despite a near linear appearance. We proposed a hypothesis to explain why excess myosin heads inhibit constriction velocity rather than enhance it. The model revealed how myosin concentration and elastic resistance tension are balanced during cytokinesis in S. pombe.

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.

Ensemble velocity of non-processive molecular motors with multiple chemical states

We study the ensemble velocity of non-processive motor proteins, described with multiple chemical states. In particular, we discuss the velocity as a function of ATP concentration. Even a simple model which neglects the strain dependence of transition rates, reverse transition rates and nonlinearities in the elasticity can show interesting functional dependencies, which deviate significantly from the frequently assumed Michaelis-Menten form. We discuss how the order of events in the duty cycle can be inferred from the measured dependence. The model also predicts the possibility of velocity reversal at a certain ATP concentration if the duty cycle contains several conformational changes of opposite directionalities.

Altered cross-bridge properties in skeletal muscle dystrophies

Force and motion generated by skeletal muscle ultimately depends on the cyclical interaction of actin with myosin. This mechanical process is regulated by intracellular Ca²⁺ through the thin filament-associated regulatory proteins i.e.; troponins and tropomyosin. Muscular dystrophies are a group of heterogeneous genetic affections characterized by progressive degeneration and weakness of the skeletal muscle as a consequence of loss of muscle tissue which directly reduces the number of potential myosin cross-bridges involved in force production. Mutations in genes responsible for skeletal muscle dystrophies (MDs) have been shown to modify the function of contractile proteins and cross-bridge interactions. Altered gene expression or RNA splicing or post-translational modifications of contractile proteins such as those related to oxidative stress, may affect cross-bridge function by modifying key proteins of the excitation-contraction coupling. Micro-architectural change in myofilament is another mechanism of altered cross-bridge performance. In this review, we provide an overview about changes in cross-bridge performance in skeletal MDs and discuss their ultimate impacts on striated muscle function.