Size and speed of the working stroke of cardiac myosin in situ

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.

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

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

SILICOFCM platform, multiscale modeling of left ventricle from echocardiographic images and drug influence for cardiomyopathy disease

BACKGROUND AND OBJECTIVE

In silico clinical trials are the future of medicine and virtual testing and simulation are the future of medical engineering. The use of a computational platform can reduce costs and time required for developing new models of medical devices and drugs. The computational platform, which is one of the main results of the SILICOFCM project, was developed using state-of-the-art finite element modeling for macro simulation of fluid-structure interaction with micro modeling at the molecular level for drug interaction with the cardiac cells. SILICOFCM platform is using for risk prediction and optimal drug therapy of familial cardiomyopathy in a specific patient.

METHODS

In order to obtain 3D image reconstruction, the U-net architecture was used to determine geometric parameters for the left ventricle which were extracted from the echocardiographic apical and M-mode views. A micro-mechanics cellular model which includes three kinetic processes of sarcomeric proteins interactions was developed. It allows simulation of the drugs which are divided into three major groups defined by the principal action of each drug. Fluid-solid coupling for the left ventricle was presented. A nonlinear material model of the heart wall that was developed by using constitutive curves which include the stress-strain relationship was used.

RESULTS

The results obtained with the parametric model of the left ventricle where pressure-volume (PV) diagrams depend on the change of Ca2+ were presented. It directly affects the ejection fraction. The presented approach with the variation of the left ventricle (LV) geometry and simulations which include the influence of different parameters on the PV diagrams are directly interlinked with drug effects on the heart function. It includes different drugs such as Entresto and Digoxin that directly affect the cardiac PV diagrams and ejection fraction.

CONCLUSIONS

Computational platforms such as the SILICOFCM platform are novel tools for risk prediction of cardiac disease in a specific patient that will certainly open a new avenue for in silico clinical trials in the future.

Cardiac ventricular myosin and slow skeletal myosin exhibit dissimilar chemomechanical properties despite bearing the same myosin heavy chain isoform

The myosin II motors are ATP-powered force-generating machines driving cardiac and muscle contraction. Myosin II heavy chain isoform-beta (β-MyHC) is primarily expressed in the ventricular myocardium and in slow-twitch muscle fibers, such as M. soleus. M. soleus-derived myosin II (SolM-II) is often used as an alternative to the ventricular β-cardiac myosin (βM-II); however, the direct assessment of biochemical and mechanical features of the native myosins is limited. By employing optical trapping, we examined the mechanochemical properties of native myosins isolated from the rabbit heart ventricle and soleus muscles at the single-molecule level. We found purified motors from the two tissue sources, despite expressing the same MyHC isoform, displayed distinct motile and ATPase kinetic properties. We demonstrate βM-II was approximately threefold faster in the actin filament-gliding assay than SolM-II. The maximum actomyosin (AM) detachment rate derived in single-molecule assays was also approximately threefold higher in βM-II, while the power stroke size and stiffness of the "AM rigor" crossbridge for both myosins were comparable. Our analysis revealed a higher AM detachment rate for βM-II, corresponding to the enhanced ADP release rates from the crossbridge, likely responsible for the observed differences in the motility driven by these myosins. Finally, we observed a distinct myosin light chain 1 isoform (MLC1sa) that associates with SolM-II, which might contribute to the observed kinetics differences between βM-II and SolM-II. These results have important implications for the choice of tissue sources and justify prerequisites for the correct myosin heavy and light chains to study cardiomyopathies.

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.

The regulatory light chain mediates inactivation of myosin motors during active shortening of cardiac muscle

The normal function of heart muscle depends on its ability to contract more strongly at longer length. Increased venous filling stretches relaxed heart muscle cells, triggering a stronger contraction in the next beat- the Frank-Starling relation. Conversely, heart muscle cells are inactivated when they shorten during ejection, accelerating relaxation to facilitate refilling before the next beat. Although both effects are essential for the efficient function of the heart, the underlying mechanisms were unknown. Using bifunctional fluorescent probes on the regulatory light chain of the myosin motor we show that its N-terminal domain may be captured in the folded OFF state of the myosin dimer at the end of the working-stroke of the actin-attached motor, whilst its C-terminal domain joins the OFF state only after motor detachment from actin. We propose that sequential folding of myosin motors onto the filament backbone may be responsible for shortening-induced de-activation in the heart.

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

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

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.

Allosteric modulation of cardiac myosin mechanics and kinetics by the conjugated omega-7,9 trans-fat rumenic acid

KEY POINTS

Direct binding of rumenic acid to the cardiac myosin-2 motor domain increases the release rate for orthophosphate and increases the Ca²⁺ responsiveness of cardiac muscle at low load. Physiological cellular concentrations of rumenic acid affect the ATP turnover rates of the super-relaxed and disordered relaxed states of β-cardiac myosin, leading to a net increase in myocardial metabolic load. In Ca²⁺ -activated trabeculae, rumenic acid exerts a direct inhibitory effect on the force-generating mechanism without affecting the number of force-generating motors. In the presence of saturating actin concentrations rumenic acid binds to the β-cardiac myosin-2 motor domain with an EC₅₀ of 200 nM. Molecular docking studies provide information about the binding site, the mode of binding, and associated allosteric communication pathways. Free rumenic acid may exceed thresholds in cardiomyocytes above which contractile efficiency is reduced and interference with small molecule therapeutics, targeting cardiac myosin, occurs.

ABSTRACT

Based on experiments using purified myosin motor domains, reconstituted actomyosin complexes and rat heart ventricular trabeculae, we demonstrate direct binding of rumenic acid, the cis-delta-9-trans-delta-11 isomer of conjugated linoleic acid, to an allosteric site located in motor domain of mammalian cardiac myosin-2 isoforms. In the case of porcine β-cardiac myosin, the EC₅₀ for rumenic acid varies from 10.5 μM in the absence of actin to 200 nM in the presence of saturating concentrations of actin. Saturating concentrations of rumenic acid increase the maximum turnover of basal and actin-activated ATPase activity of β-cardiac myosin approximately 2-fold but decrease the force output per motor by 23% during isometric contraction. The increase in ATP turnover is linked to an acceleration of the release of the hydrolysis product orthophosphate. In the presence of 5 μM rumenic acid, the difference in the rate of ATP turnover by the super-relaxed and disordered relaxed states of cardiac myosin increases from 4-fold to 20-fold. The equilibrium between the two functional myosin states is not affected by rumenic acid. Calcium responsiveness is increased under zero-load conditions but unchanged under load. Molecular docking studies provide information about the rumenic acid binding site, the mode of binding, and associated allosteric communication pathways. They show how the isoform-specific replacement of residues in the binding cleft induces a different mode of rumenic acid binding in the case of non-muscle myosin-2C and blocks binding to skeletal muscle and smooth muscle myosin-2 isoforms.

The Sliding Filament Theory Since Andrew Huxley: Multiscale and Multidisciplinary Muscle Research

Two groundbreaking papers published in 1954 laid out the theory of the mechanism of muscle contraction based on force-generating interactions between myofilaments in the sarcomere that cause filaments to slide past one another during muscle contraction. The succeeding decades of research in muscle physiology have revealed a unifying interest: to understand the multiscale processes-from atom to organ-that govern muscle function. Such an understanding would have profound consequences for a vast array of applications, from developing new biomimetic technologies to treating heart disease. However, connecting structural and functional properties that are relevant at one spatiotemporal scale to those that are relevant at other scales remains a great challenge. Through a lens of multiscale dynamics, we review in this article current and historical research in muscle physiology sparked by the sliding filament theory.

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.

Biophysically detailed mathematical models of multiscale cardiac active mechanics

We propose four novel mathematical models, describing the microscopic mechanisms of force generation in the cardiac muscle tissue, which are suitable for multiscale numerical simulations of cardiac electromechanics. Such models are based on a biophysically accurate representation of the regulatory and contractile proteins in the sarcomeres. Our models, unlike most of the sarcomere dynamics models that are available in the literature and that feature a comparable richness of detail, do not require the time-consuming Monte Carlo method for their numerical approximation. Conversely, the models that we propose only require the solution of a system of PDEs and/or ODEs (the most reduced of the four only involving 20 ODEs), thus entailing a significant computational efficiency. By focusing on the two models that feature the best trade-off between detail of description and identifiability of parameters, we propose a pipeline to calibrate such parameters starting from experimental measurements available in literature. Thanks to this pipeline, we calibrate these models for room-temperature rat and for body-temperature human cells. We show, by means of numerical simulations, that the proposed models correctly predict the main features of force generation, including the steady-state force-calcium and force-length relationships, the length-dependent prolongation of twitches and increase of peak force, the force-velocity relationship. Moreover, they correctly reproduce the Frank-Starling effect, when employed in multiscale 3D numerical simulation of cardiac electromechanics.

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

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

Hierarchical modeling of force generation in cardiac muscle

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

Orthophosphate increases the efficiency of slow muscle-myosin isoform in the presence of omecamtiv mecarbil

Omecamtiv mecarbil (OM) is a putative positive inotropic tool for treatment of systolic heart dysfunction, based on the finding that in vivo it increases the ejection fraction and in vitro it prolongs the actin-bond life time of the cardiac and slow-skeletal muscle isoforms of myosin. OM action in situ, however, is still poorly understood as the enhanced Ca²⁺-sensitivity of the myofilaments is at odds with the reduction of force and rate of force development observed at saturating Ca²⁺. Here we show, by combining fast sarcomere-level mechanics and ATPase measurements in single slow demembranated fibres from rabbit soleus, that the depressant effect of OM on the force per attached motor is reversed, without effect on the ATPase rate, by physiological concentrations of inorganic phosphate (Pi) (1-10 mM). This mechanism could underpin an energetically efficient reduction of systolic tension cost in OM-treated patients, whenever [Pi] increases with heart-beat frequency.

Omecamtiv mecarbil is a small molecule effector under clinical trial for the treatment of systolic heart failure. Here the authors define the molecular mechanisms of its inotropic action and find it can increase the efficiency of contraction in muscle fibres when the orthophosphate concentration rises with the beat frequency.

The Mechanical Power of Titin Folding

The delivery of mechanical power, a crucial component of animal motion, is constrained by the universal compromise between the force and the velocity of its constituent molecular systems. While the mechanisms of force generation have been studied at the single molecular motor level, there is little understanding of the magnitude of power that can be generated by folding proteins. Here, we use single-molecule force spectroscopy techniques to measure the force-velocity relation of folding titin domains that contain single internal disulfide bonds, a common feature throughout the titin I-band. We find that formation of the disulfide regulates the peak power output of protein folding in an all-or-none manner, providing at 6.0 pN, for example, a boost from 0 to 6,000 zW upon oxidation. This mechanism of power generation from protein folding is of great importance for muscle, where titin domains may unfold and refold with each extension and contraction of the sarcomere.

Eckels et al. use single-molecule magnetic tweezers to simultaneously probe the folding dynamics of titin Ig domains and monitor the redox status of single disulfides within the Ig fold. Oxidation of the disulfide bond greatly increases both the folding force and the magnitude of power delivered by protein folding.

Monitoring of cardiovascular physiology augmented by a patient-specific biomechanical model during general anesthesia. A proof of concept study

During general anesthesia (GA), direct analysis of arterial pressure or aortic flow waveforms may be inconclusive in complex situations. Patient-specific biomechanical models, based on data obtained during GA and capable to perform fast simulations of cardiac cycles, have the potential to augment hemodynamic monitoring. Such models allow to simulate Pressure-Volume (PV) loops and estimate functional indicators of cardiovascular (CV) system, e.g. ventricular-arterial coupling (Vva), cardiac efficiency (CE) or myocardial contractility, evolving throughout GA. In this prospective observational study, we created patient-specific biomechanical models of heart and vasculature of a reduced geometric complexity for n = 45 patients undergoing GA, while using transthoracic echocardiography and aortic pressure and flow signals acquired in the beginning of GA (baseline condition). If intraoperative hypotension (IOH) appeared, diluted norepinephrine (NOR) was administered and the model readjusted according to the measured aortic pressure and flow signals. Such patients were a posteriori assigned into a so-called hypotensive group. The accuracy of simulated mean aortic pressure (MAP) and stroke volume (SV) at baseline were in accordance with the guidelines for the validation of new devices or reference measurement methods in all patients. After NOR administration in the hypotensive group, the percentage of concordance with 10% exclusion zone between measurement and simulation was >95% for both MAP and SV. The modeling results showed a decreased Vva (0.64±0.37 vs 0.88±0.43; p = 0.039) and an increased CE (0.8±0.1 vs 0.73±0.11; p = 0.042) in hypotensive vs normotensive patients. Furthermore, Vva increased by 92±101%, CE decreased by 13±11% (p < 0.001 for both) and contractility increased by 14±11% (p = 0.002) in the hypotensive group post-NOR administration. In this work we demonstrated the application of fast-running patient-specific biophysical models to estimate PV loops and functional indicators of CV system using clinical data available during GA. The work paves the way for model-augmented hemodynamic monitoring at operating theatres or intensive care units to enhance the information on patient-specific physiology.

Dobutamine stress testing in patients with Fontan circulation augmented by biomechanical modeling

Understanding (patho)physiological phenomena and mechanisms of failure in patients with Fontan circulation-a surgically established circulation for patients born with a functionally single ventricle-remains challenging due to the complex hemodynamics and high inter-patient variations in anatomy and function. In this work, we present a biomechanical model of the heart and circulation to augment the diagnostic evaluation of Fontan patients with early-stage heart failure. The proposed framework employs a reduced-order model of heart coupled with a simplified circulation including venous return, creating a closed-loop system. We deploy this framework to augment the information from data obtained during combined cardiac catheterization and magnetic resonance exams (XMR), performed at rest and during dobutamine stress in 9 children with Fontan circulation and 2 biventricular controls. We demonstrate that our modeling framework enables patient-specific investigation of myocardial stiffness, contractility at rest, contractile reserve during stress and changes in vascular resistance. Hereby, the model allows to identify key factors underlying the pathophysiological response to stress in these patients. In addition, the rapid personalization of the model to patient data and fast simulation of cardiac cycles make our framework directly applicable in a clinical workflow. We conclude that the proposed modeling framework is a valuable addition to the current clinical diagnostic XMR exam that helps to explain patient-specific stress hemodynamics and can identify potential mechanisms of failure in patients with Fontan circulation.

The Work of Titin Protein Folding as a Major Driver in Muscle Contraction

Single-molecule atomic force microscopy and magnetic tweezers experiments have demonstrated that titin immunoglobulin (Ig) domains are capable of folding against a pulling force, generating mechanical work that exceeds that produced by a myosin motor. We hypothesize that upon muscle activation, formation of actomyosin cross bridges reduces the force on titin, causing entropic recoil of the titin polymer and triggering the folding of the titin Ig domains. In the physiological force range of 4-15 pN under which titin operates in muscle, the folding contraction of a single Ig domain can generate 200% of the work of entropic recoil and occurs at forces that exceed the maximum stalling force of single myosin motors. Thus, titin operates like a mechanical battery, storing elastic energy efficiently by unfolding Ig domains and delivering the charge back by folding when the motors are activated during a contraction. We advance the hypothesis that titin folding and myosin activation act as inextricable partners during muscle contraction.

Cardiac Disorders and Pathophysiology of Sarcomeric Proteins

The sarcomeric proteins represent the structural building blocks of heart muscle, which are essential for contraction and relaxation. During recent years, it has become evident that posttranslational modifications of sarcomeric proteins, in particular phosphorylation, tune cardiac pump function at rest and during exercise. This delicate, orchestrated interaction is also influenced by mutations, predominantly in sarcomeric proteins, which cause hypertrophic or dilated cardiomyopathy. In this review, we follow a bottom-up approach starting from a description of the basic components of cardiac muscle at the molecular level up to the various forms of cardiac disorders at the organ level. An overview is given of sarcomere changes in acquired and inherited forms of cardiac disease and the underlying disease mechanisms with particular reference to human tissue. A distinction will be made between the primary defect and maladaptive/adaptive secondary changes. Techniques used to unravel functional consequences of disease-induced protein changes are described, and an overview of current and future treatments targeted at sarcomeric proteins is given. The current evidence presented suggests that sarcomeres not only form the basis of cardiac muscle function but also represent a therapeutic target to combat cardiac disease.

Proposed mechanism for the length dependence of the force developed in maximally activated muscles

The molecular bases of the Frank-Starling law of the heart and of its cellular counterpart, the length dependent activation (LDA), are largely unknown. However, the recent discovery of the thick filament activation, a second pathway beside the well-known calcium mediated thin filament activation, is promising for elucidating these mechanisms. The thick filament activation is mediated by the tension acting on it through the mechano-sensing (MS) mechanism and can be related to the LDA via the titin passive tension. Here, we propose a mechanism to explain the higher maximum tension at longer sarcomere lengths generated by a maximally activated muscle and test it in-silico with a single fiber and a ventricle model. The active tension distribution along the thick filament generates a reservoir of inactive motors at its free-end that can be activated by passive tension on a beat-to-beat timescale. The proposed mechanism is able to quantitatively account for the observed increment in tension at the fiber level, however, the ventricle model suggests that this component of the LDA is not crucial in physiological conditions.

How myofilament strain and strain rate lead the dance of the cardiac cycle

Movement of the myocardium can modify organ-level cardiac function and its molecular (crossbridge) mechanisms. This motion, which is defined by myocardial strain and strain rate (muscle shortening, lengthening, and the speed of these movements), occurs throughout the cardiac cycle, including during isovolumic periods. This review highlights how the left ventricular myocardium moves throughout the cardiac cycle, how muscle mechanics experiments provide insight into the regulation of forces used to move blood in and out of the left ventricle, and its impact on (and regulation by) crossbridge and sarcomere kinetics. We specifically highlight how muscle mechanics experiments explain how myocardial relaxation is accelerated by lengthening (strain rate) during late systole and isovolumic relaxation, a lengthening which has been measured in human hearts. Advancing and refining both in vivo measurement and ex vivo protocols with physiologic strain and strain rates could reveal important insights into molecular (crossbridge) kinetics. These advances could provide an improvement in both diagnosis and precise treatment of cardiac dysfunction.

Inotropic interventions do not change the resting state of myosin motors during cardiac diastole

When striated (skeletal and cardiac) muscle is in its relaxed state, myosin motors are packed in helical tracks on the surface of the thick filament, folded toward the center of the sarcomere, and unable to bind actin or hydrolyze ATP (OFF state). This raises the question of whatthe mechanism is that integrates the Ca²⁺-dependent thin filament activation, making myosin heads available for interaction with actin. Here we test the interdependency of the thin and thick filament regulatory mechanisms in intact trabeculae from the rat heart. We record the x-ray diffraction signals that mark the state of the thick filament during inotropic interventions (increase in sarcomere length from 1.95 to 2.25 µm and addition of 10^-7 M isoprenaline), which potentiate the twitch force developed by an electrically paced trabecula by up to twofold. During diastole, none of the signals related to the OFF state of the thick filament are significantly affected by these interventions, except the intensity of both myosin-binding protein C- and troponin-related meridional reflections, which reduce by 20% in the presence of isoprenaline. These results indicate that recruitment of myosin motors from their OFF state occurs independently and downstream from thin filament activation. This is in agreement with the recently discovered mechanism based on thick filament mechanosensing in which the number of motors available for interaction with actin rapidly adapts to the stress on the thick filament and thus to the loading conditions of the contraction. The gain of this positive feedback may be modulated by both sarcomere length and the degree of phosphorylation of myosin-binding protein C.

The force and stiffness of myosin motors in the isometric twitch of a cardiac trabecula and the effect of the extracellular calcium concentration

KEY POINTS

Fast sarcomere-level mechanics in intact trabeculae, which allows the definition of the mechano-kinetic properties of cardiac myosin in situ, is a fundamental tool not only for understanding the molecular mechanisms of heart performance and regulation, but also for investigating the mechanisms of the cardiomyopathy-causing mutations in the myosin and testing small molecules for therapeutic interventions. The approach has been applied to measure the stiffness and force of the myosin motor and the fraction of motors attached during isometric twitches of electrically paced trabeculae under different extracellular Ca2+ concentrations. Although the average force of the cardiac myosin motor (∼6 pN) is similar to that of the fast myosin isoform of skeletal muscle, the stiffness (1.07 pN nm-1 ) is 2- to 3-fold smaller. The increase in the twitch force developed in the presence of larger extracellular Ca2+ concentrations is fully accounted for by a proportional increase in the number of attached motors.

ABSTRACT

The mechano-kinetic properties of the cardiac myosin were studied in situ, in trabeculae dissected from the right ventricle of the rat heart, by measuring the stiffness of the half-sarcomere both at the twitch force peak (Tp ) of an electrically paced intact trabecula at different extracellular Ca2+ concentrations ([Ca2+ ]o ), and in the same trabecula after skinning and induction of rigor. Taking into account the contribution of filament compliance to half-sarcomere compliance and the lattice geometry, we found that the stiffness of the cardiac myosin motor is 1.07 ± 0.09 pN nm-1 , which is slightly larger than that of the slow myosin isoform of skeletal muscle (0.6-0.8 pN nm-1 ) and 2- to 3-fold smaller than that of the fast skeletal muscle isoform. The increase in Tp from 61 ± 4 kPa to 93 ± 9 kPa, induced by raising [Ca2+ ]o from 1 to 2.5 mm at sarcomere length ∼2.2 μm, is accompanied by an increase of the half-sarcomere stiffness that is explained by an increase of the fraction of actin-attached motors from 0.08 ± 0.01 to 0.12 ± 0.02, proportional to Tp . Consequently, each myosin motor bears an average force of 6.14 ± 0.52 pN independently of Tp and [Ca2+ ]o . The application of fast sarcomere-level mechanics to intact trabeculae to define the mechano-kinetic properties of the cardiac myosin in situ represents a powerful tool for investigating cardiomyopathy-causing mutations in the myosin motor and testing specific therapeutic interventions.

Thick Filament Mechano-Sensing in Skeletal and Cardiac Muscles: A Common Mechanism Able to Adapt the Energetic Cost of the Contraction to the Task

A dual regulation of contraction operates in both skeletal and cardiac muscles. The first mechanism, based on Ca²⁺-dependent structural changes of the regulatory proteins in the thin filament, makes the actin sites available for binding of the myosin motors. The second recruits the myosin heads from the OFF state, in which they are unable to split ATP and bind to actin, in relation to the force during contraction. Comparison of the relevant X-ray diffraction signals marking the state of the thick filament demonstrates that the force feedback that controls the regulatory state of the thick filament works in the same way in skeletal as in cardiac muscles: even if in an isometric tetanus of skeletal muscle force is under the control of the firing frequency of the motor unit, while in a heartbeat force is controlled by the afterload, the stress-sensor switching the motors ON plays the same role in adapting the energetic cost of the contraction to the force. A new aspect of the Frank-Starling law of the heart emerges: independent of the diastolic filling of the ventricle, the number of myosin motors switched ON during systole, and thus the energetic cost of contraction, are tuned to the arterial pressure. Deterioration of the thick-filament regulation mechanism may explain the hyper-contractility related to hypertrophic cardiomyopathy, an inherited heart disease that in 40% of cases is due to mutations in cardiac myosin.

Successive contractile periods activate mitochondria at the onset of contractions in intact rat cardiac trabeculae

The rate of oxidative phosphorylation depends on the contractile activity of the heart. Cardiac mitochondrial oxidative phosphorylation is determined by free ADP concentration, mitochondrial Ca²⁺ accumulation, mitochondrial enzyme activities, and Krebs cycle intermediates. The purpose of the present study was to examine the factors that limit oxidative phosphorylation upon rapid changes in contractile activity in cardiac muscle. We tested the hypotheses that prior contractile performance enhances the changes in NAD(P)H and FAD concentration upon an increase in contractile activity and that this mitochondrial "priming" depends on pyruvate dehydrogenase activity. Intact rat cardiac trabeculae were electrically stimulated at 0.5 Hz for at least 30 min. Thereafter, two equal bouts at elevated stimulation frequency of 1, 2, or 3 Hz were applied for 3 min with 3 min of 0.5-Hz stimulation in between. No discernible time delay was observed in the changes in NAD(P)H and FAD fluorescence upon rapid changes in contractile activity. The amplitudes of the rapid changes in fluorescence upon an increase in stimulation frequency (the on-transients) were smaller than upon a decrease in stimulation frequency (the off-transients). A first bout in glucose-containing superfusion solution resulted, during the second bout, in an increase in the amplitudes of the on-transients, but the off-transients remained the same. No such priming effect was observed after addition of 10 mM pyruvate. These results indicate that mitochondrial priming can be observed in cardiac muscle in situ and that pyruvate dehydrogenase activity is critically involved in the mitochondrial adaptation to increases in contractile performance. NEW & NOTEWORTHY Mitochondrial respiration increases with increased cardiac contractile activity. Similar to mitochondrial "priming" in skeletal muscle, we hypothesized that cardiac mitochondrial activity is altered upon successive bouts of contractions and depends on pyruvate dehydrogenase activity. We found altered bioenergetics upon repeated contractile periods, indicative of mitochondrial priming in rat myocardium. No effect was seen when pyruvate was added to the perfusate. As such, pyruvate dehydrogenase activity is involved in the mitochondrial adaptation to increased contractile performance.

Experimental and modelling evidence of shortening heat in cardiac muscle

KEY POINTS

Heat associated with muscle shortening has been repeatedly demonstrated in skeletal muscle, but its existence in cardiac muscle remains contentious after five decades of study. By iterating between experiments and computational modelling, we show compelling evidence for the existence of shortening heat in cardiac muscle and reveal, mechanistically, the source of this excess heat. Our results clarify a long-standing uncertainty in the field of cardiac muscle energetics. We provide a revised partitioning of cardiac muscle energy expenditure to include this newly revealed thermal component.

ABSTRACT

When a muscle shortens against an afterload, the heat that it liberates is greater than that produced by the same muscle contracting isometrically at the same level of force. This excess heat is defined as 'shortening heat', and has been repeatedly demonstrated in skeletal muscle but not in cardiac muscle. Given the micro-structural similarities between these two muscle types, and since we imagine that shortening heat is the thermal accompaniment of cross-bridge cycling, we have re-examined this issue. Using our flow-through microcalorimeter, we measured force and heat generated by isolated rat trabeculae undergoing isometric contractions at different muscle lengths and work-loop (shortening) contractions at different afterloads. We simulated these experimental protocols using a thermodynamically constrained model of cross-bridge cycling and probed the mechanisms underpinning shortening heat. Predictions generated by the model were subsequently validated by a further set of experiments. Both our experimental and modelling results show convincing evidence for the existence of shortening heat in cardiac muscle. Its magnitude is inversely related to the afterload or, equivalently, directly related to the extent of shortening. Computational simulations reveal that the heat of shortening arises from the cycling of cross-bridges, and that the rate of ATP hydrolysis is more sensitive to change of muscle length than to change of afterload. Our results clarify a long-standing uncertainty in the field of cardiac muscle energetics.

Titin-mediated thick filament activation, through a mechanosensing mechanism, introduces sarcomere-length dependencies in mathematical models of rat trabecula and whole ventricle

Recent experimental evidence in skeletal muscle demonstrated the existence of a thick-filament mechanosensing mechanism, acting as a second regulatory system for muscle contraction, in addition to calcium-mediated thin filament regulation. These two systems cooperate to generate force, but the extent to which their interaction is relevant in physiologically contracting muscle was not yet assessed experimentally. Therefore, we included both regulatory mechanisms in a mathematical model of rat trabecula and whole ventricle. No additional regulatory mechanisms were considered in our model. Our simulations suggested that mechanosensing regulation is not limited to the initial phases of contraction but, instead, is crucial during physiological contraction. An important consequence of this finding is that titin mediated thick filament activation can account for several sarcomere length dependencies observed in contracting muscle. Under the hypothesis that a similar mechanism is acting on cardiac muscle, and within the limits of a finite element left ventricle model, we predict that these two regulatory mechanisms are crucial for the molecular basis of the Frank-Starling law of the heart.

Myosin filament activation in the heart is tuned to the mechanical task

The mammalian heart pumps blood through the vessels, maintaining the dynamic equilibrium in a circulatory system driven by two pumps in series. This vital function is based on the fine-tuning of cardiac performance by the Frank-Starling mechanism that relates the pressure exerted by the contracting ventricle (end systolic pressure) to its volume (end systolic volume). At the level of the sarcomere, the structural unit of the cardiac myocytes, the Frank-Starling mechanism consists of the increase in active force with the increase of sarcomere length (length-dependent activation). We combine sarcomere mechanics and micrometer-nanometer-scale X-ray diffraction from synchrotron light in intact ventricular trabeculae from the rat to measure the axial movement of the myosin motors during the diastole-systole cycle under sarcomere length control. We find that the number of myosin motors leaving the off, ATP hydrolysis-unavailable state characteristic of the diastole is adjusted to the sarcomere length-dependent systolic force. This mechanosensing-based regulation of the thick filament makes the energetic cost of the systole rapidly tuned to the mechanical task, revealing a prime aspect of the Frank-Starling mechanism. The regulation is putatively impaired by cardiomyopathy-causing mutations that affect the intramolecular and intermolecular interactions controlling the off state of the motors.

Distinct contributions of the thin and thick filaments to length-dependent activation in heart muscle

The Frank-Starling relation is a fundamental auto-regulatory property of the heart that ensures the volume of blood ejected in each heartbeat is matched to the extent of venous filling. At the cellular level, heart muscle cells generate higher force when stretched, but despite intense efforts the underlying molecular mechanism remains unknown. We applied a fluorescence-based method, which reports structural changes separately in the thick and thin filaments of rat cardiac muscle, to elucidate that mechanism. The distinct structural changes of troponin C in the thin filaments and myosin regulatory light chain in the thick filaments allowed us to identify two aspects of the Frank-Starling relation. Our results show that the enhanced force observed when heart muscle cells are maximally activated by calcium is due to a change in thick filament structure, but the increase in calcium sensitivity at lower calcium levels is due to a change in thin filament structure.

DOI:http://dx.doi.org/10.7554/eLife.24081.001

The heart needs to pump out the same volume of blood that enters it. This is not as simple as it sounds, as changes in heart rate – for example, in response to exercise – alter how hard the heart must pump.

When blood flows into the heart it stretches the heart muscle, which consists of units called sarcomeres. Sarcomeres contain two types of protein filament, known as thick filaments and thin filaments. When a heartbeat is triggered by calcium ions flowing into the heart muscle cells, the thick filaments slide over the thin filaments. This causes the heart muscle cell to contract.

The Frank–Starling mechanism helps to regulate the contraction of the heart. This mechanism has two aspects. Firstly, as the sarcomere lengthens, its protein filaments are able to contract with more force for a given high level of calcium ions. Secondly, the lengthening of the sarcomere makes the filaments more sensitive to calcium ions, which again causes the heart to contract more forcefully. However, the molecular mechanisms that underlie these effects were not clear.

Zhang et al. have now studied rat heart muscle cells using a new fluorescence-based method that can detect structural changes in the thick and thin filaments. The results show that the increased force that is generated when sarcomeres are stretched can be accounted for by changes in the structure of the thick filament. In contrast, the increase in calcium sensitivity that occurs as the sarcomere lengthens is largely due to structural alterations in the thin filament. These two processes can be controlled independently, but work together in the Frank–Starling mechanism.

Now that we better understand the molecular basis of the Frank–Starling mechanism, further work could investigate new strategies for designing and testing treatments for heart disease.

DOI:http://dx.doi.org/10.7554/eLife.24081.002

New Insights in Cardiac Calcium Handling and Excitation-Contraction Coupling

Excitation-contraction (EC) coupling denotes the conversion of electric stimulus in mechanic output in contractile cells. Several studies have demonstrated that calcium (Ca²⁺) plays a pivotal role in this process. Here we present a comprehensive and updated description of the main systems involved in cardiac Ca²⁺ handling that ensure a functional EC coupling and their pathological alterations, mainly related to heart failure.

Mitochondrial complex I dysfunction and altered NAD(P)H kinetics in rat myocardium in cardiac right ventricular hypertrophy and failure

AIMS

In cardiac hypertrophy (CH) and heart failure (HF), alterations occur in mitochondrial enzyme content and activities but the origin and implications of these changes for mitochondrial function need to be resolved.

METHODS AND RESULTS

Right ventricular CH or HF was induced by monocrotaline injection, which causes pulmonary artery hypertension, in rats. Results were compared with saline injection (CON). NAD(P)H and FAD autofluorescence were recorded in thin intact cardiac trabeculae during transitions in stimulation frequency, to assess mitochondrial complex I and complex II function, respectively. Oxygen consumption, mitochondrial morphology, protein content, and enzymatic activity were assessed. NAD(P)H autofluorescence upon an increase in stimulation frequency showed a rapid decline followed by a slow recovery. FAD autofluorescence followed a similar time course, but in opposite direction. The amplitude of the early rapid change in NAD(P)H autofluorescence was severely depressed in CH and HF compared with CON. The rapid changes in FAD autofluorescence in CH and HF were reduced to a lesser extent. Complex I-coupled respiration showed an ∼3.5-fold reduction in CH and HF; complex II-coupled respiration was depressed two-fold in HF. Western blot analyses revealed modest reductions in complex I protein content in CH and HF and in complex I activity in supercomplexes in HF. Mitochondrial volume density was similar, but mitochondrial remodelling was evident from changes in ultrastructure and fusion/fission indices in CH and HF.

CONCLUSION

These results suggest that the alterations in mitochondrial function observed in right ventricular CH and HF can be mainly attributed to complex I dysfunction.