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

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

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

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

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

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.

Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

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.

Stochastic model of translocation of knotted proteins

Knotted proteins, when forced through the pores, can get stuck if the knots in their backbone tighten under force. Alternatively, the knot can slide off the chain, making translocation possible. We construct a simple energy landscape model of this process with a time-periodic potential that mimics the action of a molecular motor. We calculate the translocation time as a function of the period of the pulling force, discuss the asymptotic limits and biological relevance of the results.

Optical tweezers across scales in cell biology

Optical tweezers (OT) provide a noninvasive approach for delivering minute physical forces to targeted objects. Controlling such forces in living cells or in vitro preparations allows for the measurement and manipulation of numerous processes relevant to the form and function of cells. As such, OT have made important contributions to our understanding of the structures of proteins and nucleic acids, the interactions that occur between microscopic structures within cells, the choreography of complex processes such as mitosis, and the ways in which cells interact with each other. In this review, we highlight recent contributions made to the field of cell biology using OT and provide basic descriptions of the physics, the methods, and the equipment that made these studies possible.

Conformational Freedom and Topological Confinement of Proteins in Biomolecular Condensates

The emergence of biomolecular condensation and liquid-liquid phase separation (LLPS) introduces a new layer of complexity into our understanding of cell and molecular biology. Evidence steadily grows indicating that condensates are not only implicated in physiology but also human disease. Macro- and mesoscale characterization of condensates as a whole have been instrumental in understanding their biological functions and dysfunctions. By contrast, the molecular level characterization of condensates and how condensates modify the properties of the molecules that constitute them thus far remain comparably scarce. In this minireview we summarize and discuss the findings of several recent studies that have focused on structure, dynamics, and interactions of proteins undergoing condensation. The mechanistic insights they provide help us identify the relevant properties nature and scientists can leverage to modulate the behavior of condensate systems. We also discuss the unique environment of the droplet surface and speculate on effects of topological constraints and physical exclusion on condensate properties.

Using Optical Tweezers to Monitor Allosteric Signals Through Changes in Folding Energy Landscapes

Signaling proteins are composed of conserved protein interaction domains that serve as allosteric regulatory elements of enzymatic or binding activities. The ubiquitous, structurally conserved cyclic nucleotide binding (CNB) domain is found covalently linked to proteins with diverse folds that perform multiple biological functions. Given that the structures of cAMP-bound CNB domains in different proteins are very similar, it remains a challenge to determine how this domain allosterically regulates such diverse protein functions and folds. Instead of a structural perspective, we focus our attention on the energy landscapes underlying the CNB domains and their responses to cAMP binding. We show that optical tweezers is an ideal tool to investigate how cAMP binding coupled to interdomain interactions remodel the energy landscape of the regulatory subunit of protein kinase A (PKA), which harbors two CNB domains. We mechanically manipulate and probe the unfolding and refolding behavior of the CNB domains as isolated structures or selectively as part of the PKA regulatory subunit, and in the presence and absence of cAMP. Optical tweezers allows us to dissect the changes in the energy landscape associated with cAMP binding, and to examine the allosteric interdomain interactions triggered by the cyclic nucleotide. This single molecule approach can be used to study other modular, multidomain signaling proteins found in nature.

Single-molecule optical tweezers reveals folding steps of the domain swapping mechanism of a protein

Domain swapping is a mechanism of protein oligomerization by which two or more subunits exchange structural elements to generate an intertwined complex. Numerous studies support a diversity of swapping mechanisms in which structural elements can be exchanged at different stages of the folding pathway of a subunit. Here, we used single-molecule optical tweezers technique to analyze the swapping mechanism of the forkhead DNA-binding domain of human transcription factor FoxP1. FoxP1 populates folded monomers in equilibrium with a swapped dimer. We generated a fusion protein linking two FoxP1 domains in tandem to obtain repetitive mechanical folding and unfolding trajectories. Thus, by stretching the same molecule several times, we detected either the independent folding of each domain or the elusive swapping step between domains. We found that a swapped dimer can be formed directly from fully or mostly folded monomer. In this situation, the interaction between the monomers in route to the domain-swapped dimer is the rate-limiting step. This approach is a useful strategy to test the different proposed domain swapping mechanisms for proteins with relevant physiological functions.

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.

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.