Joanna I. Sulkowska

Genomics-aided structure prediction

We introduce a theoretical framework that exploits the ever-increasing genomic sequence information for protein structure prediction. Structure-based models are modified to incorporate constraints by a large number of non-local contacts estimated from direct coupling analysis (DCA) of co-evolving genomic sequences. A simple hybrid method, called DCA-fold, integrating DCA contacts with an accurate knowledge of local information (e.g., the local secondary structure) is sufficient to fold proteins in the range of 1–3 Å resolution.

Mechanical stretching of proteins—a theoretical survey of the Protein Data Bank

The mechanical stretching of single proteins has been studied experimentally for about 50 proteins, yielding a variety of force patterns and peak forces. Here we perform a theoretical survey of proteins of known native structure and map out the landscape of possible dynamical behaviours under stretching at constant speed. We consider 7510 proteins comprising not more than 150 amino acids and 239 longer proteins. The model used is constructed based on the native geometry. It is solved by methods of molecular dynamics and validated by comparing the theoretical predictions to experimental results. We characterize the distribution of peak forces and investigate correlations with the system size and with the structure classification as characterized by the CATH scheme. Despite the presence of such correlations, proteins with the same CATH index may belong to different classes of dynamical behaviour. We identify proteins with the biggest forces and show that they belong to few topology classes. We determine which protein segments act as mechanical clamps and show that, in most cases, they correspond to long stretches of parallel β-strands, but other mechanisms are also possible.

Conservation of complex knotting and slipknotting patterns in proteins

While analyzing all available protein structures for the presence of knots and slipknots, we detected a strict conservation of complex knotting patterns within and between several protein families despite their large sequence divergence. Because protein folding pathways leading to knotted native protein structures are slower and less efficient than those leading to unknotted proteins with similar size and sequence, the strict conservation of the knotting patterns indicates an important physiological role of knots and slipknots in these proteins. Although little is known about the functional role of knots, recent studies have demonstrated a protein-stabilizing ability of knots and slipknots. Some of the conserved knotting patterns occur in proteins forming transmembrane channels where the slipknot loop seems to strap together the transmembrane helices forming the channel.

Dodging the crisis of folding proteins with knots

Proteins with nontrivial topology, containing knots and slipknots, have the ability to fold to their native states without any additional external forces invoked. A mechanism is suggested for folding of these proteins, such as YibK and YbeA, that involves an intermediate configuration with a slipknot. It elucidates the role of topological barriers and backtracking during the folding event. It also illustrates that native contacts are sufficient to guarantee folding in ≈1–2% of the simulations, and how slipknot intermediates are needed to reduce the topological bottlenecks. As expected, simulations of proteins with similar structure but with knot removed fold much more efficiently, clearly demonstrating the origin of these topological barriers. Although these studies are based on a simple coarse-grained model, they are already able to extract some of the underlying principles governing folding in such complex topologies.

On the remarkable mechanostability of scaffoldins and the mechanical clamp motif

Protein mechanostability is a fundamental biological property that can only be measured by single-molecule manipulation techniques. Such studies have unveiled a variety of highly mechanostable modules (mainly of the Ig-like, β-sandwich type) in modular proteins subjected to mechanical stress from the cytoskeleton and the metazoan cell–cell interface. Their mechanostability is often attributed to a “mechanical clamp” of secondary structure (a patch of backbone hydrogen bonds) fastening their ends. Here we investigate the nanomechanics of scaffoldins, an important family of scaffolding proteins that assembles a variety of cellulases into the so-called cellulosome, a microbial extracellular nanomachine for cellulose adhesion and degradation. These proteins anchor the microbial cell to cellulose substrates, which makes their connecting region likely to be subjected to mechanical stress. By using single-molecule force spectroscopy based on atomic force microscopy, polyprotein engineering, and computer simulations, here we show that the cohesin I modules from the connecting region of cellulosome scaffoldins are the most robust mechanical proteins studied experimentally or predicted from the entire Protein Data Bank. The mechanostability of the cohesin modules studied correlates well with their mechanical kinetic stability but not with their thermal stability, and it is well predicted by computer simulations, even coarse-grained. This extraordinary mechanical stability is attributed to 2 mechanical clamps in tandem. Our findings provide the current upper limit of protein mechanostability and establish shear mechanical clamps as a general structural/functional motif widespread in proteins putatively subjected to mechanical stress. These data have important implications for the scaffoldin physiology and for protein design in biotechnology and nanotechnology.

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.

Stabilizing effect of knots on proteins

Molecular dynamics studies within a coarse-grained, structure-based model were used on two similar proteins belonging to the transcarbamylase family to probe the effects of the knot in the native structure of a protein. The first protein, N-acetylornithine transcarbamylase, contains no knot, whereas human ormithine transcarbamylase contains a trefoil knot located deep within the sequence. In addition, we also analyzed a modified transferase with the knot removed by the appropriate change of a knot-making crossing of the protein chain. The studies of thermally and mechanically induced unfolding processes suggest a larger intrinsic stability of the protein with the knot.

Slipknotting upon native-like loop formation in a trefoil knot protein

Protein knots and slipknots, mostly regarded as intriguing oddities, are gradually being recognized as significant structural motifs. Recent experimental results show that knotting, starting from a fully extended polypeptide, has not yet been observed. Understanding the nucleation process of folding knots is thus a natural challenge for both experimental and theoretical investigation. In this study, we employ energy landscape theory and molecular dynamics to elucidate the entire folding mechanism. The full free energy landscape of a knotted protein is mapped using an all-atom structure-based protein model. Results show that, due to the topological constraint, the protein folds through a three-state mechanism that contains (i) a precise nucleation site that creates a correctly twisted native loop (first barrier) and (ii) a rate-limiting free energy barrier that is traversed by two parallel knot-forming routes. The main route corresponds to a slipknot conformation, a collapsed configuration where the C-terminal helix adopts a hairpin-like configuration while threading, and the minor route to an entropically limited plug motion, where the extended terminus is threaded as through a needle. Knot formation is a late transition state process and results show that random (nonspecific) knots are a very rare and unstable set of configurations both at and below folding temperature. Our study shows that a native-biased landscape is sufficient to fold complex topologies and presents a folding mechanism generalizable to all known knotted protein topologies: knotting via threading a native-like loop in a preordered intermediate.

Tightening of knots in proteins.

We perform theoretical studies of stretching of 20 proteins with knots within a coarse-grained model. The knots ends are found to jump to well defined sequential locations that are associated with sharp turns, whereas in homopolymers they diffuse around and eventually slide off. The waiting times of the jumps are increasingly stochastic as the temperature is raised. Knots typically do not return to their native locations when a protein is released after stretching.

KnotProt: a database of proteins with knots and slipknots

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knots core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.