Source author record

Marek Cieplak

Marek Cieplak appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2021arXiv

Conformational Biases of α-Synuclein and Formation of Transient Knots

We study local conformational biases in the dynamics of α-synuclein by using all-atom simulations with explicit and implicit solvents. The biases are related to the frequency of the specific contact formation. In both approaches, the protein is intrinsically disordered, and its strongest bias is to make bend and turn local structures. The explicit-solvent conformations can be substantially more extended which allows for formation of transient trefoil knots, both deep and shallow, that may last for up to 5 μs. The two-chain self-association events, both short- and long-lived, are dominated by formation of contacts in the central part of the sequence. This part tends to form helices when bound to a micelle.

preprint2016arXiv

Criteria for folding in structure-based models of proteins

In structure-based models of proteins, one often assumes that folding is accomplished when all contacts are established. This assumption may frequently lead to a conceptual problem that folding takes place in a temperature region of very low thermodynamic stability, especially when the contact map used is too sparse. We consider six different structure-based models and show that allowing for a small, but model-dependent, percentage of the native contacts not being established boosts the folding temperature substantially while affecting the time scales of folding only in a minor way. We also compare other properties of the six models. We show that the choice of the description of the backbone stiffness has a substantial effect on the values of characteristic temperatures that relate both to equilibrium and kinetic properties. Models without any backbone stiffness (like the self-organized polymer) are found to perform similar to those with the stiffness, including in the studies of stretching.

preprint2016arXiv

Stiffness of the C-terminal disordered linker affects the geometry of the active site in endoglucanase Cel8A

Cellulosomes are complex multi-enzyme machineries which efficiently degrade plant cell-wall polysaccharides. The multiple domains of the cellulosome proteins are often tethered together by intrinsically disordered regions. The properties and functions of these disordered linkers are not well understood. In this work, we study endoglucanase Cel8A, which is a relevant enzymatic component of the cellulosomes of Clostridium thermocellum. We use both all-atom and coarse-grained simulations to investigate how the equilibrium conformations of the catalytic domain of Cel8A are affected by the disordered linker at its C terminus. We find that when the endoglucanase is bound to its substrate, the effective stiffness of the linker can influence the distances between groups of amino-acid residues throughout the entire enzymatic domain. In particular, variations in the linker stiffness can lead to small changes in the geometry of the active-site cleft. We suggest that such geometrical changes may, in turn, have an effect on the catalytic activity of the enzyme.

preprint2016arXiv

Unfolding knots by proteasome-like systems: simulations of the behaviour of folded and neurotoxic proteins

Knots in proteins have been proposed to resist proteasomal degradation. Ample evidence associates proteasomal degradation with neurodegeneration. One interesting possibility is that indeed knotted conformers stall this machinery leading to toxicity. However, although the proteasome is known to unfold mechanically its substrates, at present there are no experimental methods to emulate this particular traction geometry. Here, we consider several dynamical models of the proteasome in which the complex is represented by an effective potential with an added pulling force. This force is meant to induce translocation of a protein or a polypeptide into the catalytic chamber. The force is either constant or applied periodically. The translocated proteins are modelled in a coarse-grained fashion. We do comparative analysis of several knotted globular proteins and the transiently knotted polyglutamine tracts of length 60 alone and fused in exon 1 of the huntingtin protein. Huntingtin is associated with Huntington disease, a well-known genetically-determined neurodegenerative disease. We show that the presence of a knot hinders and sometimes even jams translocation. We demonstrate that the probability to do so depends on the protein, the model of the proteasome, the magnitude of the pulling force, and the choice of the pulled terminus. In any case, the net effect would be a hindrance in the proteasomal degradation process in the cell. This would then yield toxicity \textit{via} two different mechanisms: one through toxic monomers compromising degradation and another by the formation of toxic oligomers. Our work paves the way to the mechanistic investigation of the mechanical unfolding of knotted structures by the proteasome and its relation to toxicity and disease.

preprint2015arXiv

Citrate synthase proteins in extremophilic organisms: Studies within a structure-based model

We study four citrate synthase homodimeric proteins within a structure-based coarse-grained model. Two of these proteins come from thermophilic bacteria, one from a cryophilic bacterium and one from a mesophilic organism; three are in the closed and two in the open conformations. Even though the proteins belong to the same fold, the model distinguishes the properties of these proteins in a way which is consistent with experiments. For instance, the thermophilic proteins are more stable thermodynamically than their mesophilic and cryophilic homologues, which we observe both in the magnitude of thermal fluctuations near the native state and in the kinetics of thermal unfolding. The level of stability correlates with the average coordination number for amino acids contacts and with the degree of structural compactness. The pattern of positional fluctuations along the sequence in the closed conformation is different than in the open conformation, including within the active site. The modes of correlated and anticorrelated movements of pairs of amino acids forming the active site are very different in the open and closed conformations. Taken together, our results show that the precise location of amino acid contacts in the native structure appears to be a critical element in explaining the similarities and differences in the thermodynamic properties, local flexibility and collective motions of the different forms of the enzyme.

preprint2015arXiv

Cotranslational folding of deeply knotted proteins

Proper folding of deeply knotted proteins has a very low success rate even in structure-based models which favor formation of the native contacts but have no topological bias. By employing a structure-based model, we demonstrate that cotranslational folding on a model ribosome may enhance the odds to form trefoil knots for protein YibK without any need to introduce any non-native contacts. The ribosome is represented by a repulsive wall that keeps elongating the protein. On-ribosome folding proceeds through a a slipknot conformation. We elucidate the mechanics and energetics of its formation. We show that the knotting probability in on-ribosome folding is a function of temperature and that there is an optimal temperature for the process. Our model often leads to the establishment of the native contacts without formation of the knot.

preprint2015arXiv

Multiple folding pathways of proteins with shallow knots and co-translational folding

We study the folding process in the shallowly knotted protein MJ0366 within two variants of a structure-based model. We observe that the resulting topological pathways are much richer than identified in previous studies. In addition to the single knot-loop events, we find novel, and dominant, two-loop mechanisms. We demonstrate that folding takes place in a range of temperatures and the conditions of most successful folding are at temperatures which are higher than those required for the fastest folding. We also demonstrate that nascent conditions are more favorable to knotting than off-ribosome folding.

preprint2015arXiv

Unbinding and unfolding of adhesion protein complexes through stretching: Interplay between shear and tensile mechanical clamps

Using coarse-grained molecular dynamics simulations, we analyze mechanically induced dissociation and unfolding of the protein complex CD48-2B4. This heterodimer is an indispensable component of the immunological system: 2B4 is a receptor on natural killer cells whereas CD48 is expressed on surfaces of various immune cells. So far, its mechanostability has not been assessed either experimentally or theoretically. We find that the dissociation processes strongly depend on the direction of pulling and may take place in several pathways. Interestingly, the CD48-2B4 interface can be divided into three distinct patches that act as units when resisting the pulling forces. At experimentally accessible pulling speeds, the characteristic mechanostability forces are in the range between 100 and 200 pN, depending on the pulling direction. These characteristic forces need not be associated with tensile forces involved in the act of separation of the complex because prior shear-involving unraveling within individual proteins may give rise to a higher force peak.

preprint2014arXiv

Amino acids and proteins at ZnO-water interfaces in molecular dynamics simulations

We determine potentials of the mean force for interactions of amino acids with four common surfaces of ZnO in aqueous solutions. The method involves all-atom molecular dynamics simulations combined with the umbrella sampling technique. The profiled nature of the density of water with the strongly adsorbed first layer affects the approach of amino acids to the surface and generates either repulsion or weak binding. The largest binding energy is found for tyrosine interacting with the surface in which the Zn ions are at the top. It is equal to 7 kJ/mol which is comparable to that of the hydrogen bonds in a protein. This makes the adsorption of amino acids onto the ZnO surface much weaker than onto the well studied surface of gold. Under vacuum, binding energies are more than 40 times stronger (for one of the surfaces). The precise manner in which water molecules interact with a given surface influences the binding energies in a way that depends on the surface. Among the four considered surfaces the one with Zn at the top is recognized as binding almost all amino acids with an average binding energy of 2.60 kJ/mol. Another (O at the top) is non-binding for most amino acids. For binding situations the average energy is 0.66 kJ/mol. The remaining two surfaces bind nearly as many amino acids as they do not and the average binding energies are 1.46 and 1.22 kJ/mol. For all of the surfaces the binding energies vary between amino acids significantly: the dispersion in the range of 68-154% of the mean. A small protein is shown to adsorb onto ZnO only intermittently and with only a small deformation. Various adsorption events lead to different patterns in mobilities of amino acids within the protein.

preprint2014arXiv

Aqueous Amino Acids and Proteins Near the Surface of Gold in Hydrophilic and Hydrophobic Force Fields

We calculate potentials of the mean force for twenty amino acids in the vicinity of the (111) surface of gold, for several dipeptides, and for some analogs of the side chains, using molecular dynamics simulations and the umbrella sampling method. We compare results obtained within three different force fields: one hydrophobic (for a contaminated surface) and two hydrophilic. All of these fields lead to good binding with very different specificities and different patterns in the density and polarization of water. The covalent bond with the sulfur atom on cysteine is modeled by the Morse potential. We demonstrate that binding energies of dipeptides are different than the combined binding energies of their amino-acidic components. For the hydrophobic gold, adsorption events of a small protein are driven by attraction to the strongest binding amino acids. This is not so in the hydrophilic cases - a result of smaller specificities combined with the difficulty for proteins, but not for single amino acids, to penetrate the first layer of water. The properties of water near the surface sensitively depend on the force field.

preprint2014arXiv

Interactions of aqueous amino acids and proteins with the (110) surface of ZnS in molecular dynamics simulations

The growing usage of nanoparticles of zinc sulfide as quantum dots and biosensors calls for a theoretical assessment of interactions of ZnS with biomolecules. We employ the molecular-dynamics-based umbrella sampling method to determine potentials of mean force for 20 single amino acids near the ZnS (110) surface in aqueous solutions. We find that five amino acids do not bind at all and the binding energy of the remaining amino acids does not exceed 4.3 kJ/mol. Such energies are comparable to those found for ZnO (and to hydrogen bonds in proteins) but the nature of the specificity is different. Cysteine can bind with ZnS in a covalent way, e.g. by forming the disulfide bond with S in the solid. If this effect is included within a model incorporating the Morse potential, then the potential well becomes much deeper - the binding energy is close to 98 kJ/mol. We then consider tryptophan cage, a protein of 20 residues, and characterize its events of adsorption to ZnS. We demonstrate the relevance of interactions between the amino acids in the selection of optimal adsorbed conformations and recognize the key role of cysteine in generation of lasting adsorption. We show that ZnS is more hydrophobic than ZnO and that the density profile of water is quite different than that forming near ZnO - it has only a minor articulation into layers. Furthermore, the first layer of water is disordered and mobile.

preprint2013arXiv

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.

preprint2012arXiv

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: shear-based, involving a cystine slipknot, and due to dragging of a cystine plug through a cystine ring. The latter two mechanisms are topological in nature and 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.

preprint2011arXiv

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.

preprint2010arXiv

Nanoindentation of virus capsids in a molecular model

A molecular-level model is used to study the mechanical response of empty cowpea chlorotic mottle virus (CCMV) and cowpea mosaic virus (CPMV) capsids. The model is based on the native structure of the proteins that consitute the capsids and is described in terms of the C-alpha atoms. Nanoindentation by a large tip is modeled as compression between parallel plates. Plots of the compressive force versus plate separation for CCMV are qualitatively consistent with continuum models and experiments, showing an elastic region followed by an irreversible drop in force. The mechanical response of CPMV has not been studied, but the molecular model predicts an order of magnitude higher stiffness and a much shorter elastic region than for CCMV. These large changes result from small structural changes that increase the number of bonds by only 30% and would be difficult to capture in continuum models. Direct comparison of local deformations in continuum and molecular models of CCMV shows that the molecular model undergoes a gradual symmetry breaking rotation and accommodates more strain near the walls than the continuum model. The irreversible drop in force at small separations is associated with rupturing nearly all of the bonds between capsid proteins in the molecular model while a buckling transition is observed in continuum models.

preprint2009arXiv

First-principles design of nanomachines

Learning from nature's amazing molecular machines, globular proteins, we present a framework for the predictive design of nano-machines. We show that the crucial ingredients for a chain molecule to behave as a machine are its inherent anisotropy and the coupling between the local Frenet coordinate reference frames of nearby monomers. We demonstrate that, even in the absence of heterogeneity, protein-like behavior is obtained for a simple chain molecule made up of just thirty hard spheres. This chain spontaneously switches between two distinct geometries, a single helix and a dual helix, merely due to thermal fluctuations.

preprint2008arXiv

Stretching and twisting of the DNA duplexes in coarse grained dynamical models

Three coarse-grained models of the double-stranded DNA are proposed and compared in the context of mechanical manipulation such as twisting and various schemes of stretching. The models differ in the number of effective beads (between two and five) representing each nucleotide. They all show similar behavior and, in particular, lead to a torque-force phase diagrams qualitatively consistent with experiments and all-atom simulations.

preprint2003arXiv

Folding of Proteins in Go Models with Angular Interactions

Molecular dynamics studies of Go models of proteins with the 10-12 contact potential and the bond and dihedral angle terms indicate statistical similarities to other Go models, e.g. with the Lennard-Jones contact potentials. The folding times depend on the protein size as power laws with the exponents depending on the native structural classes. There is no dependence of the folding times on the relative contact order even though the folding scenarios are governed mostly by the contact order.