Nuclear mechanosignaling in striated muscle diseases.

Mechanosignaling describes processes by which biomechanical stimuli are transduced into cellular responses. External biophysical forces can be transmitted via structural protein networks that span from the cellular membrane to the cytoskeleton and the nucleus, where they can regulate gene expression through a series of biomechanical and/or biochemical mechanosensitive mechanisms, including chromatin remodeling, translocation of transcriptional regulators, and epigenetic factors. Striated muscle cells, including cardiac and skeletal muscle myocytes, utilize these nuclear mechanosignaling mechanisms to respond to changes in their intracellular and extracellular mechanical environment and mediate gene expression and cell remodeling. In this brief review, we highlight and discuss recent experimental work focused on the pathway of biomechanical stimulus propagation at the nucleus-cytoskeleton interface of striated muscles, and the mechanisms by which these pathways regulate gene regulation, muscle structure, and function. Furthermore, we discuss nuclear protein mutations that affect mechanosignaling function in human and animal models of cardiomyopathy. Furthermore, current open questions and future challenges in investigating striated muscle nuclear mechanosignaling are further discussed.

Machine learning meets Monte Carlo methods for models of muscle's molecular machinery to classify mutations.

The timing and magnitude of force generation by a muscle depend on complex interactions in a compliant, contractile filament lattice. Perturbations in these interactions can result in cardiac muscle diseases. In this study, we address the fundamental challenge of connecting the temporal features of cardiac twitches to underlying rate constants and their perturbations associated with genetic cardiomyopathies. Current state-of-the-art metrics for characterizing the mechanical consequence of cardiac muscle disease do not utilize information embedded in the complete time course of twitch force. We pair dimension reduction techniques and machine learning methods to classify underlying perturbations that shape the timing of twitch force. To do this, we created a large twitch dataset using a spatially explicit Monte Carlo model of muscle contraction. Uniquely, we modified the rate constants of this model in line with mouse models of cardiac muscle disease and varied mutation penetrance. Ultimately, the results of this study show that machine learning models combined with biologically informed dimension reduction techniques can yield excellent classification accuracy of underlying muscle perturbations.

Subcellular Remodeling in Filamin C Deficient Mouse Hearts Impairs Myocyte Tension Development during Progression of Dilated Cardiomyopathy.

Dilated cardiomyopathy (DCM) is a life-threatening form of heart disease that is typically characterized by progressive thinning of the ventricular walls, chamber dilation, and systolic dysfunction. Multiple mutations in the gene encoding filamin C (FLNC), an actin-binding cytoskeletal protein in cardiomyocytes, have been found in patients with DCM. However, the mechanisms that lead to contractile impairment and DCM in patients with FLNC variants are poorly understood. To determine how FLNC regulates systolic force transmission and DCM remodeling, we used an inducible, cardiac-specific FLNC-knockout (icKO) model to produce a rapid onset of DCM in adult mice. Loss of FLNC reduced systolic force development in single cardiomyocytes and isolated papillary muscles but did not affect twitch kinetics or calcium transients. Electron and immunofluorescence microscopy showed significant defects in Z-disk alignment in icKO mice and altered myofilament lattice geometry. Moreover, a loss of FLNC induces a softening myocyte cortex and structural adaptations at the subcellular level that contribute to disrupted longitudinal force production during contraction. Spatially explicit computational models showed that these structural defects could be explained by a loss of inter-myofibril elastic coupling at the Z-disk. Our work identifies FLNC as a key regulator of the multiscale ultrastructure of cardiomyocytes and therefore plays an important role in maintaining systolic mechanotransmission pathways, the dysfunction of which may be key in driving progressive DCM.

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 Ca2+ 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.

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.

Modulating the tension-time integral of the cardiac twitch prevents dilated cardiomyopathy in murine hearts.

Dilated cardiomyopathy (DCM) is often associated with sarcomere protein mutations that confer reduced myofilament tension-generating capacity. We demonstrated that cardiac twitch tension-time integrals can be targeted and tuned to prevent DCM remodeling in hearts with contractile dysfunction. We employed a transgenic murine model of DCM caused by the D230N-tropomyosin (Tm) mutation and designed a sarcomere-based intervention specifically targeting the twitch tension-time integral of D230N-Tm hearts using multiscale computational models of intramolecular and intermolecular interactions in the thin filament and cell-level contractile simulations. Our models predicted that increasing the calcium sensitivity of thin filament activation using the cardiac troponin C (cTnC) variant L48Q can sufficiently augment twitch tension-time integrals of D230N-Tm hearts. Indeed, cardiac muscle isolated from double-transgenic hearts expressing D230N-Tm and L48Q cTnC had increased calcium sensitivity of tension development and increased twitch tension-time integrals compared with preparations from hearts with D230N-Tm alone. Longitudinal echocardiographic measurements revealed that DTG hearts retained normal cardiac morphology and function, whereas D230N-Tm hearts developed progressive DCM. We present a computational and experimental framework for targeting molecular mechanisms governing the twitch tension of cardiomyopathic hearts to counteract putative mechanical drivers of adverse remodeling and open possibilities for tension-based treatments of genetic cardiomyopathies.

Contracting striated muscle has a dynamic I-band spring with an undamped stiffness 100 times larger than the passive stiffness.

KEY POINTS: Fast sarcomere-level mechanics in contracting intact fibres from frog skeletal muscle reveal an I-band spring with an undamped stiffness 100 times larger than the known static stiffness. This undamped stiffness remains constant in the range of sarcomere length 2.7-3.1 µm, showing the ability of the I-band spring to adapt its length to the width of the I-band. The stiffness and tunability of the I-band spring implicate titin as a force contributor that, during contraction, allows weaker half-sarcomeres to equilibrate with in-series stronger half-sarcomeres, preventing the development of sarcomere length inhomogeneity. This work opens new possibilities for the detailed in situ description of the structural-functional basis of muscle dysfunctions related to mutations or site-directed mutagenesis in titin that alter the I-band stiffness. ABSTRACT: Force and shortening in the muscle sarcomere are due to myosin motors from thick filaments pulling nearby actin filaments toward the sarcomere centre. Thousands of serially linked sarcomeres in muscle make the shortening (and the shortening speed) macroscopic, while the intrinsic instability of in-series force generators is likely prevented by the cytoskeletal protein titin that connects the thick filament with the sarcomere end, working as an I-band spring that accounts for the rise of passive force with sarcomere length (SL). However, current estimates of titin stiffness, deduced from the passive force-SL relation and single molecule mechanics, are much smaller than what is required to avoid the development of large inhomogeneities among sarcomeres. In this work, using 4 kHz stiffness measurements on a population of sarcomeres selected along an intact fibre isolated from frog skeletal muscle contracting at different SLs (temperature 4°C), we measure the undamped stiffness of an I-band spring that at SL > 2.7 µm attains a maximum constant value of ∼6 pN nm-1 per half-thick filament, two orders of magnitude larger than expected from titin-related passive force. We conclude that a titin-like dynamic spring in the I-band, made by an undamped elastic element in-series with damped elastic elements, adapts its length to the SL with kinetics that provide force balancing among serially linked sarcomeres during contraction. In this way, the I-band spring plays a fundamental role in preventing the development of SL inhomogeneity.

Multiscale Models of Cardiac Muscle Biophysics and Tissue Remodeling in Hypertrophic Cardiomyopathies.

Myocardial hypertrophy is the result of sustained perturbations to the mechanical and/or neurohormonal homeostasis of cardiac cells and is driven by integrated, multiscale biophysical and biochemical processes that are currently not well defined. In this brief review, we highlight recent computational and experimental models of cardiac hypertrophy that span mechanisms from the molecular level to the tissue level. Specifically, we focus on: (i) molecular-level models of the structural dynamics of sarcomere proteins in hypertrophic hearts, (ii) cellular-level models of excitation-contraction coupling and mechanosensitive signaling in disease-state myocytes, and (iii) organ-level models of myocardial growth kinematics and predictors thereof. Finally, we discuss how spanning these scales and combining multiple experimental/computational models will provide new information about the processes governing hypertrophy and potential methods to prevent or reverse them.

Cardiac myosin activation with 2-deoxy-ATP via increased electrostatic interactions with actin.

The naturally occurring nucleotide 2-deoxy-adenosine 5'-triphosphate (dATP) can be used by cardiac muscle as an alternative energy substrate for myosin chemomechanical activity. We and others have previously shown that dATP increases contractile force in normal hearts and models of depressed systolic function, but the structural basis of these effects has remained unresolved. In this work, we combine multiple techniques to provide structural and functional information at the angstrom-nanometer and millisecond time scales, demonstrating the ability to make both structural measurements and quantitative kinetic estimates of weak actin-myosin interactions that underpin sarcomere dynamics. Exploiting dATP as a molecular probe, we assess how small changes in myosin structure translate to electrostatic-based changes in sarcomere function to augment contractility in cardiac muscle. Through Brownian dynamics simulation and computational structural analysis, we found that deoxy-hydrolysis products [2-deoxy-adenosine 5'-diphosphate (dADP) and inorganic phosphate (Pi)] bound to prepowerstroke myosin induce an allosteric restructuring of the actin-binding surface on myosin to increase the rate of cross-bridge formation. We then show experimentally that this predicted effect translates into increased electrostatic interactions between actin and cardiac myosin in vitro. Finally, using small-angle X-ray diffraction analysis of sarcomere structure, we demonstrate that the proposed increased electrostatic affinity of myosin for actin causes a disruption of the resting conformation of myosin motors, resulting in their repositioning toward the thin filament before activation. The dATP-mediated structural alterations in myosin reported here may provide insight into an improved criterion for the design or selection of small molecules to be developed as therapeutic agents to treat systolic 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 Ca2+-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.

A Spatially Explicit Model Shows How Titin Stiffness Modulates Muscle Mechanics and Energetics.

In striated muscle, the giant protein titin spans the entire length of a half-sarcomere and extends from the backbone of the thick filament, reversibly attaches to the thin filaments, and anchors to the dense protein network of the z-disk capping the end of the half-sarcomere. However, little is known about the relationship between the basic mechanical properties of titin and muscle contractility. Here, we build upon our previous multi-filament, spatially explicit computational model of the half-sarcomere by incorporating the nonlinear mechanics of titin filaments in the I-band. We vary parameters of the nonlinearity to understand the effects of titin stiffness on contraction dynamics and efficiency. We do so by simulating isometric contraction for a range of sarcomere lengths (SLs; 1.6-3.25 µm). Intermediate values of titin stiffness accurately reproduce the passive force-SL relation for skeletal muscle. The maximum force-SL relation is not affected by titin for SL≤2.5 µm. However, as titin stiffness increases, maximum force for the four thick filament system at SL = 3.0 µm significantly decreases from 103.2 ± 2 to 58.8 ± 1 pN. Additionally, by monitoring ATP consumption, we measure contraction efficiency as a function of titin stiffness. We find that at SL = 3.0 µm, efficiency significantly decreases from 13.9 ± 0.4 to 7.0 ± 0.3 pN/ATP when increasing titin stiffness, with little or no effect below 2.5 µm. Taken together, our results suggest that, despite an increase in the fraction of motors bound to actin along the descending limb when titin is stiffer, the force-generating capacity of the motors is reduced. These results suggest that titin stiffness has the potential to affect contractile efficiency.

Afterload promotes maturation of human induced pluripotent stem cell derived cardiomyocytes in engineered heart tissues.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) grown in engineered heart tissue (EHT) can be used for drug screening, disease modeling, and heart repair. However, the immaturity of hiPSC-CMs currently limits their use. Because mechanical loading increases during development and facilitates cardiac maturation, we hypothesized that afterload would promote maturation of EHTs. To test this we developed a system in which EHTs are suspended between a rigid post and a flexible one, whose resistance to contraction can be modulated by applying braces of varying length. These braces allow us to adjust afterload conditions over two orders of magnitude by increasing the flexible post resistance from 0.09 up to 9.2 µN/µm. After three weeks in culture, optical tracking of post deflections revealed that auxotonic twitch forces increased in correlation with the degree of afterload, whereas twitch velocities decreased with afterload. Consequently, the power and work of the EHTs were maximal under intermediate afterloads. When studied isometrically, the inotropy of EHTs increased with afterload up to an intermediate resistance (0.45 µN/µm) and then plateaued. Applied afterload increased sarcomere length, cardiomyocyte area and elongation, which are hallmarks of maturation. Furthermore, progressively increasing the level of afterload led to improved calcium handling, increased expression of several key markers of cardiac maturation, including a shift from fetal to adult ventricular myosin heavy chain isoforms. However, at the highest afterload condition, markers of pathological hypertrophy and fibrosis were also upregulated, although the bulk tissue stiffness remained the same for all levels of applied afterload tested. Together, our results indicate that application of moderate afterloads can substantially improve the maturation of hiPSC-CMs in EHTs, while high afterload conditions may mimic certain aspects of human cardiac pathology resulting from elevated mechanical overload.

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.

Minimum number of myosin motors accounting for shortening velocity under zero load in skeletal muscle.

KEY POINTS: Myosin filament mechanosensing determines the efficiency of the contraction by adapting the number of switched ON motors to the load. Accordingly, the unloaded shortening velocity (V0 ) is already set at the end of latency relaxation (LR), ∼10 ms after the start of stimulation, when the myosin filament is still in the OFF state. Here the number of actin-attached motors per half-myosin filament (n) during V0 shortening imposed either at the end of LR or at the plateau of the isometric contraction is estimated from the relation between half-sarcomere compliance and force during the force redevelopment after shortening. The value of n decreases progressively with shortening and, during V0 shortening starting at the end of LR, is 1-4. Reduction of n is accounted for by a constant duty ratio of 0.05 and a parallel switching OFF of motors, explaining the very low rate of ATP utilization found during unloaded shortening. ABSTRACT: The maximum velocity at which a skeletal muscle can shorten (i.e. the velocity of sliding between the myosin filament and the actin filament under zero load, V0 ) is already set at the end of the latency relaxation (LR) preceding isometric force generation, ∼10 ms after the start of electrical stimulation in frog muscle fibres at 4°C. At this time, Ca2+ -induced activation of the actin filament is maximal, while the myosin filament is in the OFF state characterized by most of the myosin motors lying on helical tracks on the filament surface, making them unavailable for actin binding and ATP hydrolysis. Here, the number of actin-attached motors per half-thick filament during V0 shortening (n) is estimated by imposing, on tetanized single fibres from Rana esculenta (at 4°C and sarcomere length 2.15 µm), small 4 kHz oscillations and determining the relation between half-sarcomere (hs) compliance and force during the force development following V0 shortening. When V0 shortening is superimposed on the maximum isometric force T0 , n decreases progressively with the increase of shortening (range 30-80 nm per hs) and, when V0 shortening is imposed at the end of LR, n can be as low as 1-4. Reduction of n is accounted for by a constant duty ratio of the myosin motor of ∼0.05 and a parallel switching OFF of the thick filament, providing an explanation for the very low rate of ATP utilization during extended V0 shortening.