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

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.

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.

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.

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.

A Topological Selection of Folding Pathways from Native States of Knotted Proteins

Understanding how knotted proteins fold is a challenging problem in biology. Researchers have proposed several models for their folding pathways, based on theory, simulations and experiments. The geometry of proteins with the same knot type can vary substantially and recent simulations reveal different folding behaviour for deeply and shallow knotted proteins. We analyse proteins forming open-ended trefoil knots by introducing a topologically inspired statistical metric that measures their entanglement. By looking directly at the geometry and topology of their native states, we are able to probe different folding pathways for such proteins. In particular, the folding pathway of shallow knotted carbonic anhydrases involves the creation of a double-looped structure, contrary to what has been observed for other knotted trefoil proteins. We validate this with Molecular Dynamics simulations. By leveraging the geometry and local symmetries of knotted proteins’ native states, we provide the first numerical evidence of a double-loop folding mechanism in trefoil proteins.

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.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

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.

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 entanglement of interlocked knotted-unknotted polymer rings

Topological entanglements in biopolymers could drive them to certain internal statics and dynamics with important implications for biological functions. In this study, by means of molecular dynamics simulations, we demonstrate that the minimal crossing pattern of a braid plays a major role in its structural and dynamical properties; the braid consists of a knotted ring and an interlocked entwined unknotted polymer ring. In particular, we show that depending on the bending rigidity of the chains, the conformational energy of the braid can be either lower or higher than the unlocked polymer rings. Additionally, we find that a non-identical crossing pattern in the braid could distinctly enforce concerted internal conformational fluctuations between the interlocked rings.

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.

Optimal Coarse-Grained Site Selection in Elastic Network Models of Biomolecules

Elastic network models, simple structure-based representations of biomolecules where atoms interact via short-range harmonic potentials, provide great insight into a molecule's internal dynamics and mechanical properties at extremely low computational cost. Their efficiency and effectiveness have made them a pivotal instrument in the computer-aided study of proteins and, since a few years, also of nucleic acids. In general, the coarse-grained sites, i.e. those effective force centers onto which the all-atom structure is mapped, are constructed based on intuitive rules: a typical choice for proteins is to retain only the C _α atoms of each amino acid. However, a mapping strategy relying only on the atom type and not the local properties of its embedding can be suboptimal compared to a more careful selection. Here, we present a strategy in which the subset of atoms, each of which is mapped onto a unique coarse-grained site of the model, is selected in a stochastic search aimed at optimizing a cost function. The latter is taken to be a simple measure of the consistency between the harmonic approximation of an elastic network model and the harmonic model obtained through exact integration of the discarded degrees of freedom. The method is applied to two representatives of structurally very different types of biomolecules: the protein adenylate kinase and the RNA molecule adenine riboswitch. Our analysis quantifies the substantial impact that an algorithm-driven selection of coarse-grained sites can have on a model's properties.

Topology of polymer chains under nanoscale confinement

Spatial confinement limits the conformational space accessible to biomolecules but the implications for bimolecular topology are not yet known. Folded linear biopolymers can be seen as molecular circuits formed by intramolecular contacts. The pairwise arrangement of intra-chain contacts can be categorized as parallel, series or cross, and has been identified as a topological property. Using molecular dynamics simulations, we determine the contact order distributions and topological circuits of short semi-flexible linear and ring polymer chains with a persistence length of l_p under a spherical confinement of radius R_c. At low values of l_p/R_c, the entropy of the linear chain leads to the formation of independent contacts along the chain and accordingly, increases the fraction of series topology with respect to other topologies. However, at high l_p/R_c, the fraction of cross and parallel topologies are enhanced in the chain topological circuits with cross becoming predominant. At an intermediate confining regime, we identify a critical value of l_p/R_c, at which all topological states have equal probability. Confinement thus equalizes the probability of more complex cross and parallel topologies to the level of the more simple, non-cooperative series topology. Moreover, our topology analysis reveals distinct behaviours for ring- and linear polymers under weak confinement; however, we find no difference between ring- and linear polymers under strong confinement. Under weak confinement, ring polymers adopt parallel and series topologies with equal likelihood, while linear polymers show a higher tendency for series arrangement. The radial distribution analysis of the topology reveals a non-uniform effect of confinement on the topology of polymer chains, thereby imposing more pronounced effects on the core region than on the confinement surface. Additionally, our results reveal that over a wide range of confining radii, loops arranged in parallel and cross topologies have nearly the same contact orders. Such degeneracy implies that the kinetics and transition rates between the topological states cannot be solely explained by contact order. We expect these findings to be of general importance in understanding chaperone assisted protein folding, chromosome architecture, and the evolution of molecular folds.

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.

Topology of internally constrained polymer chains

An interacting partner can provide external control over folding rates and realized topologies.

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.

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.