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

Is There a Functional Role for the Knotted Topology in Protein UCH-L1?

Knotted proteins are present in nature, but there is still an open issue regarding the existence of a universal role for these remarkable structures. To address this question, we used classical molecular dynamics (MD) simulations combined with in vitro experiments to investigate the role of the Gordian knot in the catalytic activity of UCH-L1. To create an unknotted form of UCH-L1, we modified its amino acid sequence by truncating several residues from its N-terminus. Remarkably, we find that deleting the first two N-terminal residues leads to a partial loss of enzyme activity with conservation of secondary structural content and knotted topological state. This happens because the integrity of the N-terminus is critical to ensure the correct alignment of the catalytic triad. However, the removal of five residues from the N-terminus, which significantly disrupts the native structure and the topological state, leads to a complete loss of enzymatic activity. Overall, our findings indicate that UCH-L1's catalytic activity depends critically on the integrity of the N-terminus and the secondary structure content, with the latter being strongly coupled with the knotted topological state.

Everything AlphaFold tells us about protein knots

Recent advances in Machine Learning methods in structural biology opened up new perspectives for protein analysis. Utilizing these methods allows us to go beyond the limitations of empirical research, and take advantage of the vast amount of generated data. We use a complete set of potentially knotted protein models identified in all high-quality predictions from the AlphaFold Database to search for any common trends that describe them. We show that the vast majority of knotted proteins have 31 knot and that the presence of knots is preferred in neither Bacteria, Eukaryota, or Archaea domains. On the contrary, the percentage of knotted proteins in any given proteome is around 0.4%, regardless of the taxonomical group. We also verified that the organism's living conditions do not impact the number of knotted proteins in its proteome, as previously expected. We did not encounter an organism without a single knotted protein. What is more, we found four universally present families of knotted proteins in Bacteria, consisting of SAM synthase, and TrmD, TrmH, and RsmE methyltransferases.

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

Knotted artifacts in predicted 3D RNA structures

Unlike proteins, RNAs deposited in the Protein Data Bank do not contain topological knots. Recently, admittedly, the first trefoil knot and some lasso-type conformations have been found in experimental RNA structures, but these are still exceptional cases. Meanwhile, algorithms predicting 3D RNA models have happened to form knotted structures not so rarely. Interestingly, machine learning-based predictors seem to be more prone to generate knotted RNA folds than traditional methods. A similar situation is observed for the entanglements of structural elements. In this paper, we analyze all models submitted to the CASP15 competition in the 3D RNA structure prediction category. We show what types of topological knots and structure element entanglements appear in the submitted models and highlight what methods are behind the generation of such conformations. We also study the structural aspect of susceptibility to entanglement. We suggest that predictors take care of an evaluation of RNA models to avoid publishing structures with artifacts, such as unusual entanglements, that result from hallucinations of predictive algorithms.

Tying a true topological protein knot by cyclization

Knotted proteins are fascinating to biophysicists because of their robust ability to fold into intricately defined three-dimensional structures with complex and topologically knotted arrangements. Exploring the biophysical properties of the knotted proteins is of significant interest, as they could offer enhanced chemical, thermal, and mechanostabilities. A true mathematical knot requires a closed path; in contrast, knotted protein structures have open N- and C-termini. To address the question of how a truly knotted protein differs from the naturally occurring counterpart, we enzymatically cyclized a 31 knotted YibK protein from Haemophilus influenza (HiYibK) to investigate the impact of path closure on its structure-function relationship and folding stability. Through the use of a multitude of structural and biophysical tools, including X-ray crystallography, NMR spectroscopy, small angle X-ray scattering, differential scanning calorimetry, and isothermal calorimetry, we showed that the path closure minimally perturbs the native structure and ligand binding of HiYibK. Nevertheless, the cyclization did alter the folding stability and mechanism according to chemical and thermal unfolding analysis. These molecular insights contribute to our fundamental understanding of protein folding and knotting that could have implications in the protein design with higher stabilities.

Topological Reaction Coordinate Captures the Folding Transition State Ensemble in a Pierced Lasso Protein

Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (Go̅-like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.

AlphaFold Blindness to Topological Barriers Affects Its Ability to Correctly Predict Proteins' Topology

AlphaFold is a groundbreaking deep learning tool for protein structure prediction. It achieved remarkable accuracy in modeling many 3D structures while taking as the user input only the known amino acid sequence of proteins in question. Intriguingly though, in the early steps of each individual structure prediction procedure, AlphaFold does not respect topological barriers that, in real proteins, result from the reciprocal impermeability of polypeptide chains. This study aims to investigate how this failure to respect topological barriers affects AlphaFold predictions with respect to the topology of protein chains. We focus on such classes of proteins that, during their natural folding, reproducibly form the same knot type on their linear polypeptide chain, as revealed by their crystallographic analysis. We use partially artificial test constructs in which the mutual non-permeability of polypeptide chains should not permit the formation of complex composite knots during natural protein folding. We find that despite the formal impossibility that the protein folding process could produce such knots, AlphaFold predicts these proteins to form complex composite knots. Our study underscores the necessity for cautious interpretation and further validation of topological features in protein structures predicted by AlphaFold.

Folding kinetics of an entangled protein

The possibility of the protein backbone adopting lasso-like entangled motifs has attracted increasing attention. After discovering the surprising abundance of natively entangled protein domain structures, it was shown that misfolded entangled subpopulations might become thermosensitive or escape the homeostasis network just after translation. To investigate the role of entanglement in shaping folding kinetics, we introduce a novel indicator and analyze simulations of a coarse-grained, structure-based model for two small single-domain proteins. The model recapitulates the well-known two-state folding mechanism of a non-entangled SH3 domain. However, despite its small size, a natively entangled antifreeze RD1 protein displays a rich refolding behavior, populating two distinct kinetic intermediates: a short-lived, entangled, near-unfolded state and a longer-lived, non-entangled, near-native state. The former directs refolding along a fast pathway, whereas the latter is a kinetic trap, consistently with known experimental evidence of two different characteristic times. Upon trapping, the natively entangled loop folds without being threaded by the N-terminal residues. After trapping, the native entangled structure emerges by either backtracking to the unfolded state or threading through the already formed but not yet entangled loop. Along the fast pathway, trapping does not occur because the native contacts at the closure of the lasso-like loop fold after those involved in the N-terminal thread, confirming previous predictions. Despite this, entanglement may appear already in unfolded configurations. Remarkably, a longer-lived, near-native intermediate, with non-native entanglement properties, recalls what was observed in cotranslational folding.

Temperature dependent mechanical unfolding and refolding of a protein studied by thermo-regulated optical tweezers

Temperature is a useful system variable to gather kinetic and thermodynamic information from proteins. Usually, free energy and the associated entropic and enthalpic contributions are obtained by quantifying the conformational equilibrium based on melting experiments performed in bulk conditions. Such experiments are suitable only for those small single-domain proteins whose side reactions of irreversible aggregation are unlikely to occur. Here, we avoid aggregation by pulling single-protein molecules in a thermo-regulated optical tweezers. Thus, we are able to explore the temperature dependence of the thermodynamic and kinetic parameters of MJ0366 from Methanocaldococcus jannaschii at the single-molecule level. By performing force-ramp experiments between 2°C and 40°C, we found that MJ0366 has a nonlinear dependence of free energy with temperature and a specific heat change of 2.3 ± 1.2 kcal/mol∗K. These thermodynamic parameters are compatible with a two-state unfolding/refolding mechanism for MJ0366. However, the kinetics measured as a function of the temperature show a complex behavior, suggesting a three-state folding mechanism comprising a high-energy intermediate state. The combination of two perturbations, temperature and force, reveals a high-energy species in the folding mechanism of MJ0366 not detected in force-ramp experiments at constant temperature.

A Note on the Effects of Linear Topology Preservation in Monte Carlo Simulations of Knotted Proteins

Monte Carlo simulations are a powerful technique and are widely used in different fields. When applied to complex molecular systems with long chains, such as those in synthetic polymers and proteins, they have the advantage of providing a fast and computationally efficient way to sample equilibrium ensembles and calculate thermodynamic and structural properties under desired conditions. Conformational Monte Carlo techniques employ a move set to perform the transitions in the simulation Markov chain. While accepted conformations must preserve the sequential bonding of the protein chain model and excluded volume among its units, the moves themselves may take the chain across itself. We call this a break in linear topology preservation. In this manuscript, we show, using simple protein models, that there is no difference in equilibrium properties calculated with a move set that preserves linear topology and one that does not. However, for complex structures, such as those of deeply knotted proteins, the preservation of linear topology provides correct equilibrium results but only after long relaxation. In any case, to analyze folding pathways, knotting mechanisms and folding kinetics, the preservation of linear topology may be an unavoidable requirement.

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.

Revealing Topological Barriers against Knot Untying in Thermal and Mechanical Protein Unfolding by Molecular Dynamics Simulations

The knot is one of the most remarkable topological features identified in an increasing number of proteins with important functions. However, little is known about how the knot is formed during protein folding, and untied or maintained in protein unfolding. By means of all-atom molecular dynamics simulation, here we employ methyltransferase YbeA as the knotted protein model to analyze changes of the knotted conformation coupled with protein unfolding under thermal and mechanical denaturing conditions. Our results show that the trefoil knot in YbeA is occasionally untied via knot loosening rather than sliding under enhanced thermal fluctuations. Through correlating protein unfolding with changes in the knot position and size, several aspects of barriers that jointly suppress knot untying are revealed. In particular, protein unfolding is always prior to knot untying and starts preferentially from separation of two α-helices (α1 and α5), which protect the hydrophobic core consisting of β-sheets (β1-β4) from exposure to water. These β-sheets form a loop through which α5 is threaded to form the knot. Hydrophobic and hydrogen bonding interactions inside the core stabilize the loop against loosening. In addition, residues at N-terminal of α5 define a rigid turning to impede α5 from sliding out of the loop. Site mutations are designed to specifically eliminate these barriers, and easier knot untying is achieved under the same denaturing conditions. These results provide new molecular level insights into the folding/unfolding of knotted proteins.

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.

Protein self-entanglement modulates successful folding to the native state: A multi-scale modeling study

The computer-aided investigation of protein folding has greatly benefited from coarse-grained models, that is, simplified representations at a resolution level lower than atomistic, providing access to qualitative and quantitative details of the folding process that would be hardly attainable, via all-atom descriptions, for medium to long molecules. Nonetheless, the effectiveness of low-resolution models is itself hampered by the presence, in a small but significant number of proteins, of nontrivial topological self-entanglements. Features such as native state knots or slipknots introduce conformational bottlenecks, affecting the probability to fold into the correct conformation; this limitation is particularly severe in the context of coarse-grained models. In this work, we tackle the relationship between folding probability, protein folding pathway, and protein topology in a set of proteins with a nontrivial degree of topological complexity. To avoid or mitigate the risk of incurring in kinetic traps, we make use of the elastic folder model, a coarse-grained model based on angular potentials optimized toward successful folding via a genetic procedure. This light-weight representation allows us to estimate in silico folding probabilities, which we find to anti-correlate with a measure of topological complexity as well as to correlate remarkably well with experimental measurements of the folding rate. These results strengthen the hypothesis that the topological complexity of the native state decreases the folding probability and that the force-field optimization mimics the evolutionary process these proteins have undergone to avoid kinetic traps.

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.

A Specific Set of Heterogeneous Native Interactions Yields Efficient Knotting in Protein Folding

Native interactions are crucial for folding, and non-native interactions appear to be critical for efficiently knotting proteins. Therefore, it is important to understand both their roles in the folding of knotted proteins. It has been proposed that non-native interactions drive the correct order of contact formation, which is essential to avoid backtracking and efficiently self-tie. In this study, we ask if non-native interactions are strictly necessary to tangle a protein or if the correct order of contact formation can be assured by a specific set of native, but otherwise heterogeneous (i.e., having distinct energies), interactions. In order to address this problem, we conducted extensive Monte Carlo simulations of lattice models of protein-like sequences designed to fold into a preselected knotted conformation embedding a trefoil knot. We were able to identify a specific set of heterogeneous native interactions that drives efficient knotting and is able to fold the protein when combined with the remaining native interactions modeled as homogeneous. This specific set of heterogeneous native interactions is strictly enough to efficiently self-tie. A distinctive feature of these native interactions is that they do not backtrack because their energies ensure the correct order of contact formation. Furthermore, they stabilize a knotted intermediate state, which is en route to the native structure. Our results thus show that-at least in the context of the adopted model-non-native interactions are not necessary to knot a protein. However, when they are taken into account in protein energetics, it is possible to find specific, nonlocal non-native interactions that operate as a scaffold that assists the knotting step.

Tying up the Loose Ends: A Mathematically Knotted Protein

Knots have attracted scientists in mathematics, physics, biology, and engineering. Long flexible thin strings easily knot and tangle as experienced in our daily life. Similarly, long polymer chains inevitably tend to get trapped into knots. Little is known about their formation or function in proteins despite >1,000 knotted proteins identified in nature. However, these protein knots are not mathematical knots with their backbone polypeptide chains because of their open termini, and the presence of a “knot” depends on the algorithm used to create path closure. Furthermore, it is generally not possible to control the topology of the unfolded states of proteins, therefore making it challenging to characterize functional and physicochemical properties of knotting in any polymer. Covalently linking the amino and carboxyl termini of the deeply trefoil-knotted YibK from Pseudomonas aeruginosa allowed us to create the truly backbone knotted protein by enzymatic peptide ligation. Moreover, we produced and investigated backbone cyclized YibK without any knotted structure. Thus, we could directly probe the effect of the backbone knot and the decrease in conformational entropy on protein folding. The backbone cyclization did not perturb the native structure and its cofactor binding affinity, but it substantially increased the thermal stability and reduced the aggregation propensity. The enhanced stability of a backbone knotted YibK could be mainly originated from an increased ruggedness of its free energy landscape and the destabilization of the denatured state by backbone cyclization with little contribution from a knot structure. Despite the heterogeneity in the side-chain compositions, the chemically unfolded cyclized YibK exhibited several macroscopic physico-chemical attributes that agree with theoretical predictions derived from polymer physics.

Converging experimental and computational views of the knotting mechanism of a small knotted protein

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (31) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.

Mechanically tightening, untying and retying a protein trefoil knot by single-molecule force spectroscopy

Knotted conformation is one of the most surprising topological features found in proteins, and understanding the folding mechanism of such knotted proteins remains a challenge. Here, we used optical tweezers (OT) to investigate the mechanical unfolding and folding behavior of a knotted protein Escherichia coli tRNA (guanosine-1) methyltransferase (TrmD). We found that when stretched from its N- and C-termini, TrmD can be mechanically unfolded and stretched into a tightened trefoil knot, which is composed of ca. 17 residues. Stretching of the unfolded TrmD involved a compaction process of the trefoil knot at low forces. The unfolding pathways of the TrmD were bifurcated, involving two-state and three-state pathways. Upon relaxation, the tightened trefoil knot loosened up first, leading to the expansion of the knot, and the unfolded TrmD can then fold back to its native state efficiently. By using an engineered truncation TrmD variant, we stretched TrmD along a pulling direction to allow us to mechanically unfold TrmD and untie the trefoil knot. We found that the folding of TrmD from its unfolded polypeptide without the knot is significantly slower. The knotting is the rate-limiting step of the folding of TrmD. Our results highlighted the critical importance of the knot conformation for the folding and stability of TrmD, offering a new perspective to understand the role of the trefoil knot in the biological function of TrmD.

Optical tweezers are used to stretch a knotted protein along different directions to probe its unfolding–folding behaviors, and the conformational change of its knot structure.

Mechanical unfolding of a knotted protein unveils the kinetic and thermodynamic consequences of threading a polypeptide chain

Knots are remarkable topological features in nature. The presence of knots in crystallographic structures of proteins have stimulated considerable research to determine the kinetic and thermodynamic consequences of threading a polypeptide chain. By mechanically manipulating MJ0366, a small single domain protein harboring a shallow trefoil knot, we allow the protein to refold from either the knotted or the unknotted denatured state to characterize the free energy profile associated to both folding pathways. By comparing the stability of the native state with reference to the knotted and unknotted denatured state we find that knotting the polypeptide chain of MJ0366 increase the folding energy barrier in a magnitude close to the energy cost of forming a knot randomly in the denatured state. These results support that a protein knot can be formed during a single cooperative step of folding but occurs at the expenses of a large increment on the free energy barrier.

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.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

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.

Revealing Cooperation between Knotted Conformation and Dimerization in Protein Stabilization by Molecular Dynamics Simulations

The topological knot is thought to play a stabilizing role in maintaining the global fold and nature of proteins with the underlying mechanism yet to be elucidated. Given that most proteins containing trefoil knots exist and function as homodimers with a large part of the dimer interface occupied by the knotted region, we reason that the knotted conformation cooperates with dimerization in protein stabilization. Here, we take YbeA from Escherichia coli as the knotted protein model, using molecular dynamics (MD) simulations to compare the stability of two pairs of dimeric proteins having the same sequence and secondary structures but differing in the presence or absence of a trefoil knot in each subunit. The dimer interface of YbeA is identified to involve favorable contacts among three α-helices (α1, α3, and α5), one of which (α5) is threaded through a loop connected with α3 to form the knot. Upon removal of the knot by appropriate change of the knot-making crossing of the polypeptide chain, relevant domains are less constrained and exhibit enhanced fluctuations to decrease contacts at the interface. Unknotted subunits are less compact and undergo structural changes to ease the dimer separation. Such a stabilizing effect is evidenced by steered MD simulations, showing that the mechanical force required for dimer separation is significantly reduced by removing the knot. In addition to the knotted conformation, dimerization further improves the protein stability by restricting the α1-α5 separation, which is defined as a leading step for protein unfolding. These results provide important insights into the structure-function relationship of dimerization in knotted proteins.

Hydrophobic confinement modulates thermal stability and assists knotting in the folding of tangled proteins

There is growing support for the idea that the in vivo folding process of knotted proteins is assisted by chaperonins, but the mechanism of chaperonin assisted folding remains elusive. Here, we conduct extensive Monte Carlo simulations of lattice and off-lattice models to explore the effects of confinement and hydrophobic intermolecular interactions with the chaperonin cage in the folding and knotting processes. We find that moderate to high protein-cavity interactions (which are likely to be established in the beginning of the chaperonin working cycle) cause an energetic destabilization of the protein that overcomes the entropic stabilization driven by excluded volume, and leads to a decrease of the melting temperature relative to bulk conditions. Moreover, mild-to-moderate hydrophobic interactions with the cavity (which would be established later in the cycle) lead to a significant enhancement of knotting probability in relation to bulk conditions while simultaneously moderating the effect of steric confinement in the enhancement of thermal stability. Our results thus indicate that the chaperonin may be able to assist knotting without simultaneously thermally stabilizing potential misfolded states to a point that would hamper productive folding thus compromising its functional role.

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.

Topological descriptions of protein folding

How knotted proteins fold has remained controversial since the identification of deeply knotted proteins nearly two decades ago. Both computational and experimental approaches have been used to investigate protein knot formation. Motivated by the computer simulations of Bölinger et al. [Bölinger D, et al. (2010) PLoS Comput Biol 6:e1000731] for the folding of the [Formula: see text]-knotted α-haloacid dehalogenase (DehI) protein, we introduce a topological description of knot folding that could describe pathways for the formation of all currently known protein knot types and predicts knot types that might be identified in the future. We analyze fingerprint data from crystal structures of protein knots as evidence that particular protein knots may fold according to specific pathways from our theory. Our results confirm Taylor's twisted hairpin theory of knot folding for the [Formula: see text]-knotted proteins and the [Formula: see text]-knotted ketol-acid reductoisomerases and present alternative folding mechanisms for the [Formula: see text]-knotted phytochromes and the [Formula: see text]- and [Formula: see text]-knotted 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.

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.

Protein Knotting by Active Threading of Nascent Polypeptide Chain Exiting from the Ribosome Exit Channel

The mechanism of folding of deeply knotted proteins into their native structure is still not understood. Current thinking about protein folding is dominated by the Anfinsen dogma, stating that the structure of the folded proteins is uniquely dictated by the amino acid sequence of a given protein and that the folding is driven uniquely by the energy gained from interactions between amino acids that contact each other in the native structure of the protein. The role of ribosomes in protein folding was only seen as permitting the folding to progress from the N-terminal part of nascent protein chains. We propose here that ribosomes can participate actively in the folding of knotted proteins by actively threading nascent chains emerging from the ribosome exit channels through loops formed by a synthesized earlier portion of the same protein. Our simulations of folding of deeply knotted protein Tp0624 positively verify the proposed ribosome-driven active threading mechanism leading to the formation of deeply knotted proteins.

Structural entanglements in protein complexes

We consider multi-chain protein native structures and propose a criterion that determines whether two chains in the system are entangled or not. The criterion is based on the behavior observed by pulling at both termini of each chain simultaneously in the two chains. We have identified about 900 entangled systems in the Protein Data Bank and provided a more detailed analysis for several of them. We argue that entanglement enhances the thermodynamic stability of the system but it may have other functions: burying the hydrophobic residues at the interface and increasing the DNA or RNA binding area. We also study the folding and stretching properties of the knotted dimeric proteins MJ0366, YibK, and bacteriophytochrome. These proteins have been studied theoretically in their monomeric versions so far. The dimers are seen to separate on stretching through the tensile mechanism and the characteristic unraveling force depends on the pulling direction.

Self-consistent calculation of protein folding pathways

We introduce an iterative algorithm to efficiently simulate protein folding and other conformational transitions, using state-of-the-art all-atom force fields. Starting from the Langevin equation, we obtain a self-consistent stochastic equation of motion, which directly yields the reaction pathways. From the solution of this set of equations we derive a stochastic estimate of the reaction coordinate. We validate this approach against the results of plain MD simulations of the folding of a small protein, which were performed on the Anton supercomputer. In order to explore the computational efficiency of this algorithm, we apply it to generate a folding pathway of a protein that consists of 130 amino acids and has a folding rate of the order of s^-1.

Steric confinement and enhanced local flexibility assist knotting in simple models of protein folding

The chaperonin complex GroEL-GroES is able to accelerate the folding process of knotted proteins considerably. However, the folding mechanism inside the chaperonin cage is elusive. Here we use a combination of lattice and off-lattice Monte Carlo simulations of simple Gō models to study the effect of physical confinement and local flexibility on the folding process of protein model systems embedding a trefoil knot in their native structure. This study predicts that steric confinement plays a specific role in the folding of knotted proteins by increasing the knotting probability for very high degrees of confinement. This effect is observed for protein MJ0366 even above the melting temperature for confinement sizes compatible with the size of the GroEL/GroES chaperonin cage. An enhanced local flexibility produces the same qualitative effects on the folding process. In particular, we observe that knotting probability increases up to 40% in the transition state of protein MJ0366 when flexibility is enhanced. This is underlined by a structural change in the transition state, which becomes devoid of helical content. No relation between the knotting mechanism and flexibility was found in the context of the off-lattice model adopted in this work.

All-atom calculation of protein free-energy profiles

The Bias Functional (BF) approach is a variational method which enables one to efficiently generate ensembles of reactive trajectories for complex biomolecular transitions, using ordinary computer clusters. For example, this scheme was applied to simulate in atomistic detail the folding of proteins consisting of several hundreds of amino acids and with experimental folding time of several minutes. A drawback of the BF approach is that it produces trajectories which do not satisfy microscopic reversibility. Consequently, this method cannot be used to directly compute equilibrium observables, such as free energy landscapes or equilibrium constants. In this work, we develop a statistical analysis which permits us to compute the potential of mean-force (PMF) along an arbitrary collective coordinate, by exploiting the information contained in the reactive trajectories calculated with the BF approach. We assess the accuracy and computational efficiency of this scheme by comparing its results with the PMF obtained for a small protein by means of plain molecular dynamics.

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.

Knotting and unknotting proteins in the chaperonin cage: Effects of the excluded volume

Molecular dynamics simulations are used to explore the effects of chaperonin-like cages on knotted proteins with very low sequence similarity, different depths of a knot but with a similar fold, and the same type of topology. The investigated proteins are VirC2, DndE and MJ0366 with two depths of a knot. A comprehensive picture how encapsulation influences folding rates is provided based on the analysis of different cage sizes and temperature conditions. Neither of these two effects with regard to knotted proteins has been studied by means of molecular dynamics simulations with coarse-grained structure-based models before. We show that encapsulation in a chaperonin is sufficient to self-tie and untie small knotted proteins (VirC2, DndE), for which the equilibrium process is not accessible in the bulk solvent. Furthermore, we find that encapsulation reduces backtracking that arises from the destabilisation of nucleation sites, smoothing the free energy landscape. However, this effect can also be coupled with temperature rise. Encapsulation facilitates knotting at the early stage of folding and can enhance an alternative folding route. Comparison to unknotted proteins with the same fold shows directly how encapsulation influences the free energy landscape. In addition, we find that as the size of the cage decreases, folding times increase almost exponentially in a certain range of cage sizes, in accordance with confinement theory and experimental data for unknotted proteins.

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.

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.

Crystal structure analysis of a hypothetical protein (MJ0366) from Methanocaldococcus jannaschii revealed a novel topological arrangement of the knot fold

The crystal structure of a hypothetical protein MJ0366, derived from Methanocaldococcus jannaschii was solved at 1.9 Å resolution using synchrotron radiation. MJ0366 was crystallized as a monomer and has knot structural arrangement. Intriguingly, the solved structure consists of novel 'KNOT' fold conformation. The 3₁ trefoil knot was observed in the structure. The N-terminal and C-terminal ends did not participate in knot formation.

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.

Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent

Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.

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

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

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

Knots in proteins have been proposed to resist proteasomal degradation, thought in turn to be related to neurodegenerative diseases such as Huntington.