Investigating the effects of tropomyosin mutations on its flexibility and interactions with filamentous actin using molecular dynamics simulation

Novel Mutation Glu98Lys in Cardiac Tropomyosin Alters Its Structure and Impairs Myocardial Relaxation

We characterized a novel genetic variant c.292G > A (p.E98K) in the TPM1 gene encoding cardiac tropomyosin 1.1 isoform (Tpm1.1), found in a proband with a phenotype of complex cardiomyopathy with conduction dysfunction and slow progressive neuromuscular involvement. To understand the molecular mechanism by which this mutation impairs cardiac function, we produced recombinant Tpm1.1 carrying an E98K substitution and studied how this substitution affects the structure of the Tpm1.1 molecule and its functional properties. The results showed that the E98K substitution in the N-terminal part of the Tpm molecule significantly destabilizes the C-terminal part of Tpm, thus indicating a long-distance destabilizing effect of the substitution on the Tpm coiled-coil structure. The E98K substitution did not noticeably affect Tpm's affinity for F-actin but significantly impaired Tpm's regulatory properties. It increased the Ca2+ sensitivity of the sliding velocity of regulated thin filaments over cardiac myosin in an in vitro motility assay and caused an incomplete block of the thin filament sliding at low Ca2+ concentrations. The incomplete motility block in the absence of Ca2+ can be explained by the loosening of the Tpm interaction with troponin I (TnI), thus increasing Tpm mobility on the surface of an actin filament that partially unlocks the myosin binding sites. This hypothesis is supported by the molecular dynamics (MD) simulation that showed that the E98 Tpm residue is involved in hydrogen bonding with the C-terminal part of TnI. Thus, the results allowed us to explain the mechanism by which the E98K Tpm mutation impairs sarcomeric function and myocardial relaxation.

Molecular Dynamics Assessment of Mechanical Properties of the Thin Filaments in Cardiac Muscle

Contraction of cardiac muscle is regulated by Ca2+ ions via regulatory proteins, troponin (Tn), and tropomyosin (Tpm) associated with the thin (actin) filaments in myocardial sarcomeres. The binding of Ca2+ to a Tn subunit causes mechanical and structural changes in the multiprotein regulatory complex. Recent cryo-electron microscopy (cryo-EM) models of the complex allow one to study the dynamic and mechanical properties of the complex using molecular dynamics (MD). Here we describe two refined models of the thin filament in the calcium-free state that include protein fragments unresolved by cryo-EM and reconstructed using structure prediction software. The parameters of the actin helix and the bending, longitudinal, and torsional stiffness of the filaments estimated from the MD simulations performed with these models were close to those found experimentally. However, problems revealed from the MD simulation suggest that the models require further refinement by improving the protein-protein interaction in some regions of the complex. The use of relatively long refined models of the regulatory complex of the thin filament allows one to perform MD simulation of the molecular mechanism of Ca2+ regulation of contraction without additional constraints and study the effects of cardiomyopathy-associated mutation of the thin filament proteins of cardiac muscle.

Structural analysis of M1AP variants associated with severely impaired spermatogenesis causing male infertility

Background

Impaired meiosis can result in absence of sperm in the seminal fluid. This condition, namely non-obstructive azoospermia (NOA), is one of the reasons of male infertility. Despite the low number of studies on meiosis 1-associated protein (M1AP) in the literature, M1AP is known to be crucial for spermatogenesis. Recently, seven variants (five missense, one frameshift, one splice-site) have been reported in the M1AP gene as associated with NOA, cryptozoospermia and oligozoospermia in two separate studies. However, all missense variants were evaluated as variant of uncertain significance by these studies. Therefore, we aimed to analyze their structural impacts on the M1AP protein that could lead to NOA.

Methods

We firstly performed an evolutionary conservation analysis for the variant positions. Afterwards, a comprehensive molecular modelling study was performed for the M1AP structure. By utilizing this model, protein dynamics were sampled for the wild-type and variants by performing molecular dynamics (MD) simulations.

Results

All variant positions are highly conserved, indicating that they are potentially important for function. In MD simulations, none of the variants led to a general misfolding or loss of stability in the protein structure, but they did cause severe modifications in the conformational dynamics of M1AP, particularly through changes in local interactions affecting flexibility, hinge and secondary structure.

Conclusions

Due to critical perturbations in protein dynamics, we propose that these variants may cause NOA by affecting important interactions regulating meiosis, particularly in wild-type M1AP deficiency since the variants are reported to be homozygous or bi-allelic in the infertile individuals. Our results provided reasonable insights about the M1AP structure and the effects of the variants to the structure and dynamics, which should be further investigated by experimental studies to validate.

Loss of crossbridge inhibition drives pathological cardiac hypertrophy in patients harboring the TPM1 E192K mutation

Hypertrophic cardiomyopathy (HCM) is an inherited disorder caused primarily by mutations to thick and thinfilament proteins. Although thin filament mutations are less prevalent than their oft-studied thick filament counterparts, they are frequently associated with severe patient phenotypes and can offer important insight into fundamental disease mechanisms. We have performed a detailed study of tropomyosin (TPM1) E192K, a variant of uncertain significance associated with HCM. Molecular dynamics revealed that E192K results in a more flexible TPM1 molecule, which could affect its ability to regulate crossbridges. In vitro motility assays of regulated actin filaments containing TPM1 E192K showed an overall loss of Ca2+ sensitivity. To understand these effects, we used multiscale computational models that suggested a subtle phenotype in which E192K leads to an inability to completely inhibit actin-myosin crossbridge activity at low Ca2+. To assess the physiological impact of the mutation, we generated patient-derived engineered heart tissues expressing E192K. These tissues showed disease features similar to those of the patients, including cellular hypertrophy, hypercontractility, and diastolic dysfunction. We hypothesized that excess residual crossbridge activity could be triggering cellular hypertrophy, even if the overall Ca2+ sensitivity was reduced by E192K. To test this hypothesis, the cardiac myosin-specific inhibitor mavacamten was applied to patient-derived engineered heart tissues for 4 d followed by 24 h of washout. Chronic mavacamten treatment abolished contractile differences between control and TPM1 E192K engineered heart tissues and reversed hypertrophy in cardiomyocytes. These results suggest that the TPM1 E192K mutation triggers cardiomyocyte hypertrophy by permitting excess residual crossbridge activity. These studies also provide direct evidence that myosin inhibition by mavacamten can counteract the hypertrophic effects of mutant tropomyosin.

Cardiomyopathic troponin mutations predominantly occur at its interface with actin and tropomyosin

Reversible Ca2+ binding to troponin is the primary on-off switch of the contractile apparatus of striated muscles, including the heart. Dominant missense mutations in human cardiac troponin genes are among the causes of hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy. Structural understanding of troponin action has recently advanced considerably via electron microscopy and molecular dynamics studies of the thin filament. As a result, it is now possible to examine cardiomyopathy-inducing troponin mutations in thin-filament structural context, and from that to seek new insight into pathogenesis and into the troponin regulatory mechanism. We compiled from consortium reports a representative set of troponin mutation sites whose pathogenicity was determined using standardized clinical genetics criteria. Another set of sites, apparently tolerant of amino acid substitutions, was compiled from the gnomAD v2 database. Pathogenic substitutions occurred predominantly in the areas of troponin that contact actin or tropomyosin, including, but not limited to, two regions of newly proposed structure and long-known implication in cardiomyopathy: the C-terminal third of troponin I and a part of the troponin T N terminus. The pathogenic mutations were located in troponin regions that prevent contraction under low Ca2+ concentration conditions. These regions contribute to Ca2+-regulated steric hindrance of myosin by the combined effects of troponin and tropomyosin. Loss-of-function mutations within these parts of troponin result in loss of inhibition, consistent with the hypercontractile phenotype characteristic of HCM. Notably, pathogenic mutations are absent in our dataset from the Ca2+-binding, activation-producing troponin C (TnC) N-lobe, which controls contraction by a multi-faceted mechanism. Apparently benign mutations are also diminished in the TnC N-lobe, suggesting mutations are poorly tolerated in that critical domain.

Impact of A134 and E218 Amino Acid Residues of Tropomyosin on Its Flexibility and Function

Tropomyosin (Tpm) is one of the major actin-binding proteins that play a crucial role in the regulation of muscle contraction. The flexibility of the Tpm molecule is believed to be vital for its functioning, although its role and significance are under discussion. We choose two sites of the Tpm molecule that presumably have high flexibility and stabilized them with the A134L or E218L substitutions. Applying differential scanning calorimetry (DSC), molecular dynamics (MD), co-sedimentation, trypsin digestion, and in vitro motility assay, we characterized the properties of Tpm molecules with these substitutions. The A134L mutation prevented proteolysis of Tpm molecule by trypsin, and both substitutions increased the thermal stability of Tpm and its bending stiffness estimated from MD simulation. None of these mutations affected the primary binding of Tpm to F-actin; still, both of them increased the thermal stability of the actin-Tpm complex and maximal sliding velocity of regulated thin filaments in vitro at a saturating Ca²⁺ concentration. However, the mutations differently affected the Ca²⁺ sensitivity of the sliding velocity and pulling force produced by myosin heads. The data suggest that both regions of instability are essential for correct regulation and fine-tuning of Ca²⁺-dependent interaction of myosin heads with F-actin.

Looking for Targets to Restore the Contractile Function in Congenital Myopathy Caused by Gln147Pro Tropomyosin

We have used the technique of polarized microfluorimetry to obtain new insight into the pathogenesis of skeletal muscle disease caused by the Gln¹⁴⁷Pro substitution in β-tropomyosin (Tpm2.2). The spatial rearrangements of actin, myosin and tropomyosin in the single muscle fiber containing reconstituted thin filaments were studied during simulation of several stages of ATP hydrolysis cycle. The angular orientation of the fluorescence probes bound to tropomyosin was found to be changed by the substitution and was characteristic for a shift of tropomyosin strands closer to the inner actin domains. It was observed both in the absence and in the presence of troponin, Ca²⁺ and myosin heads at all simulated stages of the ATPase cycle. The mutant showed higher flexibility. Moreover, the Gln¹⁴⁷Pro substitution disrupted the myosin-induced displacement of tropomyosin over actin. The irregular positioning of the mutant tropomyosin caused premature activation of actin monomers and a tendency to increase the number of myosin cross-bridges in a state of strong binding with actin at low Ca²⁺.

Threonine-77 Is a Determinant of the Low-Temperature Conditioning of the Most Abundant Isoform of Tropomyosin in Atlantic Salmon

The Atlantic salmon Salmo salar survives below 10 °C. The main skeletal muscle is composed of a single isoform of tropomyosin (classified as Tpm1 α-fast) that is >92% identical to the mammalian homologue. How salmon Tpm1 maintains flexibility is investigated by reversing the only full charge substitution; threonine-77(g) in salmon and lysine in other vertebrates. The mutation (Thr-77 to Lys), which falls within a known destabilizing alanine cluster, (i) yields a useful electrophoretic shift in the absence and presence of an anionic detergent, (ii) increases the T_ms of both cooperative transitions (calorimetry, 0.1 M salt, pH 7) [35 °C (minor) and 44 °C (major); ΔT_m1 = 5 °C, ΔT_m2 = 3.5 °C], (iii) increases the T_m of CN1A (residues 11-127) to 53 °C (ΔT_m = 13 °C), a value similar to that of mammalian CN1A, (iv) markedly reduces the rate of proteolysis at Leu-169, and (v) weakens the affinity of salmon Tpm1 for troponin-Sepharose. Glu-82(e), the interstrand ionic partner of Lys-77(g), is conserved. The change in ionic interactions at this locus is postulated to be "sensed" in actin period 5 (residues 166-207) and likely beyond. Wild type (acetylated) salmon Tmp1 binds more tightly to F-actin at 4 °C than at 22 °C, which is the opposite of the long-known relationship displayed by the mammalian homologue. All of the evidence indicates that the presence of a neutral 77th amino acid destabilizes a sizable portion of salmon Tpm1 that includes the midregion. Threonine-77 is a key factor in rescuing the thin filament from the peril of cold-induced rigidity.

Protein-Protein Docking Reveals Dynamic Interactions of Tropomyosin on Actin Filaments

Experimental approaches such as fiber diffraction and cryo-electron microscopy reconstruction have defined regulatory positions of tropomyosin on actin but have not, as yet, succeeded at determining key atomic-level contacts between these proteins or fully substantiated the dynamics of their interactions at a structural level. To overcome this deficiency, we have previously employed computational approaches to deduce global dynamics of thin filament components by energy landscape determination and molecular dynamics simulations. Still, these approaches remain computationally challenging for any complex and large macromolecular assembly like the thin filament. For example, tropomyosin cable wrapping around actin of thin filaments features both head-to-tail polymeric interactions and local twisting, both of which depart from strict superhelical symmetry. This produces a complex energy surface that is difficult to model and thus to evaluate globally. Therefore, at this stage of our understanding, assessing global molecular dynamics can prove to be inherently impractical. As an alternative, we adopted a "divide and conquer" protocol to investigate actin-tropomyosin interactions at an atomistic level. Here, we first employed unbiased protein-protein docking tools to identify binding specificity of individual tropomyosin pseudorepeat segments over the actin surface. Accordingly, tropomyosin "ligand" segments were rotated and translated over potential "target" binding sites on F-actin where the corresponding interaction energetics of billions of conformational poses were ranked by the programs PIPER and ClusPro. These data were used to assess favorable interactions and then to rebuild models of seamless and continuous tropomyosin cables over the F-actin substrate, which were optimized further by flexible fitting routines and molecular dynamics. The models generated azimuthally distinct regulatory positions for tropomyosin cables along thin filaments on actin dominated by stereo-specific head-to-tail overlap linkage. The outcomes are in good agreement with current cryo-electron microscopy topology and consistent with long-thought residue-to-residue interactions between actin and tropomyosin.

Multiscale Models of Cardiac Muscle Biophysics and Tissue Remodeling in Hypertrophic Cardiomyopathies

Myocardial hypertrophy is the result of sustained perturbations to the mechanical and/or neurohormonal homeostasis of cardiac cells and is driven by integrated, multiscale biophysical and biochemical processes that are currently not well defined. In this brief review, we highlight recent computational and experimental models of cardiac hypertrophy that span mechanisms from the molecular level to the tissue level. Specifically, we focus on: (i) molecular-level models of the structural dynamics of sarcomere proteins in hypertrophic hearts, (ii) cellular-level models of excitation-contraction coupling and mechanosensitive signaling in disease-state myocytes, and (iii) organ-level models of myocardial growth kinematics and predictors thereof. Finally, we discuss how spanning these scales and combining multiple experimental/computational models will provide new information about the processes governing hypertrophy and potential methods to prevent or reverse them.

A Stochastic Multiscale Model of Cardiac Thin Filament Activation Using Brownian-Langevin Dynamics

We use Brownian-Langevin dynamics principles to derive a coarse-graining multiscale myofilament model that can describe the thin-filament activation process during contraction. The model links atomistic molecular simulations of protein-protein interactions in the thin-filament regulatory unit to sarcomere-level activation dynamics. We first calculate the molecular interaction energy between tropomyosin and actin surface using Brownian dynamics simulations. This energy profile is then generalized to account for the observed tropomyosin transitions between its regulatory stable states. The generalized energy landscape then served as a basis for developing a filament-scale model using Langevin dynamics. This integrated analysis, spanning molecular to thin-filament scales, is capable of tracking the events of the tropomyosin conformational changes as it moves over the actin surface. The tropomyosin coil with flexible overlap regions between adjacent tropomyosins is represented in the model as a system of coupled stochastic ordinary differential equations. The proposed multiscale approach provides a more detailed molecular connection between tropomyosin dynamics, the trompomyosin-actin interaction-energy landscape, and the generated force by the sarcomere.

Thin filament dysfunctions caused by mutations in tropomyosin Tpm3.12 and Tpm1.1

Tropomyosin is the major regulator of the thin filament. In striated muscle its function is to bind troponin complex and control the access of myosin heads to actin in a Ca²⁺-dependent manner. It also participates in the maintenance of thin filament length by regulation of tropomodulin and leiomodin, the pointed end-binding proteins. Because the size of the overlap between actin and myosin filaments affects the number of myosin heads which interact with actin, the filament length is one of the determinants of force development. Numerous point mutations in genes encoding tropomyosin lead to single amino acid substitutions along the entire length of the coiled coil that are associated with various types of cardiomyopathy and skeletal muscle disease. Specific regions of tropomyosin interact with different binding partners; therefore, the mutations affect diverse tropomyosin functions. In this review, results of studies on mutations in the genes TPM1 and TPM3, encoding Tpm1.1 and Tpm3.12, are described. The paper is particularly focused on mutation-dependent alterations in the mechanisms of actin-myosin interactions and dynamics of the thin filament at the pointed end.

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.

Molecular dynamics simulation of tropomyosin bound to actins/myosin in the closed and open states

Tropomyosin (Tpm) is a dimeric coiled-coil protein that binds to filamentous actin, and regulates actin-myosin interaction by moving between three positions corresponding to the blocked, closed, and open states. To elucidate how Tpm undergoes transitions between these functional states, we have built structural models and conducted extensive molecular dynamics simulations of the Tpm-actins/myosin complex in the closed and open states (total simulation time >1.4 μs). Based on the simulation trajectories, we have analyzed the dynamics and energetics of a truncated Tpm interacting with actins/myosin under the physiological conditions. Our simulations have shown distinct dynamics of four Tpm periods (P3-P6), featuring pronounced biased fluctuations of P4 and P5 toward the open position in the closed state, which is consistent with a conformational selection mechanism for Tpm-regulated myosin binding. Additionally, we have identified key residues of Tpm specifically binding to actins/myosin in the closed and open state. Some of them were validated as functionally important in comparison with past functional/clinical studies, and the rest will make promising targets for future mutational experiments.

Congenital myopathy-related mutations in tropomyosin disrupt regulatory function through altered actin affinity and tropomodulin binding

Tropomyosin (Tpm) binds along actin filaments and regulates myosin binding to control muscle contraction. Tropomodulin binds to the pointed end of a filament and regulates actin dynamics, which maintains the length of a thin filament. To define the structural determinants of these Tpm functions, we examined the effects of two congenital myopathy mutations, A4V and R91C, in the Tpm gene, TPM3, which encodes the Tpm3.12 isoform, specific for slow-twitch muscle fibers. Mutation A4V is located in the tropomodulin-binding, N-terminal region of Tpm3.12. R91C is located in the actin-binding period 3 and directly interacts with actin. The A4V and R91C mutations resulted in a 2.5-fold reduced affinity of Tpm3.12 homodimers for F-actin in the absence and presence of troponin, and a two-fold decrease in actomyosin ATPase activation in the presence of Ca2+ . Actomyosin ATPase inhibition in the absence of Ca2+ was not affected. The Ca2+ sensitivity of ATPase activity was decreased by R91C, but not by A4V. In vitro, R91C altered the ability of tropomodulin 1 (Tmod1) to inhibit actin polymerization at the pointed end of the filaments, which correlated with the reduced affinity of Tpm3.12-R91C for Tmod1. Molecular dynamics simulations of Tpm3.12 in complex with F-actin suggested that both mutations reduce the affinity of Tpm3.12 for F-actin binding by perturbing the van der Waals energy, which may be attributable to two different molecular mechanisms-a reduced flexibility of Tpm3.12-R91C and an increased flexibility of Tpm3.12-A4V.

The actin 'A-triad's' role in contractile regulation in health and disease

Striated muscle contraction is regulated by Ca²⁺ -dependent modulation of myosin cross-bridge binding to F-actin by the thin filament troponin (Tn)-tropomyosin (Tm) complex. In the absence of Ca²⁺ , Tn binds to actin and constrains Tm to an azimuthal location where it sterically occludes myosin binding sites along the thin filament surface. This limits force production and promotes muscle relaxation. In addition to Tn-actin interactions, inhibitory Tm positioning requires associations between other thin filament constituents. For example, the actin 'A-triad', composed of residues K326, K328 and R147, forms numerous, highly favourable electrostatic contacts with Tm that are critical for establishing its inhibitory azimuthal binding position. Here, we review recent findings, including the identification and interrogation of modifications within and proximal to the A-triad that are associated with disease and/or altered muscle behaviour, which highlight the surface feature's role in F-actin-Tm interactions and contractile regulation.

A new twist on tropomyosin binding to actin filaments: perspectives on thin filament function, assembly and biomechanics

Tropomyosin, best known for its role in the steric regulation of muscle contraction, polymerizes head-to-tail to form cables localized along the length of both muscle and non-muscle actin-based thin filaments. In skeletal and cardiac muscles, tropomyosin, under the control of troponin and myosin, moves in a cooperative manner between blocked, closed and open positions on filaments, thereby masking and exposing actin-binding sites necessary for myosin crossbridge head interactions. While the coiled-coil signature of tropomyosin appears to be simple, closer inspection reveals surprising structural complexity required to perform its role in steric regulation. For example, component α-helices of coiled coils are typically zippered together along a continuous core hydrophobic stripe. Tropomyosin, however, contains a number of anomalous, functionally controversial, core amino acid residues. We argue that the atypical residues at this interface, including clusters of alanines and a charged aspartate, are required for preshaping tropomyosin to readily fit to the surface of the actin filament, but do so without compromising tropomyosin rigidity once the filament is assembled. Indeed, persistence length measurements of tropomyosin are characteristic of a semi-rigid cable, in this case conducive to cooperative movement on thin filaments. In addition, we also maintain that tropomyosin displays largely unrecognized and residue-specific torsional variance, which is involved in optimizing contacts between actin and tropomyosin on the assembled thin filament. Corresponding twist-induced stiffness may also enhance cooperative translocation of tropomyosin across actin filaments. We conclude that anomalous core residues of tropomyosin facilitate thin filament regulatory behavior in a multifaceted way.

Sarcomere gene variants act as a genetic trigger underlying the development of left ventricular noncompaction

BACKGROUND

Left ventricular noncompaction (LVNC) is a primary cardiomyopathy with heterogeneous genetic origins. The aim of this study was to elucidate the role of sarcomere gene variants in the pathogenesis and prognosis of LVNC.

METHODS AND RESULTS

We screened 82 Japanese patients (0-35 years old), with a diagnosis of LVNC, for mutations in seven genes encoding sarcomere proteins, by direct DNA sequencing. We identified variants in a significant proportion of cases (27%), which were associated with poor prognosis (p = 0.012), particularly variants in TPM1, TNNC1, and ACTC1 (p = 0.012). To elucidate the pathological role, we developed and studied human-induced pluripotent stem cells (hiPSCs) from a patient carrying a TPM1 p.Arg178His mutation, who underwent heart transplantation. These cells displayed pathological changes, with mislocalization of tropomyosin 1, causing disruption of the sarcomere structure in cardiomyocytes, and impaired calcium handling. Microarray analysis indicated that the TPM1 mutation resulted in the down-regulation of the expression of numerous genes involved in heart development, and positive regulation of cellular process, especially the calcium signaling pathway.

CONCLUSIONS

Sarcomere genes are implicated as genetic triggers in the development of LVNC, regulating the expression of numerous genes involved in heart development, or modifying the severity of disease.

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.

Structural and Dynamic Properties of Allergen and Non-Allergen Forms of Tropomyosin

To what extent do structural and biophysical features of food allergen proteins distinguish them from other proteins in our diet? Invertebrate tropomyosins (Tpms) as a class are considered "pan-allergens," inducing food allergy to shellfish and respiratory allergy to dust mites. Vertebrate Tpms are not known to elicit allergy or cross-reactivity, despite their high structural similarity and sequence identity to invertebrate homologs. We expect allergens are sufficiently stable against gastrointestinal proteases to survive for immune sensitization in the intestines, and that proteolytic stability will correlate with thermodynamic stability. Thermal denaturation of shrimp Tpm shows that it is more stable than non-allergen vertebrate Tpm. Shrimp Tpm is also more resistant to digestion. Molecular dynamics uncover local dynamics that select epitopes and global differences in flexibility between shrimp and pig Tpm that discriminate allergens from non-allergens. Molecular determinants of allergenicity depend not only on sequence but on contributions of protein structure and dynamics.

The primary cause of muscle disfunction associated with substitutions E240K and R244G in tropomyosin is aberrant behavior of tropomyosin and response of actin and myosin during ATPase cycle

Using the polarized photometry technique we have studied the effects of two amino acid replacements, E240K and R244G, in tropomyosin (Tpm1.1) on the position of Tpm1.1 on troponin-free actin filaments and the spatial arrangement of actin monomers and myosin heads at various mimicked stages of the ATPase cycle in the ghost muscle fibres. E240 and R244 are located in the C-terminal, seventh actin-binding period, in f and b positions of the coiled-coil heptapeptide repeat. Actin, Tpm1.1, and myosin subfragment-1 (S1) were fluorescently labeled: 1.5-IAEDANS was attached to actin and S1, 5-IAF was bound to Tpm1.1. The labeled proteins were incorporated in the ghost muscle fibres and changes in polarized fluorescence during the ATPase cycle have been measured. It was found that during the ATPase cycle both mutant tropomyosins occupied a position close to the inner domain of actin. The relative amount of the myosin heads in the strongly-bound conformations and of the switched on actin monomers increased at mimicking different stages of the ATPase cycle. This might be one of the reasons for muscle dysfunction in congenital fibre type disproportion caused by the substitutions E240K and R244G in tropomyosin.

Functional role of the core gap in the middle part of tropomyosin

Tropomyosin (Tpm) is an α-helical coiled-coil actin-binding protein playing an essential role in the regulation of muscle contraction. The middle part of the Tpm molecule has some specific features, such as the presence of noncanonical residues as well as a substantial gap at the interhelical interface, which are believed to destabilize a coiled-coil and impart structural flexibility to this part of the molecule. To study how the gap affects structural and functional properties of α-striated Tpm (the Tpm1.1 isoform that is expressed in cardiac and skeletal muscles) we replaced large conserved apolar core residues located at both sides of the gap with smaller ones by mutations M127A/I130A and M141A/Q144A. We found that in contrast with the stabilizing substitutions D137L and G126R studied earlier, these substitutions have no appreciable influence on thermal unfolding and domain structure of the Tpm molecule. They also do not affect actin-binding properties of Tpm. However, they strongly increase sliding velocity of regulated actin filaments in an in vitro motility assay and cause an oversensitivity of the velocity to Ca²⁺ similar to the stabilizing substitutions D137L and G126R. Molecular dynamics shows that the substitutions studied here increase bending stiffness of the coiled-coil structure of Tpm, like that of G126R/D137L, probably due to closure of the interhelical gap in the area of the substitutions. Our results clearly indicate that the conserved middle part of Tpm is important for the fine tuning of the Ca²⁺ regulation of actin-myosin interaction in muscle.

The Relaxation Properties of Myofibrils Are Compromised by Amino Acids that Stabilize α-Tropomyosin

We investigated the functional impact of α-tropomyosin (Tm) substituted with one (D137L) or two (D137L/G126R) stabilizing amino acid substitutions on the mechanical behavior of rabbit psoas skeletal myofibrils by replacing endogenous Tm and troponin (Tn) with recombinant Tm mutants and purified skeletal Tn. Force recordings from myofibrils (15°C) at saturating [Ca²⁺] showed that Tm-stabilizing substitutions did not significantly affect the maximal isometric tension and the rates of force activation (k_ACT) and redevelopment (k_TR). However, a clear effect was observed on force relaxation: myofibrils with D137L/G126R or D137L Tm showed prolonged durations of the slow phase of relaxation and decreased rates of the fast phase. Both Tm-stabilizing substitutions strongly decreased the slack sarcomere length (SL) at submaximal activating [Ca²⁺] and increased the steepness of the SL-passive tension relation. These effects were reversed by addition of 10 mM 2,3-butanedione 2-monoxime. Myofibrils also showed an apparent increase in Ca²⁺ sensitivity. Measurements of myofibrillar ATPase activity in the absence of Ca²⁺ showed a significant increase in the presence of these Tms, indicating that single and double stabilizing substitutions compromise the full inhibition of contraction in the relaxed state. These data can be understood with the three-state (blocked-closed-open) theory of muscle regulation, according to which the mutations increase the contribution of the active open state in the absence of Ca²⁺ (M^-). Force measurements on myofibrils substituted with C-terminal truncated TnI showed similar compromised relaxation effects, indicating the importance of TnI-Tm interactions in maintaining the blocked state. It appears that reducing the flexibility of native Tm coiled-coil structure decreases the optimum interactions of the central part of Tm with the C-terminal region of TnI. This results in a shift away from the blocked state, allowing myosin binding and activity in the absence of Ca²⁺. This work provides a basis for understanding the effects of disease-producing mutations in muscle proteins.

Investigating the inter-subunit/subdomain interactions and motions relevant to disease mutations in the N-terminal domain of ryanodine receptors by molecular dynamics simulation

The ryanodine receptors (RyR) are essential to calcium signaling in striated muscles, and numerous disease mutations have been identified in two RyR isoforms, RyR1 in skeletal muscle and RyR2 in cardiac muscle. A deep understanding of the activation/regulation mechanisms of RyRs has been hampered by the shortage of high-resolution structures and dynamic information for this giant tetrameric complex in different functional states. Toward elucidating the molecular mechanisms of disease mutations in RyRs, we performed molecular dynamics simulation of the N-terminal domain (NTD) which is not only the best-resolved structural component of RyRs, but also a hotspot of disease mutations. First, we simulated the tetrameric NTD of wild-type RyR1 and three disease mutants (K155E, R157Q, and R164Q) that perturb the inter-subunit interfaces. Our simulations identified a dynamic network of salt bridges involving charged residues at the inter-subunit/subdomain interfaces and disease-mutation sites. By perturbing this key network, the above three mutations result in greater flexibility with the highest inter-subunit opening probability for R157Q. Next, we simulated the monomeric NTD of RyR2 in the presence or absence of a central Cl^- anion which is known to stabilize the interfaces between the three NTD subdomains (A, B, and C). We found that the loss of Cl^- restructures the salt-bridge network near the Cl^- -binding site, leading to rotations of subdomain A/B relative to subdomain C and enhanced mobility between the subdomains. This finding supports a mechanism for disease mutations in the NTD of RyR2 via perturbation of the Cl^- binding. The rich structural and dynamic information gained from this study will guide future mutational and functional studies of the NTD of RyRs. Proteins 2017; 85:1633-1644. © 2017 Wiley Periodicals, Inc.

Insights and Challenges of Multi-Scale Modeling of Sarcomere Mechanics in cTn and Tm DCM Mutants-Genotype to Cellular Phenotype

Dilated Cardiomyopathy (DCM) is a leading cause of sudden cardiac death characterized by impaired pump function and dilatation of cardiac ventricles. In this review we discuss various in silico approaches to elucidating the mechanisms of genetic mutations leading to DCM. The approaches covered in this review focus on bridging the spatial and temporal gaps that exist between molecular and cellular processes. Mutations in sarcomeric regulatory thin filament proteins such as the troponin complex (cTn) and Tropomyosin (Tm) have been associated with DCM. Despite the experimentally-observed myofilament measures of contractility in the case of these mutations, the mechanisms by which the underlying molecular changes and protein interactions scale up to organ failure by these mutations remains elusive. The review highlights multi-scale modeling approaches and their applicability to study the effects of sarcomeric gene mutations in-silico. We discuss some of the insights that can be gained from computational models of cardiac biomechanics when scaling from molecular states to cellular level.

Probing the Energetics of Dynactin Filament Assembly and the Binding of Cargo Adaptor Proteins Using Molecular Dynamics Simulation and Electrostatics-Based Structural Modeling

Dynactin, a large multiprotein complex, binds with the cytoplasmic dynein-1 motor and various adaptor proteins to allow recruitment and transportation of cellular cargoes toward the minus end of microtubules. The structure of the dynactin complex is built around an actin-like minifilament with a defined length, which has been visualized in a high-resolution structure of the dynactin filament determined by cryo-electron microscopy (cryo-EM). To understand the energetic basis of dynactin filament assembly, we used molecular dynamics simulation to probe the intersubunit interactions among the actin-like proteins, various capping proteins, and four extended regions of the dynactin shoulder. Our simulations revealed stronger intersubunit interactions at the barbed and pointed ends of the filament and involving the extended regions (compared with the interactions within the filament), which may energetically drive filament termination by the capping proteins and recruitment of the actin-like proteins by the extended regions, two key features of the dynactin filament assembly process. Next, we modeled the unknown binding configuration among dynactin, dynein tails, and a number of coiled-coil adaptor proteins (including several Bicaudal-D and related proteins and three HOOK proteins), and predicted a key set of charged residues involved in their electrostatic interactions. Our modeling is consistent with previous findings of conserved regions, functional sites, and disease mutations in the adaptor proteins and will provide a structural framework for future functional and mutational studies of these adaptor proteins. In sum, this study yielded rich structural and energetic information about dynactin and associated adaptor proteins that cannot be directly obtained from the cryo-EM structures with limited resolutions.

Predicting Effects of Tropomyosin Mutations on Cardiac Muscle Contraction through Myofilament Modeling

Point mutations to the human gene TPM1 have been implicated in the development of both hypertrophic and dilated cardiomyopathies. Such observations have led to studies investigating the link between single residue changes and the biophysical behavior of the tropomyosin molecule. However, the degree to which these molecular perturbations explain the performance of intact sarcomeres containing mutant tropomyosin remains uncertain. Here, we present a modeling approach that integrates various aspects of tropomyosin's molecular properties into a cohesive paradigm representing their impact on muscle function. In particular, we considered the effects of tropomyosin mutations on (1) persistence length, (2) equilibrium between thin filament blocked and closed regulatory states, and (3) the crossbridge duty cycle. After demonstrating the ability of the new model to capture Ca-dependent myofilament responses during both dynamic and steady-state activation, we used it to capture the effects of hypertrophic cardiomyopathy (HCM) related E180G and D175N mutations on skinned myofiber mechanics. Our analysis indicates that the fiber-level effects of the two mutations can be accurately described by a combination of changes to the three tropomyosin properties represented in the model. Subsequently, we used the model to predict mutation effects on muscle twitch. Both mutations led to increased twitch contractility as a consequence of diminished cooperative inhibition between thin filament regulatory units. Overall, simulations suggest that a common twitch phenotype for HCM-linked tropomyosin mutations includes both increased contractility and elevated diastolic tension.