Dynamics of myosin-driven skeletal muscle contraction: I. Steady-state force generation

Actin Stress Fibers Response and Adaptation under Stretch

One of the many effects of soft tissues under mechanical solicitation in the cellular damage produced by highly localized strain. Here, we study the response of peripheral stress fibers (SFs) to external stretch in mammalian cells, plated onto deformable micropatterned substrates. A local fluorescence analysis reveals that an adaptation response is observed at the vicinity of the focal adhesion sites (FAs) due to its mechanosensor function. The response depends on the type of mechanical stress, from a Maxwell-type material in compression to a complex scenario in extension, where a mechanotransduction and a self-healing process takes place in order to prevent the induced severing of the SF. A model is proposed to take into account the effect of the applied stretch on the mechanics of the SF, from which relevant parameters of the healing process are obtained. In contrast, the repair of the actin bundle occurs at the weak point of the SF and depends on the amount of applied strain. As a result, the SFs display strain-softening features due to the incorporation of new actin material into the bundle. In contrast, the response under compression shows a reorganization with a constant actin material suggesting a gliding process of the SFs by the myosin II motors.

Five-Site Model for Brownian Dynamics Simulations of a Molecular Walker in Three Dimensions

Motor proteins play an important role in many biological processes and have inspired the development of synthetic analogues. Molecular walkers, such as kinesin, dynein, and myosin V, fulfill a diverse set of functions including transporting cargo along tracks, pulling molecules through membranes, and deforming fibers. The complexity of molecular motors and their environment makes it difficult to model the detailed dynamics of molecular walkers over long time scales. In this work, we present a simple, three-dimensional model for a molecular walker on a bead-spring substrate. The walker is represented by five spherically symmetric particles that interact through common intermolecular potentials and can be simulated efficiently in Brownian dynamics simulations. The movement of motor protein walkers entails energy conversion through ATP hydrolysis while artificial motors typically rely on a local conversion of energy supplied through external fields. We model energy conversion through rate equations for mechanochemical states that couple positional and chemical degrees of freedom and determine the walker conformation through interaction potential parameters. We perform Brownian dynamics simulations for two scenarios: In the first, the model walker transports cargo by walking on a substrate whose ends are fixed. In the second, a tethered motor pulls a mobile substrate chain against a variable force. We measure relative displacements and determine the effects of cargo size and retarding force on the efficiency of the walker. We find that, while the efficiency of our model walker is less than for the biological system, our simulations reproduce trends observed in single-molecule experiments on kinesin. In addition, the model and simulation method presented here can be readily adapted to biological and synthetic systems with multiple walkers.

Non-uniform distribution of myosin-mediated forces governs red blood cell membrane curvature through tension modulation

The biconcave disk shape of the mammalian red blood cell (RBC) is unique to the RBC and is vital for its circulatory function. Due to the absence of a transcellular cytoskeleton, RBC shape is determined by the membrane skeleton, a network of actin filaments cross-linked by spectrin and attached to membrane proteins. While the physical properties of a uniformly distributed actin network interacting with the lipid bilayer membrane have been assumed to control RBC shape, recent experiments reveal that RBC biconcave shape also depends on the contractile activity of nonmuscle myosin IIA (NMIIA) motor proteins. Here, we use the classical Helfrich-Canham model for the RBC membrane to test the role of heterogeneous force distributions along the membrane and mimic the contractile activity of sparsely distributed NMIIA filaments. By incorporating this additional contribution to the Helfrich-Canham energy, we find that the RBC biconcave shape depends on the ratio of forces per unit volume in the dimple and rim regions of the RBC. Experimental measurements of NMIIA densities at the dimple and rim validate our prediction that (a) membrane forces must be non-uniform along the RBC membrane and (b) the force density must be larger in the dimple than the rim to produce the observed membrane curvatures. Furthermore, we predict that RBC membrane tension and the orientation of the applied forces play important roles in regulating this force-shape landscape. Our findings of heterogeneous force distributions on the plasma membrane for RBC shape maintenance may also have implications for shape maintenance in different cell types.

The spectrin-actin network of the membrane skeleton plays an important role in controlling specialized cell membrane morphology. In the paradigmatic red blood cell (RBC), where actin filaments are present exclusively in the membrane skeleton, recent experiments reveal that nonmuscle myosin IIA (NMIIA) motor contractility maintains the RBC biconcave disk shape. In this study, we have identified criteria for micron-scale distributions of NMIIA forces at the membrane required to maintain the biconcave disk shape of an RBC in the resting condition. Supported by experimental measurements of RBC NMIIA distribution, we showed that a heterogeneous force distribution with a larger force density at the dimple is able to capture the experimentally observed biconcave morphology of an RBC with better accuracy compared to previous models that did not consider the heterogeneity in the force distribution. Furthermore, we showed that the biconcave geometry of the RBC is closely regulated by the effective membrane tension and the direction of applied forces on the membrane. These findings can be generalized to any force-mediated membrane shape, providing insight into the role of actomyosin forces in prescribing and maintaining the morphology of different cell types.

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.

Acto-myosin force organization modulates centriole separation and PLK4 recruitment to ensure centriole fidelity

The presence of aberrant number of centrioles is a recognized cause of aneuploidy and hallmark of cancer. Hence, centriole duplication needs to be tightly regulated. It has been proposed that centriole separation limits centrosome duplication. The mechanism driving centriole separation is poorly understood and little is known on how this is linked to centriole duplication. Here, we propose that actin-generated forces regulate centriole separation. By imposing geometric constraints via micropatterns, we were able to prove that precise acto-myosin force arrangements control direction, distance and time of centriole separation. Accordingly, inhibition of acto-myosin contractility impairs centriole separation. Alongside, we observed that organization of acto-myosin force modulates specifically the length of S-G2 phases of the cell cycle, PLK4 recruitment at the centrosome and centriole fidelity. These discoveries led us to suggest that acto-myosin forces might act in fundamental mechanisms of aneuploidy prevention.

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.

Sensitivity of small myosin II ensembles from different isoforms to mechanical load and ATP concentration

Based on a detailed crossbridge model for individual myosin II motors, we systematically study the influence of mechanical load and adenosine triphosphate (ATP) concentration on small myosin II ensembles made from different isoforms. For skeletal and smooth muscle myosin II, which are often used in actomyosin gels that reconstitute cell contractility, fast forward movement is restricted to a small region of phase space with low mechanical load and high ATP concentration, which is also characterized by frequent ensemble detachment. At high load, these ensembles are stalled or move backwards, but forward motion can be restored by decreasing ATP concentration. In contrast, small ensembles of nonmuscle myosin II isoforms, which are found in the cytoskeleton of nonmuscle cells, are hardly affected by ATP concentration due to the slow kinetics of the bound states. For all isoforms, the thermodynamic efficiency of ensemble movement increases with decreasing ATP concentration, but this effect is weaker for the nonmuscle myosin II isoforms.

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.

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.

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.

Dynamics of the bacterial flagellar motor: the effects of stator compliance, back steps, temperature, and rotational asymmetry

The rotation of a bacterial flagellar motor (BFM) is driven by multiple stators tethered to the cell wall. Here, we extend a recently proposed power-stroke model to study the BFM dynamics under different biophysical conditions. Our model explains several key experimental observations and reveals their underlying mechanisms. 1), The observed independence of the speed at low load on the number of stators is explained by a force-dependent stepping mechanism that is independent of the strength of the stator tethering spring. Conversely, without force-dependent stepping, an unrealistically weak stator spring is required. 2), Our model with back-stepping naturally explains the observed absence of a barrier to backward rotation. Using the same set of parameters, it also explains BFM behaviors in the high-speed negative-torque regime. 3), From the measured temperature dependence of the maximum speed, our model shows that stator-stepping is a thermally activated process with an energy barrier. 4), The recently observed asymmetry in the torque-speed curve between counterclockwise- and clockwise-rotating BFMs can be quantitatively explained by the asymmetry in the stator-rotor interaction potentials, i.e., a quasilinear form for the counterclockwise motor and a quadratic form for the clockwise motor.

Quantification of the forces driving self-assembly of three-dimensional microtissues

In a nonadhesive environment, cells will self-assemble into microtissues, a process relevant to tissue engineering. Although this has been recognized for some time, there is no basis for quantitative characterization of this complex process. Here we describe a recently developed assay designed to quantify aspects of the process and discuss its application in comparing behaviors between cell types. Cells were seeded in nonadhesive micromolded wells, each well with a circular trough at its base formed by the cylindrical sidewalls and by a central peg in the form of a right circular cone. Cells settled into the trough and coalesced into a toroid, which was then driven up the conical peg by the forces of self-assembly. The mass of the toroid and its rate of upward movement were used to calculate the cell power expended in the process against gravity. The power of the toroid was found to be 0.31 ± 0.01 pJ/h and 4.3 ± 1.7 pJ/h for hepatocyte cells and fibroblasts, respectively. Blocking Rho kinase by means of Y-27632 resulted in a 50% and greater reduction in power expended by each type of toroid, indicating that cytoskeletal-mediated contraction plays a significant role in the self-assembly of both cell types. Whereas the driving force for self-assembly has often been viewed as the binding of surface proteins, these data show that cellular contraction is important for cell-cell adhesion. The power measurement quantifies the contribution of cell contraction, and will be useful for understanding the concerted action of the mechanisms that drive self-assembly.

Two kinesins transport cargo primarily via the action of one motor: implications for intracellular transport

The number of microtubule motors attached to vesicles, organelles, and other subcellular commodities is widely believed to influence their motile properties. There is also evidence that cells regulate intracellular transport by tuning the number and/or ratio of motor types on cargos. Yet, the number of motors responsible for cargo motion is not easily characterized, and the extent to which motor copy number affects intracellular transport remains controversial. Here, we examined the load-dependent properties of structurally defined motor assemblies composed of two kinesin-1 molecules. We found that a group of kinesins can produce forces and move with velocities beyond the abilities of single kinesin molecules. However, such capabilities are not typically harnessed by the system. Instead, two-kinesin assemblies adopt a range of microtubule-bound configurations while transporting cargos against an applied load. The binding arrangement of motors on their filament dictates how loads are distributed within the two-motor system, which in turn influences motor-microtubule affinities. Most configurations promote microtubule detachment and prevent both kinesins from contributing to force production. These results imply that cargos will tend to be carried by only a fraction of the total number of kinesins that are available for transport at any given time, and provide an alternative explanation for observations that intracellular transport depends weakly on kinesin number in vivo.

A model for cell motility on soft bio-adhesive substrates

Mechanical stiffness of bio-adhesive substrates has been recognized as a major regulator of cell motility. We present a simple physical model to study the crawling locomotion of a contractile cell on a soft elastic substrate. The mechanism of rigidity sensing is accounted for using Schwarz's two-spring model Schwarz et al. (2006). The predicted dependency between the speed of motility and substrate stiffness is qualitatively consistent with experimental observations. The model demonstrates that the rigidity dependent motility of cells is rooted in the regulation of actomyosin contractile forces by substrate deformation at each anchorage point. On stiffer substrates, the traction forces required for cell translocation acquire larger magnitude but show weaker asymmetry which leads to slower cell motility. On very soft substrates, the model predicts a biphasic relationship between the substrate rigidity and the speed of locomotion, over a narrow stiffness range, which has been observed experimentally for some cell types.

A macroscopic ansatz to deduce the Hill relation

In this study, we derive the hyperbolic force-velocity relation of concentric muscular contraction, first formulated empirically by A.V. Hill in 1938, from three essential model assumptions: (1) the structural assembly of three well-known elements - i.e. active, parallel damping, and serial - fulfilling a force equilibrium, (2) the parallel damping coefficient explicitly depending on muscle force output and three parameters, and (3) the kinematic gearing ratio between active and serial element being assigned to a parameter. The energy source within the muscle represented by the force of the active element is an additional fifth parameter. As a result we find the Hill "constants" A and B as functions of our five model parameters. Using A and B values from literature on experimental data, we predict heat power release of our model. By calculating enthalpy rate and mechanical efficiency, we compare the model heat power to predictions from another Hill-type model, to Hill's original findings, and to findings from modern muscle heat measurements. We reconsider why the biggest share of heat rate during isometric contractions (maintenance heat) and the velocity-dependent heat rate during concentric contractions in addition to maintenance heat rate (shortening heat rate) may be traced back to the same mechanism represented by the kinematic gearing ratio. Namely, we suggest that the serial element transfers attachment-detachment fluctuations of actin-myosin crossbridges within one sarcomere to others in the same sarcomere and to those in parallel and in series. Numerically, in case of negligible passive muscular damping, we find the ratio between A and isometric force (relative A) to depend exclusively on the kinematic gearing ratio, whereas the maintenance heat rate scales with the square of relative A. Moreover, this mechanical coupling internal to the muscle fibres may also be behind the macroscopic force dependency of the overall parallel damping coefficient.

Cooperative molecular motors moving back and forth

We use a two-state ratchet model to study the cooperative bidirectional motion of molecular motors on cytoskeletal tracks with randomly alternating polarities. Our model is based on a previously proposed model [Badoual, Proc. Natl. Acad. Sci. U.S.A. 99, 6696 (2002)] for collective motor dynamics and, in addition, takes into account the cooperativity effect arising from the elastic tension that develops in the cytoskeletal track due to the joint action of the walking motors. We show, both computationally and analytically, that this additional cooperativity effect leads to a dramatic reduction in the characteristic reversal time of the bidirectional motion, especially in systems with a large number of motors. We also find that bidirectional motion takes place only on (almost) apolar tracks, while on even slightly polar tracks the cooperative motion is unidirectional. We argue that the origin of these observations is the sensitive dependence of the cooperative dynamics on the difference between the number of motors typically working in and against the instantaneous direction of motion.

Hysteresis in cross-bridge models of muscle

A dynamical system is said to exhibit hysteresis if its current state depends on its history. Muscle shows hysteretic properties at constant length, such as residual force enhancement after stretch. There is no generally accepted explanation for residual force enhancement. Here we examine a very simple kinetic model for the interaction between actin and myosin, the two main proteins involved in muscle contraction. We demonstrate that this model shows hysteresis at constant force. Since muscle is not a continuum but rather a group of repeating elements, called sarcomeres, arranged in series, we perform simulations of three sarcomeres. These simulations show hysteresis at constant length. This result is the first time that residual force enhancement has been demonstrated using an experimentally motivated kinetic model and multi-sarcomere simulations without passive elastic elements, damping and/or force-length relationships. We conclude by suggesting some experiments to test the model's predictions. If these experiments support the model, it becomes important to understand multiple sarcomere systems, since their behavior may be very different from most current simulations that neglect the coupling between sarcomeres.

Mechanistic role of movement and strain sensitivity in muscle contraction

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

Skeletal muscle performance determined by modulation of number of myosin motors rather than motor force or stroke size

Skeletal muscle can bear a high load at constant length, or shorten rapidly when the load is low. This force-velocity relationship is the primary determinant of muscle performance in vivo. Here we exploited the quasi-crystalline order of myosin II motors in muscle filaments to determine the molecular basis of this relationship by X-ray interference and mechanical measurements on intact single cells. We found that, during muscle shortening at a wide range of velocities, individual myosin motors maintain a force of about 6 pN while pulling an actin filament through a 6 nm stroke, then quickly detach when the motor reaches a critical conformation. Thus we show that the force-velocity relationship is primarily a result of a reduction in the number of motors attached to actin in each filament in proportion to the filament load. These results explain muscle performance and efficiency in terms of the molecular mechanism of the myosin motor.

High-frequency oscillations as a consequence of neglected serial damping in Hill-type muscle models

High-frequency vibrations e.g., induced by legs impacting with the ground during terrestrial locomotion can provoke damage within tendons even leading to ruptures. So far, macroscopic Hill-type muscle models do not account for the observed high-frequency damping at low-amplitudes. Therefore, former studies proposed that protective damping might be explained by modelling the contractile machinery of the muscles in more detail, i.e., taking the microscopic processes of the actin-myosin coupling into account. In contrast, this study formulates an alternative hypothesis: low but significant damping of the passive material in series to the contractile machinery--e.g., tendons, aponeuroses, titin--may well suffice to damp these hazardous vibrations. Thereto, we measured the contraction dynamics of a piglet muscle-tendon complex (MTC) in three contraction modes at varying loads and muscle-tendon lengths. We simulated all three respective load situations on a computer: a Hill-type muscle model including a contractile element (CE) and each an elastic element in parallel (PEE) and in series (SEE) to the CE pulled on a loading mass. By comparing the model to the measured output of the MTC, we extracted a consistent set of muscle parameters. We varied the model by introducing either linear damping in parallel or in series to the CE leading to accordant re-formulations of the contraction dynamics of the CE. The comparison of the three cases (no additional damping, parallel damping, serial damping) revealed that serial damping at a physiological magnitude suffices to explain damping of high-frequency vibrations of low amplitudes. The simulation demonstrates that any undamped serial structure within the MTC enforces SEE-load eigenoscillations. Consequently, damping must be spread all over the MTC, i.e., rather has to be de-localised than localised within just the active muscle material. Additionally, due to suppressed eigenoscillations Hill-type muscle models taking into account serial damping are numerically more efficient when used in macroscopic biomechanical neuro-musculo-skeletal models.

Mechanochemical coupling in the myosin motor domain. I. Insights from equilibrium active-site simulations

Although the major structural transitions in molecular motors are often argued to couple to the binding of Adenosine triphosphate (ATP), the recovery stroke in the conventional myosin has been shown to be dependent on the hydrolysis of ATP. To obtain a clearer mechanistic picture for such "mechanochemical coupling" in myosin, equilibrium active-site simulations with explicit solvent have been carried out to probe the behavior of the motor domain as functions of the nucleotide chemical state and conformation of the converter/relay helix. In conjunction with previous studies of ATP hydrolysis with different active-site conformations and normal mode analysis of structural flexibility, the results help establish an energetics-based framework for understanding the mechanochemical coupling. It is proposed that the activation of hydrolysis does not require the rotation of the lever arm per se, but the two processes are tightly coordinated because both strongly couple to the open/close transition of the active site. The underlying picture involves shifts in the dominant population of different structural motifs as a consequence of changes elsewhere in the motor domain. The contribution of this work and the accompanying paper [] is to propose the actual mechanism behind these "population shifts" and residues that play important roles in the process. It is suggested that structural flexibilities at both the small and large scales inherent to the motor domain make it possible to implement tight couplings between different structural motifs while maintaining small free-energy drops for processes that occur in the detached states, which is likely a feature shared among many molecular motors. The significantly different flexibility of the active site in different X-ray structures with variable level arm orientations supports the notation that external force sensed by the lever arm may transmit into the active site and influence the chemical steps (nucleotide hydrolysis and/or binding).

Mechanism of tension generation in muscle: an analysis of the forward and reverse rate constants

Tension generation in muscle occurs during the attached phase of the ATP-powered cyclic interaction of myosin heads with thin filaments. The transient nature of tension-generating intermediates and the complexity of the mechanochemical cross-bridge cycle have impeded a quantitative description of tension generation. Recent experiments performed under special conditions yielded a sigmoidal dependence of fiber tension on temperature--a unique case that simplifies the system to a two-state transition. We have applied this two-state analysis to kinetic data obtained from biexponential laser temperature-jump tension transients. Here we present the forward and reverse rate constants for de novo tension generation derived from analysis of the kinetics of the fast laser temperature-jump phase tau(2) (equivalent of the length-jump phase 2(slow)). The slow phase tau(3) is temperature-independent indicating coupling to rather than a direct role in, de novo tension generation. Increasing temperature accelerates the forward, and slows the reverse, rate constant for the creation of the tension-generating state. Arrhenius behavior of the forward and anti-Arrhenius behavior of the reverse rate constant is a kinetic signature of multistate multipathway protein-folding reactions. We conclude that locally unfolded tertiary and/or secondary structure of the actomyosin cross-bridge mediates the power stroke.

Mechanism of the Frank–Starling law—A simulation study with a novel cardiac muscle contraction model that includes titin and troponin I

A stretch-induced increase of active tension is one of the most important properties of the heart, known as the Frank-Starling law. Although a variation of myofilament Ca(2+) sensitivity with sarcomere length (SL) change was found to be involved, the underlying molecular mechanisms are not fully clarified. Some recent experimental studies indicate that a reduction of the lattice spacing between thin and thick filaments, through the increase of passive tension caused by the sarcomeric protein titin with an increase in SL within the physiological range, promotes formation of force-generating crossbridges (Xbs). However, the mechanism by which the Xb concentration determines the degree of cooperativity for a given SL has so far evaded experimental elucidation. In this simulation study, a novel, rather simple molecular-based cardiac contraction model, appropriate for integration into a ventricular cell model, was designed, being the first model to introduce experimental data on titin-based radial tension to account for the SL-dependent modulation of the interfilament lattice spacing and to include a conformational change of troponin I (TnI). Simulation results for the isometric twitch contraction time course, the length-tension and the force-[Ca(2+)] relationships are comparable to experimental data. A complete potential Frank-Starling mechanism was analyzed by this simulation study. The SL-dependent modulation of the myosin binding rate through titin's passive tension determines the Xb concentration which then alters the degree of positive cooperativity affecting the rate of the TnI conformation change and causing the Hill coefficient to be SL-dependent.

Analytical theory of the stochastic dynamics of the power stroke in nonprocessive motor proteins

Statistical distributions of the structural states of individual molecules of nonprocessive motor complexes such as actomyosins are examined theoretically by considering a two-state stochastic model coupled by chemical reactions along the reaction coordinate representing the internal conformational states of the motor. The use of a conformational reaction coordinate allows for the approximation of taking the rate constants as local in their dependence on the reaction coordinate, and yields a simple analytic solution of the stationary states. The approximation is also tested against numerical solutions with a nonlocal form of rate constants. The theory is well-suited for computational treatments based on atomic structures of protein constituents using free energy molecular dynamics simulations. With empirical sets of free energy functions, stationary distributions of forces exerted by a motor head compare well with known experimental data.

Load-dependent kinetics of myosin-V can explain its high processivity

Recent studies provide strong evidence that single myosin class V molecules transport vesicles and organelles processively along F-actin, taking several 36-nm steps, 'hand over hand', for each diffusional encounter. The mechanisms regulating myosin-V's processivity remain unknown. Here, we have used an optical-tweezers-based transducer to measure the effect of load on the mechanical interactions between rabbit skeletal F-actin and a single head of mouse brain myosin-V, which produces its working stroke in two phases. We found that the lifetimes of the first phase of the working stroke changed exponentially and about 10-fold over a range of pushing and pulling forces of +/- 1.5 pN. Stiffness measurements suggest that intramolecular forces could approach 3.6 pN when both heads are bound to F-actin, in which case extrapolation would predict the detachment kinetics of the front head to slow down 50-fold and the kinetics of the rear head to accelerate respectively. This synchronizing effect on the chemo-mechanical cycles of the heads increases the probability of the trail head detaching first and causes a strong increase in the number of forward steps per diffusional encounter over a system with no strain dependence.