Incomplete-penetrant hypertrophic cardiomyopathy MYH7 G256E mutation causes hypercontractility and elevated mitochondrial respiration.

Determining the pathogenicity of hypertrophic cardiomyopathy-associated mutations in the ß-myosin heavy chain (MYH7) can be challenging due to its variable penetrance and clinical severity. This study investigates the early pathogenic effects of the incomplete-penetrant MYH7 G256E mutation on myosin function that may trigger pathogenic adaptations and hypertrophy. We hypothesized that the G256E mutation would alter myosin biomechanical function, leading to changes in cellular functions. We developed a collaborative pipeline to characterize myosin function across protein, myofibril, cell, and tissue levels to determine the multiscale effects on structure-function of the contractile apparatus and its implications for gene regulation and metabolic state. The G256E mutation disrupts the transducer region of the S1 head and reduces the fraction of myosin in the folded-back state by 33%, resulting in more myosin heads available for contraction. Myofibrils from gene-edited MYH7WT/G256E human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) exhibited greater and faster tension development. This hypercontractile phenotype persisted in single-cell hiPSC-CMs and engineered heart tissues. We demonstrated consistent hypercontractile myosin function as a primary consequence of the MYH7 G256E mutation across scales, highlighting the pathogenicity of this gene variant. Single-cell transcriptomic and metabolic profiling demonstrated upregulated mitochondrial genes and increased mitochondrial respiration, indicating early bioenergetic alterations. This work highlights the benefit of our multiscale platform to systematically evaluate the pathogenicity of gene variants at the protein and contractile organelle level and their early consequences on cellular and tissue function. We believe this platform can help elucidate the genotype-phenotype relationships underlying other genetic cardiovascular diseases.

Structural OFF/ON transitions of myosin in relaxed porcine myocardium predict calcium-activated force.

Contraction in striated muscle is initiated by calcium binding to troponin complexes, but it is now understood that dynamic transition of myosin between resting, ordered OFF states on thick filaments and active, disordered ON states that can bind to thin filaments is critical in regulating muscle contractility. These structural OFF to ON transitions of myosin are widely assumed to correspond to transitions from the biochemically defined, energy-sparing, super-relaxed (SRX) state to the higher ATPase disordered-relaxed (DRX) state. Here we examined the effect of 2'-deoxy-ATP (dATP), a naturally occurring energy substrate for myosin, on the structural OFF to ON transitions of myosin motors in porcine cardiac muscle thick filaments. Small-angle X-ray diffraction revealed that titrating dATP in relaxation solutions progressively moves the myosin heads from ordered OFF states on the thick filament backbone to disordered ON states closer to thin filaments. Importantly, we found that the structural OFF to ON transitions are not equivalent to the biochemically defined SRX to DRX transitions and that the dATP-induced structural OFF to ON transitions of myosin motors in relaxed muscle are strongly correlated with submaximal force augmentation by dATP. These results indicate that structural OFF to ON transitions of myosin in relaxed muscle can predict the level of force attained in calcium-activated cardiac muscle. Computational modeling and stiffness measurements suggest a final step in the OFF to ON transition may involve a subset of DRX myosins that form weakly bound cross-bridges prior to becoming active force-producing cross-bridges.

α-Sheet secondary structure in amyloid ß-peptide drives aggregation and toxicity in Alzheimer's disease.

Alzheimer's disease (AD) is characterized by the deposition of ß-sheet-rich, insoluble amyloid ß-peptide (Aß) plaques; however, plaque burden is not correlated with cognitive impairment in AD patients; instead, it is correlated with the presence of toxic soluble oligomers. Here, we show, by a variety of different techniques, that these Aß oligomers adopt a nonstandard secondary structure, termed "α-sheet." These oligomers form in the lag phase of aggregation, when Aß-associated cytotoxicity peaks, en route to forming nontoxic ß-sheet fibrils. De novo-designed α-sheet peptides specifically and tightly bind the toxic oligomers over monomeric and fibrillar forms of Aß, leading to inhibition of aggregation in vitro and neurotoxicity in neuroblastoma cells. Based on this specific binding, a soluble oligomer-binding assay (SOBA) was developed as an indirect probe of α-sheet content. Combined SOBA and toxicity experiments demonstrate a strong correlation between α-sheet content and toxicity. The designed α-sheet peptides are also active in vivo where they inhibit Aß-induced paralysis in a transgenic Aß Caenorhabditis elegans model and specifically target and clear soluble, toxic oligomers in a transgenic APPsw mouse model. The α-sheet hypothesis has profound implications for further understanding the mechanism behind AD pathogenesis.

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.

Modulation of post-powerstroke dynamics in myosin II by 2'-deoxy-ADP.

Muscle myosins are molecular motors that hydrolyze ATP and generate force through coordinated interactions with actin filaments, known as cross-bridge cycling. During the cross-bridge cycle, functional sites in myosin 'sense' changes in interactions with actin filaments and the nucleotide binding region, resulting in allosteric transmission of information throughout the structure. We investigated whether the dynamics of the post-powerstroke state of the cross-bridge cycle are modulated in a nucleotide-dependent fashion. We compared molecular dynamics simulations of the myosin II motor domain (M) from Dictyostelium discoideum in the presence of ADP (M.ADP) versus 2'-deoxy-ADP bound myosin (M.dADP). We found that dADP was more flexible than ADP and the two nucleotides interacted with myosin in different ways. Replacement of ADP with dADP in the post-powerstroke state also altered the conformation of the actin binding region in myosin heads. Our results provide atomic level insights into allosteric communication networks in myosin that provide insight into the nucleotide-dependent dynamics of the cross-bridge cycle.

Edge Strand Dissociation and Conformational Changes in Transthyretin under Amyloidogenic Conditions.

During amyloidogenesis, proteins undergo conformational changes that allow them to aggregate and assemble into insoluble, fibrillar structures. Soluble oligomers that form during this process typically contain 2-24 monomeric subunits and are cytotoxic. Before the formation of these soluble oligomers, monomeric species first adopt aggregation-competent conformations. Knowledge of the structures of these intermediate states is invaluable to the development of molecular strategies to arrest pathological amyloid aggregation. However, the highly dynamic and interconverting nature of amyloidogenic species limits biophysical characterization of their structures during amyloidogenesis. Here, we use molecular dynamics simulations to probe conformations sampled by monomeric transthyretin under amyloidogenic conditions. We show that certain ß-strands in transthyretin tend to unfold and sample nonnative conformations and that the edge strands in one ß-sheet (the DAGH sheet) are particularly susceptible to conformational changes in the monomeric state. We also find that changes in the tertiary structure of transthyretin can be associated with disruptions to the secondary structure. We evaluated the conformations produced by molecular dynamics by calculating how well molecular-dynamics-derived structures reproduced NMR-derived interatomic distances. Finally, we leverage our computational results to produce experimentally testable hypotheses that may aid experimental explorations of pathological conformations of transthyretin.

Myosin dynamics during relaxation in mouse soleus muscle and modulation by 2'-deoxy-ATP.

KEY POINTS: Skeletal muscle relaxation has been primarily studied by assessing the kinetics of force decay. Little is known about the resultant dynamics of structural changes in myosin heads during relaxation. The naturally occurring nucleotide 2-deoxy-ATP (dATP) is a myosin activator that enhances cross-bridge binding and kinetics. X-ray diffraction data indicate that with elevated dATP, myosin heads were extended closer to actin in relaxed muscle and myosin heads return to an ordered, resting state after contraction more quickly. Molecular dynamics simulations of post-powerstroke myosin suggest that dATP induces structural changes in myosin heads that increase the surface area of the actin-binding regions promoting myosin interaction with actin, which could explain the observed delays in the onset of relaxation. This study of the dATP-induced changes in myosin may be instructive for determining the structural changes desired for other potential myosin-targeted molecular compounds to treat muscle diseases. ABSTRACT: Here we used time-resolved small-angle X-ray diffraction coupled with force measurements to study the structural changes in FVB mouse skeletal muscle sarcomeres during relaxation after tetanus contraction. To estimate the rate of myosin deactivation, we followed the rate of the intensity recovery of the first-order myosin layer line (MLL1) and restoration of the resting spacing of the third and sixth order of meridional reflection (SM3 and SM6 ) following tetanic contraction. A transgenic mouse model with elevated skeletal muscle 2-deoxy-ATP (dATP) was used to study how myosin activators may affect soleus muscle relaxation. X-ray diffraction evidence indicates that with elevated dATP, myosin heads were extended closer to actin in resting muscle. Following contraction, there is a slight but significant delay in the decay of force relative to WT muscle while the return of myosin heads to an ordered resting state was initially slower, then became more rapid than in WT muscle. Molecular dynamics simulations of post-powerstroke myosin suggest that dATP induces structural changes in myosin that increase the surface area of the actin-binding regions, promoting myosin interaction with actin. With dATP, myosin heads may remain in an activated state near the thin filaments following relaxation, accounting for the delay in force decay and the initial delay in recovery of resting head configuration, and this could facilitate subsequent contractions.

Drivers of α-Sheet Formation in Transthyretin under Amyloidogenic Conditions.

Amyloid diseases make up a set of fatal disorders in which proteins aggregate to form fibrils that deposit in tissues throughout the body. Amyloid-associated diseases are challenging to study because amyloid formation occurs on time scales that span several orders of magnitude and involve heterogeneous, interconverting protein conformations. The development of more effective technologies to diagnose and treat amyloid disease requires both a map of the conformations sampled during amyloidogenesis and an understanding of the molecular mechanisms that drive this process. In prior molecular dynamics simulations of amyloid proteins, we observed the formation of a nonstandard type of secondary structure, called α-sheet, that we proposed is associated with the pathogenic conformers in amyloid disease, the soluble oligomers. However, the detailed molecular interactions that drive the conversion to α-sheet remain elusive. Here we use molecular dynamics simulations to interrogate a critical event in transthyretin aggregation, the formation of aggregation-competent, monomeric species. We show that conformational changes in one of the two ß-sheets in transthyretin enable solvent molecules and polar side chains to form electrostatic interactions with main-chain peptide groups to facilitate and modulate conversion to α-sheet secondary structure. Our results shed light on the early conformational changes that drive transthyretin toward the α-sheet structure associated with toxicity. Delineation of the molecular events that lead to aggregation at atomic resolution can aid strategies to target the early, critical toxic soluble oligomers.

Identification of quinazoline based inhibitors of IRAK4 for the treatment of inflammation.

Interleukin-1 receptor associated kinase 4 (IRAK4) has been implicated in IL-1R and TLR based signaling. Therefore selective inhibition of the kinase activity of this protein represents an attractive target for the treatment of inflammatory diseases. Medicinal chemistry optimization of high throughput screening (HTS) hits with the help of structure based drug design led to the identification of orally-bioavailable quinazoline based IRAK4 inhibitors with excellent pharmacokinetic profile and kinase selectivity. These highly selective IRAK4 compounds show activity in vivo via oral dosing in a TLR7 driven model of inflammation.

Atomistic Simulations of Sarcomere Proteins.

Concerted atomic motions are requisite for sarcomere protein function and may become disrupted in HCM pathologies. Computational approaches such as molecular dynamics simulation can resolve such dynamics with unrivalled spatial and temporal resolution. This chapter describes methods to model structural and dynamical changes in biomolecules with HCM-associated perturbations.

An atomistic model of myosin interacting heads motif dynamics and their modification by 2'-deoxy-ADP.

The contraction of striated muscle is driven by cycling myosin motor proteins embedded within the thick filaments of sarcomeres. In addition to cross-bridge cycling with actin, these myosin proteins can enter an inactive, sequestered state in which the globular S1 heads rest along the thick filament surface and are unable to perform motor activities. Structurally, this state is called the interacting heads motif (IHM) and is a critical conformational state of myosin that regulates muscle contractility and energy expenditure. Structural perturbation of the sequestered state via missense mutations can pathologically disrupt the mechanical performance of muscle tissue. Thus, the IHM state has become a target for therapeutic intervention. An ATP analogue called 2'-deoxy-ATP (dATP) is a potent myosin activator which destabilizes the IHM. Here we use molecular dynamics simulations to study the molecular mechanisms by which dATP modifies the structure and dynamics of myosin in a sequestered state. Simulations with IHM containing ADP.Pi in both nucleotide binding pockets revealed residual dynamics in an otherwise 'inactive' and 'sequestered' state of a motor protein. Replacement of ADP.Pi by dADP.Pi triggered a series of structural changes that modify the protein-protein interface that stabilizes the sequestered state, and changes to this interface were accompanied by allosteric changes in remote regions of the protein complex. A comparative analysis of these dynamics predicted new structural sites that may affect IHM stability.

Myosin's powerstroke transitions define atomic scale movement of cardiac thin filament tropomyosin.

Dynamic interactions between the myosin motor head on thick filaments and the actin molecular track on thin filaments drive the myosin-crossbridge cycle that powers muscle contraction. The process is initiated by Ca2+ and the opening of troponin-tropomyosin-blocked myosin-binding sites on actin. The ensuing recruitment of myosin heads and their transformation from pre-powerstroke to post-powerstroke conformation on actin produce the force required for contraction. Cryo-EM-based atomic models confirm that during this process, tropomyosin occupies three different average positions on actin. Tropomyosin pivoting on actin away from a TnI-imposed myosin-blocking position accounts for part of the Ca2+ activation observed. However, the structure of tropomyosin on thin filaments that follows pre-powerstroke myosin binding and its translocation during myosin's pre-powerstroke to post-powerstroke transition remains unresolved. Here, we approach this transition computationally in silico. We used the myosin helix-loop-helix motif as an anchor to dock models of pre-powerstroke cardiac myosin to the cleft between neighboring actin subunits along cardiac thin filaments. We then performed targeted molecular dynamics simulations of the transition between pre- and post-powerstroke conformations on actin in the presence of cardiac troponin-tropomyosin. These simulations show Arg 369 and Glu 370 on the tip of myosin Loop-4 encountering identically charged residues on tropomyosin. The charge repulsion between residues causes tropomyosin translocation across actin, thus accounting for the final regulatory step in the activation of the thin filament, and, in turn, facilitating myosin movement along the filament. We suggest that during muscle activity, myosin-induced tropomyosin movement is likely to result in unencumbered myosin head interactions on actin at low-energy cost.

Hypertrophic cardiomyopathy mutations Y115H and E497D disrupt the folded-back state of human beta-cardiac myosin allosterically.

At the molecular level, clinical hypercontractility associated with many hypertrophic cardiomyopathy (HCM)-causing mutations in beta-cardiac myosin appears to be driven by their disruptive effect on the energy-conserving, folded-back, super relaxed (SRX) OFF-state of myosin. A pathological increase in force production results from release of heads from this OFF-state, which results in an increase in the number of heads free to interact with actin and produce force. Pathogenic mutations in myosin can conceivably disrupt the OFF-state by (1) directly affecting the intramolecular interfaces stabilizing the folded-back state, or (2) allosterically destabilizing the folded-back state via disruption of diverse conformational states of the myosin motor along its chemomechanical cycle. However, very little is understood about the mutations that fall in the latter group. Here, using recombinant human beta-cardiac myosin, we analysed the biomechanical properties of two such HCM-causing mutations, Y115H (in the transducer) and E497D (in the relay helix), neither of which falls in the regions that interact to stabilize the myosin folded-back state. We find these mutations have diverse effects on the contractility parameters of myosin, yet the primary hypercontractile change in both cases is the destabilization of the OFF-state of myosin and increased availability of active myosin heads for actin-binding. Experimental data and molecular dynamics simulations indicate that these mutations likely destabilize the pre-powerstroke state of myosin, the conformation the motor adopts in the inactive folded-back state. We propose that destabilization of the folded-back state of myosin, directly and/or allosterically, is the molecular basis of hypercontractility in HCM in a far greater number of pathogenic mutations than currently thought.

Molecular Dynamics Methods for Antibody Design.

Complex and coordinated dynamics are closely connected with protein functions, including the binding of antibodies to antigens. Knowledge of such dynamics could improve the design of antibodies. Molecular dynamics (MD) simulations provide a "computational microscope" that can resolve atomic motions and inform antibody design efforts.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects.

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7, TPM1, and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes such as MYH2, TPM2, and TNNI2 that encode parts of the skeletal muscle sarcomere cause muscle diseases affecting skeletal muscle, such as distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7) encoding sarcomeric proteins in which the same pathogenic variant affects skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain pathogenic variants that also cause cardiac abnormalities. We report five families with DA because of heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 (ACTC1). ACTC1 encodes a highly conserved actin that binds to myosin in cardiac and skeletal muscle. Pathogenic variants in ACTC1 have been found previously to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition because of variants in ACTC1 and suggests that some functions of ACTC1 are shared in cardiac and skeletal muscle.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects.

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7 , TPM1 , and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes MYH2 , TPM2 , and TNNI2 , that encode parts of the skeletal muscle sarcomere, cause muscle diseases affecting skeletal muscle, such as the distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7 ) encoding sarcomeric proteins in which the same pathogenic variant affects both skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain mutations that also cause cardiac abnormalities. We report five families with DA due to heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 ( ACTC1 ). ACTC1 encodes a highly conserved actin that binds to myosin in both cardiac and skeletal muscle. Mutations in ACTC1 have previously been found to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition due to mutations in ACTC1 and suggests that some functions of actin, alpha, cardiac muscle 1 are shared in cardiac and skeletal muscle.

Multi-scale models reveal hypertrophic cardiomyopathy MYH7 G256E mutation drives hypercontractility and elevated mitochondrial respiration.

Rationale: Over 200 mutations in the sarcomeric protein ß-myosin heavy chain (MYH7) have been linked to hypertrophic cardiomyopathy (HCM). However, different mutations in MYH7 lead to variable penetrance and clinical severity, and alter myosin function to varying degrees, making it difficult to determine genotype-phenotype relationships, especially when caused by rare gene variants such as the G256E mutation. Objective: This study aims to determine the effects of low penetrant MYH7 G256E mutation on myosin function. We hypothesize that the G256E mutation would alter myosin function, precipitating compensatory responses in cellular functions. Methods: We developed a collaborative pipeline to characterize myosin function at multiple scales (protein to myofibril to cell to tissue). We also used our previously published data on other mutations to compare the degree to which myosin function was altered. Results: At the protein level, the G256E mutation disrupts the transducer region of the S1 head and reduces the fraction of myosin in the folded-back state by 50.9%, suggesting more myosins available for contraction. Myofibrils isolated from hiPSC-CMs CRISPR-edited with G256E (MYH7 WT/G256E ) generated greater tension, had faster tension development and slower early phase relaxation, suggesting altered myosin-actin crossbridge cycling kinetics. This hypercontractile phenotype persisted in single-cell hiPSC-CMs and engineered heart tissues. Single-cell transcriptomic and metabolic profiling demonstrated upregulation of mitochondrial genes and increased mitochondrial respiration, suggesting altered bioenergetics as an early feature of HCM. Conclusions: MYH7 G256E mutation causes structural instability in the transducer region, leading to hypercontractility across scales, perhaps from increased myosin recruitment and altered crossbridge cycling. Hypercontractile function of the mutant myosin was accompanied by increased mitochondrial respiration, while cellular hypertrophy was modest in the physiological stiffness environment. We believe that this multi-scale platform will be useful to elucidate genotype-phenotype relationships underlying other genetic cardiovascular diseases.

Structure-dependent impairment of intracellular apolipoprotein E4 trafficking and its detrimental effects are rescued by small-molecule structure correctors.

Apolipoprotein (apo) E4 is the major genetic risk factor for Alzheimer disease (AD) and likely contributes to neuropathology through various pathways. Here we report that the intracellular trafficking of apoE4 is impaired in Neuro-2a cells and primary neurons, as shown by measuring fluorescence recovery after photobleaching. In Neuro-2a cells, more apoE4 than apoE3 molecules remained immobilized in the endoplasmic reticulum (ER) and the Golgi apparatus, and the lateral motility of apoE4 was significantly lower in the Golgi apparatus (but not in the ER) than that of apoE3. Likewise, the immobile fraction was larger, and the lateral motility was lower for apoE4 than apoE3 in mouse primary hippocampal neurons. ApoE4 with the R61T mutation, which abolishes apoE4 domain interaction, was less immobilized, and its lateral motility was comparable with that of apoE3. The trafficking impairment of apoE4 was also rescued by disrupting domain interaction with the small-molecule structure correctors GIND25 and PH002. PH002 also rescued apoE4-induced impairments of neurite outgrowth in Neuro-2a cells and dendritic spine development in primary neurons. ApoE4 did not affect trafficking of amyloid precursor protein, another AD-related protein, through the secretory pathway. Thus, domain interaction renders more newly synthesized apoE4 molecules immobile and slows their trafficking along the secretory pathway. Correcting the pathological structure of apoE4 by disrupting domain interaction is a potential therapeutic approach to treat or prevent AD related to apoE4.

Conformational sampling of CMT-2D associated GlyRS mutations.

During protein synthesis, aminoacyl-tRNA synthetases covalently link amino acids with their cognate tRNAs. Amino acid mutations in glycyl-tRNA synthetase can disrupt protein synthesis and lead to a neurological disorder known as Charcot-Marie-Tooth disease type 2D (CMT-2D). Several studies employing diverse techniques have identified potential disease mechanisms at the molecular level. The majority of CMT-2D mutations in glycyl-tRNA are found within its dimer interface. However, no atomic structures bearing these mutations have been solved. Consequently, the specific disease-causing structural changes that occur in glycyl-tRNA synthetase have not been definitively established. Here we use molecular dynamics simulations to probe conformational changes in glycyl-tRNA synthetase caused by one mutation within the dimer interface: G240R. Our results show that the mutation alters the number of native interactions at the dimer interface and also leads to altered dynamics of two regions of glycyl-tRNA synthetase associated with tRNA binding. Additionally, we use our simulations to make predictions about the effects of other clinically reported CMT-2D mutations. Our results identify a region of the glycyl-tRNA synthetase structure that may be disrupted in a large number of CMT-2D mutations. Structural changes in this region may be a common molecular mechanism in glycyl-tRNA synthetase CMT-2D pathologies. Statement of significance: In this study, we use molecular dynamics simulations to elucidate structural conformations accessible to glycyl-tRNA synthetase (GlyRS), an enzyme that ligates cytosolic glycine with tRNA-Gly. This protein contains multiple flexible regions with dynamics that elude in vitro structural characterization. Our computational approach provides unparalleled atomistic details of structural changes in GlyRS that contribute to its role in protein synthesis. A number of mutations in GlyRS are associated with a peripheral nerve disorder, Charcot-Marie-Tooth disease type 2D (CMT-2D). Mutation-induced structural and dynamic changes in GlyRS have similarity that elude in vitro structural characterization. Our simulations provide insights into disease mechanisms for one such mutation: G240R. Additionally, we leverage our computational data to identify regions of GlyRS critical to its function and to predict the effects of other disease-associated mutations. These results open up new directions for research into the molecular characterization of GlyRS and into hypothesis-driven studies of CMT-2D disease mechanisms.

Post-translational modification patterns on ß-myosin heavy chain are altered in ischemic and nonischemic human hearts.

Phosphorylation and acetylation of sarcomeric proteins are important for fine-tuning myocardial contractility. Here, we used bottom-up proteomics and label-free quantification to identify novel post-translational modifications (PTMs) on ß-myosin heavy chain (ß-MHC) in normal and failing human heart tissues. We report six acetylated lysines and two phosphorylated residues: K34-Ac, K58-Ac, S210-P, K213-Ac, T215-P, K429-Ac, K951-Ac, and K1195-Ac. K951-Ac was significantly reduced in both ischemic and nonischemic failing hearts compared to nondiseased hearts. Molecular dynamics (MD) simulations show that K951-Ac may impact stability of thick filament tail interactions and ultimately myosin head positioning. K58-Ac altered the solvent-exposed SH3 domain surface - known for protein-protein interactions - but did not appreciably change motor domain conformation or dynamics under conditions studied. Together, K213-Ac/T215-P altered loop 1's structure and dynamics - known to regulate ADP-release, ATPase activity, and sliding velocity. Our study suggests that ß-MHC acetylation levels may be influenced more by the PTM location than the type of heart disease since less protected acetylation sites are reduced in both heart failure groups. Additionally, these PTMs have potential to modulate interactions between ß-MHC and other regulatory sarcomeric proteins, ADP-release rate of myosin, flexibility of the S2 region, and cardiac myofilament contractility in normal and failing hearts.