Erratum to: Identification of functional differences between recombinant human α and β cardiac myosin motors

Modeling the effects of thin filament near-neighbor cooperative interactions in mammalian myocardium

The mechanisms underlying cooperative activation and inactivation of myocardial force extend from local, near-neighbor interactions involving troponin-tropomyosin regulatory units (RU) and crossbridges (XB) to more global interactions across the sarcomere. To better understand these mechanisms in the hearts of small and large mammals, we undertook a simplified mathematical approach to assess the contribution of three types of near-neighbor cooperative interactions, i.e., RU-induced, RU-activation (RU-RU), crossbridge-induced, crossbridge-binding (XB-XB), and XB-induced, RU-activation (XB-RU). We measured the Ca2+ and activation dependence of the rate constant of force redevelopment in murine- and porcine-permeabilized ventricular myocardium. Mathematical modeling of these three near-neighbor interactions yielded nonlinear expressions for the RU-RU and XB-RU rate coefficients (kon and koff) and XB-XB rate coefficients describing the attachment of force-generating crossbridges (f and f'). The derivation of single cooperative coefficient parameters (u = RU-RU, w = XB-RU, and v = XB-XB) permitted an initial assessment of the strength of each near-neighbor interaction. The parameter sets describing the effects of discrete XB-XB or XB-RU interactions failed to adequately fit the in vitro contractility data in either murine or porcine myocardium. However, the Ca2+ dependence of ktr in murine and porcine ventricular myocardium was well fit by parameter sets incorporating the RU-RU cooperative interaction. Our results indicate that a significantly stronger RU-RU interaction is present in porcine ventricular myocardium compared with murine ventricular myocardium and that the relative strength of the near-neighbor RU-RU interaction contributes to species-specific myocardial contractile dynamics in small and large mammals.

Functional control of myosin motors in the cardiac cycle

Contraction of the heart is driven by cyclical interactions between myosin and actin filaments powered by ATP hydrolysis. The modular structure of heart muscle and the organ-level synchrony of the heartbeat ensure tight reciprocal coupling between this myosin ATPase cycle and the macroscopic cardiac cycle. The myosin motors respond to the cyclical activation of the actin and myosin filaments to drive the pressure changes that control the inflow and outflow valves of the heart chambers. Opening and closing of the valves in turn switches the myosin motors between roughly isometric and roughly isotonic contraction modes. Peak filament stress in the heart is much smaller than in fully activated skeletal muscle, although the myosin filaments in the two muscle types have the same number of myosin motors. Calculations indicate that only ~5% of the myosin motors in the heart are needed to generate peak systolic pressure, although many more motors are needed to drive ejection. Tight regulation of the number of active motors is essential for the efficient functioning of the healthy heart - this control is commonly disrupted by gene variants associated with inherited heart disease, and its restoration might be a useful end point in the development of novel therapies.

Motor domain phosphorylation increases nucleotide exchange and turns MYO6 into a faster and stronger motor

Myosin motors perform many fundamental functions in eukaryotic cells by providing force generation, transport or tethering capacity. Motor activity control within the cell involves on/off switches, however, few examples are known of how myosins regulate speed or processivity and fine-tune their activity to a specific cellular task. Here, we describe a phosphorylation event for myosins of class VI (MYO6) in the motor domain, which accelerates its ATPase activity leading to a 4-fold increase in motor speed determined by actin-gliding assays, single molecule mechanics and stopped flow kinetics. We demonstrate that the serine/threonine kinase DYRK2 phosphorylates MYO6 at S267 in vitro. Single-molecule optical-tweezers studies at low load reveal that S267-phosphorylation results in faster nucleotide-exchange kinetics without change in the working stroke of the motor. The selective increase in stiffness of the acto-MYO6 complex when proceeding load-dependently into the nucleotide-free rigor state demonstrates that S267-phosphorylation turns MYO6 into a stronger motor. Finally, molecular dynamic simulations of the nucleotide-free motor reveal an alternative interaction network within insert-1 upon phosphorylation, suggesting a molecular mechanism, which regulates insert-1 positioning, turning the S267-phosphorylated MYO6 into a faster motor.

A novel kinetic model to demonstrate the independent effects of ATP and ADP/Pi concentrations on sarcomere function

Understanding muscle contraction mechanisms is a standing challenge, and one of the approaches has been to create models of the sarcomere-the basic contractile unit of striated muscle. While these models have been successful in elucidating many aspects of muscle contraction, they fall short in explaining the energetics of functional phenomena, such as rigor, and in particular, their dependence on the concentrations of the biomolecules involved in the cross-bridge cycle. Our hypothesis posits that the stochastic time delay between ATP adsorption and ADP/Pi release in the cross-bridge cycle necessitates a modeling approach where the rates of these two reaction steps are controlled by two independent parts of the total free energy change of the hydrolysis reaction. To test this hypothesis, we built a two-filament, stochastic-mechanical half-sarcomere model that separates the energetic roles of ATP and ADP/Pi in the cross-bridge cycle's free energy landscape. Our results clearly demonstrate that there is a nontrivial dependence of the cross-bridge cycle's kinetics on the independent concentrations of ATP, ADP, and Pi. The simplicity of the proposed model allows for analytical solutions of the more basic systems, which provide novel insight into the dominant mechanisms driving some of the experimentally observed contractile phenomena.

Considering impact of age and sex on cardiac cytoskeletal components

The cardiac cytoskeletal components are integral to cardiomyocyte function and are responsible for contraction, sustaining cell structure, and providing scaffolding to direct signaling. Cytoskeletal components have been implicated in cardiac pathology; however, less attention has been paid to age-related modifications of cardiac cytoskeletal components and how these contribute to dysfunction with increased age. Moreover, significant sex differences in cardiac aging have been identified, but we still lack a complete understanding to the mechanisms behind these differences. This review summarizes what is known about how key cardiomyocyte cytoskeletal components are modified because of age, as well as reported sex-specific differences. Thorough consideration of both age and sex as integral players in cytoskeletal function may reveal potential avenues for more personalized therapeutics.

Protocols for Myosin and Actin-Myosin Assays Using Rapid, Stopped-Flow Kinetics

Fast transient kinetics using stopped-flow fluorimetry is now a powerful method for defining the ATPase cycle of myosin and its subfragments and has found wide use in defining the difference between myosin isoforms, myosins carrying disease linked mutations, and the effect of small molecules on the ATPase cycle. Here the protocols for completing the classical assays of myosin and actin.myosin using the stopped-flow are described. The assays include ATP and ADP binding to myosin and actin.myosin, displacement of ADP from myosin and actin.myosin, and the cleavage of ATP to ADP and phosphate on myosin. Single and multiple turnover assays are also described.

Cooperative mechanisms underlie differences in myocardial contractile dynamics between large and small mammals

Ca2+ binding to troponin C (TnC) and myosin cross-bridge binding to actin act in a synergistic cooperative manner to modulate myocardial contraction and relaxation. The responsiveness of the myocardial thin filament to the activating effects of Ca2+ and myosin cross-bridge binding has been well-characterized in small mammals (e.g., mice). Given the nearly 10-fold difference in resting heart rates and twitch kinetics between small and large mammals, it is unlikely that the cooperative mechanisms underlying thin filament activation are identical in these two species. To test this idea, we measured the Ca2+ dependencies of steady-state force and the rate constant of force redevelopment (ktr) in murine and porcine permeabilized ventricular myocardium. While murine myocardium exhibited a steep activation-dependence of ktr, the activation-dependent profile of ktr was significantly reduced in porcine ventricular myocardium. Further insight was attained by examining force-pCa and ktr-pCa relationships. In the murine myocardium, the pCa50 for ktr was right-shifted compared with the pCa50 for force, meaning that increases in steady-state force occurred well before increases in the rate of force redevelopment were observed. In the porcine myocardium, we observed a tighter coupling of the force-pCa and ktr-pCa relationships, as evidenced by near-maximal rates of force redevelopment at low levels of Ca2+ activation. These results demonstrate that the molecular mechanisms underlying the cooperative activation of force are a dynamic property of the mammalian heart, involving, at least in part, the species- and tissue-specific expression of cardiac myosin heavy chain isoforms.

Effect of Myosin Isoforms on Cardiac Muscle Twitch of Mice, Rats and Humans

To understand how pathology-induced changes in contractile protein isoforms modulate cardiac muscle function, it is necessary to quantify the temporal-mechanical properties of contractions that occur under various conditions. Pathological responses are much easier to study in animal model systems than in humans, but extrapolation between species presents numerous challenges. Employing computational approaches can help elucidate relationships that are difficult to test experimentally by translating the observations from rats and mice, as model organisms, to the human heart. Here, we use the spatially explicit MUSICO platform to model twitch contractions from rodent and human trabeculae collected in a single laboratory. This approach allowed us to identify the variations in kinetic characteristics of α- and β-myosin isoforms across species and to quantify their effect on cardiac muscle contractile responses. The simulations showed how the twitch transient varied with the ratio of the two myosin isoforms. Particularly, the rate of tension rise was proportional to the fraction of α-myosin present, while the β-isoform dominated the rate of relaxation unless α-myosin was >50%. Moreover, both the myosin isoform and the Ca2+ transient contributed to the twitch tension transient, allowing two levels of regulation of twitch contraction.

A reverse stroke characterizes the force generation of cardiac myofilaments, leading to an understanding of heart function

Changes in the molecular properties of cardiac myosin strongly affect the interactions of myosin with actin that result in cardiac contraction and relaxation. However, it remains unclear how myosin molecules work together in cardiac myofilaments and which properties of the individual myosin molecules impact force production to drive cardiac contractility. Here, we measured the force production of cardiac myofilaments using optical tweezers. The measurements revealed that stepwise force generation was associated with a higher frequency of backward steps at lower loads and higher stall forces than those of fast skeletal myofilaments. To understand these unique collective behaviors of cardiac myosin, the dynamic responses of single cardiac and fast skeletal myosin molecules, interacting with actin filaments, were evaluated under load. The cardiac myosin molecules switched among three distinct conformational positions, ranging from pre- to post-power stroke positions, in 1 mM ADP and 0 to 10 mM phosphate solution. In contrast to cardiac myosin, fast skeletal myosin stayed primarily in the post-power stroke position, suggesting that cardiac myosin executes the reverse stroke more frequently than fast skeletal myosin. To elucidate how the reverse stroke affects the force production of myofilaments and possibly heart function, a simulation model was developed that combines the results from the single-molecule and myofilament experiments. The results of this model suggest that the reversal of the cardiac myosin power stroke may be key to characterizing the force output of cardiac myosin ensembles and possibly to facilitating heart contractions.

Identification of sequence changes in myosin II that adjust muscle contraction velocity

The speed of muscle contraction is related to body size; muscles in larger species contract at slower rates. Since contraction speed is a property of the myosin isoform expressed in a muscle, we investigated how sequence changes in a range of muscle myosin II isoforms enable this slower rate of muscle contraction. We considered 798 sequences from 13 mammalian myosin II isoforms to identify any adaptation to increasing body mass. We identified a correlation between body mass and sequence divergence for the motor domain of the 4 major adult myosin II isoforms (β/Type I, IIa, IIb, and IIx), suggesting that these isoforms have adapted to increasing body mass. In contrast, the non-muscle and developmental isoforms show no correlation of sequence divergence with body mass. Analysis of the motor domain sequence of β-myosin (predominant myosin in Type I/slow and cardiac muscle) from 67 mammals from 2 distinct clades identifies 16 sites, out of 800, associated with body mass (padj < 0.05) but not with the clade (padj > 0.05). Both clades change the same small set of amino acids, in the same order from small to large mammals, suggesting a limited number of ways in which contraction velocity can be successfully manipulated. To test this relationship, the 9 sites that differ between human and rat were mutated in the human β-myosin to match the rat sequence. Biochemical analysis revealed that the rat-human β-myosin chimera functioned like the native rat myosin with a 2-fold increase in both motility and in the rate of ADP release from the actin-myosin crossbridge (the step that limits contraction velocity). Thus, these sequence changes indicate adaptation of β-myosin as species mass increased to enable a reduced contraction velocity and heart rate.

Electrostatic interactions in the SH1-SH2 helix of human cardiac myosin modulate the time of strong actomyosin binding

Two single mutations, R694N and E45Q, were introduced in the beta isoform of human cardiac myosin to remove permanent salt bridges E45:R694 and E98:R694 in the SH1-SH2 helix of the myosin head. Beta isoform-specific bridges E45:R694 and E98:R694 were discovered in the molecular dynamics simulations of the alpha and beta myosin isoforms. Alpha and beta isoforms exhibit different kinetics, ADP dissociates slower from actomyosin containing beta myosin isoform, therefore, beta myosin stays strongly bound to actin longer. We hypothesize that the electrostatic interactions in the SH1-SH2 helix modulate the affinity of ADP to actomyosin, and therefore, the time of the strong actomyosin binding. Wild type and the mutants of the myosin head construct (1-843 amino acid residues) were expressed in differentiated C2C12 cells, and the duration of the strongly bound state of actomyosin was characterized using transient kinetics spectrophotometry. All myosin constructs exhibited a fast rate of ATP binding to actomyosin and a slow rate of ADP dissociation, showing that ADP release limits the time of the strongly bound state of actomyosin. The mutant R694N showed a faster rate of ADP release from actomyosin, compared to the wild type and the E45Q mutant, thus indicating that electrostatic interactions within the SH1-SH2 helix region of human cardiac myosin modulate ADP release and thus, the duration of the strongly bound state of actomyosin.

Multiscale modeling of twitch contractions in cardiac trabeculae

Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin-actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length-tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.

The effect of variable troponin C mutation thin filament incorporation on cardiac muscle twitch contractions

One of the complexities of understanding the pathology of familial forms of cardiac diseases is the level of mutation incorporation in sarcomeres. Computational models of the sarcomere that are spatially explicit offer an approach to study aspects of mutational incorporation into myofilaments that are more challenging to get at experimentally. We studied two well characterized mutations of cardiac TnC, L48Q and I61Q, that decrease or increase the release rate of Ca²⁺ from cTnC, k_-Ca, resulting in HCM and DCM respectively [1]. Expression of these mutations in transgenic mice was used to provide experimental data for incorporation of 30 and 50% (respectively) into sarcomeres. Here we demonstrate that fixed length twitch contractions of trabeculae from mice containing mutant differ from WT; L48Q trabeculae have slower relaxation while I61Q trabeculae have markedly reduced peak tension. Using our multiscale modelling approach [2] we were able to describe the tension transients of WT mouse myocardium. Tension transients for the mutant cTnCs were simulated with changes in k_-Ca, measured experimentally for each cTnC mutant in whole troponin complex, a change in the affinity of cTnC for cTnI, and a reduction in the number of detached crossbridges available for binding. A major advantage of the multiscale explicit 3-D model is that it predicts the effects of variable mutation incorporation, and the effects of variations in mutation distribution within thin filaments in sarcomeres. Such effects are currently impossible to explore experimentally. We explored random and clustered distributions of mutant cTnCs in thin filaments, as well as distributions of individual thin filaments with only WT or mutant cTnCs present. The effects of variable amounts of incorporation and non-random distribution of mutant cTnCs are more marked for I61Q than L48Q cTnC. We conclude that this approach can be effective for study on mutations in multiple proteins of the sarcomere. SUMMARY: A challenge in experimental studies of diseases is accounting for the effect of variable mutation incorporation into myofilaments. Here we use a spatially explicit computational approach, informed by experimental data from transgenic mice expressing one of two mutations in cardiac Troponin C that increase or decrease calcium sensitivity. We demonstrate that the model can accurately describe twitch contractions for the data and go on to explore the effect of variable mutant incorporation and localization on simulated cardiac muscle twitches.

Alpha and beta myosin isoforms and human atrial and ventricular contraction

Human atrial and ventricular contractions have distinct mechanical characteristics including speed of contraction, volume of blood delivered and the range of pressure generated. Notably, the ventricle expresses predominantly β-cardiac myosin while the atrium expresses mostly the α-isoform. In recent years exploration of the properties of pure α- & β-myosin isoforms have been possible in solution, in isolated myocytes and myofibrils. This allows us to consider the extent to which the atrial vs ventricular mechanical characteristics are defined by the myosin isoform expressed, and how the isoform properties are matched to their physiological roles. To do this we Outline the essential feature of atrial and ventricular contraction; Explore the molecular structural and functional characteristics of the two myosin isoforms; Describe the contractile behaviour of myocytes and myofibrils expressing a single myosin isoform; Finally we outline the outstanding problems in defining the differences between the atria and ventricles. This allowed us consider what features of contraction can and cannot be ascribed to the myosin isoforms present in the atria and ventricles.

Human-induced pluripotent stem cells for modelling metabolic perturbations and impaired bioenergetics underlying cardiomyopathies

Normal cardiac contractile and relaxation functions are critically dependent on a continuous energy supply. Accordingly, metabolic perturbations and impaired mitochondrial bioenergetics with subsequent disruption of ATP production underpin a wide variety of cardiac diseases, including diabetic cardiomyopathy, dilated cardiomyopathy, hypertrophic cardiomyopathy, anthracycline cardiomyopathy, peripartum cardiomyopathy, and mitochondrial cardiomyopathies. Crucially, there are no specific treatments for preventing the onset or progression of these cardiomyopathies to heart failure, one of the leading causes of death and disability worldwide. Therefore, new treatments are needed to target the metabolic disturbances and impaired mitochondrial bioenergetics underlying these cardiomyopathies in order to improve health outcomes in these patients. However, investigation of the underlying mechanisms and the identification of novel therapeutic targets have been hampered by the lack of appropriate animal disease models. Furthermore, interspecies variation precludes the use of animal models for studying certain disorders, whereas patient-derived primary cell lines have limited lifespan and availability. Fortunately, the discovery of human-induced pluripotent stem cells has provided a promising tool for modelling cardiomyopathies via human heart tissue in a dish. In this review article, we highlight the use of patient-derived iPSCs for studying the pathogenesis underlying cardiomyopathies associated with metabolic perturbations and impaired mitochondrial bioenergetics, as the ability of iPSCs for self-renewal and differentiation makes them an ideal platform for investigating disease pathogenesis in a controlled in vitro environment. Continuing progress will help elucidate novel mechanistic pathways, and discover novel therapies for preventing the onset and progression of heart failure, thereby advancing a new era of personalized therapeutics for improving health outcomes in patients with cardiomyopathy.

X-ray Crystallographic and Molecular Dynamic Analyses of Drosophila melanogaster Embryonic Muscle Myosin Define Domains Responsible for Isoform-Specific Properties

Drosophila melanogaster is a powerful system for characterizing alternative myosin isoforms and modeling muscle diseases, but high-resolution structures of fruit fly contractile proteins have not been determined. Here we report the first x-ray crystal structure of an insect myosin: the D melanogaster skeletal muscle myosin II embryonic isoform (EMB). Using our system for recombinant expression of myosin heavy chain (MHC) proteins in whole transgenic flies, we prepared and crystallized stable proteolytic S1-like fragments containing the entire EMB motor domain bound to an essential light chain. We solved the x-ray crystal structure by molecular replacement and refined the resulting model against diffraction data to 2.2 Å resolution. The protein is captured in two slightly different renditions of the rigor-like conformation with a citrate of crystallization at the nucleotide binding site and exhibits structural features common to myosins of diverse classes from all kingdoms of life. All atom molecular dynamics simulations on EMB in its nucleotide-free state and a derivative homology model containing 61 amino acid substitutions unique to the indirect flight muscle isoform (IFI) suggest that differences in the identity of residues within the relay and the converter that are encoded for by MHC alternative exons 9 and 11, respectively, directly contribute to increased mobility of these regions in IFI relative to EMB. This suggests the possibility that alternative folding or conformational stability within these regions contribute to the observed functional differences in Drosophila EMB and IFI myosins.

Molecular features of the UNC-45 chaperone critical for binding and folding muscle myosin

Myosin is a motor protein that is essential for a variety of processes ranging from intracellular transport to muscle contraction. Folding and assembly of myosin relies on a specific chaperone, UNC-45. To address its substrate-targeting mechanism, we reconstitute the interplay between Caenorhabditis elegans UNC-45 and muscle myosin MHC-B in insect cells. In addition to providing a cellular chaperone assay, the established system enabled us to produce large amounts of functional muscle myosin, as evidenced by a biochemical and structural characterization, and to directly monitor substrate binding to UNC-45. Data from in vitro and cellular chaperone assays, together with crystal structures of binding-deficient UNC-45 mutants, highlight the importance of utilizing a flexible myosin-binding domain. This so-called UCS domain can adopt discrete conformations to efficiently bind and fold substrate. Moreover, our data uncover the molecular basis of temperature-sensitive UNC-45 mutations underlying one of the most prominent motility defects in C. elegans.

Myosin, a motor protein essential for intracellular transport to muscle contraction, requires a chaperone UNC-45 for folding and assembly. Here authors use in vitro reconstitution and structural biology to characterize the interplay between UNC-45 and muscle myosin MHC-B.

CaATP prolongs strong actomyosin binding and promotes futile myosin stroke

Calcium plays an essential role in muscle contraction, regulating actomyosin interaction by binding troponin of thin filaments. There are several buffers for calcium in muscle, and those buffers play a crucial role in the formation of the transient calcium wave in sarcomere upon muscle activation. One such calcium buffer in muscle is ATP. ATP is a fuel molecule, and the important role of MgATP in muscle is to bind myosin and supply energy for the power stroke. Myosin is not a specific ATPase, and CaATP also supports myosin ATPase activity. The concentration of CaATP in sarcomeres reaches 1% of all ATP available. Since 294 myosin molecules form a thick filament, naïve estimation gives three heads per filament with CaATP bound, instead of MgATP. We found that CaATP dissociates actomyosin slower than MgATP, thus increasing the time of the strong actomyosin binding. The rate of the basal CaATPase is faster than that of MgATPase, myosin readily produces futile stroke with CaATP. When calcium is upregulated, as in malignant hyperthermia, kinetics of myosin and actomyosin interaction with CaATP suggest that myosin CaATPase activity may contribute to observed muscle rigidity and enhanced muscle thermogenesis.

A micro-volume adaptation of a stopped-flow system; use with μg quantities of muscle proteins

Stopped-flow spectroscopy is a powerful method for measuring very fast biological and chemical reactions. The technique however is often limited by the volumes of reactants needed to load the system. Here we present a simple adaptation of commercial stopped-flow system that reduces the volume needed by a factor of 4 to ≈120 μl. After evaluation the volume requirements of the system we show that many standard myosin based assays can be performed using <100 μg of myosin. This adaptation both reduces the volume and therefore mass of protein required and also produces data of similar quality to that produced using the standard set up. The 100 μg of myosin required for these assays is less than that which can be isolated from 100 mg of muscle tissue. With this reduced quantity of myosin, assays using biopsy samples become possible. This will allow assays to be used to assist diagnoses, to examine the effects of post translational modifications on muscle proteins and to test potential therapeutic drugs using patient derived samples.

The ATPase cycle of human muscle myosin II isoforms: Adaptation of a single mechanochemical cycle for different physiological roles

Striated muscle myosins are encoded by a large gene family in all mammals, including humans. These isoforms define several of the key characteristics of the different striated muscle fiber types, including maximum shortening velocity. We have previously used recombinant isoforms of the motor domains of seven different human myosin isoforms to define the actin·myosin cross-bridge cycle in solution. Here, we present data on an eighth isoform, the perinatal, which has not previously been characterized. The perinatal is distinct from the embryonic isoform, appearing to have features in common with the adult fast-muscle isoforms, including weak affinity of ADP for actin·myosin and fast ADP release. We go on to use a recently developed modeling approach, MUSICO, to explore how well the experimentally defined cross-bridge cycles for each isoform in solution can predict the characteristics of muscle fiber contraction, including duty ratio, shortening velocity, ATP economy, and load dependence of these parameters. The work shows that the parameters of the cross-bridge cycle predict many of the major characteristics of each muscle fiber type and raises the question of what sequence changes are responsible for these characteristics.

Electrostatic interactions in the force-generating region of the human cardiac myosin modulate ADP dissociation from actomyosin

Human cardiac myosin has two isoforms, alpha and beta, sharing significant sequence similarity, but different in kinetics: ADP release from actomyosin is an order of magnitude faster in the alpha myosin isoform. Apparently, small differences in the sequence are responsible for distinct local inter-residue interactions within alpha and beta isoforms, leading to such a dramatic difference in the rate of ADP release. Our analysis of structural kinetics of alpha and beta isoforms using molecular dynamics simulations revealed distinct dynamics of SH1:SH2 helix within the force-generation region of myosin head. The simulations showed that the residue R694 of the helix forms two permanent salt bridges in the beta isoform, which are not present in the alpha isoform. We hypothesized that the isoform-specific electrostatic interactions play a role in the difference of kinetic properties of myosin isoforms. We prepared R694N mutant in the beta isoform background to destabilize electrostatic interactions in the force-generating region of the myosin head. Our experimental data confirm faster ADP release from R694N actomyosin mutant, but is not as dramatic as the difference of kinetics of ADP release in the alpha and beta isoforms.

Modeling the Actin.myosin ATPase Cross-Bridge Cycle for Skeletal and Cardiac Muscle Myosin Isoforms

Modeling the complete actin.myosin ATPase cycle has always been limited by the lack of experimental data concerning key steps of the cycle, because these steps can only be defined at very low ionic strength. Here, using human β-cardiac myosin-S1, we combine published data from transient and steady-state kinetics to model a minimal eight-state ATPase cycle. The model illustrates the occupancy of each intermediate around the cycle and how the occupancy is altered by changes in actin concentration for [actin] = 1-20K_m. The cycle can be used to predict the maximal velocity of contraction (by motility assay or sarcomeric shortening) at different actin concentrations (which is consistent with experimental velocity data) and predict the effect of a 5 pN load on a single motor. The same exercise was repeated for human α-cardiac myosin S1 and rabbit fast skeletal muscle S1. The data illustrates how the motor domain properties can alter the ATPase cycle and hence the occupancy of the key states in the cycle. These in turn alter the predicted mechanical response of the myosin independent of other factors present in a sarcomere, such as filament stiffness and regulatory proteins. We also explore the potential of this modeling approach for the study of mutations in human β-cardiac myosin using the hypertrophic myopathy mutation R453C. Our modeling, using the transient kinetic data, predicts mechanical properties of the motor that are compatible with the single-molecule study. The modeling approach may therefore be of wide use for predicting the properties of myosin mutations.

Macromolecular Crowding Modulates Actomyosin Kinetics

Actomyosin kinetics is usually studied in dilute solutions, which do not reflect conditions in the cytoplasm. In cells, myosin and actin work in a dense macromolecular environment. High concentrations of macromolecules dramatically reduce the amount of free space available for all solutes, which results in an effective increase of the solutes' chemical potential and protein stabilization. Moreover, in a crowded solution, the chemical potential depends on the size of the solute, with larger molecules experiencing a larger excluded volume than smaller ones. Therefore, since myosin interacts with two ligands of different sizes (actin and ATP), macromolecular crowding can modulate the kinetics of individual steps of the actomyosin ATPase cycle. To emulate the effect of crowding in cells, we studied actomyosin cycle reactions in the presence of a high-molecular-weight polymer, Ficoll70. We observed an increase in the maximum velocity of the actomyosin ATPase cycle, and our transient-kinetics experiments showed that virtually all individual steps of the actomyosin cycle were affected by the addition of Ficoll70. The observed effects of macromolecular crowding on the myosin-ligand interaction cannot be explained by the increase of a solute's chemical potential. A time-resolved Förster resonance energy transfer experiment confirmed that the myosin head assumes a more compact conformation in the presence of Ficoll70 than in a dilute solution. We conclude that the crowding-induced myosin conformational change plays a major role in the changed kinetics of actomyosin ATPase.

Comparison of elementary steps of the cross-bridge cycle in rat papillary muscle fibers expressing α- and β-myosin heavy chain with sinusoidal analysis

In mammalian ventricles, two myosin heavy chain (MHC) isoforms have been identified. Small animals express α-MHC, whereas large animals express β-MHC, which contribute to a large difference in the heart rate. Sprague-Dawley rats possessing ~99% α-MHC were treated with propylthiouracil to result in 100% β-MHC. Papillary muscles were skinned, dissected into small fibers, and used for experiments. To understand the functional difference between α-MHC and β-MHC, skinned-fibers were activated under the intracellular ionic conditions: 5 mM MgATP, 1 mM Mg²⁺, 8 mM Pi, 200 mM ionic strength, pH 7.00 at 25 °C. Small amplitude sinusoidal length oscillations were applied in the frequency range 0.13-100 Hz (corresponding time domain: 1.6-1200 ms), and effects of Ca²⁺, Pi, and ATP were studied. The results show that Ca²⁺ sensitivity was slightly less (10-15%) in β-MHC than α-MHC containing fibers. Sinusoidal analysis at pCa 4.66 (full Ca²⁺ activation) demonstrated that, the apparent rate constants were 2-4× faster in α-MHC containing fibers. The ATP study demonstrated that, in β-MHC containing fibers, K ₁ (ATP association constant) was greater (1.7×), k ₂ and k _-2 (cross-bridge detachment and its reversal rate constants) were smaller (×0.6). The Pi study demonstrated that, in β-MHC containing fibers, k ₄ (rate constant of the force-generation step) and k _-4 were smaller (0.75× and 0.25×, respectively), resulting in greater K ₄ (3×). There were no differences in active tension, rigor stiffness, or K ₂ (equilibrium constant of the cross-bridge detachment step). Our study further demonstrated that there were no differences in parameters between fibers obtained from left and right ventricles, but with an exception in K ₅ (Pi association constant).

Myosin isoforms and the mechanochemical cross-bridge cycle

At the latest count the myosin family includes 35 distinct groups, all of which have the conserved myosin motor domain attached to a neck or lever arm, followed by a highly variable tail or cargo binding region. The motor domain has an ATPase activity that is activated by the presence of actin. One feature of the myosin ATPase cycle is that it involves an association/dissociation with actin for each ATP hydrolysed. The cycle has been described in detail for a large number of myosins from different classes. In each case the cycle is similar, but the balance between the different molecular events in the cycle has been altered to produce a range of very different mechanical activities. Myosin may spend most of the ATPase cycle attached to actin (high duty ratio), as in the processive myosin (e.g. myosin V) or the strain-sensing myosins (e.g. myosin 1c). In contrast, most muscle myosins spend 80% of their ATPase cycle detached from actin. Within the myosin IIs found in human muscle, there are 11 different sarcomeric myosin isoforms, two smooth muscle isoforms as well as three non-muscle isoforms. We have been exploring how the different myosin isoforms have adapted the cross-bridge cycle to generate different types of mechanical activity and how this goes wrong in inherited myopathies. The ideas are outlined here.

The Most Prevalent Freeman-Sheldon Syndrome Mutations in the Embryonic Myosin Motor Share Functional Defects

The embryonic myosin isoform is expressed during fetal development and rapidly down-regulated after birth. Freeman-Sheldon syndrome (FSS) is a disease associated with missense mutations in the motor domain of this myosin. It is the most severe form of distal arthrogryposis, leading to overcontraction of the hands, feet, and orofacial muscles and other joints of the body. Availability of human embryonic muscle tissue has been a limiting factor in investigating the properties of this isoform and its mutations. Using a recombinant expression system, we have studied homogeneous samples of human motors for the WT and three of the most common FSS mutants: R672H, R672C, and T178I. Our data suggest that the WT embryonic myosin motor is similar in contractile speed to the slow type I/β cardiac based on the rate constant for ADP release and ADP affinity for actin-myosin. All three FSS mutations show dramatic changes in kinetic properties, most notably the slowing of the apparent ATP hydrolysis step (reduced 5–9-fold), leading to a longer lived detached state and a slowed V_max of the ATPase (2–35-fold), indicating a slower cycling time. These mutations therefore seriously disrupt myosin function.

Establishing disease causality for a novel gene variant in familial dilated cardiomyopathy using a functional in-vitro assay of regulated thin filaments and human cardiac myosin

Background

As next generation sequencing for the genetic diagnosis of cardiovascular disorders becomes more widely used, establishing causality for putative disease causing variants becomes increasingly relevant. Diseases of the cardiac sarcomere provide a particular challenge in this regard because of the complexity of assaying the effect of genetic variants in human cardiac contractile proteins.

Results

In this study we identified a novel variant R205Q in the cardiac troponin T gene (TNNT2). Carriers of the variant allele exhibited increased chamber volumes associated with decreased left ventricular ejection fraction. To clarify the causal role of this variant, we generated recombinant variant human protein and examined its calcium kinetics as well as the maximally activated ADP release of human β-cardiac myosin with regulated thin filaments containing the mutant troponin T. We found that the R205Q mutation significantly decreased the calcium sensitivity of the thin filament by altering the effective calcium dissociation kinetics.

Conclusions

The development of moderate throughput post-genomic assays is an essential step in the realization of the potential of next generation sequencing. Although technically challenging, biochemical and functional assays of human cardiac contractile proteins of the thin filament can be achieved and provide an orthogonal source of information to inform the question of causality for individual variants.

Contractility parameters of human β-cardiac myosin with the hypertrophic cardiomyopathy mutation R403Q show loss of motor function

Hypertrophic cardiomyopathy (HCM) is the most frequently occurring inherited cardiovascular disease. It is caused by mutations in genes encoding the force-generating machinery of the cardiac sarcomere, including human β-cardiac myosin. We present a detailed characterization of the most debated HCM-causing mutation in human β-cardiac myosin, R403Q. Despite numerous studies, most performed with nonhuman or noncardiac myosin, there is no consensus about the mechanism of action of this mutation on the function of the enzyme. We use recombinant human β-cardiac myosin and new methodologies to characterize in vitro contractility parameters of the R403Q myosin compared to wild type. We extend our studies beyond pure actin filaments to include the interaction of myosin with regulated actin filaments containing tropomyosin and troponin. We find that, with pure actin, the intrinsic force generated by R403Q is ~15% lower than that generated by wild type. The unloaded velocity is, however, ~10% higher for R403Q myosin, resulting in a load-dependent velocity curve that has the characteristics of lower contractility at higher external loads compared to wild type. With regulated actin filaments, there is no increase in the unloaded velocity and the contractility of the R403Q myosin is lower than that of wild type at all loads. Unlike that with pure actin, the actin-activated adenosine triphosphatase activity for R403Q myosin with Ca(2+)-regulated actin filaments is ~30% lower than that for wild type, predicting a lower unloaded duty ratio of the motor. Overall, the contractility parameters studied fit with a loss of human β-cardiac myosin contractility as a result of the R403Q mutation.

Harmonic force spectroscopy measures load-dependent kinetics of individual human β-cardiac myosin molecules

Molecular motors are responsible for numerous cellular processes from cargo transport to heart contraction. Their interactions with other cellular components are often transient and exhibit kinetics that depend on load. Here, we measure such interactions using ‘harmonic force spectroscopy'. In this method, harmonic oscillation of the sample stage of a laser trap immediately, automatically and randomly applies sinusoidally varying loads to a single motor molecule interacting with a single track along which it moves. The experimental protocol and the data analysis are simple, fast and efficient. The protocol accumulates statistics fast enough to deliver single-molecule results from single-molecule experiments. We demonstrate the method's performance by measuring the force-dependent kinetics of individual human β-cardiac myosin molecules interacting with an actin filament at physiological ATP concentration. We show that a molecule's ADP release rate depends exponentially on the applied load, in qualitative agreement with cardiac muscle, which contracts with a velocity inversely proportional to external load.

Single molecule methods for measuring load dependence are fundamental for molecular motor research. Here, Sung et al. introduce harmonic force spectroscopy, a method that randomly applies varying loads at high frequency, allowing the determination of load dependent parameters of human β-cardiac myosin at physiological ATP concentration.

Ensemble force changes that result from human cardiac myosin mutations and a small-molecule effector

Cardiomyopathies due to mutations in human β-cardiac myosin are a significant cause of heart failure, sudden death, and arrhythmia. To understand the underlying molecular basis of changes in the contractile system's force production due to such mutations and search for potential drugs that restore force generation, an in vitro assay is necessary to evaluate cardiac myosin's ensemble force using purified proteins. Here, we characterize the ensemble force of human α- and β-cardiac myosin isoforms and those of β-cardiac myosins carrying left ventricular non-compaction (M531R) and dilated cardiomyopathy (S532P) mutations using a utrophin-based loaded in vitro motility assay and new filament-tracking software. Our results show that human α- and β-cardiac myosin, as well as the mutants, show opposite mechanical and enzymatic phenotypes with respect to each other. We also show that omecamtiv mecarbil, a previously discovered cardiac-specific myosin activator, increases β-cardiac myosin force generation.

Effects of troponin T cardiomyopathy mutations on the calcium sensitivity of the regulated thin filament and the actomyosin cross-bridge kinetics of human β-cardiac myosin

Hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM) lead to significant cardiovascular morbidity and mortality worldwide. Mutations in the genes encoding the sarcomere, the force-generating unit in the cardiomyocyte, cause familial forms of both HCM and DCM. This study examines two HCM-causing (I79N, E163K) and two DCM-causing (R141W, R173W) mutations in the troponin T subunit of the troponin complex using human β-cardiac myosin. Unlike earlier reports using various myosin constructs, we found that none of these mutations affect the maximal sliding velocities or maximal Ca²⁺-activated ADP release rates involving the thin filament human β-cardiac myosin complex. Changes in Ca²⁺ sensitivity using the human myosin isoform do, however, mimic changes seen previously with non-human myosin isoforms. Transient kinetic measurements show that these mutations alter the kinetics of Ca²⁺ induced conformational changes in the regulatory thin filament proteins. These changes in calcium sensitivity are independent of active, cycling human β-cardiac myosin.

The hypertrophic cardiomyopathy myosin mutation R453C alters ATP binding and hydrolysis of human cardiac β-myosin

The human hypertrophic cardiomyopathy mutation R453C results in one of the more severe forms of the myopathy. Arg-453 is found in a conserved surface loop of the upper 50-kDa domain of the myosin motor domain and lies between the nucleotide binding pocket and the actin binding site. It connects to the cardiomyopathy loop via a long α-helix, helix O, and to Switch-2 via the fifth strand of the central β-sheet. The mutation is, therefore, in a position to perturb a wide range of myosin molecular activities. We report here the first detailed biochemical kinetic analysis of the motor domain of the human β-cardiac myosin carrying the R453C mutation. A recent report of the same mutation (Sommese, R. F., Sung, J., Nag, S., Sutton, S., Deacon, J. C., Choe, E., Leinwand, L. A., Ruppel, K., and Spudich, J. A. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 12607-12612) found reduced ATPase and in vitro motility but increased force production using an optical trap. Surprisingly, our results show that the mutation alters few biochemical kinetic parameters significantly. The exceptions are the rate constants for ATP binding to the motor domain (reduced by 35%) and the ATP hydrolysis step/recovery stroke (slowed 3-fold), which could be the rate-limiting step for the ATPase cycle. Effects of the mutation on the recovery stroke are consistent with a perturbation of Switch-2 closure, which is required for the recovery stroke and the subsequent ATP hydrolysis.

Metal cation controls myosin and actomyosin kinetics

We have perturbed myosin nucleotide binding site with magnesium-, manganese-, or calcium-nucleotide complexes, using metal cation as a probe to examine the pathways of myosin ATPase in the presence of actin. We have used transient time-resolved FRET, myosin intrinsic fluorescence, fluorescence of pyrene labeled actin, combined with the steady state myosin ATPase activity measurements of previously characterized D.discoideum myosin construct A639C:K498C. We found that actin activation of myosin ATPase does not depend on metal cation, regardless of the cation-specific kinetics of nucleotide binding and dissociation. The rate limiting step of myosin ATPase depends on the metal cation. The rate of the recovery stroke and the reverse recovery stroke is directly proportional to the ionic radius of the cation. The rate of nucleotide release from myosin and actomyosin, and ATP binding to actomyosin depends on the cation coordination number.

Discrete effects of A57G-myosin essential light chain mutation associated with familial hypertrophic cardiomyopathy

The functional consequences of the familial hypertrophic cardiomyopathy A57G (alanine-to-glycine) mutation in the myosin ventricular essential light chain (ELC) were assessed in vitro and in vivo using previously generated transgenic (Tg) mice expressing A57G-ELC mutant vs. wild-type (WT) of human cardiac ELC and in recombinant A57G- or WT-protein-exchanged porcine cardiac muscle strips. Compared with the Tg-WT, there was a significant increase in the Ca²⁺ sensitivity of force (ΔpCa₅₀ ≅ 0.1) and an ~1.3-fold decrease in maximal force per cross section of muscle observed in the mutant preparations. In addition, a significant increase in passive tension in response to stretch was monitored in Tg-A57G vs. Tg-WT strips indicating a mutation-induced myocardial stiffness. Consistently, the hearts of Tg-A57G mice demonstrated a high level of fibrosis and hypertrophy manifested by increased heart weight-to-body weight ratios and a decreased number of nuclei indicating an increase in the two-dimensional size of Tg-A57G vs. Tg-WT myocytes. Echocardiography examination showed a phenotype of eccentric hypertrophy in Tg-A57G mice, enhanced left ventricular (LV) cavity dimension without changes in LV posterior/anterior wall thickness. Invasive hemodynamics data revealed significantly increased end-systolic elastance, defined by the slope of the pressure-volume relationship, indicating a mutation-induced increase in cardiac contractility. Our results suggest that the A57G allele causes disease by means of a discrete modulation of myofilament function, increased Ca²⁺ sensitivity, and decreased maximal tension followed by compensatory hypertrophy and enhanced contractility. These and other contributing factors such as increased myocardial stiffness and fibrosis most likely activate cardiomyopathic signaling pathways leading to pathologic cardiac remodeling.