Periodic forces trigger knot untying during translocation of knotted proteins

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.

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.

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

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

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Unfolding and Translocation of Knotted Proteins by Clp Biological Nanomachines: Synergistic Contribution of Primary Sequence and Topology Revealed by Molecular Dynamics Simulations

We use Langevin dynamics simulations to model, at an atomistic resolution, how various natively knotted proteins are unfolded in repeated allosteric translocating cycles of the ClpY ATPase. We consider proteins representative of different topologies, from the simplest knot (trefoil 3₁), to the three-twist 5₂ knot, to the most complex stevedore, 6₁, knot. We harness the atomistic detail of the simulations to address aspects that have so far remained largely unexplored, such as sequence-dependent effects on the ruggedness of the landscape traversed during knot sliding. Our simulations reveal the combined effect on translocation of the knotted protein structure, i.e., backbone topology and geometry, and primary sequence, i.e., side chain size and interactions, and show that the latter can dominate translocation hindrance. In addition, we observe that due to the interplay between the knotted topology and intramolecular contacts the transmission of tension along the polypeptide chain occurs very differently from that of homopolymers. Finally, by considering native and non-native interactions, we examine how the disruption or formation of such contacts can affect the translocation processivity and concomitantly create multiple unfolding pathways with very different activation barriers.

RNA Pore Translocation with Static and Periodic Forces: Effect of Secondary and Tertiary Elements on Process Activation and Duration

We use MD simulations to study the pore translocation properties of a pseudoknotted viral RNA. We consider the 71-nucleotide-long xrRNA from the Zika virus and establish how it responds when driven through a narrow pore by static or periodic forces applied to either of the two termini. Unlike the case of fluctuating homopolymers, the onset of translocation is significantly delayed with respect to the application of static driving forces. Because of the peculiar xrRNA architecture, activation times can differ by orders of magnitude at the two ends. Instead, translocation duration is much smaller than activation times and occurs on time scales comparable at the two ends. Periodic forces amplify significantly the differences at the two ends, for both activation times and translocation duration. Finally, we use a waiting-times analysis to examine the systematic slowing downs in xrRNA translocations and associate them to the hindrance of specific secondary and tertiary elements of xrRNA. The findings provide a useful reference to interpret and design future theoretical and experimental studies of RNA translocation.

Native-state fingerprint on the ubiquitin translocation across a nanopore

We study the translocation of the ubiquitin molecule (Ubq) across a channel with a double section which constitutes a general feature of several transmembrane nanopores such as the α-hemolysin (αHL). Our purpose is to establish the structure-dependent character of the Ubq translocation pathway. This implies to find the correspondence, if any, between the translocational unfolding steps and the Ubq native state. For this reason, it is convenient to apply a coarse-grained computational approach, where the protein is described only by the backbone and the force field only exploits the information contained in the native state (in the spirit of Gō-like models, or native-centric models). The αHL-like pore is portrayed as two coaxial confining cylinders: a larger one for the vestibule and a narrower one for the barrel (or stem). Such simplified approach allows a large number of translocation events to be collected by limited computational resources. The co-translocational unfolding of Ubq is described via a few collective variables that characterize the translocation progress. We find two translocation intermediates (stalled conformations) that can be associated with specific unfolding stages. In particular, in the earliest step, the strand S5 unfolds and enters the pore. This step splits the native conformation into two structural clusters packing against each other in the Ubq fold. A second stall occurs when the hairpin of the N terminal engages the stem region.

Protein knots provide mechano-resilience to an AAA+ protease-mediated proteolysis with profound ATP energy expenses

Knotted proteins are some of the most fascinating examples of how linear polypeptide chains can achieve intricate topological arrangements efficiently and spontaneously. The entanglements of polypeptide chains could potentially enhance their folding stabilities. We recently reported the unprecedented mechanostability of the Gordian (5₂) knotted family of human ubiquitin C-terminal hydrolases (UCHs) in the context of withstanding the mechanical unfolding of the bacterial AAA+ proteasome, ClpXP; a green fluorescence protein (GFP) was fused to the N-terminus of various UCHs as a reporter of the unfolding and degradation of these topologically knotted substrates, but it also limited the ability to examine the effect of untying the knotted topology via N-terminal truncation. In this study, we directly monitored the ClpXP-mediated degradation of UCH variants by electrophoresis and quantitative imaging analyses. We demonstrated that untying of the 5₂ knot in UCHL1 via N-terminal truncation (UCHL1_Δ11) significantly reduces its mechanostability. We further quantified the ATP expenditures of degrading different UCH variants by ClpXP. The unknotted UCHL1_Δ11 underwent accelerated ClpXP-dependent proteolysis, with a 30-fold reduction in ATP consumption compared to the knotted wild type. Unlike all other known ClpXP substrates, UCHL5, which is the most resilient substrate known to date, significantly slowed down the ATP turnover rate by ClpXP. Furthermore, UCHL5 required 1000-fold more ATP to be fully degraded by ClpXP compared to GFP. Our results underscored how the complex, knotted folding topology in UCHs may interfere with the mechano-unfolding processes of the AAA+ unfoldase, ClpX.

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.

Stokesian dynamics of sedimenting elastic rings

We consider elastic microfilaments which form closed loops. We investigate how the loops change shape and orientation while settling under gravity in a viscous fluid. Loops are circular at the equilibrium. Their dynamics are investigated numerically based on the Stokes equations for the fluid motion and the bead-spring model of the microfilament. The Rotne-Prager approximation for the bead mobility is used. We demonstrate that the relevant dimensionless parameter is the ratio of the bending resistance of the filament to the gravitation force corrected for buoyancy. The inverse of this ratio, called the elasto-gravitation number B, is widely used in the literature for sedimenting elastic linear filaments. We assume that B is of the order of 104-106, which corresponds to easily deformable loops. We find out that initially tilted circles evolve towards different sedimentation modes, depending on B. Very stiff or stiff rings attain almost planar, oval shapes, which are vertical or tilted, respectively. More flexible loops deform significantly and converge towards one of several characteristic periodic motions. These sedimentation modes are also detected when starting from various shapes, and for different loop lengths. In general, multi-stability is observed: an elastic ring converges to one of several sedimentation modes, depending on the initial conditions. This effect is pronounced for very elastic loops. The surprising diversity of long-lasting periodic motions and shapes of elastic rings found in this work gives a new perspective for the dynamics of more complex deformable objects at micrometer and nanometer scales, sedimenting under gravity or rotating in a centrifuge, such as red blood cells, ring polymers or circular DNA.

Observations of metastable states of the free swelling knots and directional motion of tensioned knots in vibrated bead chains

The free swelling of knots and the directional motion of knots under tension were studied in vertically vibrated bead chains. A metastable state of swelling was observed in the strongly vibrated two-end-free bead chains, as predicted by Grosberg and Rabin (Phys. Rev. Lett. 99, 217801 (2007)). Knots in the two-end-fixed chains were found to move directionally. The direction of motion could be changed by flipping the knot over. The velocity of motion depended on the tension in the bead chain. The effects of tension on the motion of knots were studied in one-end-fixed chains. The directional reptation might have been influenced by the random motion of the leading arc of the knot. The knots might move in a forced-reptation manner under the interaction with a simulated translocase.

Insights into protein sequencing with an α-Hemolysin nanopore by atomistic simulations

Single molecule protein sequencing would represent a disruptive burst in proteomic research with important biomedical impacts. Due to their success in DNA sequencing, nanopore based devices have been recently proposed as possible tools for the sequencing of peptide chains. One of the open questions in nanopore protein sequencing concerns the ability of such devices to provide different signals for all the 20 standard amino acids. Here, using equilibrium all-atom molecular dynamics simulations, we estimated the pore clogging in α-Hemolysin nanopore associated to 20 different homopeptides, one for each standard amino acid. Our results show that pore clogging is affected by amino acid volume, hydrophobicity and net charge. The equilibrium estimations are also supported by non-equilibrium runs for calculating the current blockades for selected homopeptides. Finally, we discuss the possibility to modify the α-Hemolysin nanopore, cutting a portion of the barrel region close to the trans side, to reduce spurious signals and, hence, to enhance the sensitivity of the nanopore.

The AAA+ protease ClpXP can easily degrade a 31 and a 52-knotted protein

Knots in proteins are hypothesized to make them resistant to enzymatic degradation by ATP-dependent proteases and recent studies have shown that whereas ClpXP can easily degrade a protein with a shallow 3₁ knot, it cannot degrade 5₂-knotted proteins if degradation is initiated at the C-terminus. Here, we present detailed studies of the degradation of both 3₁- and 5₂-knotted proteins by ClpXP using numerous constructs where proteins are tagged for degradation at both N- and C-termini. Our results confirm and extend earlier work and show that ClpXP can easily degrade a deeply 3₁-knotted protein. In contrast to recently published work on the degradation of 5₂-knotted proteins, our results show that the ClpXP machinery can also easily degrade these proteins. However, the degradation depends critically on the location of the degradation tag and the local stability near the tag. Our results are consistent with mechanisms in which either the knot simply slips along the polypeptide chain and falls off the free terminus, or one in which the tightened knot enters the translocation pore of ClpXP. Results of experiments on knotted protein fusions with a highly stable domain show partial degradation and the formation of degradation intermediates.

Periodic Motion of Sedimenting Flexible Knots

We study the dynamics of knotted deformable closed chains sedimenting in a viscous fluid. We show experimentally that trefoil and other torus knots often attain a remarkably regular horizontal toroidal structure while sedimenting, with a number of intertwined loops, oscillating periodically around each other. We then recover this motion numerically and find out that it is accompanied by a very slow rotation around the vertical symmetry axis. We analyze the dependence of the characteristic timescales on the chain flexibility and aspect ratio. It is observed in the experiments that this oscillating mode of the dynamics can spontaneously form even when starting from a qualitatively different initial configuration. In numerical simulations, the oscillating modes are usually present as transients or final stages of the evolution, depending on chain aspect ratio and flexibility, and the number of loops.

Polymer translocation under a pulling force: Scaling arguments and threshold forces

DNA translocation through nanopores is one of the most promising strategies for next-generation sequencing technologies. Most experimental and numerical works have focused on polymer translocation biased by electrophoresis, where a pulling force acts on the polymer within the nanopore. An alternative strategy, however, is emerging, which uses optical or magnetic tweezers. In this case, the pulling force is exerted directly at one end of the polymer, which strongly modifies the translocation process. In this paper, we report numerical simulations of both linear and structured (mimicking DNA) polymer models, simple enough to allow for a statistical treatment of the pore structure effects on the translocation time probability distributions. Based on extremely extended computer simulation data, we (i) propose scaling arguments for an extension of the predicted translocation times τ∼N^{2}F^{-1} over the moderate forces range and (ii) analyze the effect of pore size and polymer structuration on translocation times τ.

Knots can impair protein degradation by ATP-dependent proteases

ATP-dependent proteases translocate proteins through a narrow pore for their controlled destruction. However, how a protein substrate containing a knotted topology affects this process remains unknown. Here, we characterized the effects of the trefoil-knotted protein MJ0366 from Methanocaldococcus jannaschii on the operation of the ClpXP protease from Escherichia coli ClpXP completely degrades MJ0366 when pulling from the C-terminal ssrA-tag. However, when a GFP moiety is appended to the N terminus of MJ0366, ClpXP releases intact GFP with a 47-residue tail. The extended length of this tail suggests that ClpXP tightens the trefoil knot against GFP, which prevents GFP unfolding. Interestingly, if the linker between the knot core of MJ0366 and GFP is longer than 36 residues, ClpXP tightens and translocates the knot before it reaches GFP, enabling the complete unfolding and degradation of the substrate. These observations suggest that a knot-induced stall during degradation of multidomain proteins by AAA proteases may constitute a novel mechanism to produce partially degraded products with potentially new functions.

Frequency-control of protein translocation across an oscillating nanopore

The translocation of a lipid binding protein (LBP) is studied using a phenomenological coarse-grained computational model that simplifies both chain and pore geometry.

Polymer translocation through nano-pores in vibrating thin membranes

Polymer translocation is a promising strategy for the next-generation DNA sequencing technologies. The use of biological and synthetic nano-pores, however, still suffers from serious drawbacks. In particular, the width of the membrane layer can accommodate several bases at the same time, making difficult accurate sequencing applications. More recently, the use of graphene membranes has paved the way to new sequencing capabilities, with the possibility to measure transverse currents, among other advances. The reduced thickness of these new membranes poses new questions on the effect of deformability and vibrations of the membrane on the translocation process, two features which are not taken into account in the well established theoretical frameworks. Here, we make a first step forward in this direction. We report numerical simulation work on a model system simple enough to allow gathering significant insight on the effect of these features on the average translocation time, with appropriate statistical significance. We have found that the interplay between thermal fluctuations and the deformability properties of the nano-pore play a crucial role in determining the process. We conclude by discussing new directions for further work.

Translocation dynamics of knotted polymers under a constant or periodic external field

We perform Brownian dynamics simulations to examine how knots alter the dynamics of polymers moving through nanopores under an external field. In the first part of this paper, we study the situation when the field is constant. Here, knots halt translocation above a critical force with jamming occurring at smaller forces for twist topologies compared to non-twist topologies. Slightly below the jamming transition, the polymer's transit times exhibit large fluctuations. This phenomenon is an example of the knot's molecular individualism since the conformation of the knot plays a large role in the chain's subsequent dynamics. In the second part of the paper, we study the motion of the chain when one cycles the field on and off. If the off time is comparable to the knot's relaxation time, one can adjust the swelling of the knot at the pore and hence design strategies to ratchet the polymer in a controllable fashion. We examine how the off time affects the ratcheting dynamics. We also examine how this strategy alters the fluctuations in the polymer's transit time. We find that cycling the force field can reduce fluctuations near the knot's jamming transition, but can enhance the fluctuations at very high forces since knots get trapped in metastable states during the relaxation process. The latter effect appears to be more prominent for non-torus topologies than torus ones. We conclude by discussing the feasibility of this approach to control polymer motion in biotechnology applications such as sequencing.

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

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