AlphaKnot 2.0: a web server for the visualization of proteins' knotting and a database of knotted AlphaFold-predicted models

The availability of 3D protein models is rapidly increasing with the development of structure prediction algorithms. With the expanding availability of data, new ways of analysis, especially topological analysis, of those predictions are becoming necessary. Here, we present the updated version of the AlphaKnot service that provides a straightforward way of analyzing structure topology. It was designed specifically to determine knot types of the predicted structure models, however, it can be used for all structures, including the ones solved experimentally. AlphaKnot 2.0 provides the user's ability to obtain the knowledge necessary to assess the topological correctness of the model. Both probabilistic and deterministic knot detection methods are available, together with various visualizations (including a trajectory of simplification steps to highlight the topological complexities). Moreover, the web server provides a list of proteins similar to the queried model within AlphaKnot's database and returns their knot types for direct comparison. We pre-calculated the topology of high-quality models from the AlphaFold Database (4th version) and there are now more than 680.000 knotted models available in the AlphaKnot database. AlphaKnot 2.0 is available at https://alphaknot.cent.uw.edu.pl/.

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).

A Marcus-Type Inverted Region in the Translocation Kinetics of a Knotted Protein

Knotted proteins are rare but important species, yet how their complex topologies affect their physical properties is not fully understood. Here we combine single molecule nanopore experiments and all-atom MD simulations to study the electric-field-driven unfolding during the translocation through a model pore of individual protein knots important for methylating tRNA. One of these knots shows an unusual behavior that resembles the behavior of electrons hopping between two potential surfaces: as the electric potential driving the translocation reaction is increased, the rate eventually plateaus or slows back down in the "Marcus inverted regime". Our results shed light on the influence of topology in knotted proteins on their forced translocation through a pore connecting two electrostatic potential wells.

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.

In silico method for identifying the key residues in a knotted protein: with MJ0366 as an example

A simple in silico method for predicting the key residues for knotting and unknotting a knotted protein is put forward, with the residues ranked by the relevance to knotting and unknotting in the annealing molecular dynamics simulations.

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Factors underlying asymmetric pore dynamics of disaggregase and microtubule-severing AAA+ machines

Disaggregation and microtubule-severing nanomachines from the AAA+ (ATPases associated with various cellular activities) superfamily assemble into ring-shaped hexamers that enable protein remodeling by coupling large-scale conformational changes with application of mechanical forces within a central pore by loops protruding within the pore. We probed the asymmetric pore motions and intraring interactions that support them by performing extensive molecular dynamics simulations of single-ring severing proteins and the double-ring disaggregase ClpB. Simulations reveal that dynamic stability of hexameric pores of severing proteins and of the nucleotide-binding domain 1 (NBD1) ring of ClpB, which belong to the same clade, involves a network of salt bridges that connect conserved motifs of central pore loops. Clustering analysis of ClpB highlights correlated motions of domains of neighboring protomers supporting strong interprotomer collaboration. Severing proteins have weaker interprotomer coupling and stronger intraprotomer stabilization through salt bridges involving pore loops. Distinct mechanisms are identified in the NBD2 ring of ClpB involving weaker interprotomer coupling through salt bridges formed by noncanonical loops and stronger intraprotomer coupling. Analysis of collective motions of PL1 loops indicates that the largest amplitude motions in the spiral complex of spastin and ClpB involve axial excursions of the loops, whereas for katanin they involve opening and closing of the central pore. All three motors execute primarily axial excursions in the ring complex. These results suggest distinct substrate processing mechanisms of remodeling and translocation by ClpB and spastin compared to katanin, thus providing dynamic support for the differential action of the two severing proteins. Relaxation dynamics of the distance between the PL1 loops and the center of mass of protomers reveals observation-time-dependent dynamics, leading to predicted relaxation times of tens to hundreds of microseconds on millisecond experimental timescales. For ClpB, the predicted relaxation time is in excellent agreement with the extracted time from smFRET experiments.

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.

Transient knots in intrinsically disordered proteins and neurodegeneration

We provide a brief overview of the topological features found in structured proteins and of the dynamical processes that involve knots. We then discuss the knotted states that arise in the intrinsically disordered polyglutamine and α-synuclein. We argue that the existence of the knotted conformations stalls degradation by proteases and thus enhances aggregation. This mechanism works if the length of a peptide chain exceeds a threshold, as in the Huntington disease. We also study the cavities that form within the conformations of the disordered proteins. The volume of the cavities varies in time in a way that is different than that of the radius of gyration or the end-to-end distance. In addition, we study the traffic between the conformational basins and identify patterns associated with the deep and shallow knots. The results are obtained by molecular dynamics simulations that use coarse-grained and all-atom models (with and without the explicit solvent).

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.

Protein droplets in systems of disordered homopeptides and the amyloid glass phase

We study phase diagrams and liquid droplet formation for systems of polyglutamine and polyalanine chains in a coarse-grained model.

Protein knots provide mechano-resilience to an AAA+ protease-mediated proteolysis with profound ATP energy expenses

Knotted proteins are some of the most fascinating examples of how linear polypeptide chains can achieve intricate topological arrangements efficiently and spontaneously. The entanglements of polypeptide chains could potentially enhance their folding stabilities. We recently reported the unprecedented mechanostability of the Gordian (5₂) knotted family of human ubiquitin C-terminal hydrolases (UCHs) in the context of withstanding the mechanical unfolding of the bacterial AAA+ proteasome, ClpXP; a green fluorescence protein (GFP) was fused to the N-terminus of various UCHs as a reporter of the unfolding and degradation of these topologically knotted substrates, but it also limited the ability to examine the effect of untying the knotted topology via N-terminal truncation. In this study, we directly monitored the ClpXP-mediated degradation of UCH variants by electrophoresis and quantitative imaging analyses. We demonstrated that untying of the 5₂ knot in UCHL1 via N-terminal truncation (UCHL1_Δ11) significantly reduces its mechanostability. We further quantified the ATP expenditures of degrading different UCH variants by ClpXP. The unknotted UCHL1_Δ11 underwent accelerated ClpXP-dependent proteolysis, with a 30-fold reduction in ATP consumption compared to the knotted wild type. Unlike all other known ClpXP substrates, UCHL5, which is the most resilient substrate known to date, significantly slowed down the ATP turnover rate by ClpXP. Furthermore, UCHL5 required 1000-fold more ATP to be fully degraded by ClpXP compared to GFP. Our results underscored how the complex, knotted folding topology in UCHs may interfere with the mechano-unfolding processes of the AAA+ unfoldase, ClpX.

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.

Conformational State Hopping of Knots in Tensioned Polymer Chains

We use Brownian dynamics simulations to study the conformational states of knots on tensioned chains. Focusing specifically on the 8₁ knot, we observe knot conformational state hopping and show that the process can be described by a two-state kinetic model in the presence of an external force. The distribution of knot conformational states depends on the applied chain tension, which leads to a force-dependent distribution of knot untying pathways. We generalize our findings by considering the untying pathways of other knots and find that the way knots untie is generally governed by the force applied to the chain. From a broader perspective, being able to influence how a knot unties via external force can potentially be useful for applications of single-molecule techniques in which knots are unwanted.

The AAA+ protease ClpXP can easily degrade a 31 and a 52-knotted protein

Knots in proteins are hypothesized to make them resistant to enzymatic degradation by ATP-dependent proteases and recent studies have shown that whereas ClpXP can easily degrade a protein with a shallow 3₁ knot, it cannot degrade 5₂-knotted proteins if degradation is initiated at the C-terminus. Here, we present detailed studies of the degradation of both 3₁- and 5₂-knotted proteins by ClpXP using numerous constructs where proteins are tagged for degradation at both N- and C-termini. Our results confirm and extend earlier work and show that ClpXP can easily degrade a deeply 3₁-knotted protein. In contrast to recently published work on the degradation of 5₂-knotted proteins, our results show that the ClpXP machinery can also easily degrade these proteins. However, the degradation depends critically on the location of the degradation tag and the local stability near the tag. Our results are consistent with mechanisms in which either the knot simply slips along the polypeptide chain and falls off the free terminus, or one in which the tightened knot enters the translocation pore of ClpXP. Results of experiments on knotted protein fusions with a highly stable domain show partial degradation and the formation of degradation intermediates.

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.

Disordered peptide chains in an α-C-based coarse-grained model

We construct a one-bead-per-residue coarse-grained dynamical model to describe intrinsically disordered proteins at significantly longer timescales than in the all-atom models. In this model, inter-residue contacts form and disappear during the course of the time evolution. The contacts may arise between the sidechains, the backbones or the sidechains and backbones of the interacting residues. The model yields results that are consistent with many all-atom and experimental data on these systems. We demonstrate that the geometrical properties of various homopeptides differ substantially in this model. In particular, the average radius of gyration scales with the sequence length in a residue-dependent manner.

Topologically knotted deubiquitinases exhibit unprecedented mechanostability to withstand the proteolysis by an AAA+ protease

More than one thousand knotted protein structures have been identified so far, but the functional roles of these knots remain elusive. It has been postulated that backbone entanglement may provide additional mechanostability. Here, we employed a bacterial proteasome, ClpXP, to mechanically unfold 5₂-knotted human ubiquitin C-terminal hydrolase (UCH) paralogs from their C-termini, followed by processive translocation into the proteolytic chamber for degradation. Our results revealed unprecedentedly slow kinetics of ClpXP-mediated proteolysis for the proteasome-associated UCHL5: ten thousand times slower than that of a green fluorescence protein (GFP), which has a comparable size to the UCH domain but much higher chemical and thermal stabilities. The ClpXP-dependent mechanostability positively correlates with the intrinsic unfolding rates of the substrates, spanning over several orders of magnitude for the UCHs. The broad range of mechanostability within the same protein family may be associated with the functional requirements for their differential malleabilities.

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.

Knots can impair protein degradation by ATP-dependent proteases

ATP-dependent proteases translocate proteins through a narrow pore for their controlled destruction. However, how a protein substrate containing a knotted topology affects this process remains unknown. Here, we characterized the effects of the trefoil-knotted protein MJ0366 from Methanocaldococcus jannaschii on the operation of the ClpXP protease from Escherichia coli ClpXP completely degrades MJ0366 when pulling from the C-terminal ssrA-tag. However, when a GFP moiety is appended to the N terminus of MJ0366, ClpXP releases intact GFP with a 47-residue tail. The extended length of this tail suggests that ClpXP tightens the trefoil knot against GFP, which prevents GFP unfolding. Interestingly, if the linker between the knot core of MJ0366 and GFP is longer than 36 residues, ClpXP tightens and translocates the knot before it reaches GFP, enabling the complete unfolding and degradation of the substrate. These observations suggest that a knot-induced stall during degradation of multidomain proteins by AAA proteases may constitute a novel mechanism to produce partially degraded products with potentially new functions.

Pore translocation of knotted DNA rings

We use an accurate coarse-grained model for DNA and stochastic molecular dynamics simulations to study the pore translocation of 10-kbp-long DNA rings that are knotted. By monitoring various topological and physical observables we find that there is not one, as previously assumed, but rather two qualitatively different modes of knot translocation. For both modes the pore obstruction caused by knot passage has a brief duration and typically occurs at a late translocation stage. Both effects are well in agreement with experiments and can be rationalized with a transparent model based on the concurrent tensioning and sliding of the translocating knotted chains. We also observed that the duration of the pore obstruction event is more controlled by the knot translocation velocity than the knot size. These features should advance the interpretation and design of future experiments aimed at probing the spontaneous knotting of biopolymers.

Topological transformations in proteins: effects of heating and proximity of an interface

Using a structure-based coarse-grained model of proteins, we study the mechanism of unfolding of knotted proteins through heating. We find that the dominant mechanisms of unfolding depend on the temperature applied and are generally distinct from those identified for folding at its optimal temperature. In particular, for shallowly knotted proteins, folding usually involves formation of two loops whereas unfolding through high-temperature heating is dominated by untying of single loops. Untying the knots is found to generally precede unfolding unless the protein is deeply knotted and the heating temperature exceeds a threshold value. We then use a phenomenological model of the air-water interface to show that such an interface can untie shallow knots, but it can also make knots in proteins that are natively unknotted.

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.