Multiscale modeling shows how 2'-deoxy-ATP rescues ventricular function in heart failure.

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Comprehensive Strategies to Bicyclic Prolines: Applications in the Synthesis of Potent Arginase Inhibitors.

Comprehensive synthetic strategies afforded a diverse set of structurally unique bicyclic proline-containing arginase inhibitors with a high degree of three-dimensionality. The analogs that favored the Cγ-exo conformation of the proline improved the arginase potency over the initial lead. The novel synthetic strategies reported here not only enable access to previously unknown stereochemically complex proline derivatives but also provide a foundation for the future synthesis of bicyclic proline analogs, which incorporate inherent three-dimensional character into building blocks, medicine, and catalysts and could have a profound impact on the conformation of proline-containing peptides and macrocycles.

Structure-Based Discovery of Proline-Derived Arginase Inhibitors with Improved Oral Bioavailability for Immuno-Oncology.

Recent data suggest that the inhibition of arginase (ARG) has therapeutic potential for the treatment of a number of indications ranging from pulmonary and vascular disease to cancer. Thus, high demand exists for selective small molecule ARG inhibitors with favorable druglike properties and good oral bioavailability. In light of the significant challenges associated with the unique physicochemical properties of previously disclosed ARG inhibitors, we use structure-based drug design combined with a focused optimization strategy to discover a class of boronic acids featuring a privileged proline scaffold with superior potency and oral bioavailability. These compounds, exemplified by inhibitors 4a, 18, and 27, demonstrated a favorable overall profile, and 4a was well tolerated following multiple days of dosing at concentrations that exceed those required for serum arginase inhibition and concomitant arginine elevation in a syngeneic mouse carcinoma model.

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.

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues.

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long-range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.

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.