Abnormal myosin post-translational modifications and ATP turnover time associated with human congenital myopathy-related RYR1 mutations.

AIM: Conditions related to mutations in the gene encoding the skeletal muscle ryanodine receptor 1 (RYR1) are genetic muscle disorders and include congenital myopathies with permanent weakness, as well as episodic phenotypes such as rhabdomyolysis/myalgia. Although RYR1 dysfunction is the primary mechanism in RYR1-related disorders, other downstream pathogenic events are less well understood and may include a secondary remodeling of major contractile proteins. Hence, in the present study, we aimed to investigate whether congenital myopathy-related RYR1 mutations alter the regulation of the most abundant contractile protein, myosin. METHODS: We used skeletal muscle tissues from five patients with RYR1-related congenital myopathy and compared those with five controls and five patients with RYR1-related rhabdomyolysis/myalgia. We then defined post-translational modifications on myosin heavy chains (MyHCs) using LC/MS. In parallel, we determined myosin relaxed states using Mant-ATP chase experiments and performed molecular dynamics (MD) simulations. RESULTS: LC/MS revealed two additional phosphorylations (Thr1309-P and Ser1362-P) and one acetylation (Lys1410-Ac) on the ß/slow MyHC of patients with congenital myopathy. This method also identified six acetylations that were lacking on MyHC type IIa of these patients (Lys35-Ac, Lys663-Ac, Lys763-Ac, Lys1171-Ac, Lys1360-Ac, and Lys1733-Ac). MD simulations suggest that modifying myosin Ser1362 impacts the protein structure and dynamics. Finally, Mant-ATP chase experiments showed a faster ATP turnover time of myosin heads in the disordered-relaxed conformation. CONCLUSIONS: Altogether, our results suggest that RYR1 mutations have secondary negative consequences on myosin structure and function, likely contributing to the congenital myopathic phenotype.

Binding pocket dynamics along the recovery stroke of human ß-cardiac myosin.

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human ß-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.