GSK-3ß Localizes to the Cardiac Z-Disc to Maintain Length Dependent Activation.

BACKGROUND: Altered kinase localization is gaining appreciation as a mechanism of cardiovascular disease. Previous work suggests GSK-3ß (glycogen synthase kinase 3ß) localizes to and regulates contractile function of the myofilament. We aimed to discover GSK-3ß's in vivo role in regulating myofilament function, the mechanisms involved, and the translational relevance. METHODS: Inducible cardiomyocyte-specific GSK-3ß knockout mice and left ventricular myocardium from nonfailing and failing human hearts were studied. RESULTS: Skinned cardiomyocytes from knockout mice failed to exhibit calcium sensitization with stretch indicating a loss of length-dependent activation (LDA), the mechanism underlying the Frank-Starling Law. Titin acts as a length sensor for LDA, and knockout mice had decreased titin stiffness compared with control mice, explaining the lack of LDA. Knockout mice exhibited no changes in titin isoforms, titin phosphorylation, or other thin filament phosphorylation sites known to affect passive tension or LDA. Mass spectrometry identified several z-disc proteins as myofilament phospho-substrates of GSK-3ß. Agreeing with the localization of its targets, GSK-3ß that is phosphorylated at Y216 binds to the z-disc. We showed pY216 was necessary and sufficient for z-disc binding using adenoviruses for wild-type, Y216F, and Y216E GSK-3ß in neonatal rat ventricular cardiomyocytes. One of GSK-3ß's z-disc targets, abLIM-1 (actin-binding LIM protein 1), binds to the z-disc domains of titin that are important for maintaining passive tension. Genetic knockdown of abLIM-1 via siRNA in human engineered heart tissues resulted in enhancement of LDA, indicating abLIM-1 may act as a negative regulator that is modulated by GSK-3ß. Last, GSK-3ß myofilament localization was reduced in left ventricular myocardium from failing human hearts, which correlated with depressed LDA. CONCLUSIONS: We identified a novel mechanism by which GSK-3ß localizes to the myofilament to modulate LDA. Importantly, z-disc GSK-3ß levels were reduced in patients with heart failure, indicating z-disc localized GSK-3ß is a possible therapeutic target to restore the Frank-Starling mechanism in patients with heart failure.

Signaling network model of cardiomyocyte morphological changes in familial cardiomyopathy.

Familial cardiomyopathy is a precursor of heart failure and sudden cardiac death. Over the past several decades, researchers have discovered numerous gene mutations primarily in sarcomeric and cytoskeletal proteins causing two different disease phenotypes: hypertrophic (HCM) and dilated (DCM) cardiomyopathies. However, molecular mechanisms linking genotype to phenotype remain unclear. Here, we employ a systems approach by integrating experimental findings from preclinical studies (e.g., murine data) into a cohesive signaling network to scrutinize genotype to phenotype mechanisms. We developed an HCM/DCM signaling network model utilizing a logic-based differential equations approach and evaluated model performance in predicting experimental data from four contexts (HCM, DCM, pressure overload, and volume overload). The model has an overall prediction accuracy of 83.8%, with higher accuracy in the HCM context (90%) than DCM (75%). Global sensitivity analysis identifies key signaling reactions, with calcium-mediated myofilament force development and calcium-calmodulin kinase signaling ranking the highest. A structural revision analysis indicates potential missing interactions that primarily control calcium regulatory proteins, increasing model prediction accuracy. Combination pharmacotherapy analysis suggests that downregulation of signaling components such as calcium, titin and its associated proteins, growth factor receptors, ERK1/2, and PI3K-AKT could inhibit myocyte growth in HCM. In experiments with patient-specific iPSC-derived cardiomyocytes (MLP-W4R;MYH7-R723C iPSC-CMs), combined inhibition of ERK1/2 and PI3K-AKT rescued the HCM phenotype, as predicted by the model. In DCM, PI3K-AKT-NFAT downregulation combined with upregulation of Ras/ERK1/2 or titin or Gq protein could ameliorate cardiomyocyte morphology. The model results suggest that HCM mutations that increase active force through elevated calcium sensitivity could increase ERK activity and decrease eccentricity through parallel growth factors, Gq-mediated, and titin pathways. Moreover, the model simulated the influence of existing medications on cardiac growth in HCM and DCM contexts. This HCM/DCM signaling model demonstrates utility in investigating genotype to phenotype mechanisms in familial cardiomyopathy.

Muscle LIM Protein Force-Sensing Mediates Sarcomeric Biomechanical Signaling in Human Familial Hypertrophic Cardiomyopathy.

BACKGROUND: Familial hypertrophic cardiomyopathy (HCM) is the most common inherited cardiac disease and is typically caused by mutations in genes encoding sarcomeric proteins that regulate cardiac contractility. HCM manifestations include left ventricular hypertrophy and heart failure, arrythmias, and sudden cardiac death. How dysregulated sarcomeric force production is sensed and leads to pathological remodeling remains poorly understood in HCM, thereby inhibiting the efficient development of new therapeutics. METHODS: Our discovery was based on insights from a severe phenotype of an individual with HCM and a second genetic alteration in a sarcomeric mechanosensing protein. We derived cardiomyocytes from patient-specific induced pluripotent stem cells and developed robust engineered heart tissues by seeding induced pluripotent stem cell-derived cardiomyocytes into a laser-cut scaffold possessing native cardiac fiber alignment to study human cardiac mechanobiology at both the cellular and tissue levels. Coupled with computational modeling for muscle contraction and rescue of disease phenotype by gene editing and pharmacological interventions, we have identified a new mechanotransduction pathway in HCM, shown to be essential in modulating the phenotypic expression of HCM in 5 families bearing distinct sarcomeric mutations. RESULTS: Enhanced actomyosin crossbridge formation caused by sarcomeric mutations in cardiac myosin heavy chain (MYH7) led to increased force generation, which, when coupled with slower twitch relaxation, destabilized the MLP (muscle LIM protein) stretch-sensing complex at the Z-disc. Subsequent reduction in the sarcomeric muscle LIM protein level caused disinhibition of calcineurin-nuclear factor of activated T-cells signaling, which promoted cardiac hypertrophy. We demonstrate that the common muscle LIM protein-W4R variant is an important modifier, exacerbating the phenotypic expression of HCM, but alone may not be a disease-causing mutation. By mitigating enhanced actomyosin crossbridge formation through either genetic or pharmacological means, we alleviated stress at the Z-disc, preventing the development of hypertrophy associated with sarcomeric mutations. CONCLUSIONS: Our studies have uncovered a novel biomechanical mechanism through which dysregulated sarcomeric force production is sensed and leads to pathological signaling, remodeling, and hypertrophic responses. Together, these establish the foundation for developing innovative mechanism-based treatments for HCM that stabilize the Z-disc MLP-mechanosensory complex.

Chronic diastolic stretch unmasks conduction defects in an in vitro model of arrhythmogenic cardiomyopathy.

We seek to elucidate the precise nature of mechanical loading that precipitates conduction deficits in a concealed-phase model of arrhythmogenic cardiomyopathy (ACM). ACM is a progressive disorder often resulting from mutations in desmosomal proteins. Exercise has been shown to worsen disease progression and unmask arrhythmia vulnerability, yet the underlying pathomechanisms may depend on the type and intensity of exercise. Because exercise causes myriad changes to multiple inter-dependent hemodynamic parameters, it is difficult to isolate its effects to specific changes in mechanical load. Here, we use engineered heart tissues (EHTs) with iPSC-derived cardiomyocytes expressing R451G desmoplakin, an ACM-linked mutation, which results in a functionally null model of desmoplakin (DSP). We also use a novel bioreactor to independently perturb tissue strain at different time points during the cardiac cycle. We culture EHTs under three strain regimes: normal physiological shortening; increased diastolic stretch, simulating high preload; and isometric culture, simulating high afterload. DSPR451G EHTs that have been cultured isometrically undergo adaptation, with no change in action potential parameters, conduction velocity, or contractile function, a phenotype confirmed by global proteomic analysis. However, when DSPR451G EHTs are subjected to increased diastolic stretch, they exhibit concomitant reductions in conduction velocity and the expression of connexin-43. These effects are rescued by inhibition of both lysosome activity and ERK signaling. Our results indicate that the response of DSPR451G EHTs to mechanical stimuli depends on the strain and the timing of the applied stimulus, with increased diastolic stretch unmasking conduction deficits in a concealed-phase model of ACM.

Myofilament glycation in diabetes reduces contractility by inhibiting tropomyosin movement, is rescued by cMyBPC domains.

Diabetes doubles the risk of developing heart failure (HF). As the prevalence of diabetes grows, so will HF unless the mechanisms connecting these diseases can be identified. Methylglyoxal (MG) is a glycolysis by-product that forms irreversible modifications on lysine and arginine, called glycation. We previously found that myofilament MG glycation causes sarcomere contractile dysfunction and is increased in patients with diabetes and HF. The aim of this study was to discover the molecular mechanisms by which MG glycation of myofilament proteins cause sarcomere dysfunction and to identify therapeutic avenues to compensate. In humans with type 2 diabetes without HF, we found increased glycation of sarcomeric actin compared to non-diabetics and it correlated with decreased calcium sensitivity. Depressed calcium sensitivity is pathogenic for HF, therefore myofilament glycation represents a promising therapeutic target to inhibit the development of HF in diabetics. To identify possible therapeutic targets, we further defined the molecular actions of myofilament glycation. Skinned myocytes exposed to 100 µM MG exhibited decreased calcium sensitivity, maximal calcium-activated force, and crossbridge kinetics. Replicating MG's functional affects using a computer simulation of sarcomere function predicted simultaneous decreases in tropomyosin's blocked-to-closed rate transition and crossbridge duty cycle were consistent with all experimental findings. Stopped-flow experiments and ATPase activity confirmed MG decreased the blocked-to-closed transition rate. Currently, no therapeutics target tropomyosin, so as proof-of-principal, we used a n-terminal peptide of myosin-binding protein C, previously shown to alter tropomyosin's position on actin. C0C2 completely rescued MG-induced calcium desensitization, suggesting a possible treatment for diabetic HF.

Evidence for synergy between sarcomeres and fibroblasts in an in vitro model of myocardial reverse remodeling.

We have created a novel in-vitro platform to study reverse remodeling of engineered heart tissue (EHT) after mechanical unloading. EHTs were created by seeding decellularized porcine myocardial sections with a mixture of primary neonatal rat ventricular myocytes and cardiac fibroblasts. Each end of the ribbon-like constructs was fixed to a plastic clip, allowing the tissues to be statically stretched or slackened. Inelastic deformation was introduced by stretching tissues by 20% of their original length. EHTs were subsequently unloaded by returning tissues to their original, shorter length. Mechanical characterization of EHTs immediately after unloading and at subsequent time points confirmed the presence of a reverse-remodeling process, through which stress-free tissue length was increased after chronic stretch but gradually decreased back to its original value within 9 days. When a cardiac myosin inhibitor was applied to tissues after unloading, EHTs failed to completely recover their passive and active mechanical properties, suggesting a role for actomyosin contraction in reverse remodeling. Selectively inhibiting cardiomyocyte contraction or fibroblast activity after mechanical unloading showed that contractile activity of both cell types was required to achieve full remodeling. Similar tests with EHTs formed from human induced pluripotent stem cell-derived cardiomyocytes also showed reverse remodeling that was enhanced when treated with omecamtiv mecarbil, a myosin activator. These experiments suggest essential roles for active sarcomeric contraction and fibroblast activity in reverse remodeling of myocardium after mechanical unloading. Our findings provide a mechanistic rationale for designing potential therapies to encourage reverse remodeling in patient hearts.

Potential impacts of the cardiac troponin I mobile domain on myofilament activation and relaxation.

The cardiac thin filament is regulated in a Ca2+-dependent manner through conformational changes of troponin and tropomyosin (Tm). It has been generally understood that under conditions of low Ca2+ the inhibitory peptide domain (IP) of troponin I (TnI) binds to actin and holds Tm over the myosin binding sites on actin to prevent crossbridge formation. More recently, evidence that the C-terminal mobile domain (MD) of TnI also binds actin has made for a more complex scenario. This study uses a computational model to investigate the consequences of assuming that TnI regulates Tm movement via two actin-binding domains rather than one. First, a 16-state model of the cardiac thin filament regulatory unit was created with TnI-IP as the sole regulatory domain. Expansion of this to include TnI-MD formed a 24-state model. Comparison of these models showed that assumption of a second actin-binding site allows the individual domains to have a lower affinity for actin than would be required for IP acting alone. Indeed, setting actin affinities of the IP and MD to 25% of that assumed for the IP in the single-site model was sufficient to achieve precisely the same degree of Ca2+ regulation. We also tested the 24-state model's ability to represent steady-state experimental data in the case of disruption of either the IP or MD. We were able to capture qualitative changes in several properties that matched what was seen in the experimental data. Lastly, simulations were run to examine the effect of disruption of the IP or MD on twitch dynamics. Our results suggest that both domains are required to keep diastolic cross-bridge activity to a minimum and accelerate myofilament relaxation. Overall, our analyses support a paradigm in which two domains of TnI bind with moderate affinity to actin, working in tandem to complete Ca2+-dependent regulation of the thin filament.

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.

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.

Modelling sarcomeric cardiomyopathies with human cardiomyocytes derived from induced pluripotent stem cells.

Cardiomyocytes derived from human induced pluripotent stem cells (iPSCs) provide a unique opportunity to understand the pathophysiological effects of genetic cardiomyopathy mutations. In particular, these cells hold the potential to unmask the effects of mutations on contractile behaviour in vitro, providing new insights into genotype-phenotype relationships. With this goal in mind, several groups have established iPSC lines that contain sarcomeric gene mutations linked to cardiomyopathy in patient populations. Their studies have employed diverse systems and methods for performing mechanical measurements of contractility, ranging from single cell techniques to multicellular tissue-like constructs. Here, we review published results to date within the growing field of iPSC-based sarcomeric cardiomyopathy disease models. We devote special attention to the methods of mechanical characterization selected in each case, and how these relate to the paradigms of classical muscle mechanics. An appreciation of these somewhat subtle paradigms can inform efforts to compare the results of different studies and possibly reconcile discrepancies. Although more work remains to be done to improve and possibly standardize methods for producing, maturing, and mechanically interrogating iPSC-derived cardiomyocytes, the initial results indicate that this approach to modelling cardiomyopathies will continue to provide critical insights into these devastating diseases.

Contractile work directly modulates mitochondrial protein levels in human engineered heart tissues.

Engineered heart tissues (EHTs) have emerged as a robust in vitro model to study cardiac physiology. Although biomimetic culture environments have been developed to better approximate in vivo conditions, currently available methods do not permit full recapitulation of the four phases of the cardiac cycle. We have developed a bioreactor which allows EHTs to undergo cyclic loading sequences that mimic in vivo work loops. EHTs cultured under these working conditions exhibited enhanced concentric contractions but similar isometric contractions compared with EHTs cultured isometrically. EHTs that were allowed to shorten cyclically in culture had increased capacity for contractile work when tested acutely. Increased work production was correlated with higher levels of mitochondrial proteins and mitochondrial biogenesis; this effect was eliminated when tissues were cyclically shortened in the presence of a myosin ATPase inhibitor. Leveraging our novel in vitro method to precisely apply mechanical loads in culture, we grew EHTs under two loading regimes prescribing the same work output but with different associated afterloads. These groups showed no difference in mitochondrial protein expression. In loading regimes with the same afterload but different work output, tissues subjected to higher work demand exhibited elevated levels of mitochondrial protein. Our findings suggest that regulation of mitochondrial mass in cultured human EHTs is potently modulated by the mechanical work the tissue is permitted to perform in culture, presumably communicated through ATP demand. Precise application of mechanical loads to engineered heart tissues in culture represents a novel in vitro method for studying physiological and pathological cardiac adaptation.NEW & NOTEWORTHY In this work, we present a novel bioreactor that allows for active length control of engineered heart tissues during extended tissue culture. Specific length transients were designed so that engineered heart tissues generated complete cardiac work loops. Chronic culture with various work loops suggests that mitochondrial mass and biogenesis are directly regulated by work output.

A Microwell Cell Capture Device Reveals Variable Response to Dobutamine in Isolated Cardiomyocytes.

Isolated ventricular cardiomyocytes exhibit substantial cell-to-cell variability, even when obtained from the same small volume of myocardium. In this study, we investigated the possibility that cardiomyocyte responses to ß-adrenergic stimulus are also highly heterogeneous. To achieve the throughput and measurement duration desired for these experiments, we designed and validated a novel microwell system that immobilizes and uniformly orients isolated adult cardiomyocytes. In this configuration, detailed drug responses of dozens of cells can be followed for intervals exceeding 1 h. At the conclusion of an experiment, specific cells can also be harvested via a precision aspirator for single-cell gene expression profiling. Using this system, we followed changes in Ca2+ signaling and contractility of individual cells under sustained application of either dobutamine or omecamtiv mecarbil. Both compounds increased average cardiomyocyte contractility over the course of an hour, but responses of individual cells to dobutamine were significantly more variable. Surprisingly, some dobutamine-treated cardiomyocytes augmented Ca2+ release without increasing contractility. Other cells responded with increased contractility despite unchanged Ca2+ release. Single-cell gene expression analysis revealed significant co-expression of ß-adrenergic pathway genes PKA regulatory subunit type I, PKA regulatory subunit type II, and Ca2+/calmodulin-dependent protein kinase II across cardiomyocytes. Other data supported a connection between the effects of dobutamine on relaxation rate and the expression of protein phosphatase 2. These findings suggest that variable drug responses among cells are not merely experimental artifacts. By enabling direct comparison of the functional behavior of an individual cell and the genes it expresses, this new system constitutes a unique tool for interrogating cardiomyocyte drug responses and discovering the genes that modulate them.

The Effect of Tropomyosin Mutations on Actin-Tropomyosin Binding: In Search of Lost Time.

The initial binding of tropomyosin onto actin filaments and then its polymerization into continuous cables on the filament surface must be precisely tuned to overall thin-filament structure, function, and performance. Low-affinity interaction of tropomyosin with actin has to be sufficiently strong to localize the tropomyosin on actin, yet not so tight that regulatory movement on filaments is curtailed. Likewise, head-to-tail association of tropomyosin molecules must be favorable enough to promote tropomyosin cable formation but not so tenacious that polymerization precedes filament binding. Arguably, little molecular detail on early tropomyosin binding steps has been revealed since Wegner's seminal studies on filament assembly almost 40 years ago. Thus, interpretation of mutation-based actin-tropomyosin binding anomalies leading to cardiomyopathies cannot be described fully. In vitro, tropomyosin binding is masked by explosive tropomyosin polymerization once cable formation is initiated on actin filaments. In contrast, in silico analysis, characterizing molecular dynamics simulations of single wild-type and mutant tropomyosin molecules on F-actin, is not complicated by tropomyosin polymerization at all. In fact, molecular dynamics performed here demonstrates that a midpiece tropomyosin domain is essential for normal actin-tropomyosin interaction and that this interaction is strictly conserved in a number of tropomyosin mutant species. Elsewhere along these mutant molecules, twisting and bending corrupts the tropomyosin superhelices as they "lose their grip" on F-actin. We propose that residual interactions displayed by these mutant tropomyosin structures with actin mimic ones that occur in early stages of thin-filament generation, as if the mutants are recapitulating the assembly process but in reverse. We conclude therefore that an initial binding step in tropomyosin assembly onto actin involves interaction of the essential centrally located domain.

A short history of the development of mathematical models of cardiac mechanics.

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

Diverse relaxation rates exist among rat cardiomyocytes isolated from a single myocardial region.

KEY POINTS: Prior studies have shown variation in the functional properties of cardiomyocytes isolated from different regions of the left ventricular myocardium. We found that these region-dependent variations vanish below a tissue volume of ∼7 mm3 in the adult rat myocardium, revealing a fixed level of intrinsic relaxation rate heterogeneity that is independent of tissue volume. Within these microscopically varying cell populations, fast-relaxing cells were shown to have elevated phosphorylated troponin I compared to slow-relaxing cells. Relaxation rate was also correlated with cardiomyocyte length, in that slow-relaxing cells were longer than fast-relaxing cells. These results show a new relationship between cardiomyocyte morphology and myofilament relaxation, and suggest that functional diversity among individual myocytes at the microscale may contribute to bulk relaxation of the myocardium. ABSTRACT: The mean contractility and calcium handling properties of cardiomyocytes isolated from different regions of the ventricular myocardium are known to vary significantly. We designed experiments to quantify the variance in contractile properties among cells within the same myocardial region. Longitudinal strips of myocardial tissue were excised from the epicardial left ventricular free walls of adult Sprague-Dawley rats and then treated with collagenase to isolate individual myocytes. Cardiomyocytes were characterized by measuring sarcomere length changes and calcium transients during electrical pacing. Variance of the time from peak sarcomere shortening to 50% re-lengthening (RT50 ) was assessed in each cell population. Isolating cells from progressively shorter strips allowed an estimate of the myocardial volume below which regional variation vanished and only microscale heterogeneity remained (∼7 mm3 ). The SD of RT50 within this myocardial volume was 28% of the mean. In a series of follow-up experiments, RT50 was shown to correlate significantly with resting myocyte length, suggesting a connection between cell morphology and intrinsic relaxation behaviour. To explore the mechanistic basis of varying RT50 , a novel single-cell aspirator was employed to collect small batches of cardiomyocytes grouped according to their relaxation rates (fast or slow). Western blot analysis of the two groups revealed significantly elevated troponin I phosphorylation in fast-relaxing cells. Our observations suggest that cell-to-cell heterogeneity of active contractile properties is substantial, with implications for how we understand myocardial relaxation and design drug therapies intended to alter relaxation rate.

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

Cardiac muscle develops more force when it is activated at longer lengths. The concentration of Ca2+ required to develop half-maximal force also decreases. These effects are known as length-dependent activation and are thought to play critical roles in the Frank-Starling relationship and cardiovascular homeostasis. The molecular mechanisms underpinning length-dependent activation remain unclear, but recent experiments suggest that they may include recruitment of myosin heads from the off (sometimes called super-relaxed) state. This manuscript presents a mathematical model of muscle contraction that was developed to investigate this hypothesis. Myosin heads in the model transitioned between an off state (that could not interact with actin), an on state (that could bind to actin), and a single attached state. Simulations were fitted to experimental data using multidimensional parameter optimization. Statistical analysis showed that a model in which the rate of the off-to-on transition increased linearly with force reproduced the length-dependent behavior of chemically permeabilized myocardium better than a model with a constant off-to-on transition rate (F-test, p < 0.001). This result suggests that the thick-filament transitions are modulated by force. Additional calculations showed that the model incorporating a mechanosensitive thick filament could also reproduce twitch responses measured in a trabecula stretched to different lengths. A final set of simulations was then used to test the model. These calculations predicted how reducing passive stiffness would impact the length dependence of the calcium sensitivity of contractile force. The prediction (a 60% reduction in ΔpCa50) mimicked the 58% reduction in ΔpCa50 in myocardium from rats that expressed a giant isoform of titin and had low resting tension. Together, these computational results suggest that force-dependent recruitment of myosin heads from the thick-filament off state contributes to length-dependent activation and the Frank-Starling relationship.

Precise Binding of Tropomyosin on Actin Involves Sequence-Dependent Variance in Coiled-Coil Twisting.

Often considered an archetypal dimeric coiled coil, tropomyosin nonetheless exhibits distinctive "noncanonical" core residues located at the hydrophobic interface between its component α-helices. Notably, a charged aspartate, D137, takes the place of nonpolar residues otherwise present. Much speculation has been offered to rationalize potential local coiled-coil instability stemming from D137 and its effect on regulatory transitions of tropomyosin over actin filaments. Although experimental approaches such as electron cryomicroscopy reconstruction are optimal for defining average tropomyosin positions on actin filaments, to date, these methods have not captured the dynamics of tropomyosin residues clustered around position 137 or elsewhere. In contrast, computational biochemistry, involving molecular dynamics simulation, is a compelling choice to extend the understanding of local and global tropomyosin behavior on actin filaments at high resolution. Here, we report on molecular dynamics simulation of actin-free and actin-associated tropomyosin, showing noncanonical residue D137 as a locus for tropomyosin twist variation, with marked effects on actin-tropomyosin interactions. We conclude that D137-sponsored coiled-coil twisting is likely to optimize electrostatic side-chain contacts between tropomyosin and actin on the assembled thin filament, while offsetting disparities between tropomyosin pseudorepeat and actin subunit periodicities. We find that D137 has only minor local effects on tropomyosin coiled-coil flexibility, (i.e., on its flexural mobility). Indeed, D137-associated overtwisting may actually augment tropomyosin stiffness on actin filaments. Accordingly, such twisting-induced stiffness of tropomyosin is expected to enhance cooperative regulatory translocation of the tropomyosin cable over actin.

Decreased polycystin 2 expression alters calcium-contraction coupling and changes ß-adrenergic signaling pathways.

Cardiac disorders are the main cause of mortality in autosomal-dominant polycystic kidney disease (ADPKD). However, how mutated polycystins predispose patients with ADPKD to cardiac pathologies before development of renal dysfunction is unknown. We investigate the effect of decreased levels of polycystin 2 (PC2), a calcium channel that interacts with the ryanodine receptor, on myocardial function. We hypothesize that heterozygous PC2 mice (Pkd2(+/-)) undergo cardiac remodeling as a result of changes in calcium handling, separate from renal complications. We found that Pkd2(+/-) cardiomyocytes have altered calcium handling, independent of desensitized calcium-contraction coupling. Paradoxically, in Pkd2(+/-) mice, protein kinase A (PKA) phosphorylation of phospholamban (PLB) was decreased, whereas PKA phosphorylation of troponin I was increased, explaining the decoupling between calcium signaling and contractility. In silico modeling supported this relationship. Echocardiography measurements showed that Pkd2(+/-) mice have increased left ventricular ejection fraction after stimulation with isoproterenol (ISO), a ß-adrenergic receptor (ßAR) agonist. Blockers of ßAR-1 and ßAR-2 inhibited the ISO response in Pkd2(+/-) mice, suggesting that the dephosphorylated state of PLB is primarily by ßAR-2 signaling. Importantly, the Pkd2(+/-) mice were normotensive and had no evidence of renal cysts. Our results showed that decreased PC2 levels shifted the ßAR pathway balance and changed expression of calcium handling proteins, which resulted in altered cardiac contractility. We propose that PC2 levels in the heart may directly contribute to cardiac remodeling in patients with ADPKD in the absence of renal dysfunction.

Slowing of contractile kinetics by myosin-binding protein C can be explained by its cooperative binding to the thin filament.

Cardiac myosin binding protein C (cMyBP-C) is a thick filament-associated protein that participates in the regulation of muscle contraction. Simplified in vitro systems show that cMyBP-C binds not only to myosin, but also to the actin filament. The physiological significance of these separate binding interactions remains unclear, as does the question of whether either interaction is capable of explaining the behavior of intact muscle from which cMyBP-C has been removed. We have used a computational model to explore the characteristic effects of myosin-binding versus actin-binding by cMyBP-C. Simulations suggest that myosin-cMyBP-C interactions reduce peak force and Ca2 + sensitivity of the myofilaments, but have no appreciable effect on the rate of force redevelopment (ktr). In contrast, cMyBP-C binding to actin increases myofilament Ca2 + sensitivity and slows ktrat sub-maximal Ca2 + values. This slowing is due to cooperation between cMyBP-C 'crossbridges' and traditional myosin crossbridges as they bind to and activate the actin thin filament. We further observed that an overall recapitulation of skinned myocardial data from wild type and cMyBP-C knockout mice requires the interaction of cMyBP-C with of both of its binding targets in our model. The assumption of significant interactions with both partners was also sufficient to explain published effects of cMyBP-C ablation on twitch kinetics. These modeling results strongly support the view that both binding interactions play critical roles in the physiology of intact muscle. Furthermore, they suggest that the widely observed phenomenon of slowed force development in the presence of cMyBP-C may actually be a manifestation of cooperative binding of this protein to the thin filament.

Structural determinants of muscle thin filament cooperativity.

End-to-end connections between adjacent tropomyosin molecules along the muscle thin filament allow long-range conformational rearrangement of the multicomponent filament structure. This process is influenced by Ca(2+) and the troponin regulatory complexes, as well as by myosin crossbridge heads that bind to and activate the filament. Access of myosin crossbridges onto actin is gated by tropomyosin, and in the case of striated muscle filaments, troponin acts as a gatekeeper. The resulting tropomyosin-troponin-myosin on-off switching mechanism that controls muscle contractility is a complex cooperative and dynamic system with highly nonlinear behavior. Here, we review key information that leads us to view tropomyosin as central to the communication pathway that coordinates the multifaceted effectors that modulate and tune striated muscle contraction. We posit that an understanding of this communication pathway provides a framework for more in-depth mechanistic characterization of myopathy-associated mutational perturbations currently under investigation by many research groups.