Congenital myopathy with tremor (MYOTREM) is a recently described disorder characterized by mild myopathy and a postural and intention tremor present since early infancy. MYOTREM is associated with pathogenic variants in MYBPC1 which encodes slow myosin-binding protein C, a sarcomere protein with regulatory and structural roles. Here, we describe a family with three generations of variably affected members exhibiting a novel variant in MYBPC1 (c.656 T > C, p.Leu219Pro). Among the unique features of affected family members is the persistence of tremor in sleep. We also present the first muscle magnetic resonance images for this disorder, and report muscle atrophy and fatty infiltration.

Cardiac myosin binding protein-C (cMyBP-C) is a thick filament-associated regulatory protein frequently found mutated in patients suffering from hypertrophic cardiomyopathy (HCM). Recent in vitro experiments have highlighted the functional significance of its N-terminal region (NcMyBP-C) for heart muscle contraction, reporting regulatory interactions with both thick and thin filaments. To better understand the interactions of cMyBP-C in its native sarcomere environment, in situ Foerster resonance energy transfer-fluorescence lifetime imaging (FRET-FLIM) assays were developed to determine the spatial relationship between the NcMyBP-C and the thick and thin filaments in isolated neonatal rat cardiomyocytes (NRCs). In vitro studies showed that ligation of genetically encoded fluorophores to NcMyBP-C had no or little effect on its binding to thick and thin filament proteins. Using this assay, FRET between mTFP conjugated to NcMyBP-C and Phalloidin-iFluor 514 labeling the actin filaments in NRCs was detected by time-domain FLIM. The measured FRET efficiencies were intermediate between those observed when the donor was attached to the cardiac myosin regulatory light chain in the thick filaments and troponin T in the thin filaments. These results are consistent with the coexistence of multiple conformations of cMyBP-C, some with their N-terminal domains binding to the thin filament and others binding to the thick filament, supporting the hypothesis that the dynamic interchange between these conformations mediates interfilament signaling in the regulation of contractility. Moreover, stimulation of NRCs with β-adrenergic agonists reduces FRET between NcMyBP-C and actin-bound Phalloidin, suggesting that cMyBP-C phosphorylation reduces its interaction with the thin filament.

The pathogenesis of MYBPC3-associated hypertrophic cardiomyopathy (HCM) is still unresolved. In our HCM patient cohort, a large and well-characterized population carrying the MYBPC3:c772G>A variant (p.Glu258Lys, E258K) provides the unique opportunity to study the basic mechanisms of MYBPC3-HCM with a comprehensive translational approach.
We collected clinical and genetic data from 93 HCM patients carrying the MYBPC3:c772G>A variant. Functional perturbations were investigated using different biophysical techniques in left ventricular samples from 4 patients who underwent myectomy for refractory outflow obstruction, compared with samples from non-failing non-hypertrophic surgical patients and healthy donors. Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes and engineered heart tissues (EHTs) were also investigated.
Haplotype analysis revealed MYBPC3:c772G>A as a founder mutation in Tuscany. In ventricular myocardium, the mutation leads to reduced cMyBP-C (cardiac myosin binding protein-C) expression, supporting haploinsufficiency as the main primary disease mechanism. Mechanical studies in single myofibrils and permeabilized muscle strips highlighted faster cross-bridge cycling, and higher energy cost of tension generation. A novel approach based on tissue clearing and advanced optical microscopy supported the idea that the sarcomere energetics dysfunction is intrinsically related with the reduction in cMyBP-C. Studies in single cardiomyocytes (native and hiPSC-derived), intact trabeculae and hiPSC-EHTs revealed prolonged action potentials, slower Ca2+ transients and preserved twitch duration, suggesting that the slower excitation-contraction coupling counterbalanced the faster sarcomere kinetics. This conclusion was strengthened by in silico simulations.
HCM-related MYBPC3:c772G>A mutation invariably impairs sarcomere energetics and cross-bridge cycling. Compensatory electrophysiological changes (eg, reduced potassium channel expression) appear to preserve twitch contraction parameters, but may expose patients to greater arrhythmic propensity and disease progression. Therapeutic approaches correcting the primary sarcomeric defects may prevent secondary cardiomyocyte remodeling.

Stretch activation is defined as a delayed increase in force after rapid stretches. Although there is considerable evidence for stretch activation in isolated cardiac myofibrillar preparations, few studies have measured mechanisms of stretch activation in mammalian skeletal muscle fibers. We measured stretch activation following rapid step stretches [∼1%-4% sarcomere length (SL)] during submaximal Ca2+ activations of rat permeabilized slow-twitch skeletal muscle fibers before and after protein kinase A (PKA), which phosphorylates slow myosin binding protein-C. PKA significantly increased stretch activation during low (∼25%) Ca2+ activation and accelerated rates of delayed force development (kef) during both low and half-maximal Ca2+ activation. Following the step stretches and subsequent force development, fibers were rapidly shortened to original sarcomere length, which often elicited a shortening-induced transient force overshoot. After PKA, step shortening-induced transient force overshoot increased ∼10-fold following an ∼4% SL shortening during low Ca2+ activation levels. kdf following step shortening also increased after PKA during low and half-maximal Ca2+ activations. We next investigated thin filament regulation of stretch activation. We tested the interplay between cardiac troponin I (cTnI) phosphorylation at the canonical PKA and novel tyrosine kinase sites on stretch activation. Native slow-skeletal Tn complexes were exchanged with recombinant human cTn complex with different human cTnI N-terminal pseudo-phosphorylation molecules: 1) nonphosphorylated wild type (WT), 2) the canonical S22/23D PKA sites, 3) the tyrosine kinase Y26E site, and 4) the combinatorial S22/23D + Y26E cTnI. All three pseudo-phosphorylated cTnIs elicited greater stretch activation than WT. Following stretch activation, a new, elevated stretch-induced steady-state force was reached with pseudo-phosphorylated cTnI. Combinatorial S22/23D + Y26E pseudo-phosphorylated cTnI increased kdf. These results suggest that slow-skeletal myosin binding protein-C (sMyBP-C) phosphorylation modulates stretch activation by a combination of cross-bridge recruitment and faster cycling kinetics, whereas cTnI phosphorylation regulates stretch activation by both redundant and synergistic mechanisms; and, taken together, these sarcomere phosphoproteins offer precision targets for enhanced contractility.

Background: A 25-base pair (25bp) intronic deletion in the MYBPC3 gene enriched in South Asians (SAs) is a risk allele for late-onset left ventricular (LV) dysfunction, hypertrophy, and heart failure (HF) with several forms of cardiomyopathy. However, the effect of this variant on exercise parameters has not been evaluated. Methods: As a pilot study, 10 asymptomatic SA carriers of the MYBPC3 Δ25bp variant (52.9 ± 2.14 years) and 10 age- and gender-matched non-carriers (NCs) (50.1 ± 2.7 years) were evaluated at baseline and under exercise stress conditions using bicycle exercise echocardiography and continuous cardiac monitoring. Results: Baseline echocardiography parameters were not different between the two groups. However, in response to exercise stress, the carriers of Δ25bp had significantly higher LV ejection fraction (%) (CI: 4.57 ± 1.93; p < 0.0001), LV outflow tract peak velocity (m/s) (CI: 0.19 ± 0.07; p < 0.0001), and higher aortic valve (AV) peak velocity (m/s) (CI: 0.103 ± 0.08; p = 0.01) in comparison to NCs, and E/A ratio, a marker of diastolic compliance, was significantly lower in Δ25bp carriers (CI: 0.107 ± 0.102; p = 0.038). Interestingly, LV end-diastolic diameter (LVIDdia) was augmented in NCs in response to stress, while it did not increase in Δ25bp carriers (CI: 0.239 ± 0.125; p = 0.0002). Further, stress-induced right ventricular systolic excursion velocity s\' (m/s), as a marker of right ventricle function, increased similarly in both groups, but tricuspid annular plane systolic excursion increased more in carriers (slope: 0.008; p = 0.0001), suggesting right ventricle functional differences between the two groups. Conclusions: These data support that MYBPC3 Δ25bp is associated with LV hypercontraction under stress conditions with evidence of diastolic impairment.

Cardiac myosin binding protein-C (cMyBP-C) is a thick filament protein that influences sarcomere stiffness and modulates cardiac contraction-relaxation through its phosphorylation. Phosphorylation of cMyBP-C and ablation of cMyBP-C have been shown to increase the rate of MgADP release in the acto-myosin cross-bridge cycle in the intact sarcomere. The influence of cMyBP-C on Pi-dependent myosin kinetics has not yet been examined. We investigated the effect of cMyBP-C, and its phosphorylation, on myosin kinetics in demembranated papillary muscle strips bearing the β-cardiac myosin isoform from nontransgenic and homozygous transgenic mice lacking cMyBP-C. We used quick stretch and stochastic length-perturbation analysis to characterize rates of myosin detachment and force development over 0-12 mM Pi and at maximal (pCa 4.8) and near-half maximal (pCa 5.75) Ca2+ activation. Protein kinase A (PKA) treatment was applied to half the strips to probe the effect of cMyBP-C phosphorylation on Pi sensitivity of myosin kinetics. Increasing Pi increased myosin cross-bridge detachment rate similarly for muscles with and without cMyBP-C, although these rates were higher in muscle without cMyBP-C. Treating myocardial strips with PKA accelerated detachment rate when cMyBP-C was present over all Pi, but not when cMyBP-C was absent. The rate of force development increased with Pi in all muscles. However, Pi sensitivity of the rate force development was reduced when cMyBP-C was present versus absent, suggesting that cMyBP-C inhibits Pi-dependent reversal of the power stroke or stabilizes cross-bridge attachment to enhance the probability of completing the power stroke. These results support a functional role for cMyBP-C in slowing myosin detachment rate, possibly through a direct interaction with myosin or by altering strain-dependent myosin detachment via cMyBP-C-dependent stiffness of the thick filament and myofilament lattice. PKA treatment reduces the role for cMyBP-C to slow myosin detachment and thus effectively accelerates β-myosin detachment in the intact myofilament lattice.NEW & NOTEWORTHY Length perturbation analysis was used to demonstrate that β-cardiac myosin characteristic rates of detachment and recruitment in the intact myofilament lattice are accelerated by Pi, phosphorylation of cMyBP-C, and the absence of cMyBP-C. The results suggest that cMyBP-C normally slows myosin detachment, including Pi-dependent detachment, and that this inhibition is released with phosphorylation or absence of cMyBP-C.

The heart\'s response to varying demands of the body is regulated by signaling pathways that activate protein kinases which phosphorylate sarcomeric proteins. Although phosphorylation of cardiac myosin binding protein-C (cMyBP-C) has been recognized as a key regulator of myocardial contractility, little is known about its mechanism of action. Here, we used protein kinase A (PKA) and Cε (PKCε), as well as ribosomal S6 kinase II (RSK2), which have different specificities for cMyBP-C\'s multiple phosphorylation sites, to show that individual sites are not independent, and that phosphorylation of cMyBP-C is controlled by positive and negative regulatory coupling between those sites. PKA phosphorylation of cMyBP-C\'s N terminus on 3 conserved serine residues is hierarchical and antagonizes phosphorylation by PKCε, and vice versa. In contrast, RSK2 phosphorylation of cMyBP-C accelerates PKA phosphorylation. We used cMyBP-C\'s regulatory N-terminal domains in defined phosphorylation states for protein-protein interaction studies with isolated cardiac native thin filaments and the S2 domain of cardiac myosin to show that site-specific phosphorylation of this region of cMyBP-C controls its interaction with both the actin-containing thin and myosin-containing thick filaments. We also used fluorescence probes on the myosin-associated regulatory light chain in the thick filaments and on troponin C in the thin filaments to monitor structural changes in the myofilaments of intact heart muscle cells associated with activation of myocardial contraction by the N-terminal region of cMyBP-C in its different phosphorylation states. Our results suggest that cMyBP-C acts as a sarcomeric integrator of multiple signaling pathways that determines downstream physiological function.

Encoding the slow skeletal muscle isoform of myosin binding protein-C, MYBPC1 is associated with autosomal dominant and recessive forms of arthrogryposis. The authors describe a novel association for MYBPC1 in four patients from three independent families with skeletal muscle weakness, myogenic tremors, and hypotonia with gradual clinical improvement. The patients carried one of two de novo heterozygous variants in MYBPC1, with the p.Leu263Arg variant seen in three individuals and the p.Leu259Pro variant in one individual. Both variants are absent from controls, well conserved across vertebrate species, predicted to be damaging, and located in the M-motif. Protein modeling studies suggested that the p.Leu263Arg variant affects the stability of the M-motif, whereas the p.Leu259Pro variant alters its structure. In vitro biochemical and kinetic studies demonstrated that the p.Leu263Arg variant results in decreased binding of the M-motif to myosin, which likely impairs the formation of actomyosin cross-bridges during muscle contraction. Collectively, our data substantiate that damaging variants in MYBPC1 are associated with a new form of an early-onset myopathy with tremor, which is a defining and consistent characteristic in all affected individuals, with no contractures. Recognition of this expanded myopathic phenotype can enable identification of individuals with MYBPC1 variants without arthrogryposis.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in combination with CRISPR/Cas9 genome editing provide unparalleled opportunities to study cardiac biology and disease. However, sarcomeres, the fundamental units of myocyte contraction, are immature and nonlinear in hiPSC-CMs, which technically challenge accurate functional interrogation of contractile parameters in beating cells. Furthermore, existing analysis methods are relatively low-throughput, indirectly assess contractility, or only assess well-aligned sarcomeres found in mature cardiac tissues.
We aimed to develop an analysis platform that directly, rapidly, and automatically tracks sarcomeres in beating cardiomyocytes. The platform should assess sarcomere content, contraction and relaxation parameters, and beat rate.
We developed SarcTrack, a MatLab software that monitors fluorescently tagged sarcomeres in hiPSC-CMs. The algorithm determines sarcomere content, sarcomere length, and returns rates of sarcomere contraction and relaxation. By rapid measurement of hundreds of sarcomeres in each hiPSC-CM, SarcTrack provides large data sets for robust statistical analyses of multiple contractile parameters. We validated SarcTrack by analyzing drug-treated hiPSC-CMs, confirming the contractility effects of compounds that directly activate (CK-1827452) or inhibit (MYK-461) myosin molecules or indirectly alter contractility (verapamil and propranolol). SarcTrack analysis of hiPSC-CMs carrying a heterozygous truncation variant in the myosin-binding protein C ( MYBPC3) gene, which causes hypertrophic cardiomyopathy, recapitulated seminal disease phenotypes including cardiac hypercontractility and diminished relaxation, abnormalities that normalized with MYK-461 treatment.
SarcTrack provides a direct and efficient method to quantitatively assess sarcomere function. By improving existing contractility analysis methods and overcoming technical challenges associated with functional evaluation of hiPSC-CMs, SarcTrack enhances translational prospects for sarcomere-regulating therapeutics and accelerates interrogation of human cardiac genetic variants.

In the last decade, improvements in electron microscopy and image processing have permitted significantly higher resolutions to be achieved (sometimes <1 nm) when studying isolated actin and myosin filaments. In the case of actin filaments the changing structure when troponin binds calcium ions can be followed using electron microscopy and single particle analysis to reveal what happens on each of the seven non-equivalent pseudo-repeats of the tropomyosin α-helical coiled-coil. In the case of the known family of myosin filaments not only are the myosin head arrangements under relaxing conditions being defined, but the latest analysis, also using single particle methods, is starting to reveal the way that the α-helical coiled-coil myosin rods are packed to give the filament backbones.