Force-Dependent Recruitment from the Myosin Off State Contributes to Length-Dependent Activation

On the rate-limiting dynamics of force development in muscle

Skeletal muscles produce forces relatively slowly compared with the action potentials that excite them. The dynamics of force production are governed by multiple processes, such as calcium activation, cycling of cross-bridges between myofilaments, and contraction against elastic tissues and the body. These processes have been included piecemeal in some muscle models, but not integrated to reveal which are the most rate limiting. We therefore examined their integrative contributions to force development in two conventional types of muscle models: Hill-type and cross-bridge. We found that no combination of these processes can self-consistently reproduce classic data such as twitch and tetanus. Rather, additional dynamics are needed following calcium activation and facilitating cross-bridge cycling, such as for cooperative myofilament interaction and reconfiguration. We provisionally lump such processes into a simple first-order model of 'force facilitation dynamics' that integrate into a cross-bridge-type muscle model. The proposed model self-consistently reproduces force development for a range of excitations including twitch and tetanus and electromyography-to-force curves. The model's step response reveals relatively small timing contributions of calcium activation (3%), cross-bridge cycling (3%) and contraction (27%) to overall force development of human quadriceps, with the remainder (67%) explained by force facilitation. The same set of model parameters predicts the change in force magnitude (gain) and timing (phase delay) as a function of excitatory firing rate, or as a function of cyclic contraction frequency. Although experiments are necessary to reveal the dynamics of muscle, integrative models are useful for identifying the main rate-limiting processes.

The structural and functional effects of myosin regulatory light chain phosphorylation are amplified by increases in sarcomere length and [Ca2+]

Precise regulation of sarcomeric contraction is essential for normal cardiac function. The heart must generate sufficient force to pump blood throughout the body, but either inadequate or excessive force can lead to dysregulation and disease. Myosin regulatory light chain (RLC) is a thick-filament protein that binds to the neck of the myosin heavy chain. Post-translational phosphorylation of RLC (RLC-P) by myosin light chain kinase is known to influence acto-myosin interactions, thereby increasing force production and Ca2+-sensitivity of contraction. Here, we investigated the role of RLC-P on cardiac structure and function as sarcomere length and [Ca2+] were altered. We found that at low, non-activating levels of Ca2+, RLC-P contributed to myosin head disorder, though there were no effects on isometric stress production and viscoelastic stiffness. With increases in sarcomere length and Ca2+-activation, the structural changes due to RLC-P become greater, which translates into greater force production, greater viscoelastic stiffness, slowed myosin detachment rates and altered nucleotide handling. Altogether, these data suggest that RLC-P may alter thick-filament structure by releasing ordered, off-state myosin. These more disordered myosin heads are available to bind actin, which could result in greater force production as Ca2+ levels increase. However, prolonged cross-bridge attachment duration due to slower ADP release could delay relaxation long enough to enable cross-bridge rebinding. Together, this work further elucidates the effects of RLC-P in regulating muscle function, thereby promoting a better understanding of thick-filament regulatory contributions to cardiac function in health and disease. KEY POINTS: Myosin regulatory light chain (RLC) is a thick-filament protein in the cardiac sarcomere that can be phosphorylated (RLC-P), and changes in RLC-P are associated with cardiac dysfunction and disease. This study assesses how RLC-P alters cardiac muscle structure and function at different sarcomere lengths and calcium concentrations. At low, non-activating levels of Ca2+, RLC-P contributed to myofilament disorder, though there were no effects on isometric stress production and viscoelastic stiffness. With increases in sarcomere length and Ca2+-activation, the structural changes due to RLC-P become greater, which translates into greater force production, greater viscoelastic stiffness, slower myosin detachment rate and altered cross-bridge nucleotide handling rates. This work elucidates the role of RLC-P in regulating muscle function and facilitates understanding of thick-filament regulatory protein contributions to cardiac function in health and disease.

Multiscale fiber remodeling in the infarcted left ventricle using a stress-based reorientation law

The organization of myofibers and extra cellular matrix within the myocardium plays a significant role in defining cardiac function. When pathological events occur, such as myocardial infarction (MI), this organization can become disrupted, leading to degraded pumping performance. The current study proposes a multiscale finite element (FE) framework to determine realistic fiber distributions in the left ventricle (LV). This is achieved by implementing a stress-based fiber reorientation law, which seeks to align the fibers with local traction vectors, such that contractile force and load bearing capabilities are maximized. By utilizing the total stress (passive and active), both myofibers and collagen fibers are reoriented. Simulations are conducted to predict the baseline fiber configuration in a normal LV as well as the adverse fiber reorientation that occurs due to different size MIs. The baseline model successfully captures the transmural variation of helical fiber angles within the LV wall, as well as the transverse fiber angle variation from base to apex. In the models of MI, the patterns of fiber reorientation in the infarct, border zone, and remote regions closely align with previous experimental findings, with a significant increase in fibers oriented in a left-handed helical configuration and increased dispersion in the infarct region. Furthermore, the severity of fiber reorientation and impairment of pumping performance both showed a correlation with the size of the infarct. The proposed multiscale modeling framework allows for the effective prediction of adverse remodeling and offers the potential for assessing the effectiveness of therapeutic interventions in the future. STATEMENT OF SIGNIFICANCE: The organization of muscle and collagen fibers within the heart plays a significant role in defining cardiac function. This organization can become disrupted after a heart attack, leading to degraded pumping performance. In the current study, we implemented a stress-based fiber reorientation law into a computer model of the heart, which seeks to realign the fibers such that contractile force and load bearing capabilities are maximized. The primary goal was to evaluate the effects of different sized heart attacks. We observed substantial fiber remodeling in the heart, which matched experimental observations. The proposed computational framework allows for the effective prediction of adverse remodeling and offers the potential for assessing the effectiveness of therapeutic interventions in the future.

Cardiac length-dependent activation driven by force-dependent thick-filament dynamics

The length-dependent activation (LDA) of maximum force and calcium sensitivity are established features of cardiac muscle contraction but the dominant underlying mechanisms remain to be fully clarified. Alongside the well-documented regulation of contraction via the thin filaments, experiments have identified an additional force-dependent thick-filament activation, whereby myosin heads parked in a so-called off state become available to generate force. This process produces a feedback effect that may potentially drive LDA. Using biomechanical modeling of a human left-ventricular myocyte, this study investigates the extent to which the off-state dynamics could, by itself, plausibly account for LDA, depending on the specific mathematical formulation of the feedback. We hypothesized four different models of the off-state regulatory feedback based on (A) total force, (B) active force, (C) sarcomere strain, and (D) passive force. We tested if these models could reproduce the isometric steady-state and dynamic LDA features predicted by an earlier published model of a human left-ventricle myocyte featuring purely phenomenological length dependences. The results suggest that only total-force feedback (A) is capable of reproducing the expected behaviors, but that passive tension could provide a length-dependent signal on which to initiate the feedback. Furthermore, by attributing LDA to off-state dynamics, our proposed model also qualitatively reproduces experimentally observed effects of the off-state-stabilizing drug mavacamten. Taken together, these results support off-state dynamics as a plausible primary mechanism underlying LDA.

Multiscale modeling shows how 2'-deoxy-ATP rescues ventricular function in heart failure

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction.

Muscle short-range stiffness behaves like a maxwell element, not a spring: Implications for joint stability

INTRODUCTION

Muscles play a critical role in supporting joints during activities of daily living, owing, in part, to the phenomenon of short-range stiffness. Briefly, when an active muscle is lengthened, bound cross-bridges are stretched, yielding forces greater than what is predicted from the force length relationship. For this reason, short-range stiffness has been proposed as an attractive mechanism for providing joint stability. However, there has yet to be a forward dynamic simulation employing a cross-bridge model, that demonstrates this stabilizing role. Therefore, the purpose of this investigation was to test whether Huxley-type muscle elements, which exhibit short-range stiffness, can stabilize a joint while at constant activation.

METHODS

We analyzed the stability of an inverted pendulum (moment of inertia: 2.7 kg m2) supported by Huxley-type muscle models that reproduce the short-range stiffness phenomenon. We calculated the muscle forces that would provide sufficient short-range stiffness to stabilize the system based in minimizing the potential energy. Simulations consisted of a 50 ms long, 5 Nm square-wave perturbation, with numerical simulations carried out in ArtiSynth.

RESULTS

Despite the initial analysis predicting shared activity of antagonist and agonist muscles to maintain stable equilibrium, the inverted pendulum model was not stable, and did not maintain an upright posture even with fully activated muscles.

DISCUSSION & CONCLUSION

Our simulations suggested that short-range stiffness cannot be solely responsible for joint stability, even for modest perturbations. We argue that short-range stiffness cannot achieve stability because its dynamics do not behave like a typical spring. Instead, an alternative conceptual model for short-range stiffness is that of a Maxwell element (spring and damper in series), which can be obtained as a first-order approximation to the Huxley model. We postulate that the damping that results from short-range stiffness slows down the mechanical response and allows the central nervous system time to react and stabilize the joint. We speculate that other mechanisms, like reflexes or residual force enhancement/depression, may also play a role in joint stability. Joint stability is due to a combination of factors, and further research is needed to fully understand this complex system.

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.

Guidelines for mechanistic modeling and analysis in cardiovascular research

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

Multiscale Finite Element Modeling of Left Ventricular Growth in Simulations of Valve Disease

Multiscale models of the cardiovascular system are emerging as effective tools for investigating the mechanisms that drive ventricular growth and remodeling. These models can predict how molecular-level mechanisms impact organ-level structure and function and could provide new insights that help improve patient care. MyoFE is a multiscale computer framework that bridges molecular and organ-level mechanisms in a finite element model of the left ventricle that is coupled with the systemic circulation. In this study, we extend MyoFE to include a growth algorithm, based on volumetric growth theory, to simulate concentric growth (wall thickening/thinning) and eccentric growth (chamber dilation/constriction) in response to valvular diseases. Specifically in our model, concentric growth is controlled by time-averaged total stress along the fiber direction over a cardiac cycle while eccentric growth responds to time-averaged intracellular myofiber passive stress over a cardiac cycle. The new framework correctly predicted different forms of growth in response to two types of valvular diseases, namely aortic stenosis and mitral regurgitation. Furthermore, the model predicted that LV size and function are nearly restored (reversal of growth) when the disease-mimicking perturbation was removed in the simulations for each valvular disorder. In conclusion, the simulations suggest that time-averaged total stress along the fiber direction and time-averaged intracellular myofiber passive stress can be used to drive concentric and eccentric growth in simulations of valve disease.

Super-relaxed myosins contribute to respiratory muscle hibernation in mechanically ventilated patients

Patients receiving mechanical ventilation in the intensive care unit (ICU) frequently develop contractile weakness of the diaphragm. Consequently, they may experience difficulty weaning from mechanical ventilation, which increases mortality and poses a high economic burden. Because of a lack of knowledge regarding the molecular changes in the diaphragm, no treatment is currently available to improve diaphragm contractility. We compared diaphragm biopsies from ventilated ICU patients (N = 54) to those of non-ICU patients undergoing thoracic surgery (N = 27). By integrating data from myofiber force measurements, x-ray diffraction experiments, and biochemical assays with clinical data, we found that in myofibers isolated from the diaphragm of ventilated ICU patients, myosin is trapped in an energy-sparing, super-relaxed state, which impairs the binding of myosin to actin during diaphragm contraction. Studies on quadriceps biopsies of ICU patients and on the diaphragm of previously healthy mechanically ventilated rats suggested that the super-relaxed myosins are specific to the diaphragm and not a result of critical illness. Exposing slow- and fast-twitch myofibers isolated from the diaphragm biopsies to small-molecule compounds activating troponin restored contractile force in vitro. These findings support the continued development of drugs that target sarcomere proteins to increase the calcium sensitivity of myofibers for the treatment of ICU-acquired diaphragm weakness.

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.

Translating myosin-binding protein C and titin abnormalities to whole-heart function using a novel calcium-contraction coupling model

Mutations in cardiac myosin-binding protein C (cMyBP-C) or titin may respectively lead to hypertrophic (HCM) or dilated (DCM) cardiomyopathies. The mechanisms leading to these phenotypes remain unclear because of the challenge of translating cellular abnormalities to whole-heart and system function. We developed and validated a novel computer model of calcium-contraction coupling incorporating the role of cMyBP-C and titin based on the key assumptions: 1) tension in the thick filament promotes cross-bridge attachment mechanochemically, 2) with increasing titin tension, more myosin heads are unlocked for attachment, and 3) cMyBP-C suppresses cross-bridge attachment. Simulated stationary calcium-tension curves, isotonic and isometric contractions, and quick release agreed with experimental data. The model predicted that a loss of cMyBP-C function decreases the steepness of the calcium-tension curve, and that more compliant titin decreases the level of passive and active tension and its dependency on sarcomere length. Integrating this cellular model in the CircAdapt model of the human heart and circulation showed that a loss of cMyBP-C function resulted in HCM-like hemodynamics with higher left ventricular end-diastolic pressures and smaller volumes. More compliant titin led to higher diastolic pressures and ventricular dilation, suggesting DCM-like hemodynamics. The novel model of calcium-contraction coupling incorporates the role of cMyBP-C and titin. Its coupling to whole-heart mechanics translates changes in cellular calcium-contraction coupling to changes in cardiac pump and circulatory function and identifies potential mechanisms by which cMyBP-C and titin abnormalities may develop into HCM and DCM phenotypes. This modeling platform may help identify distinct mechanisms underlying clinical phenotypes in cardiac diseases.

Assessing Cardiac Contractility From Single Molecules to Whole Hearts

Fundamentally, the heart needs to generate sufficient force and power output to dynamically meet the needs of the body. Cardiomyocytes contain specialized structures referred to as sarcomeres that power and regulate contraction. Disruption of sarcomeric function or regulation impairs contractility and leads to cardiomyopathies and heart failure. Basic, translational, and clinical studies have adapted numerous methods to assess cardiac contraction in a variety of pathophysiological contexts. These tools measure aspects of cardiac contraction at different scales ranging from single molecules to whole organisms. Moreover, these studies have revealed new pathogenic mechanisms of heart disease leading to the development of novel therapies targeting contractility. In this review, the authors explore the breadth of tools available for studying cardiac contractile function across scales, discuss their strengths and limitations, highlight new insights into cardiac physiology and pathophysiology, and describe how these insights can be harnessed for therapeutic candidate development and translational.

Regulating Striated Muscle Contraction: Through Thick and Thin

Force generation in striated muscle is primarily controlled by structural changes in the actin-containing thin filaments triggered by an increase in intracellular calcium concentration. However, recent studies have elucidated a new class of regulatory mechanisms, based on the myosin-containing thick filament, that control the strength and speed of contraction by modulating the availability of myosin motors for the interaction with actin. This review summarizes the mechanisms of thin and thick filament activation that regulate the contractility of skeletal and cardiac muscle. A novel dual-filament paradigm of muscle regulation is emerging, in which the dynamics of force generation depends on the coordinated activation of thin and thick filaments. We highlight the interfilament signaling pathways based on titin and myosin-binding protein-C that couple thin and thick filament regulatory mechanisms. This dual-filament regulation mediates the length-dependent activation of cardiac muscle that underlies the control of the cardiac output in each heartbeat.

Design Principles and Benefits of Spatially Explicit Models of Myofilament Function

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

A multiscale finite element model of left ventricular mechanics incorporating baroreflex regulation

Cardiovascular function is regulated by a short-term hemodynamic baroreflex loop, which tries to maintain arterial pressure at a normal level. In this study, we present a new multiscale model of the cardiovascular system named MyoFE. This framework integrates a mechanistic model of contraction at the myosin level into a finite-element-based model of the left ventricle pumping blood through the systemic circulation. The model is coupled with a closed-loop feedback control of arterial pressure inspired by a baroreflex algorithm previously published by our team. The reflex loop mimics the afferent neuron pathway via a normalized signal derived from arterial pressure. The efferent pathway is represented by a kinetic model that simulates the net result of neural processing in the medulla and cell-level responses to autonomic drive. The baroreflex control algorithm modulates parameters such as heart rate and vascular tone of vessels in the lumped-parameter model of systemic circulation. In addition, it spatially modulates intracellular Ca2+ dynamics and molecular-level function of both the thick and the thin myofilaments in the left ventricle. Our study demonstrates that the baroreflex algorithm can maintain arterial pressure in the presence of perturbations such as acute cases of altered aortic resistance, mitral regurgitation, and myocardial infarction. The capabilities of this new multiscale model will be utilized in future research related to computational investigations of growth and remodeling.

Stretch Harmonizes Sarcomere Strain Across the Cardiomyocyte

BACKGROUND

Increasing cardiomyocyte contraction during myocardial stretch serves as the basis for the Frank-Starling mechanism in the heart. However, it remains unclear how this phenomenon occurs regionally within cardiomyocytes, at the level of individual sarcomeres. We investigated sarcomere contractile synchrony and how intersarcomere dynamics contribute to increasing contractility during cell lengthening.

METHODS

Sarcomere strain and Ca2+ were simultaneously recorded in isolated left ventricular cardiomyocytes during 1 Hz field stimulation at 37 °C, at resting length and following stepwise stretch.

RESULTS

We observed that in unstretched rat cardiomyocytes, differential sarcomere deformation occurred during each beat. Specifically, while most sarcomeres shortened during the stimulus, ≈10% to 20% of sarcomeres were stretched or remained stationary. This nonuniform strain was not traced to regional Ca2+ disparities but rather shorter resting lengths and lower force production in systolically stretched sarcomeres. Lengthening of the cell recruited additional shortening sarcomeres, which increased contractile efficiency as less negative, wasted work was performed by stretched sarcomeres. Given the known role of titin in setting sarcomere dimensions, we next hypothesized that modulating titin expression would alter intersarcomere dynamics. Indeed, in cardiomyocytes from mice with titin haploinsufficiency, we observed greater variability in resting sarcomere length, lower recruitment of shortening sarcomeres, and impaired work performance during cell lengthening.

CONCLUSIONS

Graded sarcomere recruitment directs cardiomyocyte work performance, and harmonization of sarcomere strain increases contractility during cell stretch. By setting sarcomere dimensions, titin controls sarcomere recruitment, and its lowered expression in haploinsufficiency mutations impairs cardiomyocyte contractility.

RLC phosphorylation amplifies Ca2+ sensitivity of force in myocardium from cMyBP-C knockout mice

Hypertrophic cardiomyopathy (HCM) is the leading genetic cause of heart disease. The heart comprises several proteins that work together to properly facilitate force production and pump blood throughout the body. Cardiac myosin binding protein-C (cMyBP-C) is a thick-filament protein, and mutations in cMyBP-C are frequently linked with clinical cases of HCM. Within the sarcomere, the N-terminus of cMyBP-C likely interacts with the myosin regulatory light chain (RLC); RLC is a subunit of myosin located within the myosin neck region that modulates contractile dynamics via its phosphorylation state. Phosphorylation of RLC is thought to influence myosin head position along the thick-filament backbone, making it more favorable to bind the thin filament of actin and facilitate force production. However, little is known about how these two proteins interact. We tested the effects of RLC phosphorylation on Ca2+-regulated contractility using biomechanical assays on skinned papillary muscle strips isolated from cMyBP-C KO mice and WT mice. RLC phosphorylation increased Ca2+ sensitivity of contraction (i.e., pCa50) from 5.80 ± 0.02 to 5.95 ± 0.03 in WT strips, whereas RLC phosphorylation increased Ca2+ sensitivity of contraction from 5.86 ± 0.02 to 6.15 ± 0.03 in cMyBP-C KO strips. These data suggest that the effects of RLC phosphorylation on Ca2+ sensitivity of contraction are amplified when cMyBP-C is absent from the sarcomere. This implies that cMyBP-C and RLC act in concert to regulate contractility in healthy hearts, and mutations to these proteins that lead to HCM (or a loss of phosphorylation with disease progression) may disrupt important interactions between these thick-filament regulatory proteins.

Muscle Mechanics and Thick Filament Activation: An Emerging Two-Way Interaction for the Vertebrate Striated Muscle Fine Regulation

Contraction in striated muscle is classically described as regulated by calcium-mediated structural changes in the actin-containing thin filaments, which release the binding sites for the interaction with myosin motors to produce force. In this view, myosin motors, arranged in the thick filaments, are basically always ready to interact with the thin filaments, which ultimately regulate the contraction. However, a new “dual-filament” activation paradigm is emerging, where both filaments must be activated to generate force. Growing evidence from the literature shows that the thick filament activation has a role on the striated muscle fine regulation, and its impairment is associated with severe pathologies. This review is focused on the proposed mechanical feedback that activates the inactive motors depending on the level of tension generated by the active ones, the so-called mechanosensing mechanism. Since the main muscle function is to generate mechanical work, the implications on muscle mechanics will be highlighted, showing: (i) how non-mechanical modulation of the thick filament activation influences the contraction, (ii) how the contraction influences the activation of the thick filament and (iii) how muscle, through the mechanical modulation of the thick filament activation, can regulate its own mechanics. This description highlights the crucial role of the emerging bi-directional feedback on muscle mechanical performance.

Increased length-dependent activation of human engineered heart tissue after chronic α1A-adrenergic agonist treatment: testing a novel heart failure therapy

Chronic stimulation of cardiac α1A-adrenergic receptors (α1A-ARs) improves symptoms in multiple preclinical models of heart failure. However, the translational significance remains unclear. Human engineered heart tissues (EHTs) provide a means of quantifying the effects of chronic α1A-AR stimulation on human cardiomyocyte physiology. EHTs were created from thin slices of decellularized pig myocardium seeded with human induced pluripotent stem cell (iPSC)-derived cardiomyocytes and fibroblasts. With a paired experimental design, EHTs were cultured for 3 wk, mechanically tested, cultured again for 2 wk with α1A-AR agonist A61603 (10 nM) or vehicle control, and retested after drug washout for 24 h. Separate control experiments determined the effects of EHT age (3-5 wk) or repeat mechanical testing. We found that chronic A61603 treatment caused a 25% increase of length-dependent activation (LDA) of contraction compared with vehicle treatment (n = 7/group, P = 0.035). EHT force was not increased after chronic A61603 treatment. However, after vehicle treatment, EHT force was increased by 35% relative to baseline testing (n = 7/group, P = 0.022), suggesting EHT maturation. Control experiments suggested that increased EHT force resulted from repeat mechanical testing, not from EHT aging. RNA-seq analysis confirmed that the α1A-AR is expressed in human EHTs and found chronic A61603 treatment affected gene expression in biological pathways known to be activated by α1A-ARs, including the MAP kinase signaling pathway. In conclusion, increased LDA in human EHT after chronic A61603 treatment raises the possibility that chronic stimulation of the α1A-AR might be beneficial for increasing LDA in human myocardium and might be beneficial for treating human heart failure by restoring LDA.NEW & NOTEWORTHY Chronic stimulation of α1A-adrenergic receptors (α1A-ARs) is known to mediate therapeutic effects in animal heart failure models. To investigate the effects of chronic α1A-AR stimulation in human cardiomyocytes, we tested engineered heart tissue (EHT) created with iPSC-derived cardiomyocytes. RNA-seq analysis confirmed human EHT expressed α1A-ARs. Chronic (2 wk) α1A-AR stimulation with A61603 (10 nM) increased length-dependent activation (LDA) of contraction. Chronic α1A-AR stimulation might be beneficial for treating human heart failure by restoring LDA.

A multiscale model of the cardiovascular system that regulates arterial pressure via closed loop baroreflex control of chronotropism, cell-level contractility, and vascular tone

Multiscale models of the cardiovascular system can provide new insights into physiological and pathological processes. PyMyoVent is a computer model that bridges from molecular- to organ-level function and which simulates a left ventricle pumping blood through the systemic circulation. Initial work with PyMyoVent focused on the end-systolic pressure volume relationship and ranked potential therapeutic strategies by their impact on contractility. This manuscript extends the PyMyoVent framework by adding closed-loop feedback control of arterial pressure. The control algorithm mimics important features of the physiological baroreflex and was developed as part of a long-term program that focuses on growth and biological remodeling. Inspired by the underlying biology, the reflex algorithm uses an afferent signal derived from arterial pressure to drive a kinetic model that mimics the net result of neural processing in the medulla and cell-level responses to autonomic drive. The kinetic model outputs control signals that are constrained between limits that represent maximum parasympathetic and maximum sympathetic drive and which modulate heart rate, intracellular Ca2+ dynamics, the molecular-level function of both the thick and the thin myofilaments, and vascular tone. Simulations show that the algorithm can regulate mean arterial pressure at user-defined setpoints as well as maintaining arterial pressure when challenged by changes in blood volume and/or valve resistance. The reflex also regulates arterial pressure when cell-level contractility is modulated to mimic the idealized impact of myotropes. These capabilities will be important for future work that uses computer modeling to investigate clinical conditions and treatments.

Dilated cardiomyopathy mutation E525K in human beta-cardiac myosin stabilizes the interacting-heads motif and super-relaxed state of myosin

The auto-inhibited, super-relaxed (SRX) state of cardiac myosin is thought to be crucial for regulating contraction, relaxation, and energy conservation in the heart. We used single ATP turnover experiments to demonstrate that a dilated cardiomyopathy (DCM) mutation (E525K) in human beta-cardiac myosin increases the fraction of myosin heads in the SRX state (with slow ATP turnover), especially in physiological ionic strength conditions. We also utilized FRET between a C-terminal GFP tag on the myosin tail and Cy3ATP bound to the active site of the motor domain to estimate the fraction of heads in the closed, interacting-heads motif (IHM); we found a strong correlation between the IHM and SRX state. Negative stain electron microscopy and 2D class averaging of the construct demonstrated that the E525K mutation increased the fraction of molecules adopting the IHM. Overall, our results demonstrate that the E525K DCM mutation may reduce muscle force and power by stabilizing the auto-inhibited SRX state. Our studies also provide direct evidence for a correlation between the SRX biochemical state and the IHM structural state in cardiac muscle myosin. Furthermore, the E525 residue may be implicated in crucial electrostatic interactions that modulate this conserved, auto-inhibited conformation of myosin.

Cardiac efficiency and Starling's Law of the Heart

The formulation by Starling of The Law of the Heart states that 'the [mechanical] energy of contraction, however measured, is a function of the length of the muscle fibre'. Starling later also stated that 'the oxygen consumption of the isolated heart … is determined by its diastolic volume, and therefore by the initial length of its muscular fibres'. This phrasing has motivated us to extend Starling's Law of the Heart to include consideration of the efficiency of contraction. In this study, we assessed both mechanical efficiency and crossbridge efficiency by studying the heat output of isolated rat ventricular trabeculae performing force-length work-loops over ranges of preload and afterload. The combination of preload and afterload allowed us, using our modelling frameworks for the end-systolic zone and the heat-force zone, to simulate cases by recreating physiologically feasible loading conditions. We found that across all cases examined, both work output and change of enthalpy increased with initial muscle length; hence it can only be that the former increases more than the latter to yield increased mechanical efficiency. In contrast, crossbridge efficiency increased with initial muscle length in cases where the extent of muscle shortening varied greatly with preload. We conclude that the efficiency of cardiac contraction increases with increasing initial muscle length and preload. An implication of our conclusion is that the length-dependent activation mechanism underlying the cellular basis of Starling's Law of the Heart is an energetically favourable process that increases the efficiency of cardiac contraction. KEY POINTS: Ernest Starling in 1914 formulated the Law of the Heart to describe the mechanical property of cardiac muscle whereby force of contraction increases with muscle length. He subsequently, in 1927, showed that the oxygen consumption of the heart is also a function of the length of the muscle fibre, but left the field unclear as to whether cardiac efficiency follows the same dependence. A century later, the field has gained an improved understanding of the factors, including the distinct effects of preload and afterload, that affect cardiac efficiency. This understanding presents an opportunity for us to investigate the elusive length-dependence of cardiac efficiency. We found that, by simulating physiologically feasible loading conditions using a mechano-energetics framework, cardiac efficiency increased with initial muscle length. A broader physiological importance of our findings is that the underlying cellular basis of Starling's Law of the Heart is an energetically favourable process that yields increased efficiency.

Reduced cardiac muscle power with low ATP simulating heart failure

For patients with heart failure, myocardial ATP level can be reduced to one-half of that observed in healthy controls. This marked reduction (from ≈8 mM in healthy controls to as low as 3-4 mM in heart failure) has been suggested to contribute to impaired myocardial contraction and to the decreased pump function characteristic of heart failure. However, in vitro measures of maximum myofilament force generation, maximum shortening velocity, and the actomyosin ATPase activity show effective KM values for MgATP ranging from ≈10 μM to 150 μM, well below the intracellular ATP level in heart failure. Thus, it is not clear that the fall of myocardial ATP observed in heart failure is sufficient to impair the function of the contractile proteins. Therefore, we tested the effect of low MgATP levels on myocardial contraction using demembranated cardiac muscle preparations that were exposed to MgATP levels typical of the range found in non-failing and failing hearts. Consistent with previous studies, we found that a 50% reduction in MgATP level (from 8 mM to 4 mM) did not reduce maximum force generation or maximum velocity of shortening. However, we found that a 50% reduction in MgATP level caused a 20%-25% reduction in maximal power generation (measured during muscle shortening against a load) and a 20% slowing of cross-bridge cycling kinetics. These results suggest that the decreased cellular ATP level occurring in heart failure contributes to the impaired pump function of the failing heart. Since the ATP-myosin ATPase dissociation constant is estimated to be submillimolar, these findings also suggest that MgATP concentration affects cross-bridge dynamics through a mechanism that is more complex than through the direct dependence of MgATP concentration on myosin ATPase activity. Finally, these studies suggest that therapies targeted to increase adenine nucleotide pool levels in cardiomyocytes might be beneficial for treating heart failure.

FiberSim: A flexible open-source model of myofilament-level contraction

FiberSim is a flexible open-source model of myofilament-level contraction. The code uses a spatially explicit technique, meaning that it tracks the position and status of each contractile molecule within the lattice framework. This allows the model to simulate some of the mechanical effects modulated by myosin-binding protein C, as well as the dose dependence of myotropes and the effects of varying isoform expression levels. This paper provides a short introduction to FiberSim and presents simulations of tension-pCa curves with and without regulation of thick and thin filament activation by myosin-binding protein C. A myotrope dose-dependent response as well as slack/re-stretch maneuvers to assess rates of tension recovery are also presented. The software was designed to be flexible (the user can define their own model and/or protocol) and computationally efficient (simulations can be performed on a regular laptop). We hope that other investigators will use FiberSim to explore myofilament level mechanisms and to accelerate research focusing on the contractile properties of sarcomeres.

Generative adversarial networks for construction of virtual populations of mechanistic models: simulations to study Omecamtiv Mecarbil action

Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM’s mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force–calcium (F–Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system.

Frank-Starling mechanism, fluid responsiveness, and length-dependent activation: Unravelling the multiscale behaviors with an in silico analysis

The Frank-Starling mechanism is a fundamental regulatory property which underlies the cardiac output adaptation to venous filling. Length-dependent activation is generally assumed to be the cellular origin of this mechanism. At the heart scale, it is commonly admitted that an increase in preload (ventricular filling) leads to an increased cellular force and an increased volume of ejected blood. This explanation also forms the basis for vascular filling therapy. It is actually difficult to unravel the exact nature of the relationship between length-dependent activation and the Frank-Starling mechanism, as three different scales (cellular, ventricular and cardiovascular) are involved. Mathematical models are powerful tools to overcome these limitations. In this study, we use a multiscale model of the cardiovascular system to untangle the three concepts (length-dependent activation, Frank-Starling, and vascular filling). We first show that length-dependent activation is required to observe both the Frank-Starling mechanism and a positive response to high vascular fillings. Our results reveal a dynamical length dependent activation-driven response to changes in preload, which involves interactions between the cellular, ventricular and cardiovascular levels and thus highlights fundamentally multiscale behaviors. We show however that the cellular force increase is not enough to explain the cardiac response to rapid changes in preload. We also show that the absence of fluid responsiveness is not related to a saturating Frank-Starling effect. As it is challenging to study those multiscale phenomena experimentally, this computational approach contributes to a more comprehensive knowledge of the sophisticated length-dependent properties of cardiac muscle.

Multiscale simulations of left ventricular growth and remodeling

Cardiomyocytes can adapt their size, shape, and orientation in response to altered biomechanical or biochemical stimuli. The process by which the heart undergoes structural changes-affecting both geometry and material properties-in response to altered ventricular loading, altered hormonal levels, or mutant sarcomeric proteins is broadly known as cardiac growth and remodeling (G&R). Although it is likely that cardiac G&R initially occurs as an adaptive response of the heart to the underlying stimuli, prolonged pathological changes can lead to increased risk of atrial fibrillation, heart failure, and sudden death. During the past few decades, computational models have been extensively used to investigate the mechanisms of cardiac G&R, as a complement to experimental measurements. These models have provided an opportunity to quantitatively study the relationships between the underlying stimuli (primarily mechanical) and the adverse outcomes of cardiac G&R, i.e., alterations in ventricular size and function. State-of-the-art computational models have shown promise in predicting the progression of cardiac G&R. However, there are still limitations that need to be addressed in future works to advance the field. In this review, we first outline the current state of computational models of cardiac growth and myofiber remodeling. Then, we discuss the potential limitations of current models of cardiac G&R that need to be addressed before they can be utilized in clinical care. Finally, we briefly discuss the next feasible steps and future directions that could advance the field of cardiac G&R.

The Super-Relaxed State and Length Dependent Activation in Porcine Myocardium

[Figure: see text].

Spontaneous myogenic fasciculation associated with the lengthening of cardiac muscle in response to static preloading

Force enhancement is one kind of myogenic spontaneous fasciculation in lengthening preload striated muscles. In cardiac muscle, the role of this biomechanical event is not well established. The physiological passive property is an essential part for maintaining normal diastole in the heart. In excessive preload heart, force enhancement relative erratic passive properties may cause muscle decompensating, implicate in the development of diastolic dysfunction. In this study, the force enhancement occurrence in mouse cardiac papillary muscle was evaluated by a microstepping stretch method. The intracellular Ca²⁺ redistribution during occurrence of force enhancement was monitored in real-time by a Flou-3 (2 mM) indicator. The force enhancement amplitude, the enhancement of the prolongation time, and the tension-time integral were analyzed by myography. The results indicated that the force enhancement occurred immediately after active stretching and was rapidly enhanced during sustained static stretch. The presence of the force and the increase in the amplitude synchronized with the acquisition and immediate transfer of Ca²⁺ to adjacent fibres. In highly preloaded fibres, the enhancement exceeded the maximum passive tension (from 4.49 ± 0.43 N/mm² to 6.20 ± 0.51 N/mm²). The occurrence of force enhancement were unstable in each static stretch. The increased enhancement amplitude combined with the reduced prolongation time to induce a reduction in the tension-time integral. We concluded that intracellular Ca²⁺-synchronized force enhancement is one kind of interruption event in excessive preload cardiac muscle. During the cardiac muscle in its passive relaxation period, the occurrence of this interruption affected the rhythmic stability of the cardiac relaxation cycle.

Hypertrophic cardiomyopathy β-cardiac myosin mutation (P710R) leads to hypercontractility by disrupting super relaxed state

Hypertrophic cardiomyopathy (HCM) is the most common inherited form of heart disease, associated with over 1,000 mutations, many in β-cardiac myosin (MYH7). Molecular studies of myosin with different HCM mutations have revealed a diversity of effects on ATPase and load-sensitive rate of detachment from actin. It has been difficult to predict how such diverse molecular effects combine to influence forces at the cellular level and further influence cellular phenotypes. This study focused on the P710R mutation that dramatically decreased in vitro motility velocity and actin-activated ATPase, in contrast to other MYH7 mutations. Optical trap measurements of single myosin molecules revealed that this mutation reduced the step size of the myosin motor and the load sensitivity of the actin detachment rate. Conversely, this mutation destabilized the super relaxed state in longer, two-headed myosin constructs, freeing more heads to generate force. Micropatterned human induced pluripotent derived stem cell (hiPSC)-cardiomyocytes CRISPR-edited with the P710R mutation produced significantly increased force (measured by traction force microscopy) compared with isogenic control cells. The P710R mutation also caused cardiomyocyte hypertrophy and cytoskeletal remodeling as measured by immunostaining and electron microscopy. Cellular hypertrophy was prevented in the P710R cells by inhibition of ERK or Akt. Finally, we used a computational model that integrated the measured molecular changes to predict the measured traction forces. These results confirm a key role for regulation of the super relaxed state in driving hypercontractility in HCM with the P710R mutation and demonstrate the value of a multiscale approach in revealing key mechanisms of disease.

Hypertrophic cardiomyopathy β-cardiac myosin mutation (P710R) leads to hypercontractility by disrupting super relaxed state

Hypertrophic cardiomyopathy (HCM) is the most common inherited form of heart disease, associated with over 1,000 mutations, many in β-cardiac myosin (MYH7). Molecular studies of myosin with different HCM mutations have revealed a diversity of effects on ATPase and load-sensitive rate of detachment from actin. It has been difficult to predict how such diverse molecular effects combine to influence forces at the cellular level and further influence cellular phenotypes. This study focused on the P710R mutation that dramatically decreased in vitro motility velocity and actin-activated ATPase, in contrast to other MYH7 mutations. Optical trap measurements of single myosin molecules revealed that this mutation reduced the step size of the myosin motor and the load sensitivity of the actin detachment rate. Conversely, this mutation destabilized the super relaxed state in longer, two-headed myosin constructs, freeing more heads to generate force. Micropatterned human induced pluripotent derived stem cell (hiPSC)-cardiomyocytes CRISPR-edited with the P710R mutation produced significantly increased force (measured by traction force microscopy) compared with isogenic control cells. The P710R mutation also caused cardiomyocyte hypertrophy and cytoskeletal remodeling as measured by immunostaining and electron microscopy. Cellular hypertrophy was prevented in the P710R cells by inhibition of ERK or Akt. Finally, we used a computational model that integrated the measured molecular changes to predict the measured traction forces. These results confirm a key role for regulation of the super relaxed state in driving hypercontractility in HCM with the P710R mutation and demonstrate the value of a multiscale approach in revealing key mechanisms of disease.

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.

Stress-dependent activation of myosin in the heart requires thin filament activation and thick filament mechanosensing

The efficiency of the heart as a pump depends on an autoregulatory mechanism, the Frank–Starling law of the heart, that potentiates the strength of contraction in response to an increase in ventricular filling. Disruption of this mechanism compromises the ability of the heart to pump blood, potentially leading to heart failure. We used fluorescent probes on myosin in heart muscle cells to investigate the molecular basis of the Frank–Starling mechanism. Our results show that the stronger contraction of heart muscle at longer lengths is due to a calcium-dependent interfilament signaling pathway that links stress sensing in the myosin-containing filaments with calcium activation of the actin-containing filaments. This pathway can potentially be targeted for treating heart failure.

Myosin-based regulation in the heart muscle modulates the number of myosin motors available for interaction with calcium-regulated thin filaments, but the signaling pathways mediating the stronger contraction triggered by stretch between heartbeats or by phosphorylation of the myosin regulatory light chain (RLC) remain unclear. Here, we used RLC probes in demembranated cardiac trabeculae to investigate the molecular structural basis of these regulatory pathways. We show that in relaxed trabeculae at near-physiological temperature and filament lattice spacing, the RLC-lobe orientations are consistent with a subset of myosin motors being folded onto the filament surface in the interacting-heads motif seen in isolated filaments. The folded conformation of myosin is disrupted by cooling relaxed trabeculae, similar to the effect induced by maximal calcium activation. Stretch or increased RLC phosphorylation in the physiological range have almost no effect on RLC conformation at a calcium concentration corresponding to that between beats. These results indicate that in near-physiological conditions, the folded myosin motors are not directly switched on by RLC phosphorylation or by the titin-based passive tension at longer sarcomere lengths in the absence of thin filament activation. However, at the higher calcium concentrations that activate the thin filaments, stretch produces a delayed activation of folded myosin motors and force increase that is potentiated by RLC phosphorylation. We conclude that the increased contractility of the heart induced by RLC phosphorylation and stretch can be explained by a calcium-dependent interfilament signaling pathway involving both thin filament sensitization and thick filament mechanosensing.

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.

Mavacamten preserves length-dependent contractility and improves diastolic function in human engineered heart tissue

Comprehensive functional characterization of cardiac tissue includes investigation of length and load dependence. Such measurements have been slow to develop in engineered heart tissues (EHTs), whose mechanical characterizations have been limited primarily to isometric and near-isometric behaviors. A more realistic assessment of myocardial function would include force-velocity curves to characterize power output and force-length loops mimicking the cardiac cycle to characterize work output. We developed a system that produces force-velocity curves and work loops in human EHTs using an adaptive iterative control scheme. We used human EHTs in this system to perform a detailed characterization of the cardiac β-myosin specific inhibitor, mavacamten. Consistent with the clinically proposed application of this drug to treat hypertrophic cardiomyopathy, our data support the premise that mavacamten improves diastolic function through reduction of diastolic stiffness and isometric relaxation time. Meanwhile, the effects of mavacamten on length- and load-dependent muscle performance were mixed. The drug attenuated the length-dependent response at small stretch values but showed normal length dependency at longer lengths. Peak power output of mavacamten-treated EHTs showed reduced power output as expected but also shifted peak power output to a lower load. Here, we demonstrate a robust method for the generation of isotonic contraction series and work loops in engineered heart tissues using an adaptive-iterative method. This approach reveals new features of mavacamten pharmacology, including previously unappreciated effects on intrinsic myosin dynamics and preservation of Frank-Starling behavior at longer muscle lengths.NEW & NOTEWORTHY We applied innovative methods to comprehensively characterize the length and load-dependent behaviors of engineered human cardiac muscle when treated with the cardiac β-myosin specific inhibitor mavacamten, a drug on the verge of clinical implementation for hypertrophic cardiomyopathy. We find mechanistic support for the role of mavacamten in improving diastolic function of cardiac tissue and note novel effects on work and power.

Effects of mavacamten on Ca2+ sensitivity of contraction as sarcomere length varied in human myocardium

BACKGROUND AND PURPOSE

Heart failure can reflect impaired contractile function at the myofilament level. In healthy hearts, myofilaments become more sensitive to Ca2+ as cells are stretched. This represents a fundamental property of the myocardium that contributes to the Frank-Starling response, although the molecular mechanisms underlying the effect remain unclear. Mavacamten, which binds to myosin, is under investigation as a potential therapy for heart disease. We investigated how mavacamten affects the sarcomere-length dependence of Ca2+ -sensitive isometric contraction to determine how mavacamten might modulate the Frank-Starling mechanism.

EXPERIMENTAL APPROACH

Multicellular preparations from the left ventricular-free wall of hearts from organ donors were chemically permeabilized and Ca2+ activated in the presence or absence of 0.5-μM mavacamten at 1.9 or 2.3-μm sarcomere length (37°C). Isometric force and frequency-dependent viscoelastic myocardial stiffness measurements were made.

KEY RESULTS

At both sarcomere lengths, mavacamten reduced maximal force and Ca2+ sensitivity of contraction. In the presence and absence of mavacamten, Ca2+ sensitivity of force increased as sarcomere length increased. This suggests that the length-dependent activation response was maintained in human myocardium, even though mavacamten reduced Ca2+ sensitivity. There were subtle effects of mavacamten reducing force values under relaxed conditions (pCa 8.0), as well as slowing myosin cross-bridge recruitment and speeding cross-bridge detachment under maximally activated conditions (pCa 4.5).

CONCLUSION AND IMPLICATIONS

Mavacamten did not eliminate sarcomere length-dependent increases in the Ca2+ sensitivity of contraction in myocardial strips from organ donors at physiological temperature. Drugs that modulate myofilament function may be useful therapies for cardiomyopathies.

Fast-relaxing cardiomyocytes exert a dominant role in the relaxation behavior of heterogeneous myocardium

Substantial variation in relaxation rate exists among cardiomyocytes within small volumes of myocardium; however, it is unknown how this variability affects the overall relaxation mechanics of heart muscle. In this study, we sought to modulate levels of cellular heterogeneity in a computational model, then validate those predictions using an engineered heart tissue platform. We formulated an in silico tissue model composed of half-sarcomeres with varied relaxation rates, incorporating single-cell cardiomyocyte experimental data. These model tissues randomly sampled relaxation parameters from two offset distributions of fast- and slow-relaxing populations of half-sarcomeres. Isometric muscle twitch simulations predicted a complex relationship between relaxation time and the proportion of fast-versus slow-relaxing cells in heterogeneous tissues. Specifically, a 50/50 mixture of fast and slow cells did not lead to relaxation time that was the mean of the relaxation times associated with the two pure cases. Rather, the mean relaxation time was achieved at a ratio of 70:30 slow:fast relaxing cells, suggesting a disproportionate impact of fast-relaxing cells on overall tissue relaxation. To examine whether this behavior persists in vitro, we constructed engineered heart tissues from two lines of fast- and slow-relaxing human iPSC-derived cardiomyocytes. Cell tracking via fluorescent nanocrystals confirmed the presence of both cell populations in the 50/50 mixed tissues at the time of mechanical characterization. Isometric muscle twitch relaxation times of these mixed-population engineered heart tissues showed agreement with the predictions from the model, namely that the measured relaxation rate of 50/50 mixed tissues more closely resembled that of tissues made with 100% fast-relaxing cells. Our observations suggest that cardiomyocyte diversity can play an important role in determining tissue-level relaxation.

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.

Multiscale Modeling of Cardiovascular Function Predicts That the End-Systolic Pressure Volume Relationship Can Be Targeted via Multiple Therapeutic Strategies

Most patients who develop heart failure are unable to elevate their cardiac output on demand due to impaired contractility and/or reduced ventricular filling. Despite decades of research, few effective therapies for heart failure have been developed. In part, this may reflect the difficulty of predicting how perturbations to molecular-level mechanisms that are induced by drugs will scale up to modulate system-level properties such as blood pressure. Computer modeling might help with this process and thereby accelerate the development of better therapies for heart failure. This manuscript presents a new multiscale model that uses a single contractile element to drive an idealized ventricle that pumps blood around a closed circulation. The contractile element was formed by linking an existing model of dynamically coupled myofilaments with a well-established model of myocyte electrophysiology. The resulting framework spans from molecular-level events (including opening of ion channels and transitions between different myosin states) to properties such as ejection fraction that can be measured in patients. Initial calculations showed that the model reproduces many aspects of normal cardiovascular physiology including, for example, pressure-volume loops. Subsequent sensitivity tests then quantified how each model parameter influenced a range of system level properties. The first key finding was that the End Systolic Pressure Volume Relationship, a classic index of cardiac contractility, was ∼50% more sensitive to parameter changes than any other system-level property. The second important result was that parameters that primarily affect ventricular filling, such as passive stiffness and Ca²⁺ reuptake via sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA), also have a major impact on systolic properties including stroke work, myosin ATPase, and maximum ventricular pressure. These results reinforce the impact of diastolic function on ventricular performance and identify the End Systolic Pressure Volume Relationship as a particularly sensitive system-level property that can be targeted using multiple therapeutic strategies.

Impaired Myocardial Energetics Causes Mechanical Dysfunction in Decompensated Failing Hearts

Cardiac mechanical function is supported by ATP hydrolysis, which provides the chemical-free energy to drive the molecular processes underlying cardiac pumping. Physiological rates of myocardial ATP consumption require the heart to resynthesize its entire ATP pool several times per minute. In the failing heart, cardiomyocyte metabolic dysfunction leads to a reduction in the capacity for ATP synthesis and associated free energy to drive cellular processes. Yet it remains unclear if and how metabolic/energetic dysfunction that occurs during heart failure affects mechanical function of the heart. We hypothesize that changes in phosphate metabolite concentrations (ATP, ADP, inorganic phosphate) that are associated with decompensation and failure have direct roles in impeding contractile function of the myocardium in heart failure, contributing to the whole-body phenotype. To test this hypothesis, a transverse aortic constriction (TAC) rat model of pressure overload, hypertrophy, and decompensation was used to assess relationships between metrics of whole-organ pump function and myocardial energetic state. A multiscale computational model of cardiac mechanoenergetic coupling was used to identify and quantify the contribution of metabolic dysfunction to observed mechanical dysfunction. Results show an overall reduction in capacity for oxidative ATP synthesis fueled by either fatty acid or carbohydrate substrates as well as a reduction in total levels of adenine nucleotides and creatine in myocardium from TAC animals compared to sham-operated controls. Changes in phosphate metabolite levels in the TAC rats are correlated with impaired mechanical function, consistent with the overall hypothesis. Furthermore, computational analysis of myocardial metabolism and contractile dynamics predicts that increased levels of inorganic phosphate in TAC compared to control animals kinetically impair the myosin ATPase crossbridge cycle in decompensated hypertrophy/heart failure.

The role of pacing rate in the modulation of mechano-induced immediate and delayed changes in the force and Ca-transient of cardiac muscle

Myocardial function is tuned by dynamic changes in length and load via mechano-calcium feedback. This regulation may be significantly affected by heart rhythm. We evaluated the mechano-induced modulation of contractility and Ca-transient (CaT) in the rat myocardium subjected to twitch-by-twitch shortening-re-lengthening (↓-↑) trains of different lengths (N = 1 … 720 cycles) at low (1 Hz) and near-physiological (3.5 Hz) pacing rates. Force/CaT characteristics were evaluated in the first post-train isometric twitch (immediate effect) and during slow changes (delayed maximal elevation/decrease) and compared with those of the pre-train twitch. The immediate inotropic effect was positive for N = 30 … 720 and negative for N = 1 … 20, while the delayed effect was always positive. The immediate and delayed inotropic effects were significantly higher at 3.5-Hz vs 1-Hz (P < 0.05). The prominent inotropism was accompanied by much smaller changes in the CaT diastolic level/amplitude. The shortening-re-lengthening train induced oscillations of the slow change in force at 3.5-Hz (always) and at 1-Hz (∼50% of muscles), which were dependent of the train length and independent of the pacing rate. We suggest that twitch-by-twitch shortening-re-lengthening of cardiac muscle decreases Ca²⁺ buffering by troponin C and elevates Ca²⁺ loading of the sarcoplasmic reticulum (SR); the latter cumulatively depends on the train length. A high pacing rate intensifies the cumulative transient shift in the SR Ca²⁺ loading, augmenting the post-train inotropic response and prolonging its recovery to the pre-train level. The pacing-dependent mechano-induced inotropic effects remain to be elucidated in the myocardium with impaired Ca handling.

Force-dependent recruitment from myosin OFF-state increases end-systolic pressure-volume relationship in left ventricle

Finite element (FE) modeling is becoming increasingly prevalent in the world of cardiac mechanics; however, many existing FE models are phenomenological and thus do not capture cellular-level mechanics. This work implements a cellular-level contraction scheme into an existing nonlinear FE code to model ventricular contraction. Specifically, this contraction model incorporates three myosin states: OFF-, ON-, and an attached force-generating state. It has been speculated that force-dependent transitions from the OFF- to ON-state may contribute to length-dependent activation at the cellular level. The current work investigates the contribution of force-dependent recruitment out of the OFF-state to ventricular-level function, specifically the Frank-Starling relationship, as seen through the end-systolic pressure-volume relationship (ESPVR). Five FE models were constructed using geometries of rat left ventricles obtained via cardiac magnetic resonance imaging. FE simulations were conducted to optimize parameters for the cellular contraction model such that the differences between FE predicted ventricular pressures for the models and experimentally measured pressures were minimized. The models were further validated by comparing FE predicted end-systolic strain to experimentally measured strain. Simulations mimicking vena cava occlusion generated descending pressure volume loops from which ESPVRs were calculated. In simulations with the inclusion of the OFF-state, using a force-dependent transition to the ON-state, the ESPVR calculated was steeper than in simulations excluding the OFF-state. Furthermore, the ESPVR was also steeper when compared to models that included the OFF-state without a force-dependent transition. This suggests that the force-dependent recruitment of thick filament heads from the OFF-state at the cellular level contributes to the Frank-Starling relationship observed at the organ level.

Length-Dependent Activation of Contractility and Ca-Transient Kinetics in Auxotonically Contracting Isolated Rat Ventricular Cardiomyocytes

Length-dependent activation (LDA) of contraction is an important mechanism of proper myocardial function that is often blunted in diseases accompanied by deficient contractility and impaired calcium homeostasis. We evaluated how the extent of LDA is related to the decreased force in healthy rat myocardium under negative inotropic conditions that affect the calcium cycle. The length-dependent effects on auxotonic twitch and Ca-transient were compared in isolated rat ventricular cardiomyocytes at room temperature (“25C”) and near-physiological temperature (“35C”) in normal Tyrode and at 25°C with thapsigargin-depleted sarcoplasmic reticulum (“25C + Thap”). At the slack length, a similar negative inotropy in “35C” and “25C + Thap” was accompanied by totally different changes in Ca-transient amplitude, time-to-peak, and time-to-decline from peak to 50% amplitude. End-systolic/end-diastolic tension-sarcomere length relationships were obtained for each individual cell, and the ratio of their slopes, the dimensionless Frank-Starling Gain index, was 2.32 ± 0.16, 1.78 ± 0.09, and 1.37 ± 0.06 in “25C,” “35C” and “25C + Thap,” respectively (mean ± S.E.M.). Ca-transient diastolic level and amplitude did not differ between “25C” and “35C” at any SL, but in “35C” it developed and declined significantly faster. In contrast, thapsigargin-induced depletion of SERCA2a significantly attenuated and retarded Ca-transient. The relative amount of Ca²⁺ utilized by troponin C, evaluated by the integral magnitude of a short-lived component of Ca-transient decline (“bump”), increased by ~25% per each 0.05 μm increase in SL in all groups. The kinetics of the Ca-TnC dissociation, evaluated by the bump time-to-peak, was significantly faster in “35C” and slower in “25C + Thap” vs. “25C” (respectively, 63.7 ± 5.3 and 253.6 ± 8.3% of the value in “25C,” mean ± S.E.M.). In conclusion, a similar inotropic effect can be observed in rat ventricular myocardium under totally different kinetics of free cytosolic calcium. The extent of LDA is not determined by actual peak systolic tension but is regulated by the level of peak systolic calcium and the kinetics of Ca-transient decline which, in turn, are governed by Ca-TnC dissociation and Ca²⁺ reuptake by the sarcoplasmic reticulum. Altogether, these findings constitute new evidence about the role of the length-dependent modulation of Ca²⁺ homeostasis in the mechanisms of calcium regulation of contraction and mechano-calcium feedback in the myocardium.

Arrhythmogenic Current Generation by Myofilament-Triggered Ca2+ Release and Sarcomere Heterogeneity

Heterogeneous mechanical dyskinesis has been implicated in many arrhythmogenic phenotypes. Strain-dependent perturbations to cardiomyocyte electrophysiology may contribute to this arrhythmogenesis through processes referred to as mechanoelectric feedback. Although the role of stretch-activated ion currents has been investigated using computational models, experimental studies suggest that mechanical strain may also promote arrhythmia by facilitating calcium wave propagation. To investigate whether strain-dependent changes in calcium affinity to the myofilament may promote arrhythmogenic intracellular calcium waves, we modified a mathematical model of rabbit excitation-contraction coupling coupled to a model of myofilament activation and force development. In a one-dimensional compartmental analysis, we bidirectionally coupled 50 sarcomere models in series to model calcium diffusion and stress transfer between adjacent sarcomeres. These considerations enabled the model to capture 1) the effects of mechanical feedback on calcium homeostasis at the sarcomeric level and 2) the combined effects of mechanical and calcium heterogeneities at the cellular level. The results suggest that in conditions of calcium overload, the vulnerable window of stretch-release to trigger suprathreshold delayed afterdepolarizations can be affected by heterogeneity in sarcomere length. Furthermore, stretch and sarcomere heterogeneity may modulate the susceptibility threshold for delayed afterdepolarizations and the aftercontraction wave propagation velocity.

The mechanism of thin filament regulation: Models in conflict?

Evidence on two- and three-state models of the calcium regulation models of muscle contractions remain in favor of three-state models.

In a recent JGP article, Heeley et al. (2019. J. Gen. Physiol. https://doi.org/10.1085/jgp.201812198) reopened the debate about two- versus three-state models of thin filament regulation. The authors review their work, which measures the rate constant of P_i release from myosin.ADP.Pi activated by actin or thin filaments under a variety of conditions. They conclude that their data can be described by a two-state model and raise doubts about the generally accepted three-state model as originally formulated by McKillop and Geeves (1993. Biophys. J.https://doi.org/10.1016/S0006-3495(93)81110-X). However, in the following article, we follow Plato’s dictum that “twice and thrice over, as they say, good it is to repeat and review what is good.” We have therefore reviewed the evidence for the three- and two-state models and present our view that the evidence is overwhelmingly in favor of three structural states of the thin filament, which regulate access of myosin to its binding sites on actin and, hence, muscle contractility.