Free energy transduction in a chemical motor model

The Problem with Inventing Molecular Mechanisms to Fit Thermodynamic Equations of Muscle

Almost every model of muscle contraction in the literature to date is a molecular power stroke model, even though this corpuscular mechanism is opposed by centuries of science, by 85 years of unrefuted evidence that muscle is a thermodynamic system, and by a quarter century of direct observations that the molecular mechanism of muscle contraction is a molecular switch, not a molecular power stroke. An ensemble of molecular switches is a binary mechanical thermodynamic system from which A.V. Hill's muscle force-velocity relationship is directly derived, where Hill's parameter a is the internal force against which unloaded muscle shortens, and Hill's parameter b is the product of the switch displacement, d, and the actin-myosin ATPase rate. Ignoring this model and the centuries of thermodynamics that preceded it, corpuscularians continue to develop molecular power stroke models, adding to their 65-year jumble of "new", "innovative", and "unconventional" molecular mechanisms for Hill's a and b parameters, none of which resemble the underlying physical chemistry. Remarkably, the corpuscularian community holds the thermodynamicist to account for these discrepancies, which, as outlined here, I have done for 25 years. It is long past time for corpuscularians to be held accountable for their mechanisms, which by all accounts have no foundation in science. The stakes are high. Molecular power stroke models are widely used in research and in clinical decision-making and have, for over half a century, muddied our understanding of the inner workings of one of the most efficient and clean-burning machines on the planet. It is problematic that corpuscularians present these models to stakeholders as science when in fact corpuscularians have been actively defending these models against science for decades. The path forward for scientists is to stop baseless rejections of muscle thermodynamics and to begin testing corpuscular and thermodynamic mechanisms with the goal of disproving one or the other of these hypotheses.

Four Phases of a Force Transient Emerge from a Binary Mechanical System

Models of muscle contraction are important for guiding drug discovery, drug validation, and clinical decision-making with the goal of improving human health. Models of muscle contraction are also key to discovering clean energy technologies from one of the most efficient and clean-burning machines on the planet. However, these important goals can only be met through muscle models that are based on science. Most every model and mechanism (e.g., a molecular power stroke) of muscle contraction described in the literature to date is based on a corpuscular mechanic philosophy that has been challenged by science for over two decades. A thermodynamic model and mechanisms (e.g., a molecular switch) of muscle contraction is supported by science but has not yet been tested against experimental data. Here, I show that following a rapid perturbation to the free energy of a thermodynamic muscle system, a transient force response emerges with four phases, each corresponding to a different clearly-defined thermodynamic (not molecular) process. I compare these four phases to those observed in two classic muscle transient experiments. The observed consistency between model and data implies that the simplest possible model of muscle contraction (a binary mechanical system) accurately describes muscle contraction.

Thermodynamics and Kinetics of a Binary Mechanical System: Mechanisms of Muscle Contraction

Biological motors function at the interface of biology, physics, and chemistry, and it remains unsettled what rules from which disciplines account for how these motors work. Myosin motors are enzymes that catalyze the hydrolysis of ATP through a mechanism involving a switch-like myosin structural change (a lever arm rotation) induced by actin binding that generates a small displacement of an actin filament. In muscle, individual myosin motors are widely assumed to function as molecular machines having mechanical properties that resemble those of muscle. In a fundamental departure from this perspective, here, I show that muscle more closely resembles a heat engine with mechanical properties that emerge from the thermodynamics of a myosin motor ensemble. The transformative impact of thermodynamics on our understanding of how a heat engine works guides a parallel transformation in our understanding of how muscle works. I consider the simplest possible model of force generation: a binary mechanical system. I develop the mechanics, energetics, and kinetics of this system and show that a single binding reaction generates force when muscle is held at a fixed length and performs work when muscle is allowed to shorten. This creates a network of thermodynamic binding pathways that resembles many of the characteristic mechanical and energetic behaviors of muscle including the muscle force-velocity relationship, heat output by shortening muscle, four phases of a muscle tension transient, spontaneous oscillatory contractions, and force redevelopment. Analogous to the thermodynamic (Carnot) cycle for a heat engine, isothermal and adiabatic binding and detachment reactions create a thermodynamic cycle for muscle that resembles cardiac pressure-volume loops (i.e., how the heart works). This paper provides an outline for how to re-interpret muscle mechanic data using thermodynamics - an ongoing effort that will continue providing novel insights into how muscle and molecular motors work.

Thermodynamic selection: mechanisms and scenarios

Thermodynamic selection is an indirect competition between agents feeding on the same energy resource and obeying the laws of thermodynamics. We examine scenarios of this selection, where the agent is modeled as a heat-engine coupled to two thermal baths and extracting work from the high-temperature bath. The agents can apply different work-extracting, game-theoretical strategies, e.g. the maximum power or the maximum efficiency. They can also have a fixed structure or be adaptive. Depending on whether the resource (i.e. the high-temperature bath) is infinite or finite, the fitness of the agent relates to the work-power or the total extracted work. These two selection scenarios lead to increasing or decreasing efficiencies of the work-extraction, respectively. The scenarios are illustrated via plant competition for sunlight, and the competition between different ATP production pathways. We also show that certain general concepts of game-theory and ecology-the prisoner's dilemma and the maximal power principle-emerge from the thermodynamics of competing agents. We emphasize the role of adaptation in developing efficient work-extraction mechanisms.

Velocity of myosin-based actin sliding depends on attachment and detachment kinetics and reaches a maximum when myosin-binding sites on actin saturate

Molecular motors such as kinesin and myosin often work in groups to generate the directed movements and forces critical for many biological processes. Although much is known about how individual motors generate force and movement, surprisingly, little is known about the mechanisms underlying the macroscopic mechanics generated by multiple motors. For example, the observation that a saturating number, N, of myosin heads move an actin filament at a rate that is influenced by actin–myosin attachment and detachment kinetics is accounted for neither experimentally nor theoretically. To better understand the emergent mechanics of actin–myosin mechanochemistry, we use an in vitro motility assay to measure and correlate the N-dependence of actin sliding velocities, actin-activated ATPase activity, force generation against a mechanical load, and the calcium sensitivity of thin filament velocities. Our results show that both velocity and ATPase activity are strain dependent and that velocity becomes maximized with the saturation of myosin-binding sites on actin at a value that is 40% dependent on attachment kinetics and 60% dependent on detachment kinetics. These results support a chemical thermodynamic model for ensemble motor mechanochemistry and imply molecularly explicit mechanisms within this framework, challenging the assumption of independent force generation.

A myosin-based mechanism for stretch activation and its possible role revealed by varying phosphate concentration in fast and slow mouse skeletal muscle fibers

Stretch activation (SA) is a delayed increase in force following a rapid muscle length increase. SA is best known for its role in asynchronous insect flight muscle, where it has replaced calcium's typical role of modulating muscle force levels during a contraction cycle. SA also occurs in mammalian skeletal muscle but has previously been thought to be too low in magnitude, relative to calcium-activated (CA) force, to be a significant contributor to force generation during locomotion. To test this supposition, we compared SA and CA force at different P_i concentrations (0-16 mM) in skinned mouse soleus (slow-twitch) and extensor digitorum longus (EDL; fast-twitch) muscle fibers. CA isometric force decreased similarly in both muscles with increasing P_i, as expected. SA force decreased with P_i in EDL (40%), leaving the SA to CA force ratio relatively constant across P_i concentrations (17-25%). In contrast, SA force increased in soleus (42%), causing a quadrupling of the SA to CA force ratio, from 11% at 0 mM P_i to 43% at 16 mM P_i, showing that SA is a significant force modulator in slow-twitch mammalian fibers. This modulation would be most prominent during prolonged muscle use, which increases P_i concentration and impairs calcium cycling. Based upon our previous Drosophila myosin isoform studies and this work, we propose that in slow-twitch fibers a rapid stretch in the presence of P_i reverses myosin's power stroke, enabling quick rebinding to actin and enhanced force production, while in fast-twitch fibers, stretch and P_i cause myosin to detach from actin.

Simple mechanics of protein machines

While belonging to the nanoscale, protein machines are so complex that tracing even a small fraction of their cycle requires weeks of calculations on supercomputers. Surprisingly, many aspects of their operation can be however already reproduced by using very simple mechanical models of elastic networks. The analysis suggests that, similar to other self-organized complex systems, functional collective dynamics in such proteins is effectively reduced to a low-dimensional attractive manifold.

An embryonic myosin isoform enables stretch activation and cyclical power in Drosophila jump muscle

The mechanism behind stretch activation (SA), a mechanical property that increases muscle force and oscillatory power generation, is not known. We used Drosophila transgenic techniques and our new muscle preparation, the jump muscle, to determine if myosin heavy chain isoforms influence the magnitude and rate of SA force generation. We found that Drosophila jump muscles show very low SA force and cannot produce positive power under oscillatory conditions at pCa 5.0. However, we transformed the jump muscle to be moderately stretch-activatable by replacing its myosin isoform with an embryonic isoform (EMB). Expressing EMB, jump muscle SA force increased by 163% and it generated net positive power. The rate of SA force development decreased by 58% with EMB expression. Power generation is Pi dependent as >4 mM Pi was required for positive power from EMB. Pi increased EMB SA force, but not wild-type SA force. Our data suggest that when muscle expressing EMB is stretched, EMB is more easily driven backward to a weakly bound state than wild-type jump muscle. This increases the number of myosin heads available to rapidly bind to actin and contribute to SA force generation. We conclude that myosin heavy chain isoforms influence both SA kinetics and SA force, which can determine if a muscle is capable of generating oscillatory power at a fixed calcium concentration.

Stochastic thermodynamics, fluctuation theorems and molecular machines

Stochastic thermodynamics as reviewed here systematically provides a framework for extending the notions of classical thermodynamics such as work, heat and entropy production to the level of individual trajectories of well-defined non-equilibrium ensembles. It applies whenever a non-equilibrium process is still coupled to one (or several) heat bath(s) of constant temperature. Paradigmatic systems are single colloidal particles in time-dependent laser traps, polymers in external flow, enzymes and molecular motors in single molecule assays, small biochemical networks and thermoelectric devices involving single electron transport. For such systems, a first-law like energy balance can be identified along fluctuating trajectories. For a basic Markovian dynamics implemented either on the continuum level with Langevin equations or on a discrete set of states as a master equation, thermodynamic consistency imposes a local-detailed balance constraint on noise and rates, respectively. Various integral and detailed fluctuation theorems, which are derived here in a unifying approach from one master theorem, constrain the probability distributions for work, heat and entropy production depending on the nature of the system and the choice of non-equilibrium conditions. For non-equilibrium steady states, particularly strong results hold like a generalized fluctuation-dissipation theorem involving entropy production. Ramifications and applications of these concepts include optimal driving between specified states in finite time, the role of measurement-based feedback processes and the relation between dissipation and irreversibility. Efficiency and, in particular, efficiency at maximum power can be discussed systematically beyond the linear response regime for two classes of molecular machines, isothermal ones such as molecular motors, and heat engines such as thermoelectric devices, using a common framework based on a cycle decomposition of entropy production.

Stochastic thermodynamics of single enzymes and molecular motors

For a single enzyme or molecular motor operating in an aqueous solution of non-equilibrated solute concentrations, a thermodynamic description is developed on the level of an individual trajectory of transitions between states. The concept of internal energy, intrinsic entropy and free energy for states follows from a microscopic description using one assumption on time scale separation. A first-law energy balance then allows the unique identification of the heat dissipated in one transition. Consistency with the second law on the ensemble level enforces both stochastic entropy as third contribution to the entropy change involved in one transition and the local detailed balance condition for the ratio between forward and backward rates for any transition. These results follow without assuming weak coupling between the enzyme and the solutes, ideal solution behavior or mass action law kinetics. The present approach highlights both the crucial role of the intrinsic entropy of each state and the physically questionable role of chemiostats for deriving the first law for molecular motors subject to an external force under realistic conditions.

The energetics of allosteric regulation of ADP release from myosin heads

Myosin molecules are involved in a wide range of transport and contractile activities in cells. A single myosin head functions through its ATPase reaction as a force generator and as a mechanosensor, and when two or more myosin heads work together in moving along an actin filament, the interplay between these mechanisms contributes to collective myosin behaviors. For example, the interplay between force-generating and force-sensing mechanisms coordinates the two heads of a myosin V molecule in its hand-over-hand processive stepping along an actin filament. In muscle, it contributes to the Fenn effect and smooth muscle latch. In both examples, a key force-sensing mechanism is the regulation of ADP release via interhead forces that are generated upon actin-myosin binding. Here we present a model describing the mechanism of allosteric regulation of ADP release from myosin heads as a change, DeltaDeltaG(-D), in the standard free energy for ADP release that results from the work, Deltamicro(mech), performed by that myosin head upon ADP release, or DeltaDeltaG(-D) = Deltamicro(mech). We show that this model is consistent with previous measurements for strain-dependent kinetics of ADP release in both myosin V and muscle myosin II. The model makes explicit the energetic cost of accelerating ADP release, showing that acceleration of ADP release during myosin V processivity requires approximately 4 kT of energy whereas the energetic cost for accelerating ADP release in a myosin II-based actin motility assay is only approximately 0.4 kT. The model also predicts that the acceleration of ADP release involves a dissipation of interhead forces. To test this prediction, we use an in vitro motility assay to show that the acceleration of ADP release from both smooth and skeletal muscle myosin II correlates with a decrease in interhead force. Our analyses provide clear energetic constraints for models of the allosteric regulation of ADP release and provide novel, testable insights into muscle and myosin V function.

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.

An accelerated state of myosin-based actin motility

We use an in vitro motility assay to determine the biochemical basis for a hypermotile state of myosin-based actin sliding. It is widely assumed that the sole biochemical determinant of actin-sliding velocities, V, is actin-myosin detachment kinetics (1/tauon), yet we recently reported that, above a critical ATP concentration of approximately 100 microM, V exceeds the detachment limit by more than 2-fold. To determine the biochemical basis for this hypermotile state, we measure the effects of ATP and inorganic phosphate, Pi, on V and observe that at low [ATP] V decreases as ln [Pi], whereas above 100 microM ATP the hypermotile V is independent of Pi. The ln [Pi] dependence of V at low [ATP] is consistent with a macroscopic model of muscle shortening, similar to Hill's contractile component, which predicts that V varies linearly with an internal force (Hill's active state) that drives actin movement against the viscous drag of myosin heads strongly bound to actin (Hill's dashpot). At high [ATP], we suggest that the hypermotile V is caused by shear thinning of the resistive population of strongly bound myosin heads. Our data and analysis indicate that, in addition to contributions from tauon and myosin's step size, d, V is influenced by the biochemistry of myosin's working step as well as resistive properties of actin and myosin.

Stochastic thermodynamics of chemical reaction networks

For chemical reaction networks in a dilute solution described by a master equation, the authors define energy and entropy on a stochastic trajectory and develop a consistent nonequilibrium thermodynamic description along a single stochastic trajectory of reaction events. A first-law like energy balance relates internal energy, applied (chemical) work, and dissipated heat for every single reaction. Entropy production along a single trajectory involves a sum over changes in the entropy of the network itself and the entropy of the medium. The latter is given by the exchanged heat identified through the first law. Total entropy production is constrained by an integral fluctuation theorem for networks arbitrarily driven by time-dependent rates and a detailed fluctuation theorem for networks in the steady state. Further exact relations such as a generalized Jarzynski relation and a generalized Clausius inequality are discussed. The authors illustrate these results for a three-species cyclic reaction network which exhibits nonequilibrium steady states as well as transitions between different steady states.

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.

Cycle kinetics, steady state thermodynamics and motors-a paradigm for living matter physics

An integration of the stochastic mathematical models for motor proteins with Hill's steady state thermodynamics yields a rather comprehensive theory for molecular motors as open systems in the nonequilibrium steady state. This theory, a natural extension of Gibbs' approach to isothermal molecular systems in equilibrium, is compared with other existing theories with dissipative structures and dynamics. The theory of molecular motors might be considered as an archetype for studying more complex open biological systems such as biochemical reaction networks inside living cells.