AlphaKnot: server to analyze entanglement in structures predicted by AlphaFold methods.

AlphaKnot is a server that measures entanglement in AlphaFold-solved protein models while considering pLDDT confidence values. AlphaKnot has two main functions: (i) providing researchers with a webserver for analyzing knotting in their own AlphaFold predictions and (ii) providing a database of knotting in AlphaFold predictions from the 21 proteomes for which models have been published prior to 2022. The knotting is defined in a probabilistic fashion. The knotting complexity of proteins is presented in the form of a matrix diagram which shows users the knot type for the entire polypeptide chain and for each of its subchains. The dominant knot types as well as the computed locations of the knot cores (i.e. minimal portions of protein backbones that form a given knot type) are shown for each protein structure. Based mainly on the pLDDT confidence values, entanglements are classified as Knots, Unsure, and Artifacts. The database portion of the server can be used, for example, to examine protein geometry and entanglement-function correlations, as a reference set for protein modeling, and for facilitating evolutional studies. The AlphaKnot server can be found at https://alphaknot.cent.uw.edu.pl/.

Mg2+-Dependent Methyl Transfer by a Knotted Protein: A Molecular Dynamics Simulation and Quantum Mechanics Study.

Mg2+ is required for the catalytic activity of TrmD, a bacteria-specific methyltransferase that is made up of a protein topological knot-fold, to synthesize methylated m1G37-tRNA to support life. However, neither the location of Mg2+ in the structure of TrmD nor its role in the catalytic mechanism is known. Using molecular dynamics (MD) simulations, we identify a plausible Mg2+ binding pocket within the active site of the enzyme, wherein the ion is coordinated by two aspartates and a glutamate. In this position, Mg2+ additionally interacts with the carboxylate of a methyl donor cofactor S-adenosylmethionine (SAM). The computational results are validated by experimental mutation studies, which demonstrate the importance of the Mg2+-binding residues for the catalytic activity. The presence of Mg2+ in the binding pocket induces SAM to adopt a unique bent shape required for the methyl transfer activity and causes a structural reorganization of the active site. Quantum mechanical calculations show that the methyl transfer is energetically feasible only when Mg2+ is bound in the position revealed by the MD simulations, demonstrating that its function is to align the active site residues within the topological knot-fold in a geometry optimal for catalysis. The obtained insights provide the opportunity for developing a strategy of antibacterial drug discovery based on targeting of Mg2+-binding to TrmD.

Restriction of S-adenosylmethionine conformational freedom by knotted protein binding sites.

S-adenosylmethionine (SAM) is one of the most important enzyme substrates. It is vital for the function of various proteins, including large group of methyltransferases (MTs). Intriguingly, some bacterial and eukaryotic MTs, while catalysing the same reaction, possess significantly different topologies, with the former being a knotted one. Here, we conducted a comprehensive analysis of SAM conformational space and factors that affect its vastness. We investigated SAM in two forms: free in water (via NMR studies and explicit solvent simulations) and bound to proteins (based on all data available in the PDB and on all-atom molecular dynamics simulations in water). We identified structural descriptors-angles which show the major differences in SAM conformation between unknotted and knotted methyltransferases. Moreover, we report that this is caused mainly by a characteristic for knotted MTs compact binding site formed by the knot and the presence of adenine-binding loop. Additionally, we elucidate conformational restrictions imposed on SAM molecules by other protein groups in comparison to conformational space in water.

KnotProt 2.0: a database of proteins with knots and other entangled structures.

The KnotProt 2.0 database (the updated version of the KnotProt database) collects information about proteins which form knots and other entangled structures. New features in KnotProt 2.0 include the characterization of both probabilistic and deterministic entanglements which can be formed by disulfide bonds and interactions via ions, a refined characterization of entanglement in terms of knotoids, the identification of the so-called cysteine knots, the possibility to analyze all or a non-redundant set of proteins, and various technical updates. The KnotProt 2.0 database classifies all entangled proteins, represents their complexity in the form of a knotting fingerprint, and presents many biological and geometrical statistics based on these results. Currently the database contains >2000 entangled structures, and it regularly self-updates based on proteins deposited in the Protein Data Bank (PDB).

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

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

TrmD: A Methyl Transferase for tRNA Methylation With m1G37.

TrmD is an S-adenosyl methionine (AdoMet)-dependent methyl transferase that synthesizes the methylated m1G37 in tRNA. TrmD is specific to and essential for bacterial growth, and it is fundamentally distinct from its eukaryotic and archaeal counterpart Trm5. TrmD is unusual by using a topological protein knot to bind AdoMet. Despite its restricted mobility, the TrmD knot has complex dynamics necessary to transmit the signal of AdoMet binding to promote tRNA binding and methyl transfer. Mutations in the TrmD knot block this intramolecular signaling and decrease the synthesis of m1G37-tRNA, prompting ribosomes to +1-frameshifts and premature termination of protein synthesis. TrmD is unique among AdoMet-dependent methyl transferases in that it requires Mg2+ in the catalytic mechanism. This Mg2+ dependence is important for regulating Mg2+ transport to Salmonella for survival of the pathogen in the host cell. The strict conservation of TrmD among bacterial species suggests that a better characterization of its enzymology and biology will have a broad impact on our understanding of bacterial pathogenesis.

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.

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.

KnotProt: a database of proteins with knots and slipknots.

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knot's core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.

Pierced Lasso Bundles are a new class of knot-like motifs.

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.

Identifying knots in proteins.

Polypeptide chains form open knots in many proteins. How these knotted proteins fold and finding the evolutionary advantage provided by these knots are among some of the key questions currently being studied in the protein folding field. The detection and identification of protein knots are substantial challenges. Different methods and many variations of them have been employed, but they can give different results for the same protein. In the present article, we review the various knot identification algorithms and compare their relative strengths when applied to the study of knots in proteins. We show that the statistical approach based on the uniform closure method is advantageous in comparison with other methods used to characterize protein knots.

The unique cysteine knot regulates the pleotropic hormone leptin.

Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.

Energy landscape of knotted protein folding.

Recent experiments have conclusively shown that proteins are able to fold from an unknotted, denatured polypeptide to the knotted, native state without the aid of chaperones. These experiments are consistent with a growing body of theoretical work showing that a funneled, minimally frustrated energy landscape is sufficient to fold small proteins with complex topologies. Here, we present a theoretical investigation of the folding of a knotted protein, 2ouf, engineered in the laboratory by a domain fusion that mimics an evolutionary pathway for knotted proteins. Unlike a previously studied knotted protein of similar length, we see reversible folding/knotting and a surprising lack of deep topological traps with a coarse-grained structure-based model. Our main interest is to investigate how evolution might further select the geometry and stiffness of the threading region of the newly fused protein. We compare the folding of the wild-type protein to several mutants. Similarly to the wild-type protein, all mutants show robust and reversible folding, and knotting coincides with the transition state ensemble. As observed experimentally, our simulations show that the knotted protein folds about ten times slower than an unknotted construct with an identical contact map. Simulated folding kinetics reflect the experimentally observed rollover in the folding limbs of chevron plots. Successful folding of the knotted protein is restricted to a narrow range of temperature as compared to the unknotted protein and fits of the kinetic folding data below folding temperature suggest slow, nondiffusive dynamics for the knotted protein.

Slipknotting upon native-like loop formation in a trefoil knot protein.

Protein knots and slipknots, mostly regarded as intriguing oddities, are gradually being recognized as significant structural motifs. Recent experimental results show that knotting, starting from a fully extended polypeptide, has not yet been observed. Understanding the nucleation process of folding knots is thus a natural challenge for both experimental and theoretical investigation. In this study, we employ energy landscape theory and molecular dynamics to elucidate the entire folding mechanism. The full free energy landscape of a knotted protein is mapped using an all-atom structure-based protein model. Results show that, due to the topological constraint, the protein folds through a three-state mechanism that contains (i) a precise nucleation site that creates a correctly twisted native loop (first barrier) and (ii) a rate-limiting free energy barrier that is traversed by two parallel knot-forming routes. The main route corresponds to a slipknot conformation, a collapsed configuration where the C-terminal helix adopts a hairpin-like configuration while threading, and the minor route to an entropically limited plug motion, where the extended terminus is threaded as through a needle. Knot formation is a late transition state process and results show that random (nonspecific) knots are a very rare and unstable set of configurations both at and below folding temperature. Our study shows that a native-biased landscape is sufficient to fold complex topologies and presents a folding mechanism generalizable to all known knotted protein topologies: knotting via threading a native-like loop in a preordered intermediate.

Mechanical strength of 17,134 model proteins and cysteine slipknots.

A new theoretical survey of proteins' resistance to constant speed stretching is performed for a set of 17,134 proteins as described by a structure-based model. The proteins selected have no gaps in their structure determination and consist of no more than 250 amino acids. Our previous studies have dealt with 7510 proteins of no more than 150 amino acids. The proteins are ranked according to the strength of the resistance. Most of the predicted top-strength proteins have not yet been studied experimentally. Architectures and folds which are likely to yield large forces are identified. New types of potent force clamps are discovered. They involve disulphide bridges and, in particular, cysteine slipknots. An effective energy parameter of the model is estimated by comparing the theoretical data on characteristic forces to the corresponding experimental values combined with an extrapolation of the theoretical data to the experimental pulling speeds. These studies provide guidance for future experiments on single molecule manipulation and should lead to selection of proteins for applications. A new class of proteins, involving cysteine slipknots, is identified as one that is expected to lead to the strongest force clamps known. This class is characterized through molecular dynamics simulations.

On the remarkable mechanostability of scaffoldins and the mechanical clamp motif.

Protein mechanostability is a fundamental biological property that can only be measured by single-molecule manipulation techniques. Such studies have unveiled a variety of highly mechanostable modules (mainly of the Ig-like, beta-sandwich type) in modular proteins subjected to mechanical stress from the cytoskeleton and the metazoan cell-cell interface. Their mechanostability is often attributed to a "mechanical clamp" of secondary structure (a patch of backbone hydrogen bonds) fastening their ends. Here we investigate the nanomechanics of scaffoldins, an important family of scaffolding proteins that assembles a variety of cellulases into the so-called cellulosome, a microbial extracellular nanomachine for cellulose adhesion and degradation. These proteins anchor the microbial cell to cellulose substrates, which makes their connecting region likely to be subjected to mechanical stress. By using single-molecule force spectroscopy based on atomic force microscopy, polyprotein engineering, and computer simulations, here we show that the cohesin I modules from the connecting region of cellulosome scaffoldins are the most robust mechanical proteins studied experimentally or predicted from the entire Protein Data Bank. The mechanostability of the cohesin modules studied correlates well with their mechanical kinetic stability but not with their thermal stability, and it is well predicted by computer simulations, even coarse-grained. This extraordinary mechanical stability is attributed to 2 mechanical clamps in tandem. Our findings provide the current upper limit of protein mechanostability and establish shear mechanical clamps as a general structural/functional motif widespread in proteins putatively subjected to mechanical stress. These data have important implications for the scaffoldin physiology and for protein design in biotechnology and nanotechnology.