Probing the flexibility of tropomyosin and its binding to filamentous actin using molecular dynamics simulations

The hypertrophic cardiomyopathy-associated A331P actin variant enhances basal contractile activity and elicits resting muscle dysfunction

Previous studies aimed at defining the mechanistic basis of hypertrophic cardiomyopathy caused by A331P cardiac actin have reported conflicting results. The mutation is located along an actin surface strand, proximal to residues that interact with tropomyosin. These F-actin-tropomyosin associations are vital for proper contractile inhibition. To help resolve disease pathogenesis, we implemented a multidisciplinary approach. Transgenic Drosophila, expressing A331P actin, displayed skeletal muscle hypercontraction and elevated basal myocardial activity. A331P thin filaments, reconstituted using recombinant human cardiac actin, exhibited higher in vitro myosin-based sliding speeds, exclusively at low Ca2+ concentrations. Cryo-EM-based reconstructions revealed no detectable A331P-related structural perturbations in F-actin. In silico, however, the P331-containing actin surface strand was less mobile and established diminished van der Waal's attractive forces with tropomyosin, which correlated with greater variability in inhibitory tropomyosin positioning. Such mutation-induced effects potentially elevate resting contractile activity among our models and may stimulate pathology in patients.

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.

Divergent Plasmodium actin residues are essential for filament localization, mosquito salivary gland invasion and malaria transmission

Actin is one of the most conserved and ubiquitous proteins in eukaryotes. Its sequence has been highly conserved for its monomers to self-assemble into filaments that mediate essential cell functions such as trafficking, cell shape and motility. The malaria-causing parasite, Plasmodium, expresses a highly sequence divergent actin that is critical for its rapid motility at different stages within its mammalian and mosquito hosts. Each of Plasmodium actin’s four subdomains have divergent regions compared to canonical vertebrate actins. We previously identified subdomains 2 and 3 as providing critical contributions for parasite actin function as these regions could not be replaced by subdomains of vertebrate actins. Here we probed the contributions of individual divergent amino acid residues in these subdomains on parasite motility and progression. Non-lethal changes in these subdomains did not affect parasite development in the mammalian host but strongly affected progression through the mosquito with striking differences in transmission to and through the insect. Live visualization of actin filaments showed that divergent amino acid residues in subdomains 2 and 4 enhanced localization associated with filaments, while those in subdomain 3 negatively affected actin filaments. This suggests that finely tuned actin dynamics are essential for efficient organ entry in the mosquito vector affecting malaria transmission. This work provides residue level insight on the fundamental requirements of actin in highly motile cells.

Actin is one of the most abundant and conserved proteins known. Actin monomers can join together to form long filaments. The malaria-causing parasite is transmitted by mosquitoes and needs actin to move very rapidly. An actin from the parasite is different to other actins: its amino acid sequence has relatively high amounts of changes compared to animal species and the actin tends to form only short filaments. We previously identified two large parts of the protein that were critical for the parasite since these large parts could not be exchanged with the equivalent regions of other species. In this study, we focused in on these regions by making more discrete mutations. Most mutations of the actin sequence were tolerated by the parasite in the blood stages. However, these mutants has striking defects in progressing through mosquitoes, especially in invading its salivary glands. We used a new filament labeler to visualize how these mutations affect the actin filaments and found surprisingly different effects. Taken together, small changes to the sequence can have large consequences for the parasite, which ultimately affects its ability to transmit to a new host.

Structural insights into the tropomodulin assembly at the pointed ends of actin filaments

Tropomodulins are a family of important regulators of actin dynamics at the pointed ends of actin filaments. Four isoforms of tropomodulin, Tmod1-Tmod4, are expressed in vertebrates. Binding of tropomodulin to the pointed end is dependent on tropomyosin, an actin binding protein that itself is represented in mammals by up to 40 isoforms. The understanding of the regulatory role of the tropomodulin/tropomyosin molecular diversity has been limited due to the lack of a three-dimensional structure of the tropomodulin/tropomyosin complex. In this study, we mapped tropomyosin residues interacting with two tropomyosin-binding sites of tropomodulin and generated a three-dimensional model of the tropomodulin/tropomyosin complex for each of these sites. The models were refined by molecular dynamics simulations and validated via building a self-consistent three-dimensional model of tropomodulin assembly at the pointed end. The model of the pointed-end Tmod assembly offers new insights in how Tmod binding ensures tight control over the pointed end dynamics.

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.

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.

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.

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.

Cardiomyopathy Mutation Alters End-to-End Junction of Tropomyosin and Reduces Calcium Sensitivity

Muscle contraction is governed by tropomyosin (Tpm) shifting azimuthally between three states on F-actin (B-, C-, and M-states) in response to calcium binding to troponin and actomyosin cross-bridge formation. The Tpm coiled coil polymerizes head to tail along the long-pitch helix of F-actin to form continuous superhelical cables that wrap around the actin filaments. The end-to-end bonds formed between the N- and C-terminus of adjacent Tpm molecules define Tpm continuity and play a critical role in the ability of Tpm to cooperatively bind to actin, thus facilitating Tpm conformational switching to cooperatively propagate along F-actin. We expect that a missense mutation in this critical overlap region associated with dilated cardiomyopathy, A277V, will alter Tpm binding and thin filament activation by altering the overlap structure. Here, we used cosedimentation assays and in vitro motility assays to determine how the mutation alters Tpm binding to actin and its ability to regulate actomyosin interactions. Analytical viscometry coupled with molecular dynamics simulations showed that the A277V mutation results in enhanced Tpm end-to-end bond strength and a reduced curvature of the Tpm overlap domain. The mutant Tpm exhibited enhanced actin-Tpm binding affinity, consistent with overlap stabilization. The observed A277V-induced decrease in cooperative activation observed with regulated thin filament motility indicates that increased overlap stabilization is not correlated with Tpm-Tpm overlap binding strength or mechanical rigidity as is often assumed. Instead, A277V-induced structural changes result in local and delocalized increases in Tpm flexibility and prominent coiled-coil twisting in pseudorepeat 4. An A277V-induced decrease in Ca²⁺ sensitivity, consistent with a mutation-induced bolstering of the B-state Tpm-actin electrostatic contacts and an increased Tpm troponin T1 binding affinity, was also observed.

Comparative dynamics of tropomyosin in vertebrates and invertebrates

Tropomyosin (Tpm) is an extended α-helical coiled-coil homodimer that regulates actinomyosin interactions in muscle. Molecular simulations of four Tpms, two from the vertebrate class Mammalia (rat and pig), and two from the invertebrate class Malacostraca (shrimp and lobster), showed that despite extensive sequence and structural homology across metazoans, dynamic behavior-particularly long-range structural fluctuations-were clearly distinct. Vertebrate Tpms were more flexible and sampled complex, multi-state conformational landscapes. Invertebrate Tpms were more rigid, sampling a highly constrained harmonic landscape. Filtering of trajectories by principle component analysis into essential subspaces showed significant overlap within but not between phyla. In vertebrate Tpms, hinge-regions decoupled long-range interhelical motions and suggested distinct domains. In contrast, crustacean Tpms did not exhibit long-range dynamic correlations-behaving more like a single rigid rod on the nanosecond time scale. These observations suggest there may be divergent mechanisms for Tpm binding to actin filaments, where conformational flexibility in mammalian Tpm allows a preorganized shape complementary to the filament surface, and where rigidity in the crustacean Tpm requires concerted bending and binding.

Tropomyosin dynamics during cardiac muscle contraction as governed by a multi-well energy landscape

The dynamic oscillations of tropomyosin molecules in the azimuthal direction over the surface of the actin filament during thin filament activation are studied here from an energy landscape perspective. A mathematical model based on principles from nonlinear dynamics and chaos theory is derived to describe these dynamical motions. In particular, an energy potential with three wells is proposed to govern the tropomyosin oscillations between the observed regulatory positions observed during muscle contraction, namely the blocked "B", closed "C" and open "M" states. Based on the variations in both the frequency and amplitude of the environmental (surrounding the thin filament system) driving tractions, such as the electrostatic, hydrophobic, and Ca2+-dependent forces, the tropomyosin movements are shown to be complex; they can change from being simple harmonic oscillations to being fully chaotic. Three cases (periodic, period-2, and chaotic patterns) are presented to showcase the different possible dynamic responses of tropomyosin sliding over the actin filament. A probability density function is used as a statistical measure to calculate the average residence time spanned out by the tropomyosin molecule when visiting each (B, C, M) equilibrium state. The results were found to depend strongly on the energy landscape profile and its featured barriers, which normally govern the transitions between the B-C-M states during striated muscle activation.

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.

Effects of an Interchain Disulfide Bond on Tropomyosin Structure: A Molecular Dynamics Study

Tropomyosin (Tpm) is a coiled-coil actin-binding dimer protein that participates in the regulation of muscle contraction. Both Tpm chains contain Cys190 residues which are normally in the reduced state, but form an interchain disulfide bond in failing heart. Changes in structural and functional properties of Tpm and its complexes with actin upon disulfide cross-linking were studied using various experimental methods. To understand the molecular mechanism underlying these changes and to reveal the possible mechanism of the involvement of the cross-linking in heart failure, molecular dynamics (MD) simulations of the middle part of Tpm were performed in cross-linked and reduced states. The cross-linking increased bending stiffness of Tpm assessed from MD trajectories at 27 °C in agreement with previous experimental observations. However, at 40 °C, the cross-linking caused a decrease in Tpm stiffness and a significant reduction in the number of main chain hydrogen bonds in the vicinity of residues 133 and 134. These data are in line with observations showing enhanced thermal unfolding of the least stable part of Tpm at 30–40 °C and accelerated trypsin cleavage at residue 133 at 40 °C (but not at 27 °C) upon cross-linking. These results allow us to speculate about the possible mechanism of involvement of Tpm cross-linking to heart failure pathogenesis.

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.

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.

Tropomyosin Structure, Function, and Interactions: A Dynamic Regulator

Tropomyosin is the archetypal-coiled coil, yet studies of its structure and function have proven it to be a dynamic regulator of actin filament function in muscle and non-muscle cells. Here we review aspects of its structure that deviate from canonical leucine zipper coiled coils that allow tropomyosin to bind to actin, regulate myosin, and interact directly and indirectly with actin-binding proteins. Four genes encode tropomyosins in vertebrates, with additional diversity that results from alternate promoters and alternatively spliced exons. At the same time that periodic motifs for binding actin and regulating myosin are conserved, isoform-specific domains allow for specific interaction with myosins and actin filament regulatory proteins, including troponin. Tropomyosin can be viewed as a universal regulator of the actin cytoskeleton that specifies actin filaments for cellular and intracellular functions.

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.

The Importance of Intrinsically Disordered Segments of Cardiac Troponin in Modulating Function by Phosphorylation and Disease-Causing Mutations

Troponin plays a central role in regulation of muscle contraction. It is the Ca²⁺ switch of striated muscles including the heart and in the cardiac muscle it is physiologically modulated by PKA-dependent phosphorylation at Ser22 and 23. Many cardiomyopathy-related mutations affect Ca²⁺ regulation and/or disrupt the relationship between Ca²⁺ binding and phosphorylation. Unlike the mechanism of heart activation, the modulation of Ca²⁺-sensitivity by phosphorylation of the cardiac specific N-terminal segment of TnI (1–30) is structurally subtle and has proven hard to investigate. The crystal structure of cardiac troponin describes only the relatively stable core of the molecule and the crucial mobile parts of the molecule are missing including TnI C-terminal region, TnI (1–30), TnI (134–149) (“inhibitory” peptide) and the C-terminal 28 amino acids of TnT that are intrinsically disordered. Recent studies have been performed to answer this matter by building structural models of cardiac troponin in phosphorylated and dephosphorylated states based on peptide NMR studies. Now these have been updated by more recent concepts derived from molecular dynamic simulations treating troponin as a dynamic structure. The emerging model confirms the stable core structure of troponin and the mobile structure of the intrinsically disordered segments. We will discuss how we can describe these segments in terms of dynamic transitions between a small number of states, with the probability distributions being altered by phosphorylation and by HCM or DCM-related mutations that can explain how Ca²⁺-sensitivity is modulated by phosphorylation and the effects of mutations.

The propensity for tropomyosin twisting in the presence and absence of F-actin

A canonical model of muscle α-tropomyosin (Tpm1.1), based on molecular-mechanics and electron microscopy of different contractile states, shows that the two-stranded coiled-coiled is pre-bent to present a specific molecular-face to the F-actin filament. This conformation is thought to facilitate both filament assembly and tropomyosin sliding across actin to modulate myosin-binding. However, to bind effectively to actin filaments, the 42 nm-long tropomyosin coiled-coil is not strictly canonical. Here, the mid-region of tropomyosin twists an additional ∼20° in order to better match the F-actin helix. In addition, the N- and C-terminal regions of tropomyosin polymerize head-to-tail to form continuous super-helical cables. In this case, 9 to 10 residue-long overlapping domains between adjacent molecules untwist relative to each other to accommodate orthogonal interactions between chains in the junctional four-helix nexus. Extensive molecular dynamics simulations show that the twisting and untwisting motions of tropomyosin vary appreciably along tropomyosin length, and in particular that substantial terminal domain winding and unwinding occurs whether tropomyosin is bound to F-actin or not. The local and regional twisting and untwisting do not appear to proceed in a concerted fashion, resembling more of a "wringing-type" behavior rather than a rotation.

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

Coarse-grained modeling and all-atom molecular dynamics simulation provide insight into the mechanism for heat activation of TRPV1 gating.

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.

Tropomyosin diffusion over actin subunits facilitates thin filament assembly

Coiled-coil tropomyosin binds to consecutive actin-subunits along actin-containing thin filaments. Tropomyosin molecules then polymerize head-to-tail to form cables that wrap helically around the filaments. Little is known about the assembly process that leads to continuous, gap-free tropomyosin cable formation. We propose that tropomyosin molecules diffuse over the actin-filament surface to connect head-to-tail to partners. This possibility is likely because (1) tropomyosin hovers loosely over the actin-filament, thus binding weakly to F-actin and (2) low energy-barriers provide tropomyosin freedom for 1D axial translation on F-actin. We consider that these unique features of the actin-tropomyosin interaction are the basis of tropomyosin cable formation.

Evolutionarily conserved sites in yeast tropomyosin function in cell polarity, transport and contractile ring formation

Tropomyosin is a coiled-coil protein that binds and regulates actin filaments. The tropomyosin gene in Schizosaccharomyces pombe, cdc8, is required for formation of actin cables, contractile rings, and polar localization of actin patches. The roles of conserved residues were investigated in gene replacement mutants. The work validates an evolution-based approach to identify tropomyosin functions in living cells and sites of potential interactions with other proteins. A cdc8 mutant with near-normal actin affinity affects patch polarization and vacuole fusion, possibly by affecting Myo52p, a class V myosin, function. The presence of labile residual cell attachments suggests a delay in completion of cell division and redistribution of cell patches following cytokinesis. Another mutant with a mild phenotype is synthetic negative with GFP-fimbrin, inferring involvement of the mutated tropomyosin sites in interaction between the two proteins. Proteins that assemble in the contractile ring region before actin do so in a mutant cdc8 strain that cannot assemble condensed actin rings, yet some cells can divide. Of general significance, LifeAct-GFP negatively affects the actin cytoskeleton, indicating caution in its use as a biomarker for actin filaments.

Tropomyosin as a Regulator of Actin Dynamics

Tropomyosin is a major regulatory protein of contractile systems and cytoskeleton, an actin-binding protein that positions laterally along actin filaments and modulates actin-myosin interaction. About 40 tropomyosin isoforms have been found in a variety of cytoskeleton systems, not necessarily connected with actin-myosin interaction and contraction. Involvement of specific tropomyosin isoforms in the regulation of key cell processes was shown, and specific features of tropomyosin genes and protein structure have been investigated with molecular biology and genetics approaches. However, the mechanisms underlying the effects of tropomyosin on cytoskeleton dynamics are still unclear. As tropomyosin is primarily an F-actin-binding protein, it is important to understand how it interacts both with actin and actin-binding proteins functioning in muscles and cytoskeleton to regulate actin dynamics. This review focuses on biochemical data on the effects of tropomyosin on actin assembly and dynamics, as well as on the modulation of these effects by actin-binding proteins. The data indicate that tropomyosin can efficiently regulate actin dynamics via allosteric conformational changes within actin filaments.

Skip residues modulate the structural properties of the myosin rod and guide thick filament assembly

The rod of sarcomeric myosins directs thick filament assembly and is characterized by the insertion of four skip residues that introduce discontinuities in the coiled-coil heptad repeats. We report here that the regions surrounding the first three skip residues share high structural similarity despite their low sequence homology. Near each of these skip residues, the coiled-coil transitions to a nonclose-packed structure inducing local relaxation of the superhelical pitch. Moreover, molecular dynamics suggest that these distorted regions can assume different conformationally stable states. In contrast, the last skip residue region constitutes a true molecular hinge, providing C-terminal rod flexibility. Assembly of myosin with mutated skip residues in cardiomyocytes shows that the functional importance of each skip residue is associated with rod position and reveals the unique role of the molecular hinge in promoting myosin antiparallel packing. By defining the biophysical properties of the rod, the structures and molecular dynamic calculations presented here provide insight into thick filament formation, and highlight the structural differences occurring between the coiled-coils of myosin and the stereotypical tropomyosin. In addition to extending our knowledge into the conformational and biological properties of coiled-coil discontinuities, the molecular characterization of the four myosin skip residues also provides a guide to modeling the effects of rod mutations causing cardiac and skeletal myopathies.

Structure and flexibility of the tropomyosin overlap junction

To be effective as a gatekeeper regulating the access of binding proteins to the actin filament, adjacent tropomyosin molecules associate head-to-tail to form a continuous super-helical cable running along the filament surface. Chimeric head-to-tail structures have been solved by NMR and X-ray crystallography for N- and C-terminal segments of smooth and striated muscle tropomyosin spliced onto non-native coiled-coil forming peptides. The resulting 4-helix complexes have a tight coiled-coil N-terminus inserted into a separated pair of C-terminal helices, with some helical unfolding of the terminal chains in the striated muscle peptides. These overlap complexes are distinctly curved, much more so than elsewhere along the superhelical tropomyosin cable. To verify whether the non-native protein adducts (needed to stabilize the coiled-coil chimeras) perturb the overlap, we carried out Molecular Dynamics simulations of head-to-tail structures having only native tropomyosin sequences. We observe that the splayed chains all refold and become helical. Significantly, the curvature of both the smooth and the striated muscle overlap domain is reduced and becomes comparable to that of the rest of the tropomyosin cable. Moreover, the measured flexibility across the junction is small. This and the reduced curvature ensure that the super-helical cable matches the contours of F-actin without manifesting localized kinking and excessive flexibility, thus enabling the high degree of cooperativity in the regulation of myosin accessibility to actin filaments.

Gly126Arg substitution causes anomalous behaviour of α-skeletal and β-smooth tropomyosins during the ATPase cycle

To investigate how TM stabilization induced by the Gly126Arg mutation in skeletal α-TM or in smooth muscle β-TM affects the flexibility of TMs and their position on troponin-free thin filaments, we labelled the recombinant wild type and mutant TMs with 5-IAF and F-actin with FITC-phalloidin, incorporated them into ghost muscle fibres and studied polarized fluorescence at different stages of the ATPase cycle. It has been shown that in the myosin- and troponin-free filaments the Gly126Arg mutation causes a shift of TM strands towards the outer domain of actin, reduces the number of switched on actin monomers and decreases the rigidity of the C-terminus of α-TM and increases the rigidity of the N-terminus of β-TMs. The binding of myosin subfragment-1 to the filaments shifted the wild type TMs towards the inner domain of actin, decreased the flexibility of both terminal parts of TMs, and increased the number of switched on actin monomers. Multistep alterations in the position of α- and β-TMs and actin monomers in the filaments and in the flexibility of TMs and F-actin during the ATPase cycle were observed. The Gly126Arg mutation uncouples a correlation between the position of TM and the number of the switched on actin monomers in the filaments.