Spatial distribution of filament elasticity determines the migratory behaviors of a cell

Any cellular response leading to morphological changes is highly tuned to balance the force generated from structural reorganization, provided by actin cytoskeleton. Actin filaments serve as the backbone of intracellular force, and transduce external mechanical signal via focal adhesion complex into the cell. During migration, cells not only undergo molecular changes but also rapid mechanical modulation. Here we focus on determining, the role of spatial distribution of mechanical changes of actin filaments in epithelial, mesenchymal, fibrotic and cancer cells with non-migration, directional migration, and non-directional migration behaviors using the atomic force microscopy. We found 1) non-migratory cells only generated one type of filament elasticity, 2) cells generating spatially distributed two types of filament elasticity showed directional migration, and 3) pathologic cells that autonomously generated two types of filament elasticity without spatial distribution were actively migrating non-directionally. The demonstration of spatial regulation of filament elasticity of different cell types at the nano-scale highlights the coupling of cytoskeletal function with physical characters at the sub-cellular level, and provides new research directions for migration related disease.     

Structure and elasticity of desmin protofibrils explored with scanning force microscopy

Desmin filaments form the intermediate filament system of muscle cells where they play important role in maintaining mechanical integrity and elasticity. Although the importance of desmin elasticity and assembly-disassembly dynamics in cellular mechanics is being increasingly recognized, the molecular basis of neither desmin's elasticity nor its disassembly pathway is well understood. In the present work, we explored the topographical structure of purified and reconstituted desmin filaments by using scanning force microscopy. With the addition of divalent cation chelators ethyleneglycoltetraacetic acid or ethylenediaminetetraacetic acid, the filaments disassembled on a time scale of hours to days into stable, thin fibrillar components with variable (up to micrometer) length, smooth surface and uniform thickness, which are identified as protofibrils. Desmin protofibrils appear as elastic structures with a persistence length of 51.5 nm, and their Young's modulus (10.6 MPa) far exceeds that of the mature filament (3.7 MPa). Protofibrillar bundling within the desmin filament results in high longitudinal tensile strength at a large bending flexibility. The stability of protofibrils suggests that desmin may disassemble along a pathway quite distinct from its assembly.     

Direct measurement of single synthetic vertebrate thick filament elasticity using nanofabricated cantilevers

Thick filaments are generally thought to be effectively inextensible. Here we use novel nanofabricated cantilevers to carry out the first direct force-elongation measurements on single vertebrate thick filaments. Cantilevers are ideal for these experiments: force ranges are from pico- to micronewtons, specimens can be visualized during the experiment, and attachment surfaces are in the same plane as the filament. Synthetic thick filaments from rabbit myosin were suspended between two cantilevers and stretched. With stretch, stiffness increased gradually and then became nearly constant after approximately 100 pN. Stretch rate had little or no effect on force-elongation behavior. Under physiological loads (approximately 240 pN axially averaged with full activation) filaments elongated by 1.1 +/- 0.3%. Previous x-ray diffraction results showed a 1.0 to 1.5% increase in myosin head spacing with activation; however, this increase in spacing has been interpreted as change in the state of the cross-bridges, not as elasticity in the thick filament backbone. Comparison with our data suggests that changes in the myosin x-ray reflections seen during activation may be due to elongation of the thick filament backbone. Recognition of thick filament elasticity is important because it affects the interpretation of mechanical experiments and inferences drawn on the molecular mechanism of contraction.     

Effect of tensile force on the mechanical behavior of actin filaments

Actin filaments are the most abundant components of the cellular cytoskeleton, and play critical roles in various cellular functions such as migration, division and shape control. In these activities, mechanical tension causes structural changes in the double-helical structure of the actin filament, which is a key modulator of cytoskeletal reorganization. This study performed large-scale molecular dynamics (MD) and steered MD simulations to quantitatively analyze the effects of tensile force on the mechanical behavior of actin filaments. The results revealed that when a tensile force of 200pN was applied to a filament consisting of 14 actin subunits, the twist angle of the filament decreased by approximately 20°, corresponding to a rotation of approximately -2° per subunit, representing a critical structural change in actin filaments. Based on these structural changes, the variance in filament length and twist angle was found to decrease, leading to increases in extensional and torsional stiffness. Torsional stiffness increased significantly under the tensile condition, and the ratio of filament stiffness under tensile force to that under no external force increased significantly on longer temporal scales. The results obtained from this study contribute to the understanding of mechano-chemical interactions concerning actin dynamics, showing that increased tensile force in the filament prevents actin regulatory proteins from binding to the filament.     

Molecular investigation of the mechanical properties of single actin filaments based on vibration analyses

The mechanical vibration properties of single actin filaments from 50 to 288 nm are investigated by the molecular dynamics simulation in this study. The natural frequencies obtained from the molecular simulations agree with those obtained from the analytical solution of the equivalent Euler–Bernoulli beam model. Through the convergence study of the mechanical properties with respect to the filament length, it was found that the Euler–Bernoulli beam model can only be reliably used when the single actin filament is of the order of hundreds of nanometre scale. This molecular investigation not only provides the evidence for the use of the continuum beam model in characterising the mechanical properties of single actin filaments, but also clarifies the criteria for the effective use of the Euler–Bernoulli beam model.     

Inhibitor discovery targeting the intermediate structure of beta-amyloid peptide on the conformational transition pathway: implications in the aggregation mechanism of beta-amyloid peptide

Abeta peptides cleaved from the amyloid precursor protein are the main components of senile plaques in Alzheimer's disease. Abeta peptides adopt a conformation mixture of random coil, beta-sheet, and alpha-helix in solution, which makes it difficult to design inhibitors based on the 3D structures of Abeta peptides. By targeting the C-terminal beta-sheet region of an Abeta intermediate structure extracted from molecular dynamics simulations of Abeta conformational transition, a new inhibitor that abolishes Abeta fibrillation was discovered using virtual screening in conjunction with thioflavin T fluorescence assay and atomic force microscopy determination. Circular dichroism spectroscopy demonstrated that the binding of the inhibitor increased the beta-sheet content of Abeta peptides either by stabilizing the C-terminal beta-sheet conformation or by inducing the intermolecular beta-sheet formation. It was proposed that the inhibitor prevented fibrillation by blocking interstrand hydrogen bond formation of the pleated beta-sheet structure commonly found in amyloid fibrils. The study not only provided a strategy for inhibitor design based on the flexible structures of amyloid peptides but also revealed some clues to understanding the molecular events involved in Abeta aggregation.     

Nanomechanical properties of desmin intermediate filaments

Desmin intermediate filaments play important role in the mechanical integrity and elasticity of muscle cells. The mechanisms of how desmin contributes to cellular mechanics are little understood. Here, we explored the nanomechanics of desmin by manipulating individual filaments with atomic force microscopy. In complex, hierarchical force responses we identified recurring features which likely correspond to distinct properties and structural transitions related to desmin's extensibility and elasticity. The most frequently observed feature is an initial unbinding transition that corresponds to the removal of approximately 45-nm-long coiled-coil dimers from the filament surface with 20-60 pN forces in usually two discrete steps. In tethers longer than 60 nm we most often observed force plateaus studded with bumps spaced approximately 16 nm apart, which are likely caused by a combination of protofilament unzipping, dimer-dimer sliding and coiled-coil-domain unfolding events. At high stresses and strains non-linear, entropic elasticity was dominant, and sometimes repetitive sawtooth force transitions were seen which might arise because of slippage within the desmin protofilament. A model is proposed in which mechanical yielding is caused by coiled-coil domain unfolding and dimer-dimer sliding/slippage, and strain hardening by the entropic elasticity of partially unfolded protofilaments.     

Amyloid fibrils of glucagon characterized by high-resolution atomic force microscopy

Glucagon solutions at pH 2.0 were subjected to mechanical agitation at 37 degrees C in the presence of a hydrophobic surface to explore the details of aggregation and fiber formation. High-resolution intermittent-contact atomic force microscopy performed in solution revealed the presence of aggregates after 0.5 h; however, longer agitation times resulted in the formation of fibrillated structures with varying levels of higher-order assembly. Height, periodicity, and amplitude measurements of these structures allowed the identification of four distinct fiber types. The most elementary fiber form, designated a filament, self-associates in a specific wound fashion to produce protofibrils composed of two filaments. Subsequent self-assembly of these filaments and protofibrils leads to two well-defined fibrillar motifs, termed Type I and Type II. Atomic force microscopy imaging of pH 2.8 glucagon solutions not agitated or exposed to elevated temperature revealed the presence of amorphous aggregates before the formation of fibrillar structures similar to those seen at pH 2.0. Time-course solution Fourier transform infrared spectroscopy and thioflavin T binding studies suggested that glucagon aggregation and fibril formation were associated with the development of beta-sheet structure. The results of these studies are used to describe a possible mechanism for glucagon aggregation and fibrillation that is consistent with a hierarchical assembly model proposed for amyloid fibril formation.     

Molecular dynamics study of interactions between polymorphic actin filaments and gelsolin segment-1

The assembly of protein actin into double-helical filaments promotes many eukaryotic cellular processes that are regulated by actin-binding proteins (ABPs). Actin filaments can adopt multiple conformations, known as structural polymorphism, which possibly influences the interaction between filaments and ABPs. Gelsolin is a Ca²⁺-regulated ABP that severs and caps actin filaments. Gelsolin binding modulates filament structure; however, it is not known how polymorphic actin filament structures influence an interaction of gelsolin S1 with the barbed-end of filament. Herein, we investigated how polymorphic structures of actin filaments affect the interactions near interfaces between the gelsolin segment 1 (S1) domain and the filament barbed-end. Using all-atom molecular dynamics simulations, we demonstrate that different tilted states of subunits modulate gelsolin S1 interactions with the barbed-end of polymorphic filaments. Hydrogen bonding and interaction energy at the filament-gelsolin S1 interface indicate distinct conformations of filament barbed ends, resulting in different interactions of gelsolin S1. This study demonstrates that filament's structural multiplicity plays important roles in the interactions of actin with ABPs.     

Assessing the flexibility of intermediate filaments by atomic force microscopy

Eukaryotic cells contain three cytoskeletal filament systems that exhibit very distinct assembly properties, supramolecular architectures, dynamic behaviour and mechanical properties. Microtubules and microfilaments are relatively stiff polar structures whose assembly is modulated by the state of hydrolysis of the bound nucleotide. In contrast, intermediate filaments (IFs) are more flexible apolar structures assembled from a approximately 45 nm long coiled-coil dimer as the elementary building block. The differences in flexibility that exist among the three filament systems have been described qualitatively by comparing electron micrographs of negatively stained dehydrated filaments and by directly measuring the persistence length of F-actin filaments (approximately 3-10 microm) and microtubules (approximately 1-8 mm) by various physical methods. However, quantitative data on the persistence length of IFs are still missing. Toward this goal, we have carried out atomic force microscopy (AFM) in physiological buffer to characterise the morphology of individual vimentin IFs adsorbed to different solid supports. In addition, we compared these images with those obtained by transmission electron microscopy (TEM) of negatively stained dehydrated filaments. For each support, we could accurately measure the apparent persistence length of the filaments, yielding values ranging between 0.3 microm and 1 microm. Making simple assumptions concerning the adsorption mechanism, we could estimate the persistence length of an IF in a dilute solution to be approximately 1 microm, indicating that the lower measured values reflect constraints induced by the adsorption process of the filaments on the corresponding support. Based on our knowledge of the structural organisation and mechanical properties of IFs, we reason that the lower persistence length of IFs compared to that of F-actin filaments is caused by the presence of flexible linker regions within the coiled-coil dimer and by postulating the occurrence of axial slipping between dimers within IFs.     

Steered molecular dynamics studies of titin I1 domain unfolding

The cardiac muscle protein titin, responsible for developing passive elasticity and extensibility of muscle, possesses about 40 immunoglobulin-like (Ig) domains in its I-band region. Atomic force microscopy (AFM) and steered molecular dynamics (SMD) have been successfully combined to investigate the reversible unfolding of individual Ig domains. However, previous SMD studies of titin I-band modules have been restricted to I27, the only structurally known Ig domain from the distal region of the titin I-band. In this paper we report SMD simulations unfolding I1, the first structurally available Ig domain from the proximal region of the titin I-band. The simulations are carried out with a view toward upcoming atomic force microscopy experiments. Both constant velocity and constant force stretching have been employed to model mechanical unfolding of oxidized I1, which has a disulfide bond bridging beta-strands C and E, as well as reduced I1, in which the disulfide bridge is absent. The simulations reveal that I1 is protected against external stress mainly through six interstrand hydrogen bonds between its A and B beta-strands. The disulfide bond enhances the mechanical stability of oxidized I1 domains by restricting the rupture of backbone hydrogen bonds between the A'- and G-strands. The disulfide bond also limits the maximum extension of I1 to approximately 220 A. Comparison of the unfolding pathways of I1 and I27 are provided and implications to AFM experiments are discussed.     

Structure of RecX protein complex with the presynaptic RecA filament: Molecular dynamics simulations and small angle neutron scattering

Using molecular modeling techniques we have built the full atomic structure and performed molecular dynamics simulations for the complexes formed by Escherichia coli RecX protein with a single-stranded oligonucleotide and with RecA presynaptic filament. Based on the modeling and SANS experimental data a sandwich-like filament structure formed two chains of RecX monomers bound to the opposite sides of the single stranded DNA is proposed for RecX::ssDNA complex. The model for RecX::RecA::ssDNA include RecX binding into the grove of RecA::ssDNA filament that occurs mainly via Coulomb interactions between RecX and ssDNA. Formation of RecX::RecA::ssDNA filaments in solution was confirmed by SANS measurements which were in agreement with the spectra computed from the molecular dynamics simulations.     

Three-dimensional reconstruction of tarantula myosin filaments suggests how phosphorylation may regulate myosin activity

Muscle contraction involves the interaction of the myosin heads of the thick filaments with actin subunits of the thin filaments. Relaxation occurs when this interaction is blocked by molecular switches on these filaments. In many muscles, myosin-linked regulation involves phosphorylation of the myosin regulatory light chains (RLCs). Electron microscopy of vertebrate smooth muscle myosin molecules (regulated by phosphorylation) has provided insight into the relaxed structure, revealing that myosin is switched off by intramolecular interactions between its two heads, the free head and the blocked head. Three-dimensional reconstruction of frozen-hydrated specimens revealed that this asymmetric head interaction is also present in native thick filaments of tarantula striated muscle. Our goal in this study was to elucidate the structural features of the tarantula filament involved in phosphorylation-based regulation. A new reconstruction revealed intra- and intermolecular myosin interactions in addition to those seen previously. To help interpret the interactions, we sequenced the tarantula RLC and fitted an atomic model of the myosin head that included the predicted RLC atomic structure and an S2 (subfragment 2) crystal structure to the reconstruction. The fitting suggests one intramolecular interaction, between the cardiomyopathy loop of the free head and its own S2, and two intermolecular interactions, between the cardiac loop of the free head and the essential light chain of the blocked head and between the Leu305-Gln327 interaction loop of the free head and the N-terminal fragment of the RLC of the blocked head. These interactions, added to those previously described, would help switch off the thick filament. Molecular dynamics simulations suggest how phosphorylation could increase the helical content of the RLC N-terminus, weakening these interactions, thus releasing both heads and activating the thick filament.     

Side chain interactions determine the amyloid organization: a single layer beta-sheet molecular structure of the calcitonin peptide segment 15-19

In this paper we present a detailed atomic model for a protofilament, the most basic organization level, of the amyloid fibre formed by the peptide DFNKF. This pentapeptide is a segment derived from the human calcitonin, a natural amyloidogenic protein. Our model, which represents the outcome of extensive explicit solvent molecular dynamics (MD) simulations of different strand/sheet organizations, is a single beta-sheet filament largely without a hydrophobic core. Nevertheless, this structure is capable of reproducing the main features of the characteristic amyloid fibril organization and provides clues to the molecular basis of its experimental aggregation behaviour. Our results show that the side chains' chemical diversity induces the formation of a complex network of interactions that finally determine the microscopic arrangement of the strands at the protofilament level. This network of interactions, consisting of both side chain-side chain and backbone-side chain interactions, confers on the final single beta-sheet arrangement an unexpected stability, both by enhancing the association of related chemical groups and, at the same time, by shielding the hydrophobic segments from the polar solvent. The chemical physical characterization of this protofilament provides hints to the possible thermodynamical basis of the supra molecular organization that allows the formation of the filaments by lateral association of the preformed protofibrils. Its regular, highly polarized structure shows how other protofilaments can assemble. In terms of structural biology, our results clearly indicate that an amyloid organization implies a degree of complexity far beyond a simple nonspecific association of peptide strands via amide hydrogen bonds.     

Coarse-Graining Provides Insights on the Essential Nature of Heterogeneity in Actin Filaments

Experiments have shown that actin is structurally polymorphic, but knowledge of the details of molecular level heterogeneity in both the dynamics of a single subunit and the interactions between subunits is still lacking. Here, using atomistic molecular dynamics simulations of the actin filament, we identify domains of atoms that move in a correlated fashion, quantify interactions between these domains using coarse-grained (CG) analysis methods, and perform CG simulations to explore the importance of filament heterogeneity. The persistence length and torsional stiffness calculated from molecular dynamics simulation data agree with experimental values. We additionally observe that distinct actin conformations coexist in actin filaments. The filaments also exhibit random twist angles that are broadly distributed. CG analysis reveals that interactions between equivalent CG pairs vary from one subunit to another. To explore the importance of heterogeneity on filament dynamics, we perform CG simulations using different methods of parameterization to show that only by including heterogeneous interactions can we reproduce the twist angles and related properties. Free energy calculations further suggest that in general the actin filament is best represented as a set of subunits with differing CG sites and interactions, and the incorporating heterogeneity into the CG interactions is more important than including that in the CG sites. Our work therefore presents a systematic method to explore molecular level detail in this large and complex biopolymer.     

Determination of the three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata: a study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata has been determined on the basis of 489 interproton and 24 hydrogen-bonding distance restraints supplemented by 23 phi backbone and 21 chi 1 side-chain torsion angle restraints derived from nuclear magnetic resonance (NMR) measurements. A total of 42 structures is calculated by a hybrid metric matrix distance geometry-dynamical simulated annealing approach. Both the backbone and side-chain atom positions are well defined. The average atomic rms difference between the 42 individual SA structures and the mean structure obtained by averaging their coordinates is 0.67 +/- 0.12 A for the backbone atoms and 0.90 +/- 0.17 A for all atoms. The core of the protein is formed by a triple-stranded antiparallel beta-sheet composed of residues 14-16 (strand 1), 30-34 (strand 2), and 37-41 (strand 3) with an additional mini-antiparallel beta-sheet at the N-terminus (residues 6-9). The first and second strands of the triple-stranded antiparallel beta-sheet are connected by a long exposed loop (residues 17-30). A number of side-chain interactions are discussed in light of the structure.     

Multiscale modeling of cellular actin filaments: from atomistic molecular to coarse-grained dynamics

In this article, we present a computational multiscale model for the characterization of subcellular proteins. The model is encoded inside a simulation tool that builds coarse-grained (CG) force fields from atomistic simulations. Equilibrium molecular dynamics simulations on an all-atom model of the actin filament are performed. Then, using the statistical distribution of the distances between pairs of selected groups of atoms at the output of the MD simulations, the force field is parameterized using the Boltzmann inversion approach. This CG force field is further used to characterize the dynamics of the protein via Brownian dynamics simulations. This combination of methods into a single computational tool flow enables the simulation of actin filaments with length up to 400 nm, extending the time and length scales compared to state-of-the-art approaches. Moreover, the proposed multiscale modeling approach allows to investigate the relationship between atomistic structure and changes on the overall dynamics and mechanics of the filament and can be easily (i) extended to the characterization of other subcellular structures and (ii) used to investigate the cellular effects of molecular alterations due to pathological conditions.     

Quantitative atomic force microscopy image analysis of unusual filaments formed by the Acanthamoeba castellanii myosin II rod domain

We describe a quantitative analysis of Acanthamoeba castellanii myosin II rod domain images collected from atomic force microscope experiments. These images reveal that the rod domain forms a novel filament structure, most likely requiring unusual head-to-tail interactions. Similar filaments are seen also in negatively stained electron microscopy images. Truncated myosins from Acanthamoeba and other model organisms have been visualized before, revealing laterally associated bipolar minifilaments. In contrast, the filament structures that we observe are dominated by axial rather than lateral polymerization. The unusually small features in this structure (1-5 nm) required the development of quantitative and statistical techniques for filament image analysis. These techniques enhance the extraction of features that hitherto have been difficult to ascertain from more qualitative imaging approaches. The heights of the filaments are observed to have a bimodal distribution consistent with the diameters of a single rod domain and a pair of close-packed rod domains. Further quantitative analysis indicates that in-plane association is limited to at most a pair of rod domains. Taken together, this implies that the filaments contain no more than four rod domains laterally associated with one another, somewhat less than that seen in bipolar minifilaments. Analysis of images of the filaments decorated with an anti-FLAG antibody reveals head-to-tail association with mean distances between the antibodies of 75 +/- 15 nm. We consider a set of molecular models to help interpret possible structures of the filaments.