Energy landscape of knotted protein folding

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

CONTEXT

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

METHODS

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

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.

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.

Folding and functions of knotted proteins

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

Long-lived metastable knots in polyampholyte chains

Knots in proteins and DNA are known to have significant effect on their equilibrium and dynamic properties as well as on their function. While knot dynamics and thermodynamics in electrically neutral and uniformly charged polymer chains are relatively well understood, proteins are generally polyampholytes, with varied charge distributions along their backbones. Here we use simulations of knotted polymer chains to show that variation in the charge distribution on a polyampholyte chain with zero net charge leads to significant variation in the resulting knot dynamics, with some charge distributions resulting in long-lived metastable knots that escape the (open-ended) chain on a timescale that is much longer than that for knots in electrically neutral chains. The knot dynamics in such systems can be described, quantitatively, using a simple one-dimensional model where the knot undergoes biased Brownian motion along a "reaction coordinate", equal to the knot size, in the presence of a potential of mean force. In this picture, long-lived knots result from charge sequences that create large electrostatic barriers to knot escape. This model allows us to predict knot lifetimes even when those times are not directly accessible by simulations.

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

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

De novo protein folding on computers. Benefits and challenges

There has been recent success in prediction of the three-dimensional folded native structures of proteins, most famously by the AlphaFold Algorithm running on Google's/Alphabet's DeepMind computer. However, this largely involves machine learning of protein structures and is not a de novo protein structure prediction method for predicting three-dimensional structures from amino acid residue sequences. A de novo approach would be based almost entirely on general principles of energy and entropy that govern protein folding energetics, and importantly do so without the use of the amino acid sequences and structural features of other proteins. Most consider that problem as still unsolved even though it has occupied leading scientists for decades. Many consider that it remains one of the major outstanding issues in modern science. There is crucial continuing help from experimental findings on protein unfolding and refolding in the laboratory, but only to a limited extent because many researchers consider that the speed by which real proteins folds themselves, often from milliseconds to minutes, is itself still not fully understood. This is unfortunate, because a practical solution to the problem would probably have a major effect on personalized medicine, the pharmaceutical industry, biotechnology, and nanotechnology, including for example "smaller" tasks such as better modeling of flexible "unfolded" regions of the SARS-COV-2 spike glycoprotein when interacting with its cell receptor, antibodies, and therapeutic agents. Some important ideas from earlier studies are given before moving on to lessons from periodic and aperiodic crystals, and a possible role for quantum phenomena. The conclusion is that better computation of entropy should be the priority, though that is presented guardedly.

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.

Topological frustration leading to backtracking in a coupled folding-binding process

Intrinsically disordered proteins (IDPs) are abundant in all species. Their discovery challenges the traditional "sequence-structure-function" paradigm of protein science because IDPs play important roles in various biological processes without preformed folded structures. Bioinformatic analysis reveals that the intrinsically conformational disorder of IDPs as well as their conformational transition upon binding to their targets is encoded by their amino acid sequences. The rRNase domain of colicin E3 and the immunity protein Im3 are a pair of proteins involved in bacterial survival. While the N-terminal segment and the central segment of E3 make comparable intermolecular contacts with Im3 in the bound state, binding of E3 with Im3 is dominantly triggered by the central segment of E3. In this work, to further investigate the binding mechanism of disordered E3 with Im3, we performed systematic free energy and transition path analysis through coarse-grained molecular dynamics simulations. We observed backtracking of the N-terminal segment of E3 in the binding process, whose occurrence depends on salt concentration. Conformational analysis revealed that initial binding of the N-terminal segment of E3 to Im3 usually leads to misorientation of a central hairpin of E3 on Im3, which generates topological frustration and results in backtracking of the N-terminal segment. Our results not only provide deeper mechanistic insights into the coupled folding-binding process of the E3/Im3 complex, but also suggest that topological frustration could be present in the coupled folding-binding process of IDPs and play an important role in regulating the binding transition pathways.

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.

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.

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.

f -distance of knotoids and protein structure

Recent studies classify the topology of proteins by analysing the distribution of their projections using knotoids. The approximation of this distribution depends on the number of projection directions that are sampled. Here, we investigate the relation between knotoids differing only by small perturbations of the direction of projection. Since such knotoids are connected by at most a single forbidden move, we characterize forbidden moves in terms of equivariant band attachment between strongly invertible knots and of strand passages between θ -curves. This allows for the determination of the optimal sample size needed to produce a well-approximated knotoid distribution. Based on that and on topological properties of the distribution, we probe the depth of knotted proteins with the trefoil as the predominant knot type without using subchain analysis.

On topology and knotty entanglement in protein folding

We investigate aspects of topology in protein folding. For this we numerically simulate the temperature driven folding and unfolding of the slipknotted archaeal virus protein AFV3-109. Due to knottiness the (un)folding is a topological process, it engages the entire backbone in a collective fashion. Accordingly we introduce a topological approach to model the process. Our simulations reveal that the (un)folding of AFV3-109 slipknot proceeds through a folding intermediate that has the topology of a trefoil knot. We observe that the final slipknot causes a slight swelling of the folded AFV3-109 structure. We disclose the relative stability of the strands and helices during both the folding and unfolding processes. We confirm results from previous studies that pointed out that it can be very demanding to simulate the formation of knotty self-entanglement, and we explain how the problems are circumvented: The slipknotted AFV3-109 protein is a very slow folder with a topologically demanding pathway, which needs to be properly accounted for in a simulation description. When we either increase the relative stiffness of bending, or when we decrease the speed of ambient cooling, the rate of slipknot formation rapidly increases.

Topological Twists in Nature

The first entangled protein was observed about 30 years ago, resulting in an increased interest for uncovering the biological functions and biophysical properties of these complex topologies. Recently, the Pierced Lasso Topology (PLT) was discovered in which a covalent bond forms an intramolecular loop, leaving one or both termini free to pierce the loop. This topology is related to knots and other entanglements. PLTs exist in many well-researched systems where the PLTs have previously been unnoticed. PLTs represents 18% of all disulfide containing proteins across all kingdoms of life. In this review, we investigate the biological implications of this specific topology in which the PLT-forming disulfide may act as a molecular switch for protein function and consequently human health.

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.

Generalized Circuit Topology of Folded Linear Chains

A wide range of physical systems can be formally mapped to a linear chain of sorted objects. Upon introduction of intrachain interactions, such a chain can "fold" to elaborate topological structures, analogous to folded linear polymer systems. Two distinct chain-topology theories, knot theory and circuit topology, have separately provided insight into the structure, dynamics, and evolution of folded linear polymers such as proteins and genomic DNA. Knot theory, however, ignores intrachain interactions (contacts), whereas chain crossings are ignored in circuit topology. Thus, there is a need for a universal approach that can provide topological description of any folded linear chain. Here, we generalize circuit topology in order to grasp particularities typically addressed by knot theory. We develop a generic approach that is simple, mathematically rigorous, and practically useful for structural classification, analysis of structural dynamics, and engineering applications.

Circuit Topology Analysis of Polymer Folding Reactions

Circuit topology is emerging as a versatile measure to classify the internal structures of folded linear polymers such as proteins and nucleic acids. The topology framework can be applied to a wide range of problems, most notably molecular folding reactions that are central to biology and molecular engineering. In this Outlook, we discuss the state-of-the art of the technology and elaborate on the opportunities and challenges that lie ahead.

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

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

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

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

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.

Stabilizing Effect of Inherent Knots on Proteins Revealed by Molecular Dynamics Simulations

A growing number of proteins have been identified as knotted in their native structures, with such entangled topological features being expected to play stabilizing roles maintaining both the global fold and the nature of proteins. However, the molecular mechanism underlying the stabilizing effect is ambiguous. Here, we combine unbiased and mechanical atomistic molecular dynamics simulations to investigate how a protein is stabilized by an inherent knot by directly comparing chemical, thermal, and mechanical denaturing properties of two proteins having the same sequence and secondary structures but differing in the presence or absence of an inherent knot. One protein is YbeA from Escherichia coli, containing a deep trefoil knot within the sequence, and the other is the modified protein with the knot of YbeA being removed. Under certain chemical denaturing conditions, the unknotted protein fully unfolds whereas the knotted protein does not, suggesting a higher intrinsic stability for the protein having a knot. Both proteins unfold under enhanced thermal fluctuations but at different rates and with distinct pathways. Opening the hydrophobic core via separation between two α-helices is identified as a crucial step initiating the protein unfolding, which, however, is restrained for the knotted protein by topological and geometrical frustrations. Energy barriers for denaturing the protein are reduced by removing the knot, as evidenced by mechanical unfolding simulations. Finally, yet importantly, no obvious change in size or location of the knot was observed during denaturing processes, indicating that YbeA may remain knotted for a relatively long time during and after denaturation.

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.

The Pathways of the iRFP713 Unfolding Induced by Different Denaturants

Near-infrared fluorescent proteins (NIR FPs) based on the complexes of bacterial phytochromes with their natural biliverdin chromophore are widely used as genetically encoded optical probes for visualization of cellular processes and deep-tissue imaging of cells and organs in living animals. In this work, we show that the steady-state and kinetic dependencies of the various spectral characteristics of iRFP713, developed from the bacterial phytochrome RpBphP2 and recorded at protein unfolding induced by guanidine hydrochloride (GdnHCl), guanidine thiocyanate (GTC), and urea, differ substantially. A study of the unfolding of three single-tryptophan mutant forms of iRFP713 expectedly revealed that protein unfolding begins with the dissociation of the native dimer, while the monomers remain compact. A further increase in the denaturant concentration leads to the formation of an intermediate state of iRFP713 having hydrophobic areas exposed on the protein surface (I). The total surface charge of iRFP713 (pI 5.86) changes from negative to positive with an increase in the concentration of GdnHCl and GTC because the negative charge of glutamic and aspartic acids is neutralized by forming salt bridges between the carboxyl groups and GdnH⁺ ions and because the guanidinium cations bind to amide groups of glutamines and asparagines. The coincidence of both the concentration of the denaturants at which the intermediate state of iRFP713 accumulates and the concentration of GdnH⁺ ions at which the neutralization of the surface charge of the protein in this state is ensured results in strong protein aggregation. This is evidently realized by iRFP713 unfolding by GTC. At the unfolding of the protein by GdnHCl, an intermediate state is populated at higher denaturant concentrations and a strong aggregation is not observed. As expected, protein aggregates are not formed in the presence of the urea. The aggregation of the protein upon neutralization of the charge on the macromolecule surface is the main indicator of the intermediate state of protein. The unfolded state of iRFP713, whose formation is accompanied by a significant decrease in the parameter A, was found to have a different residual structure in the denaturants used.

Trefoil Knot Hydrodynamic Delocalization on Sheared Ring Polymers

The behavior of unknotted and trefoil-knotted ring polymers under shear flow is here examined by means of mesoscopic simulations. In contrast to most polymers, ring polymers in a hydrodynamic solvent at high shear rates do not get shortened in the vorticity direction. This is a consequence of the backflow produced by the interaction of the sheared solvent with the end-free polymer topology. The extended structures of the ring in the vorticity-flow plane, when they are aligned in a constant velocity plane, favor ring contour fluctuations. This variety of conformations largely suppresses the tank-treading type of rotation with extended conformations in favor of the tumbling type of rotations, where stretched and collapsed conformations alternate. The extension of trefoil knots is also enhanced, so that the knots become delocalized. We anticipate that these effects, which disappear in the absence of hydrodynamic interactions, will have a crucial impact on the rheological properties of concentrated ring solutions, and will also influence the behavior of more complicated systems such as mixtures of polymers with different topologies.

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.

Exploring the selectivity of inhibitor complexes with Bcl-2 and Bcl-XL: A molecular dynamics simulation approach

B-cell lymphoma 2 (Bcl-2) family proteins are potential drug targets in cancer and have a relatively flat and flexible binding site. ABT-199 is one of the most promising selective Bcl-2 inhibitors, and A-1155463 selectively inhibits Bcl-XL. Although the amino acid sequences of the binding sites of these two inhibitors are similar, the inhibitors selectively bind the target protein. In order to determine the origin of the selectivity of these inhibitors, we conducted molecular dynamics simulations using protein-inhibitor modeling. We confirmed that ASP103 of Bcl-2 is a key residue and that hydrogen bonding between ASP103 and ABT-199 confers the Bcl-2 selectivity of this inhibitor. For Bcl-XL selectivity, the secondary structure of α-helix 3 is a key factor. PHE105, SER106, and LEU108 in the loose α-helix 3 interact with A-1155463 to confer Bcl-XL selectivity. These findings provide important insights into the molecular mechanisms of selective inhibitors of Bcl-2 family proteins.

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.

Studies of global and local entanglements of individual protein chains using the concept of knotoids

We study here global and local entanglements of open protein chains by implementing the concept of knotoids. Knotoids have been introduced in 2012 by Vladimir Turaev as a generalization of knots in 3-dimensional space. More precisely, knotoids are diagrams representing projections of open curves in 3D space, in contrast to knot diagrams which represent projections of closed curves in 3D space. The intrinsic difference with classical knot theory is that the generalization provided by knotoids admits non-trivial topological entanglement of the open curves provided that their geometry is frozen as it is the case for crystallized proteins. Consequently, our approach doesn’t require the closure of chains into loops which implies that the geometry of analysed chains does not need to be changed by closure in order to characterize their topology. Our study revealed that the knotoid approach detects protein regions that were classified earlier as knotted and also new, topologically interesting regions that we classify as pre-knotted.

Accelerating physical simulations of proteins by leveraging external knowledge

It is challenging to compute structure-function relationships of proteins using molecular physics. The problem arises from the exponential scaling of the computational searching and sampling of large conformational spaces. This scaling challenge is not met by today's methods, such as Monte Carlo, simulated annealing, genetic algorithms, or molecular dynamics (MD) or its variants such as replica exchange. Such methods of searching for optimal states on complex probabalistic landscapes are referred to more broadly as Explore-and-Exploit (EE), including in contexts such as computational learning, games, industrial planning and modeling military strategies. Here we describe a Bayesian method, called MELD, that 'melds' together explore-and-exploit approaches with externally added information that can be vague, combinatoric, noisy, intuitive, heuristic, or from experimental data. MELD is shown to accelerate physical MD simulations when using experimental data to determine protein structures; for predicting protein structures by using heuristic directives; and when predicting binding affinities of proteins from limited information about the binding site. Such Guided Explore-and-Exploit approaches might also be useful beyond proteins and beyond molecular science.

Topological knots and links in proteins

Twenty years after their discovery, knots in proteins are now quite well understood. They are believed to be functionally advantageous and provide extra stability to protein chains. In this work, we go one step further and search for links-entangled structures, more complex than knots, which consist of several components. We derive conditions that proteins need to meet to be able to form links. We search through the entire Protein Data Bank and identify several sequentially nonhomologous chains that form a Hopf link and a Solomon link. We relate topological properties of these proteins to their function and stability and show that the link topology is characteristic of eukaryotes only. We also explain how the presence of links affects the folding pathways of proteins. Finally, we define necessary conditions to form Borromean rings in proteins and show that no structure in the Protein Data Bank forms a link of this type.

Pierced Lasso Topology Controls Function in Leptin

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

Protein knotting through concatenation significantly reduces folding stability

Concatenation by covalent linkage of two protomers of an intertwined all-helical HP0242 homodimer from Helicobacter pylori results in the first example of an engineered knotted protein. While concatenation does not affect the native structure according to X-ray crystallography, the folding kinetics is substantially slower compared to the parent homodimer. Using NMR hydrogen-deuterium exchange analysis, we showed here that concatenation destabilises significantly the knotted structure in solution, with some regions close to the covalent linkage being destabilised by as much as 5 kcal mol^-1. Structural mapping of chemical shift perturbations induced by concatenation revealed a pattern that is similar to the effect induced by concentrated chaotrophic agent. Our results suggested that the design strategy of protein knotting by concatenation may be thermodynamically unfavourable due to covalent constrains imposed on the flexible fraying ends of the template structure, leading to rugged free energy landscape with increased propensity to form off-pathway folding intermediates.

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.

In Search of Functional Advantages of Knots in Proteins

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

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

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

Methyl transfer by substrate signaling from a knotted protein fold

Proteins with knotted configurations, in comparison with unknotted proteins, are restricted in conformational space. Little is known regarding whether knotted proteins have sufficient dynamics to communicate between spatially separated substrate-binding sites. TrmD is a bacterial methyltransferase that uses a knotted protein fold to catalyze methyl transfer from S-adenosyl methionine (AdoMet) to G37-tRNA. The product, m1G37-tRNA, is essential for life and maintains protein-synthesis reading frames. Using an integrated approach of structural, kinetic, and computational analysis, we show that the structurally constrained TrmD knot is required for its catalytic activity. Unexpectedly, the TrmD knot undergoes complex internal movements that respond to AdoMet binding and signaling. Most of the signaling propagates the free energy of AdoMet binding, thereby stabilizing tRNA binding and allowing assembly of the active site. This work demonstrates new principles of knots as organized structures that capture the free energies of substrate binding and facilitate catalysis.

Folding analysis of the most complex Stevedore's protein knot

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

Stochastic Resonance in Protein Folding Dynamics

Although protein folding reactions are usually studied under static external conditions, it is likely that proteins fold in a locally fluctuating cellular environment in vivo. To mimic such behavior in in vitro experiments, the local temperature of the solvent can be modulated either harmonically or using correlated noise. In this study, coarse-grained molecular simulations are used to investigate these possibilities, and it is found that both periodic and correlated random fluctuations of the environment can indeed accelerate folding kinetics if the characteristic frequencies of the applied fluctuations are commensurate with the internal timescale of the folding reaction; this is consistent with the phenomenon of stochastic resonance observed in many other condensed-matter processes. To test this theoretical prediction, the folding dynamics of phosphoglycerate kinase under harmonic temperature fluctuations are experimentally probed using Förster resonance energy transfer fluorescence measurements. To analyze these experiments, a combination of theoretical approaches is developed, including stochastic simulations of folding kinetics and an analytical mean-field kinetic theory. The experimental observations are consistent with the theoretical predictions of stochastic resonance in phosphoglycerate kinase folding. When combined with an alternative experiment on the protein VlsE using a power spectrum analysis, elaborated in Dave et al., ChemPhysChem 2016, 10.1002/cphc.201501041, the overall data overwhelmingly point to the experimental confirmation of stochastic resonance in protein folding dynamics.

The Knotted Protein UCH-L1 Exhibits Partially Unfolded Forms under Native Conditions that Share Common Structural Features with Its Kinetic Folding Intermediates

The human ubiquitin C-terminal hydrolase, UCH-L1, is an abundant neuronal deubiquitinase that is associated with Parkinson's disease. It contains a complex Gordian knot topology formed by the polypeptide chain alone. Using a combination of fluorescence-based kinetic measurements, we show that UCH-L1 has two distinct kinetic folding intermediates that are transiently populated on parallel pathways between the denatured and native states. NMR hydrogen-deuterium exchange (HDX) experiments indicate the presence of partially unfolded forms (PUFs) of UCH-L1 under native conditions. HDX measurements as a function of urea concentration were used to establish the structure of the PUFs and pulse-labelled HDX NMR was used to show that the PUFs and the folding intermediates are likely the same species. In both cases, a similar stable core encompassing most of the central β-sheet is highly structured and α-helix 3, which is partially formed, packs against it. In contrast to the stable β-sheet core, the peripheral α-helices display significant local fluctuations leading to rapid exchange. The results also suggest that the main difference between the two kinetic intermediates is structure and packing of α-helices 3 and 7 and the degree of structure in β-strand 5. Together, the fluorescence and NMR results establish that UCH-L1 neither folds through a continuum of pathways nor by a single discrete pathway. Its folding is complex, the β-sheet core forms early and is present in both intermediate states, and the rate-limiting step which is likely to involve the threading of the chain to form the 52-knot occurs late on the folding pathway.