Mechanically untying a protein slipknot: multiple pathways revealed by force spectroscopy and steered molecular dynamics simulations

Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein

Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.

Recent advances in sensing the inter-biomolecular interactions at the nanoscale - A comprehensive review of AFM-based force spectroscopy

Biomolecular interactions underpin most processes inside the cell. Hence, a precise and quantitative understanding of molecular association and dissociation events is crucial, not only from a fundamental perspective, but also for the rational design of biomolecular platforms for state-of-the-art biomedical and industrial applications. In this context, atomic force microscopy (AFM) appears as an invaluable experimental technique, allowing the measurement of the mechanical strength of biomolecular complexes to provide a quantitative characterization of their interaction properties from a single molecule perspective. In the present review, the most recent methodological advances in this field are presented with special focus on bioconjugation, immobilization and AFM tip functionalization, dynamic force spectroscopy measurements, molecular recognition imaging and theoretical modeling. We expect this work to significantly aid in grasping the principles of AFM-based force spectroscopy (AFM-FS) technique and provide the necessary tools to acquaint the type of data that can be achieved from this type of experiments. Furthermore, a critical assessment is done with other nanotechnology techniques to better visualize the future prospects of AFM-FS.

Single-Molecule Force Spectroscopy Reveals Stability of mitoNEET and its [2Fe2Se] Cluster in Weakly Acidic and Basic Solutions

The outer mitochondrial membrane protein mitoNEET (mNT) is a recently identified iron-sulfur protein containing a unique Fe₂ S₂ (His)₁ (Cys)₃ metal cluster with a single Fe-N(His87) coordinating bond. This labile Fe-N bond led to multiple unfolding/rupture pathways of mNT and its cluster by atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS), one of most common tools for characterizing the molecular mechanics. Although previous ensemble studies showed that this labile Fe-N(His) bond is essential for protein function, they also indicated that the protein and its [2Fe2S] cluster are stable under acidic conditions. Thus, we applied AFM-SMFS to measure the stability of mNT and its cluster at pH values of 6, 7, and 8. Indeed, all previous multiple unfolding pathways of mNT were still observed. Moreover, single-molecule measurements revealed that the stabilities of the protein and the [2Fe2S] cluster are consistent at these pH values with only ≈20 pN force differences. Thus, we found that the behavior of the protein is consistent in both weakly acidic and basic solutions despite a labile Fe-N bond.

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.

Interpretation of Single-Molecule Force Experiments on Proteins Using Normal Mode Analysis

Single-molecule force spectroscopy experiments allow protein folding and unfolding to be explored using mechanical force. Probably the most informative technique for interpreting the results of these experiments at the structural level makes use of steered molecular dynamics (MD) simulations, which can explicitly model the protein under load. Unfortunately, this technique is computationally expensive for many of the most interesting biological molecules. Here, we find that normal mode analysis (NMA), a significantly cheaper technique from a computational perspective, allows at least some of the insights provided by MD simulation to be gathered. We apply this technique to three non-homologous proteins that were previously studied by force spectroscopy: T4 lysozyme (T4L), Hsp70 and the glucocorticoid receptor domain (GCR). The NMA results for T4L and Hsp70 are compared with steered MD simulations conducted previously, and we find that we can recover the main results. For the GCR, which did not undergo MD simulation, our approach identifies substructures that correlate with experimentally identified unfolding intermediates. Overall, we find that NMA can make a valuable addition to the analysis toolkit for the structural analysis of single-molecule force experiments on proteins.

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.

Reconstruction of ARNT PAS-B Unfolding Pathways by Steered Molecular Dynamics and Artificial Neural Networks

Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains.

Type III secretion system effector proteins are mechanically labile

Multiple gram-negative bacteria encode type III secretion systems (T3SS) that allow them to inject effector proteins directly into host cells to facilitate colonization. To be secreted, effector proteins must be at least partially unfolded to pass through the narrow needle-like channel (diameter <2 nm) of the T3SS. Fusion of effector proteins to tightly packed proteins-such as GFP, ubiquitin, or dihydrofolate reductase (DHFR)-impairs secretion and results in obstruction of the T3SS. Prior observation that unfolding can become rate-limiting for secretion has led to the model that T3SS effector proteins have low thermodynamic stability, facilitating their secretion. Here, we first show that the unfolding free energy ([Formula: see text]) of two Salmonella effector proteins, SptP and SopE2, are 6.9 and 6.0 kcal/mol, respectively, typical for globular proteins and similar to published [Formula: see text] for GFP, ubiquitin, and DHFR. Next, we mechanically unfolded individual SptP and SopE2 molecules by atomic force microscopy (AFM)-based force spectroscopy. SptP and SopE2 unfolded at low force (F _unfold ≤ 17 pN at 100 nm/s), making them among the most mechanically labile proteins studied to date by AFM. Moreover, their mechanical compliance is large, as measured by the distance to the transition state (Δx ^‡ = 1.6 and 1.5 nm for SptP and SopE2, respectively). In contrast, prior measurements of GFP, ubiquitin, and DHFR show them to be mechanically robust (F _unfold > 80 pN) and brittle (Δx ^‡ < 0.4 nm). These results suggest that effector protein unfolding by T3SS is a mechanical process and that mechanical lability facilitates efficient effector protein secretion.

Ultrahigh Adhesion Force Between Silica-Binding Peptide SB7 and Glass Substrate Studied by Single-Molecule Force Spectroscopy and Molecular Dynamic Simulation

Many proteins and peptides have been identified to effectively and specifically bind on certain surfaces such as silica, polystyrene and titanium dioxide. It is of great interest, in many areas such as enzyme immobilization, surface functionalization and nanotechnology, to understand how these proteins/peptides bind to solid surfaces. Here we use single-molecule force spectroscopy (SMFS) based on atomic force microscopy to directly measure the adhesion force between a silica-binding peptide SB7 and glass surface at single molecule level. SMFS results show that the adhesion force of a single SB7 detaching from the glass surface distributes in two populations at ~220 pN and 610 pN, which is higher than the unfolding forces of most mechanically stable proteins and the unbinding forces of most stable protein-protein interactions. Molecular dynamics simulation reveals that the electrostatic interactions between positively charged arginine residues and the silica surface dominates the binding of SB7 on silica. Our study provides experimental evidence and molecular mechanism at the single-molecule level for the SB7-based immobilization of proteins on silica-based surface, which is able to withstand high mechanical forces, making it an ideal fusion tag for silica surface immobilization or peptide-base adhesive materials.

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.

Enzymatic Construction of Protein Polymer/Polyprotein Using OaAEP1 and TEV Protease

The development of chemical and biological coupling technologies in recent years has made possible of protein polymers engineering. We have developed an enzymatic method for building polyproteins using a protein ligase OaAEP1 (asparagine endopeptidase 1) and protease TEV (tobacco etching virus). Using a mobile TEV protease site compatible with the OaAEP1 ligation, we achieved a stepwise polymerization of the protein on the surface. The produced polyprotein can be verified by protein unfolding scenario using atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS). Thus, this study provides an alternative method for polyprotein engineering and immobilization.

The conformational phase space of the photoprotective switch in the major light harvesting complex II

Extensive conformational sampling of the major light harvesting complex II defines the collective variables of the photoprotective switch.

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.

Fine tuning of the photosystem II major antenna mobility within the thylakoid membrane of higher plants

Depending on the amount of light, the photosystem II (PSII) antennae or Light Harvesting Complexes (LHCII) switch between two states within the thylakoid membranes of higher plants, i.e., a light-harvesting and a photoprotective mode. This switch is co-regulated by a pH gradient (ΔpH) across the membrane and the interaction with the PSII subunit S (PsbS) that is proposed to induce LHCII aggregation. Herein we employ all-atom and coarse-grained molecular simulations of the major LHCII trimer at low and excess ΔpH, as well as in complexation with PsbS within a native thylakoid membrane model. Our results demonstrate the aggregation potential of LHCII and, consistent with the experimental literature, reveal the role of PsbS at atomic resolution. PsbS alters the LHCII-thylakoid lipid interactions and restores the LHCII mobility that is lost in the transition to photoprotective conditions (low lumenal pH). In agreement with this finding, diffusion of the integral membrane protein LHCII is dependent on both, electrostatic interactions and hydrophobic mismatch, while it does not obey the Saffman-Delbrück diffusion model.

Multistep Protein Unfolding Scenarios from the Rupture of a Complex Metal Cluster Cd3S9

Protein (un)folding is a complex and essential process. With the rapid development of single-molecule techniques, we can detect multiple and transient proteins (un)folding pathways/intermediates. However, the observation of multiple multistep (>2) unfolding scenarios for a single protein domain remains limited. Here, we chose metalloprotein with relatively stable and multiple metal-ligand coordination bonds as a system for such a purpose. Using AFM-based single-molecule force spectroscopy (SMFS), we successfully demonstrated the complex and multistep protein unfolding scenarios of the β-domain of a human protein metallothionein-3 (MT). MT is a protein of ~60 amino acids (aa) in length with 20 cysteines for various metal binding, and the β-domain (βMT) is of ~30 aa with an M3S9 metal cluster. We detected four different types of three-step protein unfolding scenarios from the Cd-βMT, which can be possibly explained by the rupture of Cd-S bonds in the complex Cd3S9 metal cluster. In addition, complex unfolding scenarios with four rupture peaks were observed. The Cd-S bonds ruptured in both single bond and multiple bonds modes. Our results provide not only evidence for multistep protein unfolding phenomena but also reveal unique properties of metalloprotein system using single-molecule AFM.

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.

The Photosystem II Subunit S under Stress

Nonphotochemical quenching is the protective mechanism against overexcitation of photosystem II, triggered by excess ΔpH in photosynthetic membranes. The light-harvesting complexes (LHCs), the de-epoxidation of violaxanthin to zeaxanthin, and the photosystem II subunit S (PsbS) work in synergy for an optimized multilevel response. Understanding the fine details of this synergy has proven challenging to scientific research. Here, we employ large-scale, all-atom molecular simulations and beyond experimental insight, we proceed a step further in identifying the PsbS dynamics that could possibly be associated with this synergy. For the first time, to our knowledge, we probe the distinct behavior of PsbS under ΔpH that probes the details of the potential dimer-to-monomer transition, and in a violaxanthin/zeaxanthin-rich membrane, at an all-atom resolution. We propose that the lumen-exposed residues, threonine 162 and glutamic acid 173, form stabilizing hydrogen bonds between the PsbS monomers only at high lumen pH, whereas at low pH (excess ΔpH) this interaction is lost, and leads to higher flexibility of the protein and potentially to the dimer-to-monomer transition. Lastly, we discuss how conformational changes under the presence of ΔpH/zeaxanthin are related to the PsbS role in the current nonphotochemical quenching model in the literature. For the latter, we probe a PsbS-monomeric LHCII association. The association is proposed to potentially alter the monomeric LHCII sensitivity to ΔpH by changing the pKa values of interacting LHCII residues. This serves as an example where protonation-ligation events enhance protein-protein interactions fundamental to many life processes.

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.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.

Effect of Chain Conformation on the Single-Molecule Melting Force in Polymer Single Crystals: Steered Molecular Dynamics Simulations Study

Understanding the relationship between polymer chain conformation as well as the chain composition within the single crystal and the mechanical properties of the corresponding single polymer chain will facilitate the rational design of high performance polymer materials. Here three model systems of polymer single crystals, namely poly(ethylene oxide) (PEO), polyethylene (PE), and nylon-66 (PA66) have been chosen to study the effects of chain conformation, helical (PEO) versus planar zigzag conformation (PE, PA66), and chain composition (PE versus PA66) on the mechanical properties of a single polymer chain. To do that, steered molecular dynamics simulations were performed on those polymer single crystals by pulling individual polymer chains out of the crystals. Our results show that the patterns of force-extension curve as well as the chain moving mode are closely related to the conformation of the polymer chain in the single crystal. In addition, hydrogen bonds can enhance greatly the force required to stretch the polymer chain out of the single crystal. The dynamic breaking and reformation of multivalent hydrogen bonds have been observed for the first time in PA66 at the single molecule level.

Thermal Unthreading of the Lasso Peptides Astexin-2 and Astexin-3

Lasso peptides are a class of knot-like polypeptides in which the C-terminal tail of the peptide threads through a ring formed by an isopeptide bond between the N-terminal amine group and a side chain carboxylic acid. The small size (∼20 amino acids) and simple topology of lasso peptides make them a good model system for studying the unthreading of entangled polypeptides, both with experiments and atomistic simulation. Here, we present an in-depth study of the thermal unthreading behavior of two lasso peptides astexin-2 and astexin-3. Quantitative kinetics and energetics of the unthreading process were determined for variants of these peptides using a series of chromatography and mass spectrometry experiments and biased molecular dynamics (MD) simulations. In addition, we show that the Tyr15Phe variant of astexin-3 unthreads via an unprecedented "tail pulling" mechanism. MD simulations on a model ring-thread system coupled with machine learning approaches also led to the discovery of physicochemical descriptors most important for peptide unthreading.

A simple DNA handle attachment method for single molecule mechanical manipulation experiments

Manipulating single molecules and systems of molecules with mechanical force is a powerful technique to examine their physical properties. Applying force requires attachment of the target molecule to larger objects using some sort of molecular tether, such as a strand of DNA. DNA handle attachment often requires difficult manipulations of the target molecule, which can preclude attachment to unstable, hard to obtain, and/or large, complex targets. Here we describe a method for covalent DNA handle attachment to proteins that simply requires the addition of a preprepared reagent to the protein and a short incubation. The handle attachment method developed here provides a facile approach for studying the biomechanics of biological systems.

Rapid and Robust Polyprotein Production Facilitates Single-Molecule Mechanical Characterization of β-Barrel Assembly Machinery Polypeptide Transport Associated Domains

Single-molecule force spectroscopy by atomic force microscopy exploits the use of multimeric protein constructs, namely, polyproteins, to decrease the impact of nonspecific interactions, to improve data accumulation, and to allow the accommodation of benchmarking reference domains within the construct. However, methods to generate such constructs are either time- and labor-intensive or lack control over the length or the domain sequence of the obtained construct. Here, we describe an approach that addresses both of these shortcomings that uses Gibson assembly (GA) to generate a defined recombinant polyprotein rapidly using linker sequences. To demonstrate the feasibility of this approach, we used GA to make a polyprotein composed of alternating domains of I27 and TmCsp, (I27-TmCsp)3-I27)(GA), and showed the mechanical fingerprint, mechanical strength, and pulling speed dependence are the same as an analogous polyprotein constructed using the classical approach. After this benchmarking, we exploited this approach to facilitiate the mechanical characterization of POTRA domain 2 of BamA from E. coli (EcPOTRA2) by assembling the polyprotein (I27-EcPOTRA2)3-I27(GA). We show that, as predicted from the α + β topology, EcPOTRA2 domains are mechanically robust over a wide range of pulling speeds. Furthermore, we identify a clear correlation between mechanical robustness and brittleness for a range of other α + β proteins that contain the structural feature of proximal terminal β-strands in parallel geometry. We thus demonstrate that the GA approach is a powerful tool, as it circumvents the usual time- and labor-intensive polyprotein production process and allows for rapid production of new constructs for single-molecule studies. As shown for EcPOTRA2, this approach allows the exploration of the mechanical properties of a greater number of proteins and their variants. This improves our understanding of the relationship between structure and mechanical strength, increasing our ability to design proteins with tailored mechanical properties.

Molecular knots in biology and chemistry

Knots and entanglements are ubiquitous. Beyond their aesthetic appeal, these fascinating topological entities can be either useful or cumbersome. In recent decades, the importance and prevalence of molecular knots have been increasingly recognised by scientists from different disciplines. In this review, we provide an overview on the various molecular knots found in naturally occurring biological systems (DNA, RNA and proteins), and those created by synthetic chemists. We discuss the current knowledge in these fields, including recent developments in experimental and, in some cases, computational studies which are beginning to shed light into the complex interplay between the structure, formation and properties of these topologically intricate molecules.

Engineering covalent loops in proteins can serve as an on/off switch to regulate threaded topologies

Knots in proteins are under active investigation motivating refinements of current techniques and the development of tools to better understand the knotted topology. A strong focus is to identify new knots and expand upon the current understanding of their complex topology. Previous work has shown that the knotted topology, even in the simplest case of knots, encompasses a variety of unique challenges in folding and tying a chain. To bypass many of the in vitro experimental complications involved in working with knots, it is useful to apply methodologies to a more simplified system. The pierced lasso bundles (PLB), we discovered where a single disulphide bridge holds the threaded topology together, presents a simpler system to study knots in vitro. Having a disulphide bridge as an on/off switch between the threaded/unthreaded topology is advantageous because a covalent loop allows manipulation of the knot without directly altering affecting secondary and tertiary structure. Because disulphide bridges are commonly used in protein engineering, a pierced lasso (PL) topology can be easily introduced into a protein of interest to form a knotted topology within a given secondary structure. It is also important to take into account that if formed, disulphides can inadvertently introduce an unwanted PL. This was found upon determination of the crystal structure (PDB code 2YHG) of the recently de novo designed nucleoside hydrolase. Our detailed investigations of the PL presented here will allow researchers to look at the introduction of disulphide bridges in a larger context with respect to potential geometrical consequences on the structure and functional properties of proteins.

The mechanochemistry of copper reports on the directionality of unfolding in model cupredoxin proteins

Understanding the directionality and sequence of protein unfolding is crucial to elucidate the underlying folding free energy landscape. An extra layer of complexity is added in metalloproteins, where a metal cofactor participates in the correct, functional fold of the protein. However, the precise mechanisms by which organometallic interactions are dynamically broken and reformed on (un)folding are largely unknown. Here we use single molecule force spectroscopy AFM combined with protein engineering and MD simulations to study the individual unfolding pathways of the blue-copper proteins azurin and plastocyanin. Using the nanomechanical properties of the native copper centre as a structurally embedded molecular reporter, we demonstrate that both proteins unfold via two independent, competing pathways. Our results provide experimental evidence of a novel kinetic partitioning scenario whereby the protein can stochastically unfold through two distinct main transition states placed at the N and C termini that dictate the direction in which unfolding occurs.

Polyproteins in structural biology

•

Structures have been determined for natural and recombinant polyproteins.

•

Native HIV Gag polyprotein architecture was revealed by cryo-EM of immature capsids.

•

Recombinant polyprotein technology has resolved sample preparation bottlenecks.

•

The high-resolution structure of influenza polymerase has been solved.

•

Single-molecule analysis of polyproteins revealed their folding characteristics.

Polyproteins are chains of covalently conjoined smaller proteins that occur in nature as versatile means to organize the proteome of viruses including HIV. During maturation, viral polyproteins are typically cleaved into the constituent proteins with different biological functions by highly specific proteases, and structural analyses at defined stages of this maturation process can provide clues for antiviral intervention strategies. Recombinant polyproteins that use similar mechanisms are emerging as powerful tools for producing hitherto inaccessible protein targets such as the influenza polymerase, for high-resolution structure determination by X-ray crystallography. Conversely, covalent linking of individual protein subunits into single polypeptide chains are exploited to overcome sample preparation bottlenecks. Moreover, synthetic polyproteins provide a promising tool to dissect dynamic folding of polypeptide chains into three-dimensional architectures in single-molecule structure analysis by atomic force microscopy (AFM). The recent use of natural and synthetic polyproteins in structural biology and major achievements are highlighted in this contribution.

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.

Mechanical unfolding of ribose binding protein and its comparison with other periplasmic binding proteins

Folding and unfolding studies on large, multidomain proteins are still rare despite their high abundance in genomes of prokaryotes and eukaryotes. Here, we investigate the unfolding properties of a 271 residue, two-domain ribose binding protein (RBP) from the bacterial periplasm using single-molecule force spectroscopy. We observe that RBP predominately unfolds via a two-state pathway with an unfolding force of ∼80 pN and an unfolding contour length of ∼95 nm. Only a small population (∼15%) of RBP follows three-state pathways. The ligand binding neither increases the mechanical stability nor influences the unfolding flux of RBP through different pathways. The kinetic partitioning between two-state and three-state pathways, which has been reported earlier for other periplasmic proteins, is also observed in RBP, albeit to a lesser extent. These results provide important insights into the mechanical stability and unfolding processes of large two-domain proteins and highlight the contrasting features upon ligand binding. Protein structural topology diagrams are used to explain the differences in the mechanical unfolding behavior of RBP with other periplasmic binding proteins.

Knotting pathways in proteins

Most proteins, in order to perform their biological function, have to fold to a compact native state. The increasing number of knotted and slipknotted proteins identified suggests that proteins are able to manoeuvre around topological barriers during folding. In the present article, we review the current progress in elucidating the knotting process in proteins. Although we concentrate on theoretical approaches, where a knotted topology can be unambiguously detected, comparison with experiments is also reviewed. Numerical simulations suggest that the folding process for small knotted proteins is composed of twisted loop formation and then threading by either slipknot geometries or flipping. As the size of the knotted proteins increases, particularly for more deeply threaded termini, the prevalence of traps in the free energy landscape also increases. Thus, in the case of longer knotted and slipknotted proteins, the folding mechanism is probably supported by chaperones. Overall, results imply that knotted proteins can be folded efficiently and survive evolutionary pressure in order to perform their biological functions.

Tightening and Untying the Knot in Human Carbonic Anhydrase III

The forced mechanical unfolding of the knotted protein Human Carbonic Anhydrase (HCA) III is examined by steered, explicit-water molecular dynamics computer simulations. In agreement with previous indications from experiments and coarse-grained simulations, knot tightening by pulling near-terminal amino acids (4 and 267) leads to a much higher resistance to unfolding than for knot untying, where pulling amino acids 4 and 253 untangles the knot by threading the C-terminal end out of the knotting loop. In particular, the resistance during knot tightening is observed to diverge due to a tightly tied-up enzymatic core of the HCA if it is coordinated by the catalytically important zinc ion. The underlying structural pictures are presented and discussed.

Steered molecular dynamics studies of the potential of mean force for peptide amphiphile self-assembly into cylindrical nanofibers

Steered molecular dynamics (SMD) simulations were applied to determine the potential of mean force for the self-assembly of peptide amphiphile (PA) nanofibers, specifically considering a single PA adding to a growing cylindrical nanofiber at 310 K. It is found that the free energy, enthalpy, and entropy differences for this assembly process are -67 kcal/mol, -71.5 kcal/ml, and -14.5 cal/(mol K), respectively, and therefore that enthalpy provides the driving force for self-assembly to form a fiber. A pairwise interaction analysis shows that both electrostatic and van der Waals interactions play important roles in the self-assembly process, with the van der Waals interaction being the larger effect. The mechanistic picture that emerges from this work is that as the PA is pulled from the fiber, the interaction evolves through three stages: (1) initially electrostatic interactions between the charged head of the pulled PA and other PAs, and between the pulled PA and solvent are dominant, (2) after the charged head emerges, the rest of the peptide comes out, with both PA-solvent electrostatic interactions and van der Waals interactions being significant, and (3) in the last step, the alkane tail emerges, dominated by van der Waals interactions with either peptide or solvent.

Formation of cystine slipknots in dimeric proteins

We consider mechanical stability of dimeric and monomeric proteins with the cystine knot motif. A structure based dynamical model is used to demonstrate that all dimeric and some monomeric proteins of this kind should have considerable resistance to stretching that is significantly larger than that of titin. The mechanisms of the large mechanostability are elucidated. In most cases, it originates from the induced formation of one or two cystine slipknots. Since there are four termini in a dimer, there are several ways of selecting two of them to pull by. We show that in the cystine knot systems, there is strong anisotropy in mechanostability and force patterns related to the selection. We show that the thermodynamic stability of the dimers is enhanced compared to the constituting monomers whereas machanostability is either lower or higher.

Force as a single molecule probe of multidimensional protein energy landscapes

Force spectroscopy has developed into an indispensable tool for studying folding and binding of proteins on a single molecule level in real time. Design of the pulling geometry allows tuning the reaction coordinate in a very precise manner. Many recent experiments have taken advantage of this possibility and have provided detailed insight the folding pathways on the complex high dimensional energy landscape. Beyond its potential to provide control over the reaction coordinate, force is also an important physiological parameter that affects protein conformation under in vivo conditions. Single molecule force spectroscopy studies have started to unravel the response and adaptation of force bearing protein structures to mechanical loads.

Using simulations to provide the framework for experimental protein folding studies

Molecular dynamics simulations are a powerful theoretical tool to model the protein folding process in atomistic details under realistic conditions. Combined with a number of experimental techniques, simulations provide a detailed picture of how a protein folds or unfolds in the presence of explicit solvent and other molecular species, such as cosolvents, osmolytes, cofactors, active binding partners or inert crowding agents. The denaturing effects of temperature, pressure and external mechanical forces can also be probed. Qualitative and quantitative agreement with experiment contributes to a comprehensive molecular picture of protein states along the folding/unfolding pathway. The variety of systems examined reveals key features of the protein folding process.