Jamming proteins with slipknots and their free energy landscape

Computational Design of Extraordinarily Stable Peptide Structures through Side-Chain-Locked Knots

Extraordinarily stable protein and peptide structures are critically demanded in many applications. Typical approaches to enhance protein and peptide stability are strengthening certain interactions. Here, we develop a very different approach: stabilizing peptide structures through side-chain-locked knots. More specifically, a peptide core consists of a knot, which is prevented from unknotting and unfolding by large side chains of amino acids at knot boundaries. These side chains impose free energy barriers for unknotting. The free energy barriers are quantified using all-atom and coarse-grained simulations. The barriers become infinitely high for large side chains and tight knot cores, resulting in stable peptide structures, which never unfold unless one chemical bond is broken. The extraordinary stability is essentially kinetic stability. Our new approach lifts the thermodynamic restriction in designing peptide structures, provides extra freedom in selecting sequence and structural motifs that are thermodynamically unstable, and should expand the functionality of peptides. This work also provides a bottom-up understanding of how knotting enhances protein stability.

Slipknotted and unknotted monovalent cation-proton antiporters evolved from a common ancestor

While the slipknot topology in proteins has been known for over a decade, its evolutionary origin is still a mystery. We have identified a previously overlooked slipknot motif in a family of two-domain membrane transporters. Moreover, we found that these proteins are homologous to several families of unknotted membrane proteins. This allows us to directly investigate the evolution of the slipknot motif. Based on our comprehensive analysis of 17 distantly related protein families, we have found that slipknotted and unknotted proteins share a common structural motif. Furthermore, this motif is conserved on the sequential level as well. Our results suggest that, regardless of topology, the proteins we studied evolved from a common unknotted ancestor single domain protein. Our phylogenetic analysis suggests the presence of at least seven parallel evolutionary scenarios that led to the current diversity of proteins in question. The tools we have developed in the process can now be used to investigate the evolution of other repeated-domain proteins.

In proteins with the slipknot topology, the polypeptide chain forms a slipknot—a structure that is not necessarily manifest to a naked eye, but it can be detected using mathematical methods. Slipknots are conserved motifs often found at catalytic sites and are directly involved in molecular transport. Although the first proteins with slipknots were found in 2007, many questions remain unanswered, e.g. how these proteins appeared, or whether the slipknotted proteins evolved from unknotted ones or vice versa. Here we provide the first analysis of homologous slipknotted and unknotted transmembrane proteins in order to elucidate their evolutionary relationship. We show that two-domain slipknotted and unknotted membrane transporters share the same one-domain unknotted protein as an ancestor. The ancestor gene duplicated and underwent various diversification and fusion events during the evolution, which have led to the appearance of a large superfamily of secondary active transporters. The slipknot motif seems to have been created by chance after a fusion of two single domain genes. Therefore, we show here that the slipknotted transporter evolved from an unknotted one-domain protein and that there are at least seven different evolutionary scenarios that gave rise to this large superfamily of transporters.

On folding of entangled proteins: knots, lassos, links and θ-curves

Around 6% of protein structures deposited in the PDB are entangled, forming knots, slipknots, lassos, links, and θ-curves. In each of these cases, the protein backbone weaves through itself in a complex way, and at some point passes through a closed loop, formed by other regions of the protein structure. Such a passing can be interpreted as crossing a topological barrier. How proteins overcome such barriers, and therefore different degrees of frustration, challenged scientists and has shed new light on the field of protein folding. In this review, we summarize the current knowledge about the free energy landscape of proteins with non-trivial topology. We describe identified mechanisms which lead proteins to self-tying. We discuss the influence of excluded volume, such as crowding and chaperones, on tying, based on available data. We briefly discuss the diversity of topological complexity of proteins and their evolution. We also list available tools to investigate non-trivial topology. Finally, we formulate intriguing and challenging questions at the boundary of biophysics, bioinformatics, biology, and mathematics, which arise from the discovery of entangled proteins.

Computational methods in the study of self-entangled proteins: a critical appraisal

The existence of self-entangled proteins, the native structure of which features a complex topology, unveils puzzling, and thus fascinating, aspects of protein biology and evolution. The discovery that a polypeptide chain can encode the capability to self-entangle in an efficient and reproducible way during folding, has raised many questions, regarding the possible function of these knots, their conservation along evolution, and their role in the folding paradigm. Understanding the function and origin of these entanglements would lead to deep implications in protein science, and this has stimulated the scientific community to investigate self-entangled proteins for decades by now. In this endeavour, advanced experimental techniques are more and more supported by computational approaches, that can provide theoretical guidelines for the interpretation of experimental results, and for the effective design of new experiments. In this review we provide an introduction to the computational study of self-entangled proteins, focusing in particular on the methodological developments related to this research field. A comprehensive collection of techniques is gathered, ranging from knot theory algorithms, that allow detection and classification of protein topology, to Monte Carlo or molecular dynamics strategies, that constitute crucial instruments for investigating thermodynamics and kinetics of this class of proteins.

Direct observation of the fast and robust folding of a slipknotted protein by optical tweezers

Understanding the folding mechanism of knotted and slipknotted proteins has attracted considerable interest. Due to their topological complexity, knotted and slipknotted proteins are predicted to fold slowly and involve large topological barriers. Molecular dynamics simulation studies suggest that a slipknotted conformation can serve as an important intermediate to help greatly reduce the topological difficulty during the folding of some knotted proteins. Here we use a single molecule optical tweezers technique to directly probe the folding of a small slipknotted protein AFV3-109. We found that stretching AFV3-109 can lead to the untying of the slipknot and complete unfolding of AFV3-109. Upon relaxation, AFV3-109 can readily refold back to its native slipknot conformation with high fidelity when the stretching force is lower than 6 pN. The refolding of AFV3-109 occurs in a sharp two-state like transition. Our results indicate that, different from knotted proteins, the folding of a slipknotted protein like AFV3-109 can be fast, and may not necessarily involve a high topological barrier.

KnotProt 2.0: a database of proteins with knots and other entangled structures

The KnotProt 2.0 database (the updated version of the KnotProt database) collects information about proteins which form knots and other entangled structures. New features in KnotProt 2.0 include the characterization of both probabilistic and deterministic entanglements which can be formed by disulfide bonds and interactions via ions, a refined characterization of entanglement in terms of knotoids, the identification of the so-called cysteine knots, the possibility to analyze all or a non-redundant set of proteins, and various technical updates. The KnotProt 2.0 database classifies all entangled proteins, represents their complexity in the form of a knotting fingerprint, and presents many biological and geometrical statistics based on these results. Currently the database contains >2000 entangled structures, and it regularly self-updates based on proteins deposited in the Protein Data Bank (PDB).

Proteins' Knotty Problems

Knots in proteins are increasingly being recognized as an important structural concept, and the folding of these peculiar structures still poses considerable challenges. From a functional point of view, most protein knots discovered so far are either enzymes or DNA-binding proteins. Our comprehensive topological analysis of the Protein Data Bank reveals several novel structures including knotted mitochondrial proteins and the most deeply embedded protein knot discovered so far. For the latter, we propose a novel folding pathway based on the idea that a loose knot forms at a terminus and slides to its native position. For the mitochondrial proteins, we discuss the folding problem from the perspective of transport and suggest that they fold inside the mitochondria. We also discuss the evolutionary origin of a novel class of knotted membrane proteins and argue that a novel knotted DNA-binding protein constitutes a new fold. Finally, we have also discovered a knot in an artificially designed protein structure.

Unraveling the mechanism of the cadherin-catenin-actin catch bond

The adherens junctions between epithelial cells involve a protein complex formed by E-cadherin, β-catenin, α-catenin and F-actin. The stability of this complex was a puzzle for many years, since in vitro studies could reconstitute various stable subsets of the individual proteins, but never the entirety. The missing ingredient turned out to be mechanical tension: a recent experiment that applied physiological forces to the complex with an optical tweezer dramatically increased its lifetime, a phenomenon known as catch bonding. However, in the absence of a crystal structure for the full complex, the microscopic details of the catch bond mechanism remain mysterious. Building on structural clues that point to α-catenin as the force transducer, we present a quantitative theoretical model for how the catch bond arises, fully accounting for the experimental lifetime distributions. The underlying hypothesis is that force induces a rotational transition between two conformations of α-catenin, overcoming a significant energy barrier due to a network of salt bridges. This transition allosterically regulates the energies at the interface between α-catenin and F-actin. The model allows us to predict these energetic changes, as well as highlighting the importance of the salt bridge rotational barrier. By stabilizing one of the α-catenin states, this barrier could play a role in how the complex responds to additional in vivo binding partners like vinculin. Since significant conformational energy barriers are a common feature of other adhesion systems that exhibit catch bonds, our model can be adapted into a general theoretical framework for integrating structure and function in a variety of force-regulated protein complexes.

KnotGenome: a server to analyze entanglements of chromosomes

The KnotGenome server enables the topological analysis of chromosome model data using three-dimensional coordinate files of chromosomes as input. In particular, it detects prime and composite knots in single chromosomes, and links between chromosomes. The knotting complexity of the chromosome is presented in the form of a matrix diagram that reveals the knot type of the entire polynucleotide chain and of each of its subchains. Links are determined by means of the Gaussian linking integral and the HOMFLY-PT polynomial. Entangled chromosomes are presented graphically in an intuitive way. It is also possible to relax structure with short molecular dynamics runs before the analysis. KnotGenome is freely available at http://knotgenom.cent.uw.edu.pl/.

To Tie or Not to Tie? That Is the Question

In this review, we provide an overview of entangled proteins. Around 6% of protein structures deposited in the PBD are entangled, forming knots, slipknots, lassos and links. We present theoretical methods and tools that enabled discovering and classifying such structures. We discuss the advantages and disadvantages of the non-trivial topology in proteins, based on available data about folding, stability, biological properties and evolutionary conservation. We also formulate intriguing and challenging questions on the border of biophysics, bioinformatics, biology and mathematics, which arise from the discovery of an entanglement in proteins. Finally, we discuss possible applications of entangled proteins in medicine and nanotechnology, such as the chance to design super stable proteins, whose stability could be controlled by chemical potential.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Folding analysis of the most complex Stevedore's protein knot

DehI is a homodimeric haloacid dehalogenase from Pseudomonas putida that contains the most complex 61 Stevedore's protein knot within its folding topology. To examine how DehI attains such an intricate knotted topology we combined far-UV circular dichroism (CD), intrinsic fluorescence spectroscopy and small angle X-ray scattering (SAXS) to investigate its folding mechanism. Equilibrium unfolding of DehI by chemical denaturation indicated the presence of two highly populated folding intermediates, I and I'. While the two intermediates vary in secondary structure contents and tertiary packing according to CD and intrinsic fluorescence, respectively, their overall dimension and compactness are similar according to SAXS. Three single-tryptophan variants (W34, W53, and W196) were generated to probe non-cooperative unfolding events localized around the three fluorophores. Kinetic fluorescence measurements indicated that the transition from the intermediate I' to the unfolded state is rate limiting. Our multiparametric folding analyses suggest that DehI unfolds through a linear folding pathway with two distinct folding intermediates by initial hydrophobic collapse followed by nucleation condensation, and that knotting precedes the formation of secondary structures.

Decoding the mechanical fingerprints of biomolecules

The capacity of biological macromolecules to act as exceedingly sophisticated and highly efficient cellular machines – switches, assembly factors, pumps, or motors – is realized through their conformational transitions, that is, their folding into distinct shapes and selective binding to other molecules. Conformational transitions can be induced, monitored, and manipulated by pulling individual macromolecules apart with an applied force. Pulling experiments reveal, for a given biomolecule, the relationship between applied force and molecular extension. Distinct signatures in the force–extension relationship identify a given biomolecule and thus serve as the molecule's ‘mechanical fingerprints’. But, how can these fingerprints be decoded to uncover the energy barriers crossed by the molecule in the course of its conformational transition, as well as the associated timescales? This review summarizes a powerful class of approaches to interpreting single-molecule force spectroscopy measurements – namely, analytically tractable approaches. On the fundamental side, analytical theories have the power to reveal the unifying principles underneath the bewildering diversity of biomolecules and their behaviors. On the practical side, analytical expressions that result from these theories are particularly well suited for a direct fit to experimental data, yielding the important parameters that govern biological processes at the molecular level.

Prediction of the optimal set of contacts to fold the smallest knotted protein

Knotted protein chains represent a new motif in protein folds. They have been linked to various diseases, and recent extensive analysis of the Protein Data Bank shows that they constitute 1.5% of all deposited protein structures. Despite thorough theoretical and experimental investigations, the role of knots in proteins still remains elusive. Nonetheless, it is believed that knots play an important role in mechanical and thermal stability of proteins. Here, we perform a comprehensive analysis of native, shadow-specific and non-native interactions which describe free energy landscape of the smallest knotted protein (PDB id 2efv). We show that the addition of shadow-specific contacts in the loop region greatly enhances folding kinetics, while the addition of shadow-specific contacts along the C-terminal region (H3 or H4) results in a new folding route with slower kinetics. By means of direct coupling analysis (DCA) we predict non-native contacts which also can accelerate kinetics. Next, we show that the length of the C-terminal knot tail is responsible for the shape of the free energy barrier, while the influence of the elongation of the N-terminus is not significant. Finally, we develop a concept of a minimal contact map sufficient for 2efv protein to fold and analyze properties of this protein using this map.

Molecular knots in biology and chemistry

Knots and entanglements are ubiquitous. Beyond their aesthetic appeal, these fascinating topological entities can be either useful or cumbersome. In recent decades, the importance and prevalence of molecular knots have been increasingly recognised by scientists from different disciplines. In this review, we provide an overview on the various molecular knots found in naturally occurring biological systems (DNA, RNA and proteins), and those created by synthetic chemists. We discuss the current knowledge in these fields, including recent developments in experimental and, in some cases, computational studies which are beginning to shed light into the complex interplay between the structure, formation and properties of these topologically intricate molecules.

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 (knot's 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.

Pierced Lasso Bundles are a new class of knot-like motifs

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.

Plasticity of hydrogen bond networks regulates mechanochemistry of cell adhesion complexes

Mechanical forces acting on cell adhesion receptor proteins regulate a range of cellular functions by formation and rupture of noncovalent interactions with ligands. Typically, force decreases the lifetimes of intact complexes ("slip bonds"), making the discovery that these lifetimes can also be prolonged ("catch bonds") a surprise. We created a microscopic analytic theory by incorporating the structures of selectin and integrin receptors into a conceptual framework based on the theory of stochastic equations, which quantitatively explains a wide range of experimental data (including catch bonds at low forces and slip bonds at high forces). Catch bonds arise due to force-induced remodeling of hydrogen bond networks, a finding that also accounts for unbinding in structurally unrelated integrin-fibronectin and actomyosin complexes. For the selectin family, remodeling of hydrogen bond networks drives an allosteric transition resulting in the formation of the maximum number of hydrogen bonds determined only by the structure of the receptor and independent of the ligand. A similar transition allows us to predict the increase in the number of hydrogen bonds in a particular allosteric state of α5β1 integrin-fibronectin complex, a conformation which is yet to be crystallized. We also make a testable prediction that a single point mutation (Tyr51Phe) in the ligand associated with selectin should dramatically alter the nature of the catch bond compared with the wild type. Our work suggests that nature uses a ductile network of hydrogen bonds to engineer function over a broad range of forces.

Knotting pathways in proteins

Most proteins, in order to perform their biological function, have to fold to a compact native state. The increasing number of knotted and slipknotted proteins identified suggests that proteins are able to manoeuvre around topological barriers during folding. In the present article, we review the current progress in elucidating the knotting process in proteins. Although we concentrate on theoretical approaches, where a knotted topology can be unambiguously detected, comparison with experiments is also reviewed. Numerical simulations suggest that the folding process for small knotted proteins is composed of twisted loop formation and then threading by either slipknot geometries or flipping. As the size of the knotted proteins increases, particularly for more deeply threaded termini, the prevalence of traps in the free energy landscape also increases. Thus, in the case of longer knotted and slipknotted proteins, the folding mechanism is probably supported by chaperones. Overall, results imply that knotted proteins can be folded efficiently and survive evolutionary pressure in order to perform their biological functions.

Knot localization in proteins

The backbones of proteins form linear chains. In the case of some proteins, these chains can be characterized as forming linear open knots. The knot type of the full chain reveals only limited information about the entanglement of the chain since, for example, subchains of an unknotted protein can form knots and subchains of a knotted protein can form different types of knots than the entire protein. To understand fully the entanglement within the backbone of a given protein, a complete analysis of the knotting within all of the subchains of that protein is necessary. In the present article, we review efforts to characterize the full knotting complexity within individual proteins and present a matrix that conveys information about various aspects of protein knotting. For a given protein, this matrix identifies the precise localization of knotted regions and shows the knot types formed by all subchains. The pattern in the matrix can be considered as a knotting fingerprint of that protein. We observe that knotting fingerprints of distantly related knotted proteins are strongly conserved during evolution and discuss how some characteristic motifs in the knotting fingerprints are related to the structure of the knotted regions and their possible biological role.

The unique cysteine knot regulates the pleotropic hormone leptin

Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.

Mechanically untying a protein slipknot: multiple pathways revealed by force spectroscopy and steered molecular dynamics simulations

Protein structure is highly diverse when considering a wide range of protein types, helping to give rise to the multitude of functions that proteins perform. In particular, certain proteins are known to adopt a knotted or slipknotted fold. How such proteins undergo mechanical unfolding was investigated utilizing a combination of single molecule atomic force microscopy (AFM), protein engineering, and steered molecular dynamics (SMD) simulations to show the mechanical unfolding mechanism of the slipknotted protein AFV3-109. Our results reveal that the mechanical unfolding of AFV3-109 can proceed via multiple parallel unfolding pathways that all cause the protein slipknot to untie and the polypeptide chain to completely extend. These distinct unfolding pathways proceed via either a two- or three-state unfolding process involving the formation of a well-defined, stable intermediate state. SMD simulations predict the same contour length increments for different unfolding pathways as single molecule AFM results, thus providing a plausible molecular mechanism for the mechanical unfolding of AFV3-109. These SMD simulations also reveal that two-state unfolding is initiated from both the N- and C-termini, while three-state unfolding is initiated only from the C-terminus. In both pathways, the protein slipknot was untied during unfolding, and no tightened slipknot conformation was observed. Detailed analysis revealed that interactions between key structural elements lock the knotting loop in place, preventing it from shrinking and the formation of a tightened slipknot conformation. Our results demonstrate the bifurcation of the mechanical unfolding pathway of AFV3-109 and point to the generality of a kinetic partitioning mechanism for protein folding/unfolding.

Structures and folding pathways of topologically knotted proteins

In the last decade, a new class of proteins has emerged that contain a topological knot in their backbone. Although these structures are rare, they nevertheless challenge our understanding of protein folding. In this review, we provide a short overview of topologically knotted proteins with an emphasis on newly discovered structures. We discuss the current knowledge in the field, including recent developments in both experimental and computational studies that have shed light on how these intricate structures fold.

Structure and folding of a designed knotted protein

A very small number of natural proteins have folded configurations in which the polypeptide backbone is knotted. Relatively little is known about the folding energy landscapes of such proteins, or how they have evolved. We explore those questions here by designing a unique knotted protein structure. Biophysical characterization and X-ray crystal structure determination show that the designed protein folds to the intended configuration, tying itself in a knot in the process, and that it folds reversibly. The protein folds to its native, knotted configuration approximately 20 times more slowly than a control protein, which was designed to have a similar tertiary structure but to be unknotted. Preliminary kinetic experiments suggest a complicated folding mechanism, providing opportunities for further characterization. The findings illustrate a situation where a protein is able to successfully traverse a complex folding energy landscape, though the amino acid sequence of the protein has not been subjected to evolutionary pressure for that ability. The success of the design strategy--connecting two monomers of an intertwined homodimer into a single protein chain--supports a model for evolution of knotted structures via gene duplication.

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 protein's resistance to stretching.

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.

A Stevedore's protein knot

Protein knots, mostly regarded as intriguing oddities, are gradually being recognized as significant structural motifs. Seven distinctly knotted folds have already been identified. It is by and large unclear how these exceptional structures actually fold, and only recently, experiments and simulations have begun to shed some light on this issue. In checking the new protein structures submitted to the Protein Data Bank, we encountered the most complex and the smallest knots to date: A recently uncovered alpha-haloacid dehalogenase structure contains a knot with six crossings, a so-called Stevedore knot, in a projection onto a plane. The smallest protein knot is present in an as yet unclassified protein fragment that consists of only 92 amino acids. The topological complexity of the Stevedore knot presents a puzzle as to how it could possibly fold. To unravel this enigma, we performed folding simulations with a structure-based coarse-grained model and uncovered a possible mechanism by which the knot forms in a single loop flip.

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.