Knotted vs. unknotted proteins: evidence of knot-promoting loops

Exploring the electronic structure of knotted proteins: the case of two ornithine transcarbamylase family

CONTEXT

Geometrical knots are rare structural arrangements in proteins in which the polypeptide chain ties itself into a knot, which is very intriguing due to the uncertainty of their impact on the protein properties. Presently, classical molecular dynamics is the most employed technique in the few studies found on this topic, so any information on how the presence of knots affects the reactivity and electronic properties of proteins is even scarcer. Using the electronic structure methods and quantum chemical descriptors analysis, we found that the same amino-acid residues in the knot core have statistically larger values for the unknotted protein, for both hard-hard and soft-soft interaction descriptors. In addition, we present a computationally feasible protocol, where we show it is possible to separate the contribution of the geometrical knot to the reactivity and other electronic structure properties.

METHODS

In order to investigate these systems, we used PRIMoRDiA, a new software developed by our research group, to explore the electronic structure of biological macromolecules. We evaluated several local quantum chemical descriptors to unveil relevant patterns potentially originating from the presence of the geometrical knot in two proteins, belonging to the ornithine transcarbamylase family. We compared several sampled structures from these two enzymes that are highly similar in both tertiary structure and function, but one of them has a knot whereas the other does not. The sampling was carried out through molecular dynamics simulations using ff14SB force field along 50 ns, and the semiempirical convergence was performed with PM7 Hamiltonian.

Knot or not? Identifying unknotted proteins in knotted families with sequence-based Machine Learning model

Knotted proteins, although scarce, are crucial structural components of certain protein families, and their roles continue to be a topic of intense research. Capitalizing on the vast collection of protein structure predictions offered by AlphaFold (AF), this study computationally examines the entire UniProt database to create a robust dataset of knotted and unknotted proteins. Utilizing this dataset, we develop a machine learning (ML) model capable of accurately predicting the presence of knots in protein structures solely from their amino acid sequences. We tested the model's capabilities on 100 proteins whose structures had not yet been predicted by AF and found agreement with our local prediction in 92% cases. From the point of view of structural biology, we found that all potentially knotted proteins predicted by AF can be classified only into 17 families. This allows us to discover the presence of unknotted proteins in families with a highly conserved knot. We found only three new protein families: UCH, DUF4253, and DUF2254, that contain both knotted and unknotted proteins, and demonstrate that deletions within the knot core could potentially account for the observed unknotted (trivial) topology. Finally, we have shown that in the majority of knotted families (11 out of 15), the knotted topology is strictly conserved in functional proteins with very low sequence similarity. We have conclusively demonstrated that proteins AF predicts as unknotted are structurally accurate in their unknotted configurations. However, these proteins often represent nonfunctional fragments, lacking significant portions of the knot core (amino acid sequence).

Are there double knots in proteins? Prediction and in vitro verification based on TrmD-Tm1570 fusion from C. nitroreducens

We have been aware of the existence of knotted proteins for over 30 years-but it is hard to predict what is the most complicated knot that can be formed in proteins. Here, we show new and the most complex knotted topologies recorded to date-double trefoil knots (31 #31). We found five domain arrangements (architectures) that result in a doubly knotted structure in almost a thousand proteins. The double knot topology is found in knotted membrane proteins from the CaCA family, that function as ion transporters, in the group of carbonic anhydrases that catalyze the hydration of carbon dioxide, and in the proteins from the SPOUT superfamily that gathers 31 knotted methyltransferases with the active site-forming knot. For each family, we predict the presence of a double knot using AlphaFold and RoseTTaFold structure prediction. In the case of the TrmD-Tm1570 protein, which is a member of SPOUT superfamily, we show that it folds in vitro and is biologically active. Our results show that this protein forms a homodimeric structure and retains the ability to modify tRNA, which is the function of the single-domain TrmD protein. However, how the protein folds and is degraded remains unknown.

Folding and functions of knotted proteins

Topologically knotted proteins have entangled structural elements within their native structures that cannot be disentangled simply by pulling from the N- and C-termini. Systematic surveys have identified different types of knotted protein structures, constituting as much as 1% of the total entries within the Protein Data Bank. Many knotted proteins rely on their knotted structural elements to carry out evolutionarily conserved biological functions. Being knotted may also provide mechanical stability to withstand unfolding-coupled proteolysis. Reconfiguring a knotted protein topology by circular permutation or cyclization provides insights into the importance of being knotted in the context of folding and functions. With the explosion of predicted protein structures by artificial intelligence, we are now entering a new era of exploring the entangled protein universe.

Homology of homologous knotted proteins

Quantification and classification of protein structures, such as knotted proteins, often requires noise-free and complete data. Here, we develop a mathematical pipeline that systematically analyses protein structures. We showcase this geometric framework on proteins forming open-ended trefoil knots, and we demonstrate that the mathematical tool, persistent homology, faithfully represents their structural homology. This topological pipeline identifies important geometric features of protein entanglement and clusters the space of trefoil proteins according to their depth. Persistence landscapes quantify the topological difference between a family of knotted and unknotted proteins in the same structural homology class. This difference is localized and interpreted geometrically with recent advancements in systematic computation of homology generators. The topological and geometric quantification we find is robust to noisy input data, which demonstrates the potential of this approach in contexts where standard knot theoretic tools fail.

AlphaFold predicts the most complex protein knot and composite protein knots

The computer artificial intelligence system AlphaFold has recently predicted previously unknown three-dimensional structures of thousands of proteins. Focusing on the subset with high-confidence scores, we algorithmically analyze these predictions for cases where the protein backbone exhibits rare topological complexity, that is, knotting. Amongst others, we discovered a 71 -knot, the most topologically complex knot ever found in a protein, as well several six-crossing composite knots comprised of two methyltransferase or carbonic anhydrase domains, each containing a simple trefoil knot. These deeply embedded composite knots occur evidently by gene duplication and interconnection of knotted dimers. Finally, we report two new five-crossing knots including the first 51 -knot. Our list of analyzed structures forms the basis for future experimental studies to confirm these novel-knotted topologies and to explore their complex folding mechanisms.

Investigation of the structural dynamics of a knotted protein and its unknotted analog using molecular dynamics

The role of knots in proteins remains elusive. Some studies suggest an impact on stability; the difficulty in comparing systems to assess this effect, however, has been a significant challenge. In this study, we produced and analyzed molecular dynamic trajectories considering three different temperatures of two variants of ornithine transcarbamylase (OTC), only one of which has a 3₁ knot, in order to evaluate the relative stability of the two molecules. RMSD showed equilibrated structures for the produced trajectories, and RMSF showed subtle differences in flexibility. In the knot moiety, the knotted protein did not show a great deal of fluctuation at any temperature. For the unknotted protein, the residue GLY243 showed a high fluctuation in the corresponding moiety. The fraction of native contacts (Q) showed a similar profile at all temperatures, with the greatest decrease by 436 K. The investigation of conformational behavior with principal component analysis (PCA) and dynamic cross-correlation map (DCCM) showed that knotted protein is less likely to undergo changes in its conformation under the conditions employed compared to unknotted. PCA data showed that the unknotted protein had greater dispersion in its conformations, which suggests that it has a greater capacity for conformation transitions in response to thermal changes. DCCM graphs comparing the 310 K and 436 K temperatures showed that the knotted protein had less change in its correlation and anti-correlation movements, indicating stability compared to the unknotted.

Efficient compressed database of equilibrated configurations of ring-linear polymer blends for MD simulations

To effectively archive configuration data during molecular dynamics (MD) simulations of polymer systems, we present an efficient compression method with good numerical accuracy that preserves the topology of ring-linear polymer blends. To compress the fraction of floating-point data, we used the Jointed Hierarchical Precision Compression Number - Data Format (JHPCN-DF) method to apply zero padding for the tailing fraction bits, which did not affect the numerical accuracy, then compressed the data with Huffman coding. We also provided a dataset of well-equilibrated configurations of MD simulations for ring-linear polymer blends with various lengths of linear and ring polymers, including ring complexes composed of multiple rings such as polycatenane. We executed 109 MD steps to obtain 150 equilibrated configurations. The combination of JHPCN-DF and SZ compression achieved the best compression ratio for all cases. Therefore, the proposed method enables efficient archiving of MD trajectories. Moreover, the publicly available dataset of ring-linear polymer blends can be employed for studies of mathematical methods, including topology analysis and data compression, as well as MD simulations.

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.

Unfolding and Translocation of Knotted Proteins by Clp Biological Nanomachines: Synergistic Contribution of Primary Sequence and Topology Revealed by Molecular Dynamics Simulations

We use Langevin dynamics simulations to model, at an atomistic resolution, how various natively knotted proteins are unfolded in repeated allosteric translocating cycles of the ClpY ATPase. We consider proteins representative of different topologies, from the simplest knot (trefoil 3₁), to the three-twist 5₂ knot, to the most complex stevedore, 6₁, knot. We harness the atomistic detail of the simulations to address aspects that have so far remained largely unexplored, such as sequence-dependent effects on the ruggedness of the landscape traversed during knot sliding. Our simulations reveal the combined effect on translocation of the knotted protein structure, i.e., backbone topology and geometry, and primary sequence, i.e., side chain size and interactions, and show that the latter can dominate translocation hindrance. In addition, we observe that due to the interplay between the knotted topology and intramolecular contacts the transmission of tension along the polypeptide chain occurs very differently from that of homopolymers. Finally, by considering native and non-native interactions, we examine how the disruption or formation of such contacts can affect the translocation processivity and concomitantly create multiple unfolding pathways with very different activation barriers.

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.

Searching the Optimal Folding Routes of a Complex Lasso Protein

Understanding how polypeptides can efficiently and reproducibly attain a self-entangled conformation is a compelling biophysical challenge that might shed new light on our general knowledge of protein folding. Complex lassos, namely self-entangled protein structures characterized by a covalent loop sealed by a cysteine bridge, represent an ideal test system in the framework of entangled folding. Indeed, because cysteine bridges form in oxidizing conditions, they can be used as on/off switches of the structure topology to investigate the role played by the backbone entanglement in the process. In this work, we have used molecular dynamics to simulate the folding of a complex lasso glycoprotein, granulocyte-macrophage colony-stimulating factor, modeling both reducing and oxidizing conditions. Together with a well-established Gō-like description, we have employed the elastic folder model, a coarse-grained, minimalistic representation of the polypeptide chain driven by a structure-based angular potential. The purpose of this study is to assess the kinetically optimal pathways in relation to the formation of the native topology. To this end, we have implemented an evolutionary strategy that tunes the elastic folder model potentials to maximize the folding probability within the early stages of the dynamics. The resulting protein model is capable of folding with high success rate, avoiding the kinetic traps that hamper the efficient folding in the other tested models. Employing specifically designed topological descriptors, we could observe that the selected folding routes avoid the topological bottleneck by locking the cysteine bridge after the topology is formed. These results provide valuable insights on the selection of mechanisms in self-entangled protein folding while, at the same time, the proposed methodology can complement the usage of established minimalistic models and draw useful guidelines for more detailed simulations.

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.

KymoKnot: A web server and software package to identify and locate knots in trajectories of linear or circular polymers

The KymoKnot software package and web server identifies and locates physical knots or proper knots in a series of polymer conformations. It is mainly intended as an analysis tool for trajectories of linear or circular polymers, but it can be used on single instances too, e.g. protein structures in PDB format. A key element of the software package is the so-called minimally interfering chain closure algorithm that is used to detect physical knots in open chains and to locate the knotted region in both open and closed chains. The web server offers a user-friendly graphical interface that identifies the knot type and highlights the knotted region on each frame of the trajectory, which the user can visualize interactively from various viewpoints. The dynamical evolution of the knotted region along the chain contour is presented as a kymograph. All data can be downloaded in text format. The KymoKnot package is licensed under the BSD 3-Clause licence. The server is publicly available at http://kymoknot.sissa.it/kymoknot/interactive.php .

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.

Molecular Knots

The first synthetic molecular trefoil knot was prepared in the late 1980s. However, it is only in the last few years that more complex small-molecule knot topologies have been realized through chemical synthesis. The steric restrictions imposed on molecular strands by knotting can impart significant physical and chemical properties, including chirality, strong and selective ion binding, and catalytic activity. As the number and complexity of accessible molecular knot topologies increases, it will become increasingly useful for chemists to adopt the knot terminology employed by other disciplines. Here we give an overview of synthetic strategies towards molecular knots and outline the principles of knot, braid, and tangle theory appropriate to chemistry and molecular structure.

Are There Knots in Chromosomes?

Recent developments have for the first time allowed the determination of three-dimensional structures of individual chromosomes and genomes in nuclei of single haploid mouse embryonic stem (ES) cells based on Hi⁻C chromosome conformation contact data. Although these first structures have a relatively low resolution, they provide the first experimental data that can be used to study chromosome and intact genome folding. Here we further analyze these structures and provide the first evidence that G1 phase chromosomes are knotted, consistent with the fact that plots of contact probability vs sequence separation show a power law dependence that is intermediate between that of a fractal globule and an equilibrium structure.

Pierced Lasso Topology Controls Function in Leptin

Protein engineering is a powerful tool in drug design and therapeutics, where disulphide bridges are commonly introduced to stabilize proteins. However, these bonds also introduce covalent loops, which are often neglected. These loops may entrap the protein backbone on opposite sides, leading to a "knotted" topology, forming a so-called Pierced Lasso (PL). In this elegant system, the "knot" is held together with a single disulphide bridge where part of the polypeptide chain is threaded through. The size and position of these covalent loops can be manipulated through protein design in vitro, whereas nature uses polymorphism to switch the PL topology. The PL protein leptin shows genetic modification of an N-terminal residue, adding a third cysteine to the same sequence. In an effort to understand the mechanism of threading of these diverse topologies, we designed three loop variants to mimic the polymorphic sequence. This adds elegance to the system under study, as it allows the generation of three possible covalent loops; they are the original wild-type C-terminal loop protein, the fully circularized unthreaded protein, and the N-terminal loop protein, responsible for different lasso topologies. The size of the loop changes the threading mechanism from a slipknotting to a plugging mechanism, with increasing loop size. Interestingly, the ground state of the native protein structure is largely unaffected, but biological assays show that the activity is maximized by properly controlled dynamics in the threaded state. A threaded topology with proper conformational dynamics is important for receptor interaction and activation of the signaling pathways in vivo.

Untangling the Influence of a Protein Knot on Folding

Entanglement and knots occur across all aspects of the physical world. Despite the common belief that knots are too complicated for incorporation into proteins, knots have been identified in the native fold of a growing number of proteins. The discovery of proteins with this unique backbone characteristic has challenged the preconceptions about the complexity of biological structures, as well as current folding theories. Given the intricacies of the knotted geometry, the interplay between a protein's fold, structure, and function is of particular interest. Interestingly, for most of these proteins, the knotted region appears critical both in folding and function, although full understanding of these contributions is still incomplete. Here, we experimentally reveal the impact of the knot on the landscape, the origin of the bistable nature of the knotted protein, and broaden the view of knot formation as uniquely decoupled from folding.

The elusive quest for RNA knots

Physical entanglement, and particularly knots arise spontaneously in equilibrated polymers that are sufficiently long and densely packed. Biopolymers are no exceptions: knots have long been known to occur in proteins as well as in encapsidated viral DNA. The rapidly growing number of RNA structures has recently made it possible to investigate the incidence of physical knots in this type of biomolecule, too. Strikingly, no knots have been found to date in the known RNA structures. In this Point of View Article we discuss the absence of knots in currently available RNAs and consider the reasons why knots in RNA have not yet been found, despite the expectation that they should exist in Nature. We conclude by singling out a number of RNA sequences that, based on the properties of their predicted secondary structures, are good candidates for knotted RNAs.

In Search of Functional Advantages of Knots in Proteins

We analysed the structure of deeply knotted proteins representing three unrelated families of knotted proteins. We looked at the correlation between positions of knotted cores in these proteins and such local structural characteristics as the number of intra-chain contacts, structural stability and solvent accessibility. We observed that the knotted cores and especially their borders showed strong enrichment in the number of contacts. These regions showed also increased thermal stability, whereas their solvent accessibility was decreased. Interestingly, the active sites within these knotted proteins preferentially located in the regions with increased number of contacts that also have increased thermal stability and decreased solvent accessibility. Our results suggest that knotting of polypeptide chains provides a favourable environment for the active sites observed in knotted proteins. Some knotted proteins have homologues without a knot. Interestingly, these unknotted homologues form local entanglements that retain structural characteristics of the knotted cores.

How to fold intricately: using theory and experiments to unravel the properties of knotted proteins

Over the years, advances in experimental and computational methods have helped us to understand the role of thermodynamic, kinetic and active (chaperone-aided) effects in coordinating the folding steps required to achieving a knotted native state. Here, we review such developments by paying particular attention to the complementarity of experimental and computational studies. Key open issues that could be tackled with either or both approaches are finally pointed out.

Folding of small knotted proteins: Insights from a mean field coarse-grained model

A small but relevant number of proteins whose native structure is known features nontrivial topology, i.e., they are knotted. Understanding the process of folding from a swollen unknotted state to the biologically relevant native conformation is, for these proteins, particularly difficult, due to their rate-limiting topological entanglement. To shed some light into this conundrum, we introduced a structure-based coarse-grained model of the protein, where the information about the folded conformation is encoded in bonded angular interactions only, which do not favor the formation of native contacts. A stochastic search scheme in parameter space is employed to identify a set of interactions that maximizes the probability to attain the knotted state. The optimal knotting pathways of the two smallest knotted proteins, obtained through this approach, are consistent with the results derived by means of coarse-grained as well as full atomistic simulations.

Using Local States To Drive the Sampling of Global Conformations in Proteins

Conformational changes associated with protein function often occur beyond the time scale currently accessible to unbiased molecular dynamics (MD) simulations, so that different approaches have been developed to accelerate their sampling. Here we investigate how the knowledge of backbone conformations preferentially adopted by protein fragments, as contained in precalculated libraries known as structural alphabets (SA), can be used to explore the landscape of protein conformations in MD simulations. We find that (a) enhancing the sampling of native local states in both metadynamics and steered MD simulations allows the recovery of global folded states in small proteins; (b) folded states can still be recovered when the amount of information on the native local states is reduced by using a low-resolution version of the SA, where states are clustered into macrostates; and (c) sequences of SA states derived from collections of structural motifs can be used to sample alternative conformations of preselected protein regions. The present findings have potential impact on several applications, ranging from protein model refinement to protein folding and design.

Discretized knot motion on a tensioned fiber induced by transverse waves

Topological entanglement is a ubiquitous feature of many biological as well as artificial polymers and fibers. While the equilibrium properties of entangled chains have been the subject of several studies, little is known about their out-of-equilibrium behavior. Here, we address the problem of a stretched knotted fiber driven by a periodic force applied to one of its termini. We show that the onset of standing waves kinetically traps the knot in spatially localized states where the amplitude of the oscillations is maximal, while the knot normal diffusive dynamics is replaced by a discrete jump dynamics.

Scaling behavior of topologically constrained polymer rings in a melt

Large scale molecular dynamics simulations on graphic processing units (GPUs) are employed to study the scaling behavior of ring polymers with various topological constraints in melts. Typical sizes of rings containing 3(1), 5(1) knots and catenanes made up of two unknotted rings scale like N(1/3) in the limit of large ring sizes N. This is consistent with the crumpled globule model and similar findings for unknotted rings. For small ring lengths knots occupy a significant fraction of the ring. The scaling of typical ring sizes for small N thus depends on the particular knot type and the exponent is generally larger than 0.4.

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.

Knotted proteins: A tangled tale of Structural Biology

Knotted proteins have their native structures arranged in the form of an open knot. In the last ten years researchers have been making significant efforts to reveal their folding mechanism and understand which functional advantage(s) knots convey to their carriers. Molecular simulations have been playing a fundamental role in this endeavor, and early computational predictions about the knotting mechanism have just been confirmed in wet lab experiments. Here we review a collection of simulation results that allow outlining the current status of the field of knotted proteins, and discuss directions for future research.

Two Adhesive Sites Can Enhance the Knotting Probability of DNA

Self-entanglement, or knotting, is entropically favored in long polymers. Relatively short polymers such as proteins can knot as well, but in this case the entanglement is mainly driven by fine-tuned, sequence-specific interactions. The relation between the sequence of a long polymer and its topological state is here investigated by means of a coarse-grained model of DNA. We demonstrate that the introduction of two adhesive regions along the sequence of a self-avoiding chain substantially increases the probability of forming a knot.

Absence of knots in known RNA structures

The ongoing effort to detect and characterize physical entanglement in biopolymers has so far established that knots are present in many globular proteins and also, abound in viral DNA packaged inside bacteriophages. RNA molecules, however, have not yet been systematically screened for the occurrence of physical knots. We have accordingly undertaken the systematic profiling of the several thousand RNA structures present in the Protein Data Bank (PDB). The search identified no more than three deeply knotted RNA molecules. These entries are rRNAs of about 3,000 nt solved by cryo-EM. Their genuine knotted state is, however, doubtful based on the detailed structural comparison with homologs of higher resolution, which are all unknotted. Compared with the case of proteins and viral DNA, the observed incidence of knots in available RNA structures is, therefore, practically negligible. This fact suggests that either evolutionary selection or thermodynamic and kinetic folding mechanisms act toward minimizing the entanglement of RNA to an extent that is unparalleled by other types of biomolecules. A possible general strategy for designing synthetic RNA sequences capable of self-tying in a twist-knot fold is finally proposed.

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.

An intramolecular lock facilitates folding and stabilizes the tertiary structure of Streptococcus mutans adhesin P1

The cariogenic bacterium Streptococcus mutans uses adhesin P1 to adhere to tooth surfaces, extracellular matrix components, and other bacteria. A composite model of P1 based on partial crystal structures revealed an unusual complex architecture in which the protein forms an elongated hybrid alpha/polyproline type II helical stalk by folding back on itself to display a globular head at the apex and a globular C-terminal region at the base. The structure of P1's N terminus and the nature of its critical interaction with the C-terminal region remained unknown, however. We have cocrystallized a stable complex of recombinant N- and C-terminal fragments and here describe a previously unidentified topological fold in which these widely discontinuous domains are intimately associated. The structure reveals that the N terminus forms a stabilizing scaffold by wrapping behind the base of P1's elongated stalk and physically "locking" it into place. The structure is stabilized through a highly favorable ΔG(solvation) on complex formation, along with extensive hydrogen bonding. We confirm the functional relevance of this intramolecular interaction using differential scanning calorimetry and circular dichroism to show that disruption of the proper spacing of residues 989-1001 impedes folding and diminishes stability of the full-length molecule, including the stalk. Our findings clarify previously unexplained functional and antigenic properties of P1.

How molecular knots can pass through each other

We propose a mechanism in which two molecular knots pass through each other and swap positions along a polymer strand. Associated free energy barriers in our simulations only amount to a few kBT, which may enable the interchange of knots on a single DNA strand.

Computational study on the progressive factorization of composite polymer knots into separated prime components

Using Monte Carlo simulations and advanced knot localization methods, we analyze the length and distribution of prime components in composite knots tied on freely jointed rings. For increasing contour length, we observe the progressive factorization of composite knots into separated prime components. However, we observe that a complete factorization, equivalent to the "decorated ring" picture, is not obtained even for rings of contour lengths N ≃ 3 N(0), about tens of times the most probable length of the prime knots tied on the rings. The decorated ring hypothesis has been used in the literature to justify the factorization of composite knot probabilities into the knotting probabilities of their prime components. Following our results, we suggest that such a hypothesis may not be necessary to explain the factorization of the knotting probabilities, at least when polymers excluding volume is not relevant. We rationalize the behavior of the system through a simple one-dimensional model in which prime knots are replaced by slip links randomly placed on a circle, with the only constraint being that the length of the loops has the same distribution as that of the length of the corresponding prime knots.

The role of non-native interactions in the folding of knotted proteins: insights from molecular dynamics simulations

For several decades, the presence of knots in naturally-occurring proteins was largely ruled out a priori for its supposed incompatibility with the efficiency and robustness of folding processes. For this very same reason, the later discovery of several unrelated families of knotted proteins motivated researchers to look into the physico-chemical mechanisms governing the concerted sequence of folding steps leading to the consistent formation of the same knot type in the same protein location. Besides experiments, computational studies are providing considerable insight into these mechanisms. Here, we revisit a number of such recent investigations within a common conceptual and methodological framework. By considering studies employing protein models with different structural resolution (coarse-grained or atomistic) and various force fields (from pure native-centric to realistic atomistic ones), we focus on the role of native and non-native interactions. For various unrelated instances of knotted proteins, non-native interactions are shown to be very important for favoring the emergence of conformations primed for successful self-knotting events.

Protein folding on the ribosome studied using NMR spectroscopy

NMR spectroscopy is a powerful tool for the investigation of protein folding and misfolding, providing a characterization of molecular structure, dynamics and exchange processes, across a very wide range of timescales and with near atomic resolution. In recent years NMR methods have also been developed to study protein folding as it might occur within the cell, in a de novo manner, by observing the folding of nascent polypeptides in the process of emerging from the ribosome during synthesis. Despite the 2.3 MDa molecular weight of the bacterial 70S ribosome, many nascent polypeptides, and some ribosomal proteins, have sufficient local flexibility that sharp resonances may be observed in solution-state NMR spectra. In providing information on dynamic regions of the structure, NMR spectroscopy is therefore highly complementary to alternative methods such as X-ray crystallography and cryo-electron microscopy, which have successfully characterized the rigid core of the ribosome particle. However, the low working concentrations and limited sample stability associated with ribosome-nascent chain complexes means that such studies still present significant technical challenges to the NMR spectroscopist. This review will discuss the progress that has been made in this area, surveying all NMR studies that have been published to date, and with a particular focus on strategies for improving experimental sensitivity.

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.

Influence of chain stiffness on knottedness in single polymers

In the present article, we investigate and review the influence of chain stiffness on self-entanglements and knots in a single polymer chain with Monte Carlo simulations spanning good solvent, theta and globular phases. The last-named are of particular importance as a model system for DNA in viral capsids. Intriguingly, the dependence of knot occurrence and complexity with increasing stiffness is non-trivial, but can be understood with a few simple concepts outlined in the present article.

Folding pathways of a knotted protein with a realistic atomistic force field

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native -sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting.

It has been recently observed that the native structure of proteins can contain knots. These are formed during the folding process and are tightened in a specific (i.e. native) location, along the poly-peptide chain. The existence of knots hence implies a high degree coordination of local and global conformational changes, during the folding reaction. In this work we investigate how the knot is formed and what are the dynamical mechanisms which drive the self-entanglement process. To this end, we report on the first atomistically detailed numerical simulation of the folding of a knotted protein, based on a realistic description of the inter-atomic forces. These simulations show that the knot is formed by following a specific sequence of contacts. The comparison of the findings with those based on simplified folding models suggest that the productive succession of contacts is aided by a concerted interplay of amino acid interactions, arguably including non-native ones.

Polymer uncrossing and knotting in protein folding, and their role in minimal folding pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold [Formula: see text], [Formula: see text], [Formula: see text], and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation "alignment". The consensus minimal pathway is constructed and shown schematically for representative cases of an [Formula: see text], [Formula: see text], and knotted protein. An overlap parameter is defined between pathways; we find that [Formula: see text] proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and [Formula: see text] proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.

The role of non-native interactions in the folding of knotted proteins

Stochastic simulations of coarse-grained protein models are used to investigate the propensity to form knots in early stages of protein folding. The study is carried out comparatively for two homologous carbamoyltransferases, a natively-knotted N-acetylornithine carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase (OTCase). In addition, two different sets of pairwise amino acid interactions are considered: one promoting exclusively native interactions, and the other additionally including non-native quasi-chemical and electrostatic interactions. With the former model neither protein shows a propensity to form knots. With the additional non-native interactions, knotting propensity remains negligible for the natively-unknotted OTCase while for AOTCase it is much enhanced. Analysis of the trajectories suggests that the different entanglement of the two transcarbamylases follows from the tendency of the C-terminal to point away from (for OTCase) or approach and eventually thread (for AOTCase) other regions of partly-folded protein. The analysis of the OTCase/AOTCase pair clarifies that natively-knotted proteins can spontaneously knot during early folding stages and that non-native sequence-dependent interactions are important for promoting and disfavouring early knotting events.

Knotted proteins provide an ideal ground for examining how amino acid interactions (which are local) can favor their folding into a native state of non-trivial topology (which is a global property). Some of the mechanisms that can aid knot formation are investigated here by comparing coarse-grained folding simulations of two enzymes that are structurally similar, and yet have natively knotted and unknotted states, respectively. In folding simulations that exclusively promote the formation of native contacts, neither protein forms knots. Strikingly, when sequence-dependent non-native interactions between amino acids are introduced, one observes knotting events but only for the natively-knotted protein. The results support the importance of non-native interactions in favoring or disfavoring knotting events in the early stages of folding.