States and transitions during forced unfolding of a single spectrin repeat

Conditional nmy-1 and nmy-2 alleles establish that nonmuscle myosins are required for late Caenorhabditis elegans embryonic elongation

The elongation of Caenorhabditis elegans embryos allows examination of mechanical interactions between adjacent tissues. Muscle contractions during late elongation induce the remodeling of epidermal circumferential actin filaments through mechanotransduction. Force inputs from the muscles deform circumferential epidermal actin filament, which causes them to be severed, eventually reformed, and shortened. This squeezing force drives embryonic elongation. We investigated the possible role of the nonmuscle myosins NMY-1 and NMY-2 in this process using nmy-1 and nmy-2 thermosensitive alleles. Our findings show these myosins act redundantly in late elongation, since double nmy-2(ts); nmy-1(ts) mutants immediately stop elongation when raised to 25°C. Their inactivation does not reduce muscle activity, as measured from epidermis deformation, suggesting that they are directly involved in the multistep process of epidermal remodeling. Furthermore, NMY-1 and NMY-2 inactivation is reversible when embryos are kept at the nonpermissive temperature for a few hours. However, after longer exposure to 25°C double mutant embryos fail to resume elongation, presumably because NMY-1 was seen to form protein aggregates. We propose that the two C. elegans nonmuscle myosin II act during actin remodeling either to bring severed ends or hold them.

Folding-unfolding asymmetry and a RetroFold computational algorithm

We treat protein folding as molecular self-assembly, while unfolding is viewed as disassembly. Fracture is typically a much faster process than self-assembly. Self-assembly is often an exponentially decaying process, since energy relaxes due to dissipation, while fracture is a constant-rate process as the driving force is opposed by damping. Protein folding takes two orders of magnitude longer than unfolding. We suggest a mathematical transformation of variables, which makes it possible to view self-assembly as time-reversed disassembly, thus folding can be studied as reversed unfolding. We investigate the molecular dynamics modelling of folding and unfolding of the short Trp-cage protein. Folding time constitutes about 800 ns, while unfolding (denaturation) takes only about 5.0 ns and, therefore, fewer computational resources are needed for its simulation. This RetroFold approach can be used for the design of a novel computation algorithm, which, while approximate, is less time-consuming than traditional folding algorithms.

The spectrin cytoskeleton integrates endothelial mechanoresponses

Physiological blood flow induces the secretion of vasoactive compounds, notably nitric oxide, and promotes endothelial cell elongation and reorientation parallel to the direction of applied shear. How shear is sensed and relayed to intracellular effectors is incompletely understood. Here, we demonstrate that an apical spectrin network is essential to convey the force imposed by shear to endothelial mechanosensors. By anchoring CD44, spectrins modulate the cell surface density of hyaluronan and sense and translate shear into changes in plasma membrane tension. Spectrins also regulate the stability of apical caveolae, where the mechanosensitive PIEZO1 channels are thought to reside. Accordingly, shear-induced PIEZO1 activation and the associated calcium influx were absent in spectrin-deficient cells. As a result, cell realignment and flow-induced endothelial nitric oxide synthase stimulation were similarly dependent on spectrin. We conclude that the apical spectrin network is not only required for shear sensing but also transmits and distributes the resulting tensile forces to mechanosensors that elicit protective and vasoactive responses.

Mechanical behavior of actin and spectrin subjected to high strain rate: A molecular dynamics simulation study

Recent nanoscopy and super-resolution microscopy studies have substantiated the structural contribution of periodic actin-spectrin lattice to the axonal cytoskeleton of neuron. However, sufficient mechanical insight is not present for spectrin and actin-spectrin network, especially in high strain rate scenario. To quantify the mechanical behavior of actin-spectrin cytoskeleton in such conditions, this study determines individual stretching characteristics of actin and spectrin at high strain rate by molecular dynamics (MD) simulation. The actin-spectrin separation criteria are also determined. It is found that both actin and spectrin have high stiffness when susceptible to high strain rate and show strong dependence on applied strain rate. The stretching stiffness of actin and forced unfolding mechanism of spectrin are in harmony with the current literature. Actin-spectrin model provides novel insight into their interaction and separation stretch. It is shown that the region vulnerable to failure is the actin-spectrin interface at lower strain rate, while it is the inter-repeat region of spectrin at higher strain rate.

Manipulation of Focal Adhesion Signaling by Pathogenic Microbes

Focal adhesions (FAs) serve as dynamic signaling hubs within the cell. They connect intracellular actin to the extracellular matrix (ECM) and respond to environmental cues. In doing so, these structures facilitate important processes such as cell-ECM adhesion and migration. Pathogenic microbes often modify the host cell actin cytoskeleton in their pursuit of an ideal replicative niche or during invasion to facilitate uptake. As actin-interfacing structures, FA dynamics are also intimately tied to actin cytoskeletal organization. Indeed, exploitation of FAs is another avenue by which pathogenic microbes ensure their uptake, survival and dissemination. This is often achieved through the secretion of effector proteins which target specific protein components within the FA. Molecular mimicry of the leucine-aspartic acid (LD) motif or vinculin-binding domains (VBDs) commonly found within FA proteins is a common microbial strategy. Other effectors may induce post-translational modifications to FA proteins through the regulation of phosphorylation sites or proteolytic cleavage. In this review, we present an overview of the regulatory mechanisms governing host cell FAs, and provide examples of how pathogenic microbes have evolved to co-opt them to their own advantage. Recent technological advances pose exciting opportunities for delving deeper into the mechanistic details by which pathogenic microbes modify FAs.

Biological physics by high-speed atomic force microscopy

While many fields have contributed to biological physics, nanotechnology offers a new scale of observation. High-speed atomic force microscopy (HS-AFM) provides nanometre structural information and dynamics with subsecond resolution of biological systems. Moreover, HS-AFM allows us to measure piconewton forces within microseconds giving access to unexplored, fast biophysical processes. Thus, HS-AFM provides a tool to nourish biological physics through the observation of emergent physical phenomena in biological systems. In this review, we present an overview of the contribution of HS-AFM, both in imaging and force spectroscopy modes, to the field of biological physics. We focus on examples in which HS-AFM observations on membrane remodelling, molecular motors or the unfolding of proteins have stimulated the development of novel theories or the emergence of new concepts. We finally provide expected applications and developments of HS-AFM that we believe will continue contributing to our understanding of nature, by serving to the dialogue between biology and physics. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.

Cotranslational folding cooperativity of contiguous domains of α-spectrin

Proteins synthesized in the cell can begin to fold during translation before the entire polypeptide has been produced, which may be particularly relevant to the folding of multidomain proteins. Here, we study the cotranslational folding of adjacent domains from the cytoskeletal protein α-spectrin using force profile analysis (FPA). Specifically, we investigate how the cotranslational folding behavior of the R15 and R16 domains are affected by their neighboring R14 and R16, and R15 and R17 domains, respectively. Our results show that the domains impact each other's folding in distinct ways that may be important for the efficient assembly of α-spectrin, and may reduce its dependence on chaperones. Furthermore, we directly relate the experimentally observed yield of full-length protein in the FPA assay to the force exerted by the folding protein in piconewtons. By combining pulse-chase experiments to measure the rate at which the arrested protein is converted into full-length protein with a Bell model of force-induced rupture, we estimate that the R16 domain exerts a maximal force on the nascent chain of ∼15 pN during cotranslational folding.

The axonal actin-spectrin lattice acts as a tension buffering shock absorber

Axons span extreme distances and are subject to significant stretch deformations during limb movements or sudden head movements, especially during impacts. Yet, axon biomechanics, and its relation to the ultrastructure that allows axons to withstand mechanical stress, is poorly understood. Using a custom developed force apparatus, we demonstrate that chick dorsal root ganglion axons exhibit a tension buffering or strain-softening response, where its steady state elastic modulus decreases with increasing strain. We then explore the contributions from the various cytoskeletal components of the axon to show that the recently discovered membrane-associated actin-spectrin scaffold plays a prominent mechanical role. Finally, using a theoretical model, we argue that the actin-spectrin skeleton acts as an axonal tension buffer by reversibly unfolding repeat domains of the spectrin tetramers to release excess mechanical stress. Our results revise the current viewpoint that microtubules and their associated proteins are the only significant load-bearing elements in axons.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

The plakin domain of C. elegans VAB-10/plectin acts as a hub in a mechanotransduction pathway to promote morphogenesis

Mechanical forces can elicit a mechanotransduction response through junction-associated proteins. In contrast to the wealth of knowledge available for focal adhesions and adherens junctions, much less is known about mechanotransduction at hemidesmosomes. Here, we focus on the C. elegans plectin homolog VAB-10A, the only evolutionary conserved hemidesmosome component. In C. elegans, muscle contractions induce a mechanotransduction pathway in the epidermis through hemidesmosomes. We used CRISPR to precisely remove spectrin repeats (SRs) or a partially hidden Src homology 3 (SH3) domain within the VAB-10 plakin domain. Deleting the SH3 or SR8 domains in combination with mutations affecting mechanotransduction, or just the part of SR5 shielding the SH3 domain, induced embryonic elongation arrest because hemidesmosomes collapse. Notably, recruitment of GIT-1, the first mechanotransduction player, requires the SR5 domain and the hemidesmosome transmembrane receptor LET-805. Furthermore, molecular dynamics simulations confirmed that forces acting on VAB-10 could make the central SH3 domain, otherwise in contact with SR4, available for interaction. Collectively, our data strongly indicate that the plakin domain plays a central role in mechanotransduction and raise the possibility that VAB-10/plectin might act as a mechanosensor.

Distinct mechanical properties in homologous spectrin-like repeats of utrophin

Patients with Duchenne muscular dystrophy (DMD) lack the protein dystrophin, which is a critical molecular component of the dystrophin-glycoprotein complex (DGC). Dystrophin is hypothesized to function as a molecular shock absorber that mechanically stabilizes the sarcolemma of striated muscle through interaction with the cortical actin cytoskeleton via its N-terminal half and with the transmembrane protein β-dystroglycan via its C-terminal region. Utrophin is a fetal homologue of dystrophin that can subserve many dystrophin functions and is therefore under active investigation as a dystrophin replacement therapy for DMD. Here, we report the first mechanical characterization of utrophin using atomic force microscopy (AFM). Our data indicate that the mechanical properties of spectrin-like repeats in utrophin are more in line with the PEVK and Ig-like repeats of titin rather than those reported for repeats in spectrin or dystrophin. Moreover, we measured markedly different unfolding characteristics for spectrin repeats within the N-terminal actin-binding half of utrophin compared to those in the C-terminal dystroglycan-binding half, even though they exhibit identical thermal denaturation profiles. Our results demonstrate dramatic differences in the mechanical properties of structurally homologous utrophin constructs and suggest that utrophin may function as a stiff elastic element in series with titin at the myotendinous junction.

α-Helix Unwinding as Force Buffer in Spectrins

Spectrins are cytoskeletal proteins located at the inner face of the plasma membrane, making connections between membrane anchors and the actin cortex, and between actin filaments. Spectrins share a common structure forming a bundle of 3 α-helices and play a major role during cell deformation. Here, we used high-speed force spectroscopy and steered molecular dynamics simulations to understand the mechanical stability of spectrin, revealing a molecular force buffering function. We find that spectrin acts as a soft spring at short extensions (70-100 Å). Under continuous external stretching, its α-helices unwind, leading to a viscous mechanical response over larger extensions (100-300 Å), represented by a constant-force plateau in force/extension curves. This viscous force buffering emerges from a quasi-equilibrium competition between disruption and re-formation of α-helical hydrogen bonds. Our results suggest that, in contrast to β-sheet proteins, which unfold in a catastrophic event, α-helical spectrins dominantly unwind, providing a viscous force buffer over extensions about 5 times their folded length.

Nuclear Envelope-Associated Chromosome Dynamics during Meiotic Prophase I

Chromosome dynamics during meiotic prophase I are associated with a series of major events such as chromosomal reorganization and condensation, pairing/synapsis and recombination of the homologs, and chromosome movements at the nuclear envelope (NE). The NE is the barrier separating the nucleus from the cytoplasm and thus plays a central role in NE-associated chromosomal movements during meiosis. Previous studies have shown in various species that NE-linked chromosome dynamics are actually driven by the cytoskeleton. The linker of nucleoskeleton and cytoskeleton (LINC) complexes are important constituents of the NE that facilitate in the transfer of cytoskeletal forces across the NE to individual chromosomes. The LINCs consist of the inner and outer NE proteins Sad1/UNC-84 (SUN), and Klarsicht/Anc-1/Syne (KASH) domain proteins. Meiosis-specific adaptations of the LINC components and unique modifications of the NE are required during chromosomal movements. Nonetheless, the actual role of the NE in chromosomic dynamic movements in plants remains elusive. This review summarizes the findings of recent studies on meiosis-specific constituents and modifications of the NE and corresponding nucleoplasmic/cytoplasmic adaptors being involved in NE-associated movement of meiotic chromosomes, as well as describes the potential molecular network of transferring cytoplasm-derived forces into meiotic chromosomes in model organisms. It helps to gain a better understanding of the NE-associated meiotic chromosomal movements in plants.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 Å·ps^-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

Protein unfolding mechanisms and their effects on folding experiments

In this review, I discuss the various methods researchers use to unfold proteins in the lab in order to understand protein folding both in vitro and in vivo. The four main techniques, chemical-, heat-, pressure- and force-denaturation, produce distinctly different unfolded conformational ensembles. Recent measurements have revealed different folding kinetics from different unfolding mechanisms. Thus, comparing these distinct unfolded ensembles sheds light on the underlying free energy landscape of folding.

The nesprin-cytoskeleton interface probed directly on single nuclei is a mechanically rich system

The cytoskeleton provides structure and plays an important role in cellular function such as migration, resisting compression forces, and transport. The cytoskeleton also reacts to physical cues such as fluid shear stress or extracellular matrix remodeling by reorganizing filament associations, most commonly focal adhesions and cell-cell cadherin junctions. These mechanical stimuli can result in genome-level changes, and the physical connection of the cytoskeleton to the nucleus provides an optimal conduit for signal transduction by interfacing with nuclear envelope proteins, called nesprins, within the LINC (linker of the nucleus to the cytoskeleton) complex. Using single-molecule on single nuclei assays, we report that the interactions between the nucleus and the cytoskeleton, thought to be nesprin-cytoskeleton interactions, are highly sensitive to force magnitude and direction depending on whether cells are historically interfaced with the matrix or with cell aggregates. Application of ∼10-30 pN forces to these nesprin linkages yielded structural transitions, with a base transition size of 5-6 nm, which are speculated to be associated with partial unfoldings of the spectrin domains of the nesprins and/or structural changes of histones within the nucleus.

Determination of contact maps in proteins: A combination of structural and chemical approaches

Contact map selection is a crucial step in structure-based molecular dynamics modelling of proteins. The map can be determined in many different ways. We focus on the methods in which residues are represented as clusters of effective spheres. One contact map, denoted as overlap (OV), is based on the overlap of such spheres. Another contact map, named Contacts of Structural Units (CSU), involves the geometry in a different way and, in addition, brings chemical considerations into account. We develop a variant of the CSU approach in which we also incorporate Coulombic effects such as formation of the ionic bridges and destabilization of possible links through repulsion. In this way, the most essential and well defined contacts are identified. The resulting residue-residue contact map, dubbed repulsive CSU (rCSU), is more sound in its physico-chemical justification than CSU. It also provides a clear prescription for validity of an inter-residual contact: the number of attractive atomic contacts should be larger than the number of repulsive ones - a feature that is not present in CSU. However, both of these maps do not correlate well with the experimental data on protein stretching. Thus, we propose to use rCSU together with the OV map. We find that the combined map, denoted as OV+rCSU, performs better than OV. In most situations, OV and OV+rCSU yield comparable folding properties but for some proteins rCSU provides contacts which improve folding in a substantial way. We discuss the likely residue-specificity of the rCSU contacts. Finally, we make comparisons to the recently proposed shadow contact map, which is derived from different principles.

Spectroscopic Monitoring of Mechanical Forces during Protein Folding by using Molecular Force Probes

Detailed folding pathways of proteins are still largely unknown. Real-time monitoring of mechanical forces acting in proteins during structural transitions would provide deep insights into these highly complex processes. Here, we propose two molecular force probes that can be incorporated into the protein backbone to gain insight into the magnitude and direction of mechanical forces acting in proteins during natural folding and unfolding through their optical spectroscopic response. In fact, changes in the infrared and Raman spectra are proportional to the mechanical force deforming the force probes, and the relevant bands can be intensified and shifted to a transparent window in the protein spectrum by isotopic substitution. As a result, the proposed molecular force probes can act as "force rulers", allowing the spectroscopic observation and measurement of mechanical forces acting within the proteins under natural conditions without external perturbation.

Synaptopodin couples epithelial contractility to α-actinin-4-dependent junction maturation

A novel tension-sensitive junctional protein, synaptopodin, can relay biophysical input from cellular actomyosin contractility to induce biochemical changes at cell–cell contacts, resulting in structural reorganization of the junctional complex and epithelial barrier maturation.

The epithelial junction experiences mechanical force exerted by endogenous actomyosin activities and from interactions with neighboring cells. We hypothesize that tension generated at cell–cell adhesive contacts contributes to the maturation and assembly of the junctional complex. To test our hypothesis, we used a hydraulic apparatus that can apply mechanical force to intercellular junction in a confluent monolayer of cells. We found that mechanical force induces α-actinin-4 and actin accumulation at the cell junction in a time- and tension-dependent manner during junction development. Intercellular tension also induces α-actinin-4–dependent recruitment of vinculin to the cell junction. In addition, we have identified a tension-sensitive upstream regulator of α-actinin-4 as synaptopodin. Synaptopodin forms a complex containing α-actinin-4 and β-catenin and interacts with myosin II, indicating that it can physically link adhesion molecules to the cellular contractile apparatus. Synaptopodin depletion prevents junctional accumulation of α-actinin-4, vinculin, and actin. Knockdown of synaptopodin and α-actinin-4 decreases the strength of cell–cell adhesion, reduces the monolayer permeability barrier, and compromises cellular contractility. Our findings underscore the complexity of junction development and implicate a control process via tension-induced sequential incorporation of junctional components.

Mechanical stability analysis of the protein L immunoglobulin-binding domain by full alanine screening using molecular dynamics simulations

This article is the first to study the mechanical properties of the immunoglobulin-binding domain of protein L (referred to as protein L) and its mutants at the atomic level. In the structure of protein L, each amino acid residue (except for alanines and glycines) was replaced sequentially by alanine. Thus, 49 mutants of protein L were obtained. The proteins were stretched at their termini at constant velocity using molecular dynamics simulations in water, i.e. by forced unfolding. 19 out of 49 mutations resulted in a large decrease of mechanical protein stability. These amino acids were affecting either the secondary structure (11 mutations) or loop structures (8 mutations) of protein L. Analysis of mechanical unfolding of the generated protein that has the same topology as protein L but consists of only alanines and glycines allows us to suggest that the mechanical stability of proteins, and specifically protein L, is determined by interactions between certain amino acid residues, although the unfolding pathway depends on the protein topology. This insight can now be used to modulate the mechanical properties of proteins and their unfolding pathways in the desired direction for using them in various biochips, biosensors and biomaterials for medicine, industry, and household purposes.

Experimental and theoretical studies of mechanical unfolding of different proteins

Mechanical properties of proteins are important for a wide range of biological processes including cell adhesion, muscle contraction, and protein translocation across biological membranes. It is necessary to reveal how proteins achieve their required mechanical stability under natural conditions in order to understand the biological processes and also to use the knowledge for constructing new biomaterials for medical and industrial purposes. In this connection, it is important to know how a protein will behave in response to various impacts. Theoretical and experimental works on mechanical unfolding of globular proteins will be considered in detail in this review.

Protein interactions at the higher plant nuclear envelope: evidence for a linker of nucleoskeleton and cytoskeleton complex

Following the description of SAD1/UNC84 (SUN) domain proteins in higher plants, evidence has rapidly increased that plants contain a functional linker of nucleoskeleton and cytoskeleton (LINC) complex bridging the nuclear envelope (NE). While the SUN domain proteins appear to be highly conserved across kingdoms, other elements of the complex are not and some key components and interactions remain to be identified. This mini review examines components of the LINC complex, including proteins of the SUN domain family and recently identified plant Klarsicht/Anc/Syne-1 homology (KASH) domain proteins. First of these to be described were WIPs (WPP domain interacting proteins), which act as protein anchors in the outer NE. The plant KASH homologs are C-terminally anchored membrane proteins with the extreme C-terminus located in the nuclear periplasm; AtWIPs contain a highly conserved X-VPT motif at the C-terminus in contrast to PPPX in opisthokonts. The role of the LINC complex in organisms with a cell wall, and description of further LINC complex components will be considered, together with other potential plant-specific functions.

Demonstration of specific binding of heparin to Plasmodium falciparum-infected vs. non-infected red blood cells by single-molecule force spectroscopy

Glycosaminoglycans (GAGs) play an important role in the sequestration of Plasmodium falciparum-infected red blood cells (pRBCs) in the microvascular endothelium of different tissues, as well as in the formation of small clusters (rosettes) between infected and non-infected red blood cells (RBCs). Both sequestration and rosetting have been recognized as characteristic events in severe malaria. Here we have used heparin and pRBCs infected by the 3D7 strain of P. falciparum as a model to study GAG-pRBC interactions. Fluorescence microscopy and fluorescence-assisted cell sorting assays have shown that exogenously added heparin has binding specificity for pRBCs (preferentially for those infected with late forms of the parasite) vs. RBCs. Heparin-pRBC adhesion has been probed by single-molecule force spectroscopy, obtaining an average binding force ranging between 28 and 46 pN depending on the loading rate. No significant binding of heparin to non-infected RBCs has been observed in control experiments. This work represents the first approach to quantitatively evaluate GAG-pRBC molecular interactions at the individual molecule level.

Mechanical strength of 17,134 model proteins and cysteine slipknots

A new theoretical survey of proteins' resistance to constant speed stretching is performed for a set of 17 134 proteins as described by a structure-based model. The proteins selected have no gaps in their structure determination and consist of no more than 250 amino acids. Our previous studies have dealt with 7510 proteins of no more than 150 amino acids. The proteins are ranked according to the strength of the resistance. Most of the predicted top-strength proteins have not yet been studied experimentally. Architectures and folds which are likely to yield large forces are identified. New types of potent force clamps are discovered. They involve disulphide bridges and, in particular, cysteine slipknots. An effective energy parameter of the model is estimated by comparing the theoretical data on characteristic forces to the corresponding experimental values combined with an extrapolation of the theoretical data to the experimental pulling speeds. These studies provide guidance for future experiments on single molecule manipulation and should lead to selection of proteins for applications. A new class of proteins, involving cystein slipknots, is identified as one that is expected to lead to the strongest force clamps known. This class is characterized through molecular dynamics simulations.

The advances in nanotechnology have allowed for manipulation of single biomolecules and determination of their elastic properties. Titin was among the first proteins studied in this way. Its unravelling by stretching requires a 204 pN force. The resistance to stretching comes mostly from a localized region known as a force clamp. In titin, the force clamp is simple as it is formed by two parallel β-strands that are sheared on pulling. Studies of a set of under a hundred proteins accomplished in the last decade have revealed a variety of the force clamps that lead to forces ranging from under 20 pN to about 500 pN. This set comprises only a tiny fraction of proteins known. Thus one needs guidance as to what proteins should be considered for specific mechanical properties. Such a guidance is provided here through simulations within simplified coarse-grained models on 17 134 proteins that are stretched at constant speed. We correlate their unravelling forces with two structure classification schemes. We identify proteins with large resistance to unravelling and characterize their force clamps. Quite a few top strength proteins owe their sturdiness to a new type of the force clamp: the cystein slipknot in which the force peak is due to dragging of a piece of the backbone through a closed ring formed by two other pieces of the backbone and two connecting disulphide bonds.

Bringing KASH under the SUN: the many faces of nucleo-cytoskeletal connections

The nucleus is the most prominent cellular organelle, and its sharp boundaries suggest the compartmentalization of the nucleoplasm from the cytoplasm. However, the recent identification of evolutionarily conserved linkers of the nucleoskeleton to the cytoskeleton (LINC) complexes, a family of macromolecular assemblies that span the double membrane of the nuclear envelope, reveals tight physical connections between the two compartments. Here, we review the structure and evolutionary conservation of SUN and KASH domain–containing proteins, whose interaction within the perinuclear space forms the “nuts and bolts” of LINC complexes. Moreover, we discuss the function of these complexes in nuclear, centrosomal, and chromosome dynamics, and their connection to human disease.

Mechanical signaling on the single protein level studied using steered molecular dynamics

Efficient communication between the cell and its external environment is of the utmost importance to the function of multicellular organisms. While signaling events can be generally characterized as information exchange by means of controlled energy conversion, research efforts have hitherto mainly been concerned with mechanisms involving chemical and electrical energy transfer. Here, we review recent computational efforts addressing the function of mechanical force in signal transduction. Specifically, we focus on the role of steered molecular dynamics (SMD) simulations in providing details at the atomic level on a group of protein domains, which play a fundamental role in signal exchange by responding properly to mechanical strain. We start by giving a brief introduction to the SMD technique and general properties of mechanically stable protein folds, followed by specific examples illustrating three general regimes of signal transfer utilizing mechanical energy: purely mechanical, mechanical to chemical, and chemical to mechanical. Whenever possible the physiological importance of the example at hand is stressed to highlight the diversity of the processes in which mechanical signaling plays a key role. We also provide an overview of future challenges and perspectives for this rapidly developing field.

A multi-timescale strength model of alpha-helical protein domains

Here we report a constitutive model that characterizes the strength of an alpha-helical protein domain subjected to tensile deformation, covering more than ten orders of magnitude in timescales. The model elucidates multiple physical mechanisms of failure in dependence on the deformation rate, quantitatively linking atomistic simulation results with experimental strength measurements of alpha-helical protein domains. The model provides a description of the strength of alpha-helices based on fundamental physical parameters such as the H-bond energy and the polypeptide's persistence length, showing that strength is controlled by energetic, nonequilibrium processes at high rates and by thermodynamical, equilibrium processes at low rates. Our model provides a novel perspective on the strength of protein domains at ultra-slow pulling speeds relevant under physiologic and experimental conditions.

Actomyosin interaction: mechanical and energetic properties in different nucleotide binding states

The mechanics of the actomyosin interaction is central in muscle contraction and intracellular trafficking. A better understanding of the events occurring in the actomyosin complex requires the examination of all nucleotide-dependent states and of the energetic features associated with the dynamics of the cross-bridge cycle. The aim of the present study is to estimate the interaction strength between myosin in nucleotide-free, ATP, ADP·Pi and ADP states and actin monomer. The molecular models of the complexes were constructed based on cryo-electron microscopy maps and the interaction properties were estimated by means of a molecular dynamics approach, which simulate the unbinding of the complex applying a virtual spring to the core of myosin protein. Our results suggest that during an ATP hydrolysis cycle the affinity of myosin for actin is modulated by the presence and nature of the nucleotide in the active site of the myosin motor domain. When performing unbinding simulations with a pulling rate of 0.001 nm/ps, the maximum pulling force applied to the myosin during the experiment is about 1nN. Under these conditions the interaction force between myosin and actin monomer decreases from 0.83 nN in the nucleotide-free state to 0.27 nN in the ATP state, and increases to 0.60 nN after ATP hydrolysis and Pi release from the complex (ADP state).

Interaction forces and interface properties of KIF1A kinesin-alphabeta tubulin complex assessed by molecular dynamics

Kinesin is a microtubule-based motor protein that generates motion involved in intracellular trafficking and cell division. Even if the force-generating and enzymatic properties of kinesin were extensively studied, the molecular basis of its interaction with the microtubule is still not well understood. The aim of the present study is to provide a detailed description, in terms of conformational changes and interaction properties, of the kinesin-alphabeta tubulin complex during a cycle of ATP hydrolysis. Four different nucleotide-dependent conformations (nucleotide-free, ATP, ADP.Pi and ADP) of the kinesin-alphabeta tubulin were constructed and investigated by performing molecular dynamics simulations. Computational results show that small conformational changes, in the order of few Angstrom, occurring in the kinesin structure reflect on its affinity for the filament substrate. Indeed the rotation of the alpha4 helix due to the transition from the bound (ADP.Pi) to the unbound (ADP) state, when the Pi is released from the complex, coupled with the modification occurred in the loop L9 of switch I domain are associated to a marked decrease (approximately 45%) of the maximum interaction force between the kinesin motor and the tubulin dimer.

Selection of optimal variants of Gō-like models of proteins through studies of stretching

The Gō-like models of proteins are constructed based on the knowledge of the native conformation. However, there are many possible choices of a Hamiltonian for which the ground state coincides with the native state. Here, we propose to use experimental data on protein stretching to determine what choices are most adequate physically. This criterion is motivated by the fact that stretching processes usually start with the native structure, in the vicinity of which the Gō-like models should work the best. Our selection procedure is applied to 62 different versions of the Gō model and is based on 28 proteins. We consider different potentials, contact maps, local stiffness energies, and energy scales--uniform and nonuniform. In the latter case, the strength of the nonuniformity was governed either by specificity or by properties related to positioning of the side groups. Among them is the simplest variant: uniform couplings with no i, i + 2 contacts. This choice also leads to good folding properties in most cases. We elucidate relationship between the local stiffness described by a potential which involves local chirality and the one which involves dihedral and bond angles. The latter stiffness improves folding but there is little difference between them when it comes to stretching.

Conformational equilibria in monomeric alpha-synuclein at the single-molecule level

Human α-Synuclein (αSyn) is a natively unfolded protein whose aggregation into amyloid fibrils is involved in the pathology of Parkinson disease. A full comprehension of the structure and dynamics of early intermediates leading to the aggregated states is an unsolved problem of essential importance to researchers attempting to decipher the molecular mechanisms of αSyn aggregation and formation of fibrils. Traditional bulk techniques used so far to solve this problem point to a direct correlation between αSyn's unique conformational properties and its propensity to aggregate, but these techniques can only provide ensemble-averaged information for monomers and oligomers alike. They therefore cannot characterize the full complexity of the conformational equilibria that trigger the aggregation process. We applied atomic force microscopy–based single-molecule mechanical unfolding methodology to study the conformational equilibrium of human wild-type and mutant αSyn. The conformational heterogeneity of monomeric αSyn was characterized at the single-molecule level. Three main classes of conformations, including disordered and “β-like” structures, were directly observed and quantified without any interference from oligomeric soluble forms. The relative abundance of the “β-like” structures significantly increased in different conditions promoting the aggregation of αSyn: the presence of Cu²⁺, the pathogenic A30P mutation, and high ionic strength. This methodology can explore the full conformational space of a protein at the single-molecule level, detecting even poorly populated conformers and measuring their distribution in a variety of biologically important conditions. To the best of our knowledge, we present for the first time evidence of a conformational equilibrium that controls the population of a specific class of monomeric αSyn conformers, positively correlated with conditions known to promote the formation of aggregates. A new tool is thus made available to test directly the influence of mutations and pharmacological strategies on the conformational equilibrium of monomeric αSyn.

Natively unstructured proteins defy the classical “one sequence–one structure” paradigm of protein science. In pathological conditions, monomers of these proteins can aggregate in the cell, a process that underlies neurodegenerative diseases such as Alzheimer and Parkinson. A key step in the aggregation process—the formation of misfolded intermediates—remains obscure. To shed light on this process, we characterized the folding and conformational diversity of αSyn, a natively unstructured protein involved in Parkinson disease, by mechanically stretching single molecules of this protein and recording their mechanical properties. These experiments permitted us to observe directly and quantify three main classes of conformations that, under in vitro physiological conditions, exist simultaneously in the αSyn sample. We found that one class of conformations, “β-like” structures, is directly related to αSyn aggregation. In fact, their relative abundance increases drastically in three different conditions known to promote the formation of αSyn fibrils. We expect that a critical concentration of αSyn with a “β-like” structure must be reached to trigger fibril formation. This critical concentration is therefore controlled by a chemical equilibrium. Novel pharmacological strategies can now be tailored to act upstream, before the aggregation process ensues, by targeting this equilibrium. To this end, single-molecule force spectroscopy can be an effective tool to tailor and test new pharmacological agents.

A single-molecule study detects structured and unstructured conformers in equilibrium in monomeric α-synuclein. The β-like conformers increase with pathological mutations and under other conditions known to promote aggregation.

Transient violations of the second law of thermodynamics in protein unfolding examined using synthetic atomic force microscopy and the fluctuation theorem

The synthetic atomic force microscopy (AFM) method is developed to simulate a periodically replicated atomistic system subject to force and length fluctuations characteristic of an AFM experiment. This new method is used to examine the forced-extension and subsequent rupture of the alpha-helical linker connecting periodic images of a spectrin protein repeat unit. A two-dimensional potential of mean force (PMF) along the length and a reaction coordinate describing the state of the linker was calculated. This PMF reveals that the basic material properties of the spectrin repeat unit are sensitive to the state of linker, an important feature that cannot be accounted for in a one-dimensional PMF. Furthermore, nonequilibrium simulations were generated to examine the rupture event in the context of the fluctuation theorem. These atomistic simulations demonstrate that trajectories which are in apparent violation of the second law can overcome unfolding barriers at significantly reduced rupture forces.

Stretching to understand proteins - a survey of the protein data bank

We make a survey of resistance of 7510 proteins to mechanical stretching at constant speed as studied within a coarse-grained molecular dynamics model. We correlate the maximum force of resistance with the native structure, predict proteins which should be especially strong, and identify the nature of their force clamps.

The effect of protein complexation on the mechanical stability of Im9

Force mode microscopy can be used to examine the effect of mechanical manipulation on the noncovalent interactions that stabilize proteins and their complexes. Here we describe the effect of complexation by the high affinity protein ligand E9 on the mechanical resistance of the simple four-helical protein, Im9. When concatenated into a construct of alternating I27 domains, Im9 unfolded below the thermal noise limit of the instrument ( approximately 20 pN). Complexation of E9 had little effect on the mechanical resistance of Im9 (unfolding force approximately 30 pN) despite the high avidity of this complex (K(d) approximately 10 fM).

Mechanical unfolding of proteins: insights into biology, structure and folding

Since user-friendly atomic force microscopes came onto the market a few years ago, scientists have explored the response of many proteins to applied force. This field has now matured beyond the phenomenological with exciting recent developments, particularly with regards to research into biological questions. For example, detailed mechanistic studies have suggested how mechanically active proteins perform their functions. Also, in vitro forced unfolding has been compared with in vivo protein import and degradation. Additionally, investigations have been carried out that probe the relationship between protein structure and response to applied force, an area that has benefited significantly from synergy between experiments and simulations. Finally, recent technological developments offer exciting new avenues for experimental studies.

Spectrin domains lose cooperativity in forced unfolding

Spectrin is a multidomain cytoskeletal protein, the component three-helix bundle domains are expected to experience mechanical force in vivo. In thermodynamic and kinetic studies, neighboring domains of chicken brain alpha-spectrin R16 and R17 have been shown to behave cooperatively. Is this cooperativity maintained under force? The effect of force on these spectrin domains was investigated using atomic force microscopy. The response of the individual domains to force was compared to that of the tandem repeat R1617. Importantly, nonhelical linkers (all-beta immunoglobulin domains) were used to avoid formation of nonnative helical linkers. We show that, in contrast to previous studies on spectrin repeats, only 3% of R1617 unfolding events gave an increase in contour length consistent with cooperative two-domain unfolding events. Furthermore, the unfolding forces for R1617 were the same as those for the unfolding of R16 or R17 alone. This is a strong indication that the cooperative unfolding behavior observed in the stopped-flow studies is absent between these spectrin domains when force is acting as a denaturant. Our evidence suggests that the rare double unfolding events result from misfolding between adjacent repeats. We suggest that this switch from cooperative to independent behavior allows multidomain proteins to maintain integrity under applied force.

Mechanical response and conformational changes of alpha-actinin domains during unfolding: a molecular dynamics study

Alpha-actinin is a cytoskeleton-binding protein involved in the assembly and regulation of the actin filaments. In this work molecular dynamics method was applied to investigate the mechanical behaviour of the human skeletal muscle alpha-actinin. Five configurations were unfolded at an elongation speed of 0.1 nm/ps in order to investigate the conformational changes occurring during the extension process. Moreover, a sensitivity analysis at different velocities was performed for one of the R2-R3 spectrin-like repeat configuration extracted in order to evaluate the effect of the pulling speed on the mechanical behaviour of the molecule. Two different behaviours were recognized with respect to the pulling speed. In particular, at speed higher than 0.025 nm/ps a continuous rearrangement without evident force peaks was obtained, on the contrary at lower speed evident peaks in the range 500-750 pN were detected. R3 repeat resulted more stable than R2 during mechanical unfolding, due to the lower hydrophobic surface available to the solvent. The characterization of the R2-R3 units can be useful for the development of cytoskeleton network models based on stiffness values obtained by analyses performed at the molecular level.

Examining the influence of linkers and tertiary structure in the forced unfolding of multiple-repeat spectrin molecules

The unfolding pathways of multiple-repeat spectrin molecules were examined using steered molecular dynamics (SMD) simulations to forcibly unfold double- and triple-repeat spectrin molecules. Although SMD has previously been used to study other repeating-domain proteins, spectrin offers a unique challenge in that the linker connecting repeat units has a definite secondary structure, that of an alpha-helix. Therefore, the boundary conditions imposed on a double- or triple-repeat spectrin must be carefully considered if any relationship to the real system is to be deduced. This was accomplished by imposing additional forces on the system which ensure that the terminal alpha-helices behave as if there were no free noncontiguous helical ends. The results of the SMD simulations highlight the importance of the rupture of the alpha-helical linker on the subsequent unfolding events. Rupture of the linker propagates unfolding in the adjacent repeat units by destabilizing the tertiary structure, ultimately resulting in complete unfolding of the affected repeat unit. Two dominant classes of unfolding pathways are observed after the initial rupture of a linker which involve either rupture of another linker (possibly adjacent) or rupture of the basic tertiary structure of a repeat unit. The relationship between the force response observed on simulation timescales and those of experiment or physiological conditions is also discussed.

Extending a spectrin repeat unit. I: linear force-extension response

Nonequilibrium molecular dynamics simulations were used to calculate the elastic properties of a spectrin repeat unit. A contiguous alpha-helical linker was constructed by employing periodic boundary conditions, allowing a novel scheme for evaluating the thermodynamic force as a function of extension. By measuring the force-extension response under small extensions, spectrin was observed to behave primarily as an elastic material with a spring constant of 1700 +/- 100 pN/nm. The implications of this spring constant, in terms of the properties of the spectrin tetramer, are also discussed.

Extending a spectrin repeat unit. II: rupture behavior

A spectrin repeat unit was subject to extension using cyclic expansion nonequilibrium molecular dynamics. Periodic boundary conditions were used to examine the effects of the contiguous alpha-helical linker on the force response. The measured force-extension curve shows a linear increase in the force response when the spectrin repeat unit is extended by approximately 0.4 nm. After that point, the force response peaks and subsequently declines. The peak in the force response marks the point where the spectrin repeat unit undergoes a change in its material properties from a strongly elastic material to a mostly viscous one, on the timescales of the simulations. The force peak is also correlated with rupture of the alpha-helical linker, and is likely the event responsible for the peaks in the sawtooth-pattern force-extension curves measured by atomic force microscopy experiments. Rupture of the linker involves simultaneously breaking approximately four hydrogen bonds that maintain the alpha-helical linker. After this initial rupture, the linker undergoes simple helix-to-coil transitions as the spectrin repeat unit continues to be extended. The implications of linker rupture in the interpretation of unfolding and atomic force microscopy experiments are also discussed.

Mechanical resistance of proteins explained using simple molecular models

Recent experiments have demonstrated that proteins unfold when two atoms are mechanically pulled apart, and that this process is different to when heated or when a chemical denaturant is added to the solution. Experiments have also shown that the response of proteins to external forces is very diverse, some of them being "hard," and others "soft." Mechanical resistance originates from the presence of barriers on the energy landscape; together, experiment and simulation have demonstrated that unfolding occurs through alternative pathways when different pairs of atoms undergo mechanical extension. Here we use simulation to probe the mechanical resistance of six structurally diverse proteins when pulled in different directions. For this, we use two very different models: a detailed, transferable one, and a coarse-grained, structure-based one. The coarse-grained model gives results that are surprisingly similar to the detailed one and qualitatively agree with experiment; i.e., the mechanical resistance of different proteins or of a single protein pulled in different directions can be predicted by simulation. The results demonstrate the importance of pulling direction relative to the local topology in determining mechanical stability, and rationalize the effect of the location of importation/degradation tags on the rates of mitochondrial import or protein degradation in vivo.

Single adhesive nanofibers from a live diatom have the signature fingerprint of modular proteins

The adhesive and mechanical properties of a cell-substratum adhesive secreted by live diatom cells were examined in situ using atomic force microscopy. The resulting force curves have a regular saw-tooth pattern, the characteristic fingerprint of modular proteins, and when bridged between tip and surface can repeatedly be stretched and relaxed resulting in precisely overlaying saw-tooth curves (up to approximately 600 successive cycles). The average rupture force of the peaks is 0.794 +/- 0.007 (mean +/- SE) nN at a loading rate of 0.8 microm/s and the average persistence length is 0.026 +/- <0.001 (mean +/- SE) nm (fit using the worm-like chain model). We propose that we are pulling on single adhesive nanofibers, each a cohesive unit composed of a set number of modular proteins aligned in register. Furthermore, we can observe and differentiate when up to three adhesive nanofibers are pulled based upon multimodal distributions of force and persistence length. The high force required for bond rupture, high extensibility (approximately 1.2 microm), and the accurate and rapid refolding upon relaxation, together provide strong and flexible properties ideally suited for the cell-substratum adhesion of this fouling diatom and allow us to understand the mechanism responsible for the strength of adhesion.

Mechanically unfolding the small, topologically simple protein L

beta-sheet proteins are generally more able to resist mechanical deformation than alpha-helical proteins. Experiments measuring the mechanical resistance of beta-sheet proteins extended by their termini led to the hypothesis that parallel, directly hydrogen-bonded terminal beta-strands provide the greatest mechanical strength. Here we test this hypothesis by measuring the mechanical properties of protein L, a domain with a topology predicted to be mechanically strong, but with no known mechanical function. A pentamer of this small, topologically simple protein is resistant to mechanical deformation over a wide range of extension rates. Molecular dynamics simulations show the energy landscape for protein L is highly restricted for mechanical unfolding and that this protein unfolds by the shearing apart of two structural units in a mechanism similar to that proposed for ubiquitin, which belongs to the same structural class as protein L, but unfolds at a significantly higher force. These data suggest that the mechanism of mechanical unfolding is conserved in proteins within the same fold family and demonstrate that although the topology and presence of a hydrogen-bonded clamp are of central importance in determining mechanical strength, hydrophobic interactions also play an important role in modulating the mechanical resistance of these similar proteins.