Direct observation of a force-induced switch in the anisotropic mechanical unfolding pathway of a protein

Anomalous Force-Dependent Transition Rates Unveil Dual Pathways in Folding and Unfolding Dynamics of Acyl-coenzyme A Binding Protein

All-α proteins typically fold rapidly and are unable to withstand high forces. Acyl-coenzyme A binding protein (ACBP), a four-α-helix bundle protein, serves as a model protein for studying the folding dynamics of all-α proteins. In previous biochemistry and single molecule force spectroscopy experiments, a controversy exists for the folding pathway and the conformation of the transition state. In this article, we investigate the folding and unfolding dynamics of ACBP in a force range of 4-10 pN using magnetic tweezers, revealing anomalous force-dependent transition rates. The unfolding rate of ACBP remains nearly constant when force is below 6 pN, and it increases sharply when the force exceeds this threshold, while the logarithm of its folding rate is almost a linear function of force. Detailed analysis combined with molecular dynamics simulations indicates that ACBP has two transition pathways: one dominating at zero or low force and the other dominating at high force. Our results provide strong evidence that stretching force not only modulates the folding and unfolding rates but also switches the transition pathways, leading to complex force response behaviors.

Capturing the impact of protein unfolding on the dynamic assembly of protein networks

The rapid assembly of molecular or nanoscale building blocks into extended arrays is crucial to the construction of functional networks in vivo and in vitro and depends on various factors. One factor seldom considered is the dynamic changes of the building block shape. Folded protein building blocks offer a unique system to investigate dynamic shape changes due to their intrinsic ability to change from a compact and specific folded structure to an extended unfolded structure in response to a perturbation such as force. Here, we use photochemically crosslinked folded protein hydrogels constructed from force labile protein building blocks as a model dynamic shape-changing network system and characterise them by combining time-resolved rheology and small-angle X-ray scattering (SAXS). This approach probes both the load-bearing network structures, using rheology, and network architectures, using SAXS, thereby providing a crosslength scale understanding of the network formation. We propose a triple assembly model for the structural evolution of networks constructed from force labile protein building block consisting of: primary formation where monomeric folded proteins create the preliminary protein network scaffold; a subsequent secondary formation phase, where larger oligomers of protein diffuse to join the preliminary network scaffold; and finally in situ unfolding and relaxation which leads to the mature network structure of connected larger and denser fractal-like clusters. The time-resolved SAXS data provides evidence that protein unfolding occurs on the edges of the fractal-like clusters, resulting in a population of unfolded proteins in the space between clusters. Identifying the key stages of assembly in protein networks constructed from force labile proteins provides a greater understanding of the importance of protein unfolding in hierarchical biomechanics in vivo and creates future opportunities to develop bespoke biomaterials for novel biomedical applications.

Pulling Forces Differentially Affect Refolding Pathways Due to Entangled Misfolded States in SARS-CoV-1 and SARS-CoV-2 Receptor Binding Domain

Single-molecule force spectroscopy (SMFS) experiments can monitor protein refolding by applying a small force of a few piconewtons (pN) and slowing down the folding process. Bell theory predicts that in the narrow force regime where refolding can occur, the folding time should increase exponentially with increased external force. In this work, using coarse-grained molecular dynamics simulations, we compared the refolding pathways of SARS-CoV-1 RBD and SARS-CoV-2 RBD (RBD refers to the receptor binding domain) starting from unfolded conformations with and without a force applied to the protein termini. For SARS-CoV-2 RBD, the number of trajectories that fold is significantly reduced with the application of a 5 pN force, indicating that, qualitatively consistent with Bell theory, refolding is slowed down when a pulling force is applied to the termini. In contrast, the refolding times of SARS-CoV-1 RBD do not change meaningfully when a force of 5 pN is applied. How this lack of a Bell response could arise at the molecular level is unknown. Analysis of the entanglement changes of the folded conformations revealed that in the case of SARS-CoV-1 RBD, an external force minimizes misfolding into kinetically trapped states, thereby promoting efficient folding and offsetting any potential slowdown due to the external force. These misfolded states contain non-native entanglements that do not exist in the native state of either SARS-CoV-1-RBD or SARS-CoV-2-RBD. These results indicate that non-Bell behavior can arise from this class of misfolding and, hence, may be a means of experimentally detecting these elusive, theoretically predicted states.

Theory and Examples of Catch Bonds

We discuss slip bonds, catch bonds, and the tug-of-war mechanism using mathematical arguments. The aim is to explain the theoretical tool of molecular potential energy surfaces (PESs). For this, we propose simple 2-dimensional surface models to demonstrate how a molecule under an external force behaves. Examples are selectins. Catch bonds, in particular, are explained in more detail, and they are contrasted to slip bonds. We can support special two-dimensional molecular PESs for E- and L-selectin which allow the catch bond property. We demonstrate that Newton trajectories (NT) are powerful tools to describe these phenomena. NTs form the theoretical background of mechanochemistry.

Exploring the free energy landscape of proteins using magnetic tweezers

Proteins fold to their native states by searching through the free energy landscapes. As single-domain proteins are the basic building block of multiple-domain proteins or protein complexes composed of subunits, the free energy landscapes of single-domain proteins are of critical importance to understand the folding and unfolding processes of proteins. To explore the free energy landscapes of proteins over large conformational space, the stability of native structure is perturbed by biochemical or mechanical means, and the conformational transition process is measured. In single molecular manipulation experiments, stretching force is applied to proteins, and the folding and unfolding transitions are recorded by the extension time course. Due to the broad force range and long-time stability of magnetic tweezers, the free energy landscape over large conformational space can be obtained. In this article, we describe the magnetic tweezers instrument design, protein construct design and preparation, fluid chamber preparation, common-used measuring protocols including force-ramp and force-jump measurements, and data analysis methods to construct the free energy landscape. Single-domain cold shock protein is introduced as an example to build its free energy landscape by magnetic tweezers measurements.

Mechanical forces control the valency of the malaria adhesin VAR2CSA by exposing cryptic glycan binding sites

Plasmodium falciparum (Pf) is responsible for the most lethal form of malaria. VAR2CSA is an adhesin protein expressed by this parasite at the membrane of infected erythrocytes for attachment to the placenta, leading to pregnancy-associated malaria. VAR2CSA is a large 355 kDa multidomain protein composed of nine extracellular domains, a transmembrane helix, and an intracellular domain. VAR2CSA binds to Chondroitin Sulphate A (CSA) of the proteoglycan matrix of the placenta. Shear flow, as the one occurring in blood, has been shown to enhance the (VAR2CSA-mediated) adhesion of Pf-infected erythrocytes on the CSA-matrix. However, the underlying molecular mechanism governing this enhancement has remained elusive. Here, we address this question by using equilibrium, force-probe, and docking-based molecular dynamics simulations. We subjected the VAR2CSA protein-CSA sugar complex to a force mimicking the tensile force exerted on this system due to the shear of the flowing blood. We show that upon this force exertion, VAR2CSA undergoes a large opening conformational transition before the CSA sugar chain dissociates from its main binding site. This preferential order of events is caused by the orientation of the molecule during elongation, as well as the strong electrostatic attraction of the sugar to the main protein binding site. Upon opening, two additional cryptic CSA binding sites get exposed and a functional dodecameric CSA molecule can be stably accommodated at these force-exposed positions. Thus, our results suggest that mechanical forces increase the avidity of VAR2CSA by turning it from a monovalent to a multivalent state. We propose this to be the molecular cause of the observed shear-enhanced adherence. Mechanical control of the valency of VAR2CSA is an intriguing hypothesis that can be tested experimentally and which is of relevance for the understanding of the malaria infection and for the development of anti placental-malaria vaccines targeting VAR2CSA.

Mechanical Unfolding and Refolding of NanoLuc via Single-Molecule Force Spectroscopy and Computer Simulations

A highly bioluminescent protein, NanoLuc (Nluc), has seen numerous applications in biological assays since its creation. We recently engineered a NanoLuc polyprotein that showed high bioluminescence but displayed a strong misfolding propensity after mechanical unfolding. Here, we present our single-molecule force spectroscopy (SMFS) studies by atomic force microscopy (AFM) and steered molecular dynamics (SMD) simulations on two new hybrid protein constructs comprised of Nluc and I91 titin domains, I91-I91-Nluc-I91-I91-I91-I91 (I91₂-Nluc-I91₄) and I91-Nluc-I91-Nluc-I91-Nluc-I91, to characterize the unfolding behavior of Nluc in detail and to further investigate its misfolding properties that we observed earlier for the I91₂-Nluc₃-I91₂ construct. Our SMFS results confirm that Nluc's unfolding proceeds similarly in all constructs; however, Nluc's refolding differs in these constructs, and its misfolding is minimized when Nluc is monomeric or separated by I91 domains. Our simulations on monomeric Nluc, Nluc dyads, and Nluc triads pinpointed the origin of its mechanical stability and captured interesting unfolding intermediates, which we also observed experimentally.

Histidine-Specific Bioconjugation for Single-Molecule Force Spectroscopy

Atomic force microscopy (AFM) based single-molecule force spectroscopy (SMFS) is a powerful tool to study the mechanical properties of proteins. In these experiments, site-specific immobilization of proteins is critical, as the tether determines the direction and amplitude of forces applied to the protein of interest. However, existing methods are mainly based on thiol chemistry or specific protein tags, which cannot meet the need of many challenging experiments. Here, we developed a histidine-specific phosphorylation strategy to covalently anchor proteins to an AFM cantilever tip or the substrate via their histidine tag or surface-exposed histidine residues. The formed covalent linkage was mechanically stable with rupture forces of over 1.3 nN. This protein immobilization method considerably improved the pickup rate and data quality of SMFS experiments. We further demonstrated the use of this method to explore the pulling-direction-dependent mechanical stability of green fluorescent protein and the unfolding of the membrane protein archaerhodopsin-3.

Probing and Manipulating Noncovalent Interactions in Functional Polymeric Systems

Noncovalent interactions, which usually feature tunable strength, reversibility, and environmental adaptability, have been recognized as driving forces in a variety of biological and chemical processes, contributing to the recognition between molecules, the formation of molecule clusters, and the establishment of complex structures of macromolecules. The marriage of noncovalent interactions and conventional covalent polymers offers the systems novel mechanical, physicochemical, and biological properties, which are highly dependent on the binding mechanisms of the noncovalent interactions that can be illuminated via quantification. This review systematically discusses the nanomechanical characterization of typical noncovalent interactions in polymeric systems, mainly through direct force measurements at microscopic, nanoscopic, and molecular levels, which provide quantitative information (e.g., ranges, strengths, and dynamics) on the binding behaviors. The fundamental understandings of intermolecular and interfacial interactions are then correlated to the macroscopic performances of a series of noncovalently bonded polymers, whose functions (e.g., stimuli-responsiveness, self-healing capacity, universal adhesiveness) can be customized through the manipulation of the noncovalent interactions, providing insights into the rational design of advanced materials with applications in biomedical, energy, environmental, and other engineering fields.

Free-energy landscape of two-state protein acylphosphatase with large contact order revealed by force-dependent folding and unfolding dynamics

Acylphosphatase (AcP) is a small protein with 98 amino acid residues that catalyzes the hydrolysis of carboxyl-phosphate bonds. AcP is a typical two-state protein with slow folding rate due to its relatively large contact order in the native structure. The mechanical properties and unfolding behavior of AcP has been studied by atomic force microscope. Here using stable magnetic tweezers, we measured the force-dependent folding rates within a force range 1-3 pN, and unfolding rates 15-40 pN. The obtained unfolding rates show different force sensitivities at forces below and above ∼27 pN, which determines a free-energy landscape with two energy barriers. Our results indicate that the free-energy landscape of small globule proteins have general Bactrian camel shape, and large contact order of the native state produces a high barrier dominate at low forces.

Exploring the Effect of Mechanical Anisotropy of Protein Structures in the Unfoldase Mechanism of AAA+ Molecular Machines

Essential cellular processes of microtubule disassembly and protein degradation, which span lengths from tens of μm to nm, are mediated by specialized molecular machines with similar hexameric structure and function. Our molecular simulations at atomistic and coarse-grained scales show that both the microtubule-severing protein spastin and the caseinolytic protease ClpY, accomplish spectacular unfolding of their diverse substrates, a microtubule lattice and dihydrofolate reductase (DHFR), by taking advantage of mechanical anisotropy in these proteins. Unfolding of wild-type DHFR requires disruption of mechanically strong β-sheet interfaces near each terminal, which yields branched pathways associated with unzipping along soft directions and shearing along strong directions. By contrast, unfolding of circular permutant DHFR variants involves single pathways due to softer mechanical interfaces near terminals, but translocation hindrance can arise from mechanical resistance of partially unfolded intermediates stabilized by β-sheets. For spastin, optimal severing action initiated by pulling on a tubulin subunit is achieved through specific orientation of the machine versus the substrate (microtubule lattice). Moreover, changes in the strength of the interactions between spastin and a microtubule filament, which can be driven by the tubulin code, lead to drastically different outcomes for the integrity of the hexameric structure of the machine.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

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

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

Single-Molecule Force Spectroscopy of Protein Folding

The use of force probes to induce unfolding and refolding of single molecules through the application of mechanical tension, known as single-molecule force spectroscopy (SMFS), has proven to be a powerful tool for studying the dynamics of protein folding. Here we provide an overview of what has been learned about protein folding using SMFS, from small, single-domain proteins to large, multi-domain proteins. We highlight the ability of SMFS to measure the energy landscapes underlying folding, to map complex pathways for native and non-native folding, to probe the mechanisms of chaperones that assist with native folding, to elucidate the effects of the ribosome on co-translational folding, and to monitor the folding of membrane proteins.

Anisotropy in mechanical unfolding of protein upon partner-assisted pulling and handle-assisted pulling

Proteins as force-sensors respond to mechanical cues and regulate signaling in physiology. Proteins commonly connect the source and response points of mechanical cues in two conformations, independent proteins in end-to-end geometry and protein complexes in handshake geometry. The force-responsive property of independent proteins in end-to-end geometry is studied extensively using single-molecule force spectroscopy (SMFS). The physiological significance of the complex conformations in force-sensing is often disregarded as mere surge protectors. However, with the potential of force-steering, protein complexes possess a distinct mechano-responsive property over individual force-sensors. To decipher, we choose a force-sensing protein, cadherin-23, from tip-link complex and perform SMFS using end-to-end geometry and handshake complex geometry. We measure higher force-resilience of cadherin-23 with preferential shorter extensions in handshake mode of pulling over the direct mode. The handshake geometry drives the force-response of cadherin-23 through different potential-energy landscapes than direct pulling. Analysis of the dynamic network structure of cadherin-23 under tension indicates narrow force-distributions among residues in cadherin-23 in direct pulling, resulting in low force-dissipation paths and low resilience to force. Overall, the distinct and superior mechanical responses of cadherin-23 in handshake geometry than single protein geometry highlight a probable evolutionary drive of protein-protein complexes as force-conveyors over independent ones.

Energy Landscape of Ubiquitin Is Weakly Multidimensional

Single molecule pulling experiments report time-dependent changes in the extension (X) of a biomolecule as a function of the applied force (f). By fitting the data to one-dimensional analytical models of the energy landscape, we can extract the hopping rates between the folded and unfolded states in two-state folders as well as the height and the location of the transition state (TS). Although this approach is remarkably insightful, there are cases for which the energy landscape is multidimensional (catch bonds being the most prominent). To assess if the unfolding energy landscape in small single domain proteins could be one-dimensional, we simulated force-induced unfolding of ubiquitin (Ub) using the coarse-grained self-organized polymer-side chain (SOP-SC) model. Brownian dynamics simulations using the SOP-SC model reveal that the Ub energy landscape is weakly multidimensional (WMD), governed predominantly by a single barrier. The unfolding pathway is confined to a narrow reaction pathway that could be described as diffusion in a quasi-1D X-dependent free energy profile. However, a granular analysis using the P_fold analysis, which does not assume any form for the reaction coordinate, shows that X alone does not account for the height and, more importantly, the location of the TS. The f-dependent TS location moves toward the folded state as f increases, in accord with the Hammond postulate. Our study shows that, in addition to analyzing the f-dependent hopping rates, the transition state ensemble must also be determined without resorting to X as a reaction coordinate to describe the unfolding energy landscapes of single domain proteins, especially if they are only WMD.

SlyD Accelerates trans-to-cis Prolyl Isomerization in a Mechanosignaling Protein under Load

Prolyl isomerization is recognized as one of the key regulatory mechanisms, which plays a crucial role in cell signaling, ion channel gating, phage virus infection, and molecular timing. This isomerization is usually slow but often accelerated by an enzyme, called peptidyl-prolyl isomerase (PPIase). In the current project, we investigate using single-molecule force spectroscopy (SMFS) the impact of a bacterial PPIase, SlyD, on the cis-trans isomerization of the proline 2225 (P2225) in an isolated 20th domain of a cytoskeletal mechanosensing protein filamin-A (FlnA20). To explore the FlnA20-PPIase interaction, we have used multiple SMFS modes, like constant velocity, constant distance, and jumping trap experiments. In our previous study, we reported the unique nature of the P2225, which is conserved in all naturally occurring filamins and can slowly (minutes) interconvert between cis-trans isomers, in absence of any PPIase. Our current results show a staggering 25-fold acceleration of the trans-to-cis isomerization rate in the presence of saturating SlyD concentration (7.25 μM) compared to the unenzymatic condition. A SlyD concentration-dependent depletion of the trans isomeric lifetime was also observed. Additionally, we observed that SlyD stabilizes the cis-isomer in the native state of FlnA20 by ∼2 k_BT. This is the first single-molecule observation of the cis-trans isomerization catalysis by a PPIase in a mechanosensing protein.

SpyTag/SpyCatcher tether as a fingerprint and force marker in single-molecule force spectroscopy experiments

Single molecule force spectroscopy has emerged as a powerful tool to study protein folding dynamics, ligand-receptor interactions, and various mechanobiological processes. High force precision does not necessarily lead to high force accuracy, as the uncertainties in calibration can bring serious systematic errors. In the case of magnetic tweezers, accurate determination of the applied forces for short biomolecular tethers, by measuring thermal fluctuations of inhomogeneous magnetic beads, remains difficult. Here we address this challenge by showing that the SpyTag/SpyCatcher complex is not only a convenient and genetically encodable covalent linker but also an ideal molecular fingerprint and force marker in single molecule force spectroscopy experiments. By stretching the N-termini of both SpyCatcher and SpyTag, the complex unfolds locally up to the isopeptide bond position in an unzipping geometry, resulting in equilibrium transitions at ∼30 pN with step sizes of ∼3.4 nm. This mechanical feature can be used as the fingerprint to identify single-molecular events. Moreover, the transitions occur with a fast exchange rate and in a narrow force range. Therefore, the real applied forces can be determined accurately based on the force-dependent transitions. The equilibrium forces are insensitive to buffer conditions and temperature, making the calibration applicable to many complicated experimental systems. We provide an example to calibrate protein unfolding forces using this force marker and expect that this method can greatly simplify force calibration in single-molecule force spectroscopy experiments and improve the force accuracy.

Protein Hydrogels: The Swiss Army Knife for Enhanced Mechanical and Bioactive Properties of Biomaterials

Biomaterials are dynamic tools with many applications: from the primitive use of bone and wood in the replacement of lost limbs and body parts, to the refined involvement of smart and responsive biomaterials in modern medicine and biomedical sciences. Hydrogels constitute a subtype of biomaterials built from water-swollen polymer networks. Their large water content and soft mechanical properties are highly similar to most biological tissues, making them ideal for tissue engineering and biomedical applications. The mechanical properties of hydrogels and their modulation have attracted a lot of attention from the field of mechanobiology. Protein-based hydrogels are becoming increasingly attractive due to their endless design options and array of functionalities, as well as their responsiveness to stimuli. Furthermore, just like the extracellular matrix, they are inherently viscoelastic in part due to mechanical unfolding/refolding transitions of folded protein domains. This review summarizes different natural and engineered protein hydrogels focusing on different strategies followed to modulate their mechanical properties. Applications of mechanically tunable protein-based hydrogels in drug delivery, tissue engineering and mechanobiology are discussed.

Reconstruction of mechanical unfolding and refolding pathways of proteins with atomic force spectroscopy and computer simulations

Most proteins in proteomes are large, typically consist of more than one domain and are structurally complex. This often makes studying their mechanical unfolding pathways challenging. Proteins composed of tandem repeat domains are a subgroup of multi-domain proteins that, when stretched, display a saw-tooth pattern in their mechanical unfolding force extension profiles due to their repetitive structure. However, the assignment of force peaks to specific repeats undergoing mechanical unraveling is complicated because all repeats are similar and they interact with their neighbors and form a contiguous tertiary structure. Here, we describe in detail a combination of experimental and computational single-molecule force spectroscopy methods that proved useful for examining the mechanical unfolding and refolding pathways of ankyrin repeat proteins. Specifically, we explain and delineate the use of atomic force microscope-based single molecule force spectroscopy (SMFS) to record the mechanical unfolding behavior of ankyrin repeat proteins and capture their unusually strong refolding propensity that is responsible for generating impressive refolding force peaks. We also describe Coarse Grain Steered Molecular Dynamic (CG-SMD) simulations which complement the experimental observations and provide insights in understanding the unfolding and refolding of these proteins. In addition, we advocate the use of novel coiled-coils-based mechanical polypeptide probes which we developed to demonstrate the vectorial character of folding and refolding of these repeat proteins. The combination of AFM-based SMFS on native and CC-equipped proteins with CG-SMD simulations is powerful not only for ankyrin repeat polypeptides, but also for other repeat proteins and more generally to various multidomain, non-repetitive proteins with complex topologies.

Co-translational folding of nascent polypeptides: Multi-layered mechanisms for the efficient biogenesis of functional proteins

The coupling of protein synthesis and folding is a crucial yet poorly understood aspect of cellular protein folding. Over the past few years, it has become possible to experimentally follow and define protein folding on the ribosome, revealing principles that shape co-translational folding and distinguish it from refolding in solution. Here, we highlight some of these recent findings from biochemical and biophysical studies and their potential significance for cellular protein biogenesis. In particular, we focus on nascent chain interactions with the ribosome, interactions within the nascent protein, modulation of translation elongation rates, and the role of mechanical force that accompanies nascent protein folding. The ability to obtain mechanistic insight in molecular detail has set the stage for exploring the intricate process of nascent protein folding. We believe that the aspects discussed here will be generally important for understanding how protein synthesis and folding are coupled and regulated.

Persistent repression of tau in the brain using engineered zinc finger protein transcription factors

Neuronal tau reduction confers resilience against β-amyloid and tau-related neurotoxicity in vitro and in vivo. Here, we introduce a novel translational approach to lower expression of the tau gene MAPT at the transcriptional level using gene-silencing zinc finger protein transcription factors (ZFP-TFs). Following a single administration of adeno-associated virus (AAV), either locally into the hippocampus or intravenously to enable whole-brain transduction, we selectively reduced tau messenger RNA and protein by 50 to 80% out to 11 months, the longest time point studied. Sustained tau lowering was achieved without detectable off-target effects, overt histopathological changes, or molecular alterations. Tau reduction with AAV ZFP-TFs was able to rescue neuronal damage around amyloid plaques in a mouse model of Alzheimer's disease (APP/PS1 line). The highly specific, durable, and controlled knockdown of endogenous tau makes AAV-delivered ZFP-TFs a promising approach for the treatment of tau-related human brain diseases.

Protein folding modulates the chemical reactivity of a Gram-positive adhesin

Gram-positive bacteria colonize mucosal tissues, withstanding large mechanical perturbations such as coughing, which generate shear forces that exceed the ability of non-covalent bonds to remain attached. To overcome these challenges, the pathogen Streptococcus pyogenes utilizes the protein Cpa, a pilus tip-end adhesin equipped with a Cys-Gln thioester bond. The reactivity of this bond towards host surface ligands enables covalent anchoring; however, colonization also requires cell migration and spreading over surfaces. The molecular mechanisms underlying these seemingly incompatible requirements remain unknown. Here we demonstrate a magnetic tweezers force spectroscopy assay that resolves the dynamics of the Cpa thioester bond under force. When folded at forces <6 pN, the Cpa thioester bond reacts reversibly with amine ligands, which are common in inflammation sites; however, mechanical unfolding and exposure to forces >6 pN block thioester reformation. We hypothesize that this folding-coupled reactivity switch (termed a smart covalent bond) could allow the adhesin to undergo binding and unbinding to surface ligands under low force and remain covalently attached under mechanical stress.

Mapping Distinct Sequences of Structure Formation Differentiating Multiple Folding Pathways of a Small Protein

To determine experimentally how the multiple folding pathways of a protein differ, in the order in which the structural parts are assembled, has been a long-standing challenge. To resolve whether structure formation during folding can progress in multiple ways, the complex folding landscape of monellin has been characterized, structurally and temporally, using the multisite time-resolved FRET methodology. After an initial heterogeneous polypeptide chain collapse, structure formation proceeds on parallel pathways. Kinetic analysis of the population evolution data across various protein segments provides a clear structural distinction between the parallel pathways. The analysis leads to a phenomenological model that describes how and when discrete segments acquire structure independently of each other in different subensembles of protein molecules. When averaged over all molecules, structure formation is seen to progress as α-helix formation, followed by core consolidation, then β-sheet formation, and last end-to-end distance compaction. Parts of the protein that are closer in the primary sequence acquire structure before parts separated by longer sequence.

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long-range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.

Hidden Intermediate State and Second Pathway Determining Folding and Unfolding Dynamics of GB1 Protein at Low Forces

Atomic force microscopy experiments found that GB1, a typical two-state model protein used for study of folding and unfolding dynamics, can sustain forces of more than 100 pN, but its response to low forces still remains unclear. Using ultrastable magnetic tweezers, we discovered that GB1 has an unexpected nonmonotonic force-dependent unfolding rate at 5-160 pN, from which a free energy landscape with two main barriers and a hidden intermediate state was constructed. A model combining two separate models by Dudko et al. with two pathways between the native state and this intermediate state is proposed to rebuild the unfolding dynamics over the full experimental force range. One candidate of this transient intermediate state is the theoretically proposed molten globule state with a loosely collapsed conformation, which might exist universally in the folding and unfolding processes of two-state proteins.

The SecA motor generates mechanical force during protein translocation

The Sec translocon moves proteins across lipid bilayers in all cells. The Sec channel enables passage of unfolded proteins through the bacterial plasma membrane, driven by the cytosolic ATPase SecA. Whether SecA generates mechanical force to overcome barriers to translocation posed by structured substrate proteins is unknown. Here, we kinetically dissect Sec-dependent translocation by monitoring translocation of a folded substrate protein with tunable stability at high time resolution. We find that substrate unfolding constitutes the rate-limiting step during translocation. Using single-molecule force spectroscopy, we also define the response of the protein to mechanical force. Relating the kinetic and force measurements reveals that SecA generates at least 10 piconewtons of mechanical force to actively unfold translocating proteins, comparable to cellular unfoldases. Combining biochemical and single-molecule measurements thus allows us to define how the SecA motor ensures efficient and robust export of proteins that contain stable structure.

The ATPase SecA drives Sec-dependent protein translocation across the bacterial plasma membrane. Here, the authors combine kinetic translocation measurements with single-molecule force spectroscopy and demonstrate that the SecA motor generates mechanical force to unfold and translocate preproteins.

Single-Molecule Studies of Protein Folding with Optical Tweezers

Manipulation of individual molecules with optical tweezers provides a powerful means of interrogating the structure and folding of proteins. Mechanical force is not only a relevant quantity in cellular protein folding and function, but also a convenient parameter for biophysical folding studies. Optical tweezers offer precise control in the force range relevant for protein folding and unfolding, from which single-molecule kinetic and thermodynamic information about these processes can be extracted. In this review, we describe both physical principles and practical aspects of optical tweezers measurements and discuss recent advances in the use of this technique for the study of protein folding. In particular, we describe the characterization of folding energy landscapes at high resolution, studies of structurally complex multidomain proteins, folding in the presence of chaperones, and the ability to investigate real-time cotranslational folding of a polypeptide.

The generic unfolding of a biomimetic polymer during force spectroscopy

With the help of force spectroscopy, several analytical theories aim at estimating the rate coefficient of folding for various proteins. Nevertheless, a chief bottleneck lies in the fact that there is still no perfect consensus on how does a force generally perturb the crystal-coil transition. Consequently, the goal of our work is in clarifying the generic behavior of most proteins in force spectroscopy; in other words, what general signature does an arbitrary protein exhibit for its rate coefficient as a function of the applied force? By employing a biomimetic polymer in molecular simulations, we focus on evaluating its respective activation energy for unfolding, while pulling on various pairs of its monomers. Above all, we find that in the vicinity of the force-free scenario, this activation energy possesses a negative slope and a negative curvature as a function of the applied force. Our work is in line with the most recent theories for unfolding, which suggest that such a signature is expected for most proteins, and thus, we further reiterate that many of the classical formulae, that estimate the rate coefficient of the crystal-coil transition, are inadequate. Besides, we also present here an analytical expression which experimentalists can use for approximating the activation energy for unfolding; importantly, it is based on measurements for the mean and variance of the distance between the beads which are being pulled. In summary, our work presents an interesting view for protein folding in force spectroscopy.

Principles for enhancing virus capsid capacity and stability from a thermophilic virus capsid structure

The capsids of double-stranded DNA viruses protect the viral genome from the harsh extracellular environment, while maintaining stability against the high internal pressure of packaged DNA. To elucidate how capsids maintain stability in an extreme environment, we use cryoelectron microscopy to determine the capsid structure of thermostable phage P74-26 to 2.8-Å resolution. We find P74-26 capsids exhibit an overall architecture very similar to those of other tailed bacteriophages, allowing us to directly compare structures to derive the structural basis for enhanced stability. Our structure reveals lasso-like interactions that appear to function like catch bonds. This architecture allows the capsid to expand during genome packaging, yet maintain structural stability. The P74-26 capsid has T = 7 geometry despite being twice as large as mesophilic homologs. Capsid capacity is increased with a larger, flatter major capsid protein. Given these results, we predict decreased icosahedral complexity (i.e. T ≤ 7) leads to a more stable capsid assembly.

Viral capsids need to protect the genome against harsh environmental conditions and cope with high internal pressure from the packaged genome. Here, the authors determine the structure of the thermostable phage P74-26 capsid at 2.8-Å resolution and identify features underlying enhanced capsid capacity and structural stability.

A modular DNA scaffold to study protein-protein interactions at single-molecule resolution

The residence time of a drug on its target has been suggested as a more pertinent metric of therapeutic efficacy than the traditionally used affinity constant. Here, we introduce junctured-DNA tweezers as a generic platform that enables real-time observation, at the single-molecule level, of biomolecular interactions. This tool corresponds to a double-strand DNA scaffold that can be nanomanipulated and on which proteins of interest can be engrafted thanks to widely used genetic tagging strategies. Thus, junctured-DNA tweezers allow a straightforward and robust access to single-molecule force spectroscopy in drug discovery, and more generally in biophysics. Proof-of-principle experiments are provided for the rapamycin-mediated association between FKBP12 and FRB, a system relevant in both medicine and chemical biology. Individual interactions were monitored under a range of applied forces and temperatures, yielding after analysis the characteristic features of the energy profile along the dissociation landscape.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

Multiplexed protein force spectroscopy reveals equilibrium protein folding dynamics and the low-force response of von Willebrand factor

Single-molecule force spectroscopy has provided unprecedented insights into protein folding, force regulation, and function. So far, the field has relied primarily on atomic force microscope and optical tweezers assays that, while powerful, are limited in force resolution, throughput, and require feedback for constant force measurements. Here, we present a modular approach based on magnetic tweezers (MT) for highly multiplexed protein force spectroscopy. Our approach uses elastin-like polypeptide linkers for the specific attachment of proteins, requiring only short peptide tags on the protein of interest. The assay extends protein force spectroscopy into the low force (<1 pN) regime and enables parallel and ultra-stable measurements at constant forces. We present unfolding and refolding data for the small, single-domain protein ddFLN4, commonly used as a molecular fingerprint in force spectroscopy, and for the large, multidomain dimeric protein von Willebrand factor (VWF) that is critically involved in primary hemostasis. For both proteins, our measurements reveal exponential force dependencies of unfolding and refolding rates. We directly resolve the stabilization of the VWF A2 domain by Ca²⁺ and discover transitions in the VWF C domain stem at low forces that likely constitute the first steps of VWF's mechano-activation. Probing the force-dependent lifetime of biotin-streptavidin bonds, we find that monovalent streptavidin constructs with specific attachment geometry are significantly more force stable than commercial, multivalent streptavidin. We expect our modular approach to enable multiplexed force-spectroscopy measurements for a wide range of proteins, in particular in the physiologically relevant low-force regime.

Mechanical unfolding of spectrin reveals a super-exponential dependence of unfolding rate on force

We investigated the mechanical unfolding of single spectrin molecules over a broad range of loading rates and thus unfolding forces by combining magnetic tweezers with atomic force microscopy. We find that the mean unfolding force increases logarithmically with loading rate at low loading rates, but the increase slows at loading rates above 1pN/s. This behavior indicates an unfolding rate that increases exponentially with the applied force at low forces, as expected on the basis of one-dimensional models of protein unfolding. At higher forces, however, the increase of the unfolding rate with the force becomes faster than exponential, which may indicate anti-Hammond behavior where the structures of the folded and transition states become more different as their free energies become more similar. Such behavior is rarely observed and can be explained by either a change in the unfolding pathway or as a reflection of a multidimensional energy landscape of proteins under force.

Nonnative structure in a peptide model of the unfolded state of superoxide dismutase 1 (SOD1): Implications for ALS-linked aggregation

Dozens of mutations throughout the sequence of the gene encoding superoxide dismutase 1 (SOD1) have been linked to toxic protein aggregation in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). A parsimonious explanation for numerous genotypes resulting in a common phenotype would be mutation-induced perturbation of the folding free-energy surface that increases the populations of high-energy states prone to aggregation. The absence of intermediates in the folding of monomeric SOD1 suggests that the unfolded ensemble is a potential source of aggregation. To test this hypothesis, here we dissected SOD1 into a set of peptides end-labeled with FRET probes to model the local behavior of the corresponding sequences in the unfolded ensemble. Using time-resolved FRET, we observed that the peptide corresponding to the Loop VII-β8 sequence at the SOD1 C terminus was uniquely sensitive to denaturant. Utilizing a two-dimensional form of maximum entropy modeling, we demonstrate that the sensitivity to denaturant is the surprising result of a two-state-like transition from a compact to an expanded state. Variations of the peptide sequence revealed that the compact state involves a nonnative interaction between the disordered N terminus and the hydrophobic C terminus of the peptide. This nonnative intramolecular structure could serve as a precursor for intermolecular association and result in aggregation associated with ALS. We propose that this precursor would provide a common molecular target for therapeutic intervention in the dozens of ALS-linked SOD1 mutations.

Using Single-Molecule Chemo-Mechanical Unfolding to Simultaneously Probe Multiple Structural Parameters in Protein Folding

While single-molecule force spectroscopy has greatly advanced the study of protein folding, there are limitations to what can be learned from studying the effect of force alone. We developed a novel technique, chemo-mechanical unfolding, that combines multiple perturbants—force and chemical denaturant—to more fully characterize the folding process by simultaneously probing multiple structural parameters—the change in end-to-end distance, and solvent accessible surface area. Here, we describe the theoretical background, experimental design, and data analysis for chemo-mechanical unfolding experiments probing protein folding thermodynamics and kinetics. This technique has been applied to characterize parallel protein folding pathways, the protein denatured state, protein folding on the ribosome, and protein folding intermediates.

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.

Methylation Status of the Adeno-Associated Virus Type 2 (AAV2)

To analyze the methylation status of wild-type adeno-associated virus type 2 (AAV2), bisulfite PCR sequencing (BPS) of the packaged viral genome and its integrated form was performed and 262 of the total 266 CG dinucleotides (CpG) were mapped. In virion-packaged DNA, the ratio of the methylated cytosines ranged between 0⁻1.7%. In contrast, the chromosomally integrated AAV2 genome was hypermethylated with an average of 76% methylation per CpG site. The methylation level showed local minimums around the four known AAV2 promoters. To study the effect of methylation on viral rescue and replication, the replication initiation capability of CpG methylated and non-CpG methylated AAV DNA was compared. The in vitro hypermethylation of the viral genome does not inhibit its rescue and replication from a plasmid transfected into cells. This insensitivity of the viral replicative machinery to methylation may permit the rescue of the integrated heavily methylated AAV genome from the host's chromosomes.

Three Weaknesses for Three Perturbations: Comparing Protein Unfolding Under Shear, Force, and Thermal Stresses

The perturbation of a protein conformation by a physiological fluid flow is crucial in various biological processes including blood clotting and bacterial adhesion to human tissues. Investigating such mechanisms by computer simulations is thus of great interest, but it requires development of ad hoc strategies to mimic the complex hydrodynamic interactions acting on the protein from the surrounding flow. In this study, we apply the Lattice Boltzmann Molecular Dynamics (LBMD) technique built on the implicit solvent coarse-grained model for protein Optimized Potential for Efficient peptide structure Prediction (OPEP) and a mesoscopic representation of the fluid solvent, to simulate the unfolding of a small globular cold-shock protein in shear flow and to compare it to the unfolding mechanisms caused either by mechanical or thermal perturbations. We show that each perturbation probes a specific weakness of the protein and causes the disruption of the native fold along different unfolding pathways. Notably, the shear flow and the thermal unfolding exhibit very similar pathways, while because of the directionality of the perturbation, the unfolding under force is quite different. For force and thermal disruption of the native state, the coarse-grained simulations are compared to all-atom simulations in explicit solvent, showing an excellent agreement in the explored unfolding mechanisms. These findings encourage the use of LBMD based on the OPEP model to investigate how a flow can affect the function of larger proteins, for example, in catch-bond systems.

Structural-elastic determination of the force-dependent transition rate of biomolecules

The force-dependent unfolding/refolding of protein domains and ligand-receptor association/dissociation are crucial for mechanosensitive functions, while many aspects of how force affects the transition rate still remain poorly understood. Here, we report a new analytical expression of the force-dependent rate of molecules for transitions overcoming a single barrier. Unlike previous models derived in the framework of Kramers theory that requires a presumed one-dimensional free energy landscape, our model is derived based on the structural-elastic properties of molecules which are not restricted by the shape and dimensionality of the underlying free energy landscape. Importantly, the parameters of this model provide direct information on the structural-elastic features of the molecules between their transition and initial states. We demonstrate the applications of this model by applying it to explain force-dependent transition kinetics for several molecules and predict the structural-elastic properties of the transition states of these molecules.

Rationally designed synthetic protein hydrogels with predictable mechanical properties

Designing synthetic protein hydrogels with tailored mechanical properties similar to naturally occurring tissues is an eternal pursuit in tissue engineering and stem cell and cancer research. However, it remains challenging to correlate the mechanical properties of protein hydrogels with the nanomechanics of individual building blocks. Here we use single-molecule force spectroscopy, protein engineering and theoretical modeling to prove that the mechanical properties of protein hydrogels are predictable based on the mechanical hierarchy of the cross-linkers and the load-bearing modules at the molecular level. These findings provide a framework for rationally designing protein hydrogels with independently tunable elasticity, extensibility, toughness and self-healing. Using this principle, we demonstrate the engineering of self-healable muscle-mimicking hydrogels that can significantly dissipate energy through protein unfolding. We expect that this principle can be generalized for the construction of protein hydrogels with customized mechanical properties for biomedical applications.

Denaturants Alter the Flux through Multiple Pathways in the Folding of PDZ Domain

Although we understand many aspects of how small proteins (number of residues less than about hundred) fold, it is a major challenge to quantitatively describe how large proteins self-assemble. To partially overcome this challenge, we performed simulations using the self-organized polymer model with side chains (SOP-SC) in guanidinium chloride (GdmCl), using the molecular transfer model (MTM), to describe the folding of the 110-residue PDZ3 domain. The simulations reproduce the folding thermodynamics accurately including the melting temperature (T_m), the stability of the folded state with respect to the unfolded state. We show that the calculated dependence of ln k_obs (k_obs is the relaxation rate) has the characteristic chevron shape. The slopes of the chevron plots are in good agreement with experiments. We show that PDZ3 folds by four major pathways populating two metastable intermediates, in accord with the kinetic partitioning mechanism. The structure of one of the intermediates, populated after polypeptide chain collapse, is structurally similar to an equilibrium intermediate. Surprisingly, the connectivities between the intermediates and hence, the fluxes through the pathways depend on the concentration of GdmCl. The results are used to predict possible outcomes for unfolding of PDZ domain subject to mechanical forces. Our study demonstrates that, irrespective of the size or topology, simulations based on MTM and SOP-SC offer a theoretical framework for describing the folding of proteins, mimicking precisely the conditions used in experiments.

Investigations on the elasticity of functional gold nanoparticles using single-molecule force spectroscopy

In this focused review we turn our attention towards several approaches for detecting the elasticity of NPs, systematically summarizing the divergent elasticity values of distinct gold nanoparticles (AuNPs) with different surfaces.