Mateusz Sikora

Mechanical stability of multidomain proteins and novel mechanical clamps

We estimate the size of mechanostability for 318 multidomain proteins which are single‐chain and contain up to 1021 amino acids. We predict existence of novel types of mechanical clamps in which interdomain contacts play an essential role. Mechanical clamps are structural regions which are the primary source of a proteins resistance to pulling. Among these clamps there is one that opposes tensile stress due to two domains swinging apart. This movement strains and then ruptures the contacts that hold the two domains together. Another clamp also involves tensile stress but it originates from an immobilization of a structural region by a surrounding knot‐loop (without involving any disulfide bonds). Still another mechanism involves shear between helical regions belonging to two domains. We also consider the amyloid‐prone cystatin C which provides an example of a two‐chain 3D domain‐swapped protein. We predict that this protein should withstand remarkably large stress, perhaps of order 800 pN, when inducing a shearing strain. The survey is generated through molecular dynamics simulations performed within a structure‐based coarse grained model. Proteins 2011;

Actin Rings of Power.

Circular or ring-like actin structures play important roles in various developmental and physiological processes. Commonly, these rings are composed of actin filaments and myosin motors (actomyosin) that, upon activation, trigger ring constriction. Actomyosin ring constriction, in turn, has been implicated in key cellular processes ranging from cytokinesis to wound closure. Non-constricting actin ring-like structures also form at cell-cell contacts, where they exert a stabilizing function. Here, we review recent studies on the formation and function of actin ring-like structures in various morphogenetic processes, shedding light on how those different rings have been adapted to fulfill their specific roles.

BSDB: the biomolecule stretching database

We describe the Biomolecule Stretching Data Base that has been recently set up at http://www.ifpan.edu.pl/BSDB/. It provides information about mechanostability of proteins. Its core is based on simulations of stretching of 17 134 proteins within a structure-based model. The primary information is about the heights of the maximal force peaks, the force–displacement patterns, and the sequencing of the contact-rupturing events. We also summarize the possible types of the mechanical clamps, i.e. the motifs which are responsible for a proteins resistance to stretching.

Linker-mediated assembly of gold nanoparticles into multimeric motifs

We present a theoretical description of linker-mediated self-assembly of gold nanoparticles (Au-NP). Using mesoscale simulations with a coarse-grained model for the Au NPs and dirhenium-based linker molecules, we investigate the conditions under which large clusters can grow and construct a phase diagram that identifies favorable growth conditions in terms of floating and bound linker concentrations. The findings can be considered as generic, as we expect other NP-linker systems to behave in a qualitatively similar way. In particular, we also discuss the case of antibody-functionalised Au NPs connected by the C-reactive proteins (CRPs). We extract some general rules for NP linking that may aid the production of size- and shape-specific NP clusters for technology applications.

Molecular jamming—The cystine slipknot mechanical clamp in all-atom simulations

A recent survey of 17 134 proteins has identified a new class of proteins which are expected to yield stretching induced force peaks in the range of 1 nN. Such high force peaks should be due to forcing of a slip-loop through a cystine ring, i.e., by generating a cystine slipknot. The survey has been performed in a simple coarse grained model. Here, we perform all-atom steered molecular dynamics simulations on 15 cystine knot proteins and determine their resistance to stretching. In agreement with previous studies within a coarse grained structure based model, the level of resistance is found to be substantially higher than in proteins in which the mechanical clamp operates through shear. The large stretching forces arise through formation of the cystine slipknot mechanical clamp and the resulting steric jamming. We elucidate the workings of such a clamp in an atomic detail. We also study the behavior of five top strength proteins with the shear-based mechanostability in which no jamming is involved. We show that in the atomic model, the jamming state is relieved by moving one amino acid at a time and there is a choice in the selection of the amino acid that advances the first. In contrast, the coarse grained model also allows for a simultaneous passage of two amino acids.

Cystine plug and other novel mechanisms of large mechanical stability in dimeric proteins.

We identify three dimeric proteins whose mechanostability is anisotropic and should exceed 1 nN along some directions. They come with distinct mechanical clamps: either shear-based, or involving a cystine slipknot, or due to dragging of a cystine plug through a cystine ring. The latter two mechanisms are topological in nature; the cystine plug mechanism has not yet been discussed but it turns out to provide the largest resistance to stretching. Its possible applications in elastomers are discussed.

An Exploration of the Universe of Polyglutamine Structures.

Deposits of misfolded proteins in the human brain are associated with the development of many neurodegenerative diseases. Recent studies show that these proteins have common traits even at the monomer level. Among them, a polyglutamine region that is present in huntingtin is known to exhibit a correlation between the length of the chain and the severity as well as the earliness of the onset of Huntington disease. Here, we apply bias exchange molecular dynamics to generate structures of polyglutamine expansions of several lengths and characterize the resulting independent conformations. We compare the properties of these conformations to those of the standard proteins, as well as to other homopolymeric tracts. We find that, similar to the previously studied polyvaline chains, the set of possible transient folds is much broader than the set of known-to-date folds, although the conformations have different structures. We show that the mechanical stability is not related to any simple geometrical characteristics of the structures. We demonstrate that long polyglutamine expansions result in higher mechanical stability than the shorter ones. They also have a longer life span and are substantially more prone to form knotted structures. The knotted region has an average length of 35 residues, similar to the typical threshold for most polyglutamine-related diseases. Similarly, changes in shape and mechanical stability appear once the total length of the peptide exceeds this threshold of 35 glutamine residues. We suggest that knotted conformers may also harm the cellular machinery and thus lead to disease.

Theoretical tests of the mechanical protection strategy in protein nanomechanics

We provide theoretical tests of a novel experimental technique to determine mechanostability of proteins based on stretching a mechanically protected protein by single‐molecule force spectroscopy. This technique involves stretching a homogeneous or heterogeneous chain of reference proteins (single‐molecule markers) in which one of them acts as host to the guest protein under study. The guest protein is grafted into the host through genetic engineering. It is expected that unraveling of the host precedes the unraveling of the guest removing ambiguities in the reading of the force‐extension patterns of the guest protein. We study examples of such systems within a coarse‐grained structure‐based model. We consider systems with various ratios of mechanostability for the host and guest molecules and compare them to experimental results involving cohesin I as the guest molecule. For a comparison, we also study the force‐displacement patterns in proteins that are linked in a serial fashion. We find that the mechanostability of the guest is similar to that of the isolated or serially linked protein. We also demonstrate that the ideal configuration of this strategy would be one in which the host is much more mechanostable than the single‐molecule markers. We finally show that it is troublesome to use the highly stable cystine knot proteins as a host to graft a guest in stretching studies because this would involve a cleaving procedure. Proteins 2014; 82:717–726.

Formation of Cystine Slipknots in Dimeric Proteins

We consider mechanical stability of dimeric and monomeric proteins with the cystine knot motif. A structure based dynamical model is used to demonstrate that all dimeric and some monomeric proteins of this kind should have considerable resistance to stretching that is significantly larger than that of titin. The mechanisms of the large mechanostability are elucidated. In most cases, it originates from the induced formation of one or two cystine slipknots. Since there are four termini in a dimer, there are several ways of selecting two of them to pull by. We show that in the cystine knot systems, there is strong anisotropy in mechanostability and force patterns related to the selection. We show that the thermodynamic stability of the dimers is enhanced compared to the constituting monomers whereas machanostability is either lower or higher.

Overcoming the limitations of the MARTINI force field in simulations of polysaccharides

Polysaccharides (carbohydrates) are key regulators of a large number of cell biological processes. However, precise biochemical or genetic manipulation of these often complex structures is laborious and hampers experimental structure-function studies. Molecular Dynamics (MD) simulations provide a valuable alternative tool to generate and test hypotheses on saccharide function. Yet, currently used MD force fields often overestimate the aggregation propensity of polysaccharides, affecting the usability of those simulations. Here we tested MARTINI, a popular coarse-grained (CG) force field for biological macromolecules, for its ability to accurately represent molecular forces between saccharides. To this end, we calculated a thermodynamic solution property, the second virial coefficient of the osmotic pressure (B22). Comparison with light scattering experiments revealed a nonphysical aggregation of a prototypical polysaccharide in MARTINI, pointing at an imbalance of the nonbonded solute-solute, solute-water, and water-water interactions. This finding also applies to smaller oligosaccharides which were all found to aggregate in simulations even at moderate concentrations, well below their solubility limit. Finally, we explored the influence of the Lennard-Jones (LJ) interaction between saccharide molecules and propose a simple scaling of the LJ interaction strength that makes MARTINI more reliable for the simulation of saccharides.