Nano Letters
Home Molecular fluctuations as a ruler of force-induced protein conformations
Molecular fluctuations as a ruler of force-induced protein conformations
Molecular fluctuations as a ruler of force-induced protein conformations

Article Type: Research Article Article History
Publisher: None
Abstract

Molecular fluctuations directly reflect the underlying energy landscape. Variance analysis examines protein dynamics in several biochemistry-driven approaches, yet measurement of probe-independent fluctuations in proteins exposed to mechanical forces remains only accessible through steered molecular dynamics simulations. Using single molecule magnetic tweezers, here we conduct variance analysis to show that individual unfolding and refolding transitions occurring in dynamic equilibrium in a single protein under force are hallmarked by a change in the protein’s end-to-end fluctuations, revealing a change in protein stiffness. By unfolding and refolding three structurally distinct proteins under a wide range of constant forces, we demonstrate that the associated change in protein compliance to reach force-induced thermodynamically-stable states scales with the protein’s contour length, in agreement with the sequence-independent freely-jointed chain model of polymer physics. Our findings will help elucidate the conformational dynamics of proteins exposed to mechanical force at high resolution, of central importance in mechanosensing and mechanotransduction.

Keywords
Stannard,Mora,Beedle,Castro-López,Board,and Garcia-Manyes: Molecular fluctuations as a ruler of force-induced protein conformations

Mechanical unfolding of proteins regulates a variety of biological processes, encompassing muscle elasticity1, focal adhesion mechanosensing2, 3 and nuclear translocation4. Whether working as monomers or in form of polyproteins5, mechanical unfolding provides an efficient shock absorbing strategy that works as a safety mechanism to preserve the mechanical integrity of adhesive and mechanotransducing protein complexes6. Several strategies —including the protein’s intrinsic topology, the introduction of point mutations, ligand binding, the pulling direction or the dynamic rupture/reformation of covalent bonds— have been reported to govern the mechanical resistance of proteins in their folded form7. Once unfolded, mechanically stretched proteins exhibit previously cryptic sites to the solvent, which become suddenly available to protein binding or post-translational modifications, ultimately affecting protein elasticity8. Given the emerging evidence of conformation-dependent functionality9, it is enticing to experimentally capture the mechanical properties of the different protein forms, collectively defining the rich conformational repertoire of proteins under force. For example, the so far elusive measurement of protein stiffness could provide an additional tool, further to the unfolding force and contour length increment measurement7, to characterise the diverse conformational dynamics of proteins, using the spatial magnitude of the fluctuations of each force-induced equilibrium conformation as its structural fingerprint.

Single molecule nanomechanical experiments have provided a detailed in-vitro insight of the molecular determinants governing the conformational dynamics of proteins under force8, 10. Collectively, these measurements have concentrated on providing an extensive characterization of the mechanical stability of the native state of proteins11, and on capturing the kinetics and out-of-equilibrium transition paths underpinning individual unfolding and refolding trajectories12. By contrast, a thorough experimental description of the extended states —often perceived as a single, featureless, compliant conformation —has comparatively received much less attention. This largely stems from the lack of experimental approaches able to capture the evolution of protein fluctuations with force, which emerge naturally from a general and rather simplified picture of the folding energy landscape that defines the behaviour of an entropic ideal polymer under force (predicted by models of protein elasticity13 e.g. the freely-jointed chain (FJC) or worm-like chain (WLC)) added to an enthalpic, short-range attractive contribution (e.g. Morse or Lennard-Jones potential)14. In the absence of force, proteins populate this latter energy basin, which defines the native state of the protein (Fig. 1). The application of a constant force tilts the 1D energy landscape and creates a new energy minimum that corresponds to the protein’s unfolded state. Increasing the stretching force lowers the energy of the unfolded well and displaces its position toward increasingly extended conformations, as predicted by the FJC/WLC model. Consequently, the width of the unfolded minimum decreases when increasing the stretching force. Hence, the measurement of the length fluctuations as a function of force should provide an accurate and direct measurement of protein stiffness, which is a direct, natural consequence of the curvature of the energy well (Fig. 1). Despite its theoretical simplicity, a systematic experimental quantification of the dependency of protein stiffness with protein extension has evaded characterization, mainly due to the technical inability to capture the natural fluctuation modes of the protein when exposed to the wide range of stretching forces that fingerprint its entire conformational space, encompassing both the native and unfolded conformations.

Here we employ single molecule magnetic tweezers1518, which enable (i) capturing of the unfolding/refolding trajectories of the same single protein over extended (hour-long) periods of time with high stability16, (ii) under a wide range of constant forces (2-100 pN), (iii) in passive clamp mode (hence avoiding potential masking of the protein fluctuations due to the active electronic feedback19), (iv) in the absence of long tethers (which necessarily add additional inherent fluctuations to the protein construct20), and, most importantly, (v) in an experimental set-up where the measured stiffness is dominated by the intrinsic stiffness of the protein-of-interest (and not the pulling device, such as e.g. the AFM cantilever2123, Fig. S1 and SI text) to capture the change in compliance upon mechanical unfolding/refolding of individual proteins, in both monomeric and polymeric forms. Our experiments enable direct observation of the change in protein stiffness upon mechanical unfolding and over a continuum range of force-induced conformations in thermodynamic equilibrium, and demonstrate that molecular fluctuations are a reliable reporter of force-induced protein conformations.

Results

Using single molecule magnetic tweezers, we monitored the unfolding and refolding trajectories of an individual protein L (PL) monomer bracketed between two mechanically rigid Ig32 titin domains either side. The resulting (Ig32)2-PL-(Ig32)2 polyprotein is covalently attached to the glass substrate through a Halo-Tag on the N-terminus, and to the paramagnetic bead via biotin-streptavidin binding on the C-terminus (Fig. 2A and Methods). Using a force quench protocol, we first applied a constant force of 38 pN, which elicited the unfolding of the individual PL monomer, fingerprinted by a step elongation of the protein by Δ<x> = 16.8 nm 16. We note that, at this force, we did not observe the unfolding of the Ig32 domains, ensuring that the folding dynamics of the PL monomer are not masked by the dynamics of the stiff Ig32 molecular handles (Fig. S2). After 60 seconds, the force was quenched down to 8.1 pN –the force, F 0.5, at which PL exhibits a 50% folded probability– and held constant for 60 minutes. During this period of time, we observed the stochastic unfolding and refolding of the PL monomer (hallmarked by an average step-change in extension of Δ<x> = 10.3 ± 0.3 nm) occurring with approximate equal rates16, Fig. 2B. The protein was finally stretched back again at 38 pN to ensure that, over this long pulling time, the protein had not misfolded.

To experimentally measure changes in the mechanical properties of PL each time the monomer hops between the folded and unfolded states, we first consider that the measured stiffness, k M, of the entire series construct (combining the stiffness of PL, k P, and the rest of the construct, k C) is given by 1/k M = 1/k P + 1/k C. As such, upon un/folding, the change in measured stiffness does not correspond to the change in PL stiffness (Δk M ≠ Δk P). If instead we consider compliance –the inverse of stiffness, c = 1/k– we have the simple situation where c M = c P + c C and, crucially, the change in measured compliance is the change in PL compliance (Δc M = Δc P). In the simple case of the classic Hookean spring, equipartition theory implies that constant (force-independent) compliance is given by (1)

c=βσ2 (1)
where β = 1/k B T is the inverse of thermal energy, and σ 2 = <x 2> – <x>2 is the variance in end-to-end length, x. The trivial consequence is that the length fluctuations are force independent. By contrast, in the case of an entropic spring (e.g. in the FJC and WLC models) where the average length <x> is non-linear with force, compliance is force-dependent since c ≡ d<x>/dF. Even in this case, Eq. 1 still holds true (SI text), implying that fluctuations are force-dependent, being lower at higher forces, and can be used to determine compliance of non-linear springs such as proteins. Combining these two arguments, we see that the change in PL compliance can be found from the change in measured variance, Δc P = βΔσ M 2.

We next processed the hopping trajectory of Fig. 2B by flattening individual sections to remove the effects of slow timescale drift of the magnetic bead relative to the reference bead (<10 nm/hour average). Applying Eq. 1 to the resulting trajectory (Fig. 2C) enables direct measurement of the protein’s compliance (Fig. 2D) corresponding to the unfolded and folded states before and after each individual (un)folding transition24 (SI text and Fig. S3). By averaging the compliance for both, the natively-folded (blue) and unfolded (red) states underpinning the 13 individual transitions in Fig. 2, we obtained an increase in protein compliance Δc = 0.86 ± 0.09 nm/pN upon PL unfolding (see also Fig. S4). Of note, the fluctuations underpinning the (un)folding dynamics of the PL monomer are independent of the rest of the polyprotein construct as observed when replacing the stiff Ig32 handles by the inextensible Spy0128 protein (Fig. S5). Conceptually, this increase in compliance corresponds to the transition of the protein, at a constant force (8.1 pN), from the narrow basin defining the native state to the much-broader well defining the unfolded state of the protein (Fig. 1).

Having obtained a measurable change in protein compliance upon mechanical unfolding of an individual protein monomer, we conjectured that, in the more complex case of a polyprotein (Fig. 3A), we might be able to single out the (un)folding of each individual domain by directly measuring the fluctuation variance of the entire construct at a given time and under constant force. To test this hypothesis, we repeated our force-quench protocol on a PL8 polyprotein (Fig. 3B). The long quench pulse at 8.1 pN exhibited complex dynamics, underpinning the stochastic unfolding/refolding transitions of each distinct protein monomer forming the polyprotein chain (Fig. 3C). Measuring the variance of the overall fluctuations of the polyprotein construct at each state (number of unfolded domains) level revealed a linear correlation between the construct’s compliance and the number of unfolded domains (Fig. 3D and Fig S6). Notably, the slope of the trend (Δc/domain = 0.81 ± 0.15 nm/pN) is in very good agreement with that obtained for the monomeric protein (Fig. 2D). These results further highlight that the different monomers forming a polyprotein chain fold independently, implying that the folding dynamics of an individual monomer is independent of the composition of the rest of the multidomain construct25, 26. Most importantly, our results reveal that the dynamic change in compliance of a multi-protein construct can be used as a direct proxy to identify the number of folded (and unfolded) modules within an actively folding homopolyprotein at a given time.

In the equilibrium experiments reported thus far, the pulling force was held constant to acquire sufficient statistics of the two protein states as the protein hopped between the unfolded/folded conformations. To systematically study the effect of force on protein compliance, we expanded the range of probed forces by conducting successive unfolding/refolding cycles on the same single protein monomer (Fig. 2A), where the unfolding/refolding force was varied in the range 4.3-22 pN (Figs. 4A&B). We then measured the relative change in the fluctuation variance (and the related change in protein compliance) for each particular force (Fig. 4C&D). These results clearly demonstrate that, when PL is mechanically unfolded, the greater the change in average extension (due to higher force being applied), the lesser the change in compliance. Finally, our resulting data on the extension-dependent change in compliance was compared to the prediction by the FJC model – which has analytical solutions for Δ<x> and Δc (Supplementary Information). Our experimental data, including both the monomer and polymer measurements, quantitatively agree with the FJC model (Fig. 4D), confirming that the change in compliance (ultimately measured through the protein’s end-to-end length fluctuations) is reduced as the pulling forces are increased, suggesting that the findings might be general and independent of the studied protein since the FJC has no consideration for protein structure.

A closer look to the analytical expressions obtained from the FJC model (Supplementary Information) reveal that both the force-dependent extension change, and the force-dependent compliance change, are predicted to be linear with contour length increment. To experimentally test this hypothesis, we studied the change in compliance as an individual monomer IVVI mutant of the R3 rod domain of talin – a focal adhesion sensor3 featuring alpha-helix topology (Fig. 5A) and an associated increment in contour length of ΔL ≈ 37 nm27, 28 (which is significantly larger than that of PL, ΔL = 18.6 nm29) – hops in equilibrium between the folded and unfolded states at a constant force of 8.35 pN30 (Fig. 5B). Analysis of the hopping trajectories (following the same procedure used in Fig. 2B,C) reveals a significantly larger average increase in protein compliance Δc = 1.28 ± 0.08 nm/pN upon unfolding.

A particular characteristics of the R3 domain of talin is its fast (un)folding/dynamics and its steep force-dependency for both the unfolding and refolding reactions30, which implies that at forces only slightly higher [lower] than the hopping (F 0.5) force, the lifetime of the folded [unfolded] state becomes too short-lived to make a thorough measurement of compliance through the length variance (Fig. S7). To circumvent this limitation, we moved to conducting folding measurements on the spectrin 73 repeat (SR73) of nesprin, a nuclear envelope mechanosensor protein that has been suggested to play a key role in propagating extracellular mechanical forces to the cell nucleus31 likely through mechanical unfolding32 and which exhibits a comparable contour length increment (ΔL = 34.2 nm, Fig. S8C). Stretching an individual [Ig27-SR73]4 polyprotein (Fig. 5C) in a force-ramp (initial pulse) elicits the stepwise unfolding of the 4 individual SR73 nesprin domains, concomitant to an average increase in the protein length of Δ<x>= 28.6 nm (Fig. 5D, Fig. S8A,C). After complete unfolding of the nesprin SR73 domains with force ramped up to 38 pN—at this force we did not observe any signature of unfolding of the Ig27 marker (Fig. S2)—, the force was then quenched down to a lower value of 7.4 pN to trigger protein folding. We chose this force value because, according to the force-dependent folding probability that we measured for the nesprin SR73 domains (Fig. S8A,B), nesprin folds slow enough to directly capture the independent folding transitions corresponding to the sequential refolding of each individual SR73 monomer. Figure 5D shows how, akin to PL16 and to several domains of titin33, the previously mechanically uncharacterised nesprin mechanosensor effectively refolds against a constant force (7.4 pN), marked by the sequential stepwise shrinking (Δ<x> = 16.8 ± 0.2 nm) of the protein construct. To elucidate whether the sequential stepwise recoil of the protein (Fig. 5E) is mirrored by a step decrease in the end-to-end fluctuations of the polyprotein as each individual monomer refolds, we conducted variance analysis of the refolding trajectory over time and found that, on average the folding of each individual domain results in a concomitant change of the construct’s compliance by Δc = 1.99 ± 0.12 nm/pN (Fig. 5F and Fig. S8D,E). Normalization of the measured change in compliance with the normalised change in extension for proteins with markedly different contour lengths and over a wide range of forces collapses data onto a FJC master curve, further demonstrating that the change in protein compliance is force- and length-dependent (through the FJC relationship) but—at least for the different proteins tested here—protein(sequence)- independent (Fig. 5G).

Discussion

Analysis of molecular fluctuations has enabled capturing the conformational dynamics of proteins using a number of complementary experimental techniques, such as nuclear magnetic resonance (NMR) and small-angle X-ray scattering (SAXS)34, 35. At the single molecule level, sm-FRET has provided invaluable insight into the broad conformational heterogeneity and dynamics of the unfolded state of proteins, as well as of intrinsically disordered motifs36. More recently, interferometric scattering microscopy (iSCAT) has shown its ability to capture conformational transitions, provided they are associated to a change in the protein’s center of mass37. However, none of these experimental approaches can capture the effect of mechanical load on the conformational dynamics of individual proteins.

When denatured by force, proteins undergo a conformational change that brings the mechanically stiff, natively folded form into a stretched, compliant and unfolded conformation that cannot be populated in classical biochemical experiments using temperature or chemical denaturants7, 38. Single molecule nanomechanical techniques, mostly using AFM in the high force-regime7 and optical tweezers to probe proteins with lower mechanical stabilities3941, have enabled sampling the regions of the protein free energy landscape that cannot be accessed through any other experimental means, and have rationalized the effect of force on modulating the conformational dynamics of an individual protein during its folding journey.

In this vein, earlier experiments using force-clamp AFM showed that, during collapse from highly-extended states, proteins exhibit increased end-to-end length fluctuations, reminiscent of the temporal occupancy of locally broad energy minima, before reaching their natively folded length 14, 25, 42, 43. Conformational heterogeneity was also captured at the molten globule44 and native states45 levels, however these experiments required using the mechanical stability of each conformation as its indirect structural fingerprint. Nonetheless, the high applied forces, the intrinsic stiffness (and mass) of the AFM probe and its associated force uncertainty (and decreased diffusion coefficient21), in addition to the active PID feedback, have precluded direct capture of the intrinsic fluctuations and underlying stiffness defining each force-dependent protein conformation. Our experiments using magnetic tweezers (where the stiffness of the protein is not masked by the pulling probe) demonstrate that the end-to-end fluctuations of a single protein under force can be used as a direct reporter of its precise conformation and folding status, hence introducing an additional experimentally measurable dimension to the sudden changes in protein length accompanying protein unfolding/refolding events. Noteworthy, these measurements are independent of the mechanical stability or the fold of the studied protein. The approach presented here does not require knowledge of the exact protein conformation of the protein before and after a transition, or even whether that transition has fleetingly-occupied intermediate states, and it is applicable to any conformational change that involves a change in the protein’s contour length. This methodology may be of particular relevance to identify time-dependent events occurring without obvious changes in the protein’s end-to-end length, such as ligand/drug binding46 and post-translational modifications4749 that place dihedral angle constrictions to the force-induced protein conformation50, 51. For example, binding of vinculin to unfolded talin52 and α-catenin17 is a central step in cellular mechanotransduction. Importantly, all these functionalities can be finely regulated by the pulling force through precise modulation of protein stiffness and conformation, as revealed, with atomic detail, by steered molecular dynamics (SMD) simulations53. It will be enticing to compare, for distinct proteins, how our experimental results agree with the equilibrium fluctuations measured for equivalent force-induced protein conformations in SMD simulations53. More broadly, our experiments focusing on the size of natural molecular fluctuations defining the continuum force-induced thermodynamic equilibrium states complement recent elegant single molecule experiments concentrating on the non-equilibrium transition paths12, 54, 55 between those defined protein states.

Methods

Polyprotein engineering

The different protein constructs were designed, cloned and purified as reported elsewhere16. Details for each construct are described in detail in the Supplementary Information section.

Single-molecule magnetic tweezers experiments

Sample preparation and the single-molecule magnetic tweezers setup are based on those reported elsewhere16, with slight variations detailed in the the Supplementary Information section.

Acknowledgements

We thank Julio Fernández for his help in building the magnetic tweezers set-up and for sharing the plasmid for the talin polyprotein. We thank Ionel Popa for help in assemblying the instrument, and Ainhoa Lezamiz, Palma Rico-Lastres and Aisling Williams for the engineering of the nesprin polyprotein. This work was supported in part by the Francis Crick Institute which receives its core funding from Cancer Research UK (FC001002), the UK Medical Research Council (FC001002), and the Wellcome Trust (FC001002). This research was funded in part, by the Wellcome Trust [212218/Z/18/Z]. For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission. A.E.M.B. is recipient of a Sir Henry Wellcome fellowship (210887/Z/18/Z). This work was supported by the European Commission (Mechanocontrol, grant agreement 731957), EPSRC Fellowship K00641X/1, Leverhulme Trust Research Leadership Award RL-2016-015, Wellcome Trust Investigator Award 212218/Z/18/Z and Royal Society Wolfson Fellowship RSWF/R3/183006 to S.G.M.

Author contributions:

A.S. and S.G.M. designed research. A.E.M.B., M.M. and A.S. conducted experiments, M.L-C. built and A.E.M.B. calibrated the magnetic tweezers spectrometer, A.S. devised and performed data analysis, S.B. conducted molecular biology and S.G-M. wrote the paper. All authors participated in the discussion and editing of the manuscript.

The Supporting Information is available free of charge on the ACS Publications website at DOI: XXXX

Supporting Information contains details of polyprotein engineering, single-molecule magnetic tweezers instrumentation and set-up, calculations of the magnetic bead potential, derivations of compliance of chain models under constant force, analytical expressions for the freely-jointed chain model, the data processing workflow, and additional data sets.

References

    Rief M, Gautel M, Oesterhelt F, Fernandez JM, Gaub HE. Reversible unfolding of individual titin immunoglobulin domains by AFM. Science1997. 276: 53151109-12

    Elosegui-Artola A, Oria R, Chen Y, Kosmalska A, Perez-Gonzalez C, Castro N, Zhu C, Trepat X, Roca-Cusachs P. Mechanical regulation of a molecular clutch defines force transmission and transduction in response to matrix rigidity. Nat Cell Biol2016. 18: 5540-8

    Yao M, Goult BT, Klapholz B, Hu X, Toseland CP, Guo Y, Cong P, Sheetz MP, Yan J. The mechanical response of talin. Nat Commun2016. 7: 11966

    Infante E, Stannard A, Board SJ, Rico-Lastres P, Rostkova E, Beedle AEM, Lezamiz A, Wang YJ, Gulaidi Breen S, Panagaki F, Sundar Rajan V, . The mechanical stability of proteins regulates their translocation rate into the cell nucleus. Nature Physics2019. 15: 9973-981

    Vogel V. Mechanotransduction involving multimodular proteins: converting force into biochemical signals. Annu Rev Biophys Biomol Struct2006. 35: 459-88

    Cao Y, Li H. Engineered elastomeric proteins with dual elasticity can be controlled by a molecular regulator. Nat Nanotechnol2008. 3: 8512-6

    Mora M, Stannard A, Garcia-Manyes S. The nanomechanics of individual proteins. Chem Soc Rev2020. 49: 6816-6832

    Garcia-Manyes S, Beedle AEM. Steering chemical reactions with force. Nature Reviews Chemistry2017. 1: 110083

    Schonfelder J, Alonso-Caballero A, De Sancho D, Perez-Jimenez R. The life of proteins under mechanical force. Chem Soc Rev2018. 47: 103558-3573

10 

    Puchner EM, Gaub HE. Force and function: probing proteins with AFM-based force s pectroscopy. Curr Opin Struct Biol2009. 19: 5605-14

11 

    Hoffmann T, Dougan L. Single molecule force spectroscopy using polyproteins. Chem Soc Rev2012. 41: 144781-96

12 

    Cossio P, Hummer G, Szabo A. Transition paths in single-molecule force spectroscopy. J Chem Phys2018. 148: 12123309

13 

    Smith SB, Cui Y, Bustamante C. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science1996. 271: 5250795-9

14 

    Berkovich R, Garcia-Manyes S, Urbakh M, Klafter J, Fernandez JM. Collapse dynamics of single proteins extended by force. Biophys J2010. 98: 112692-701

15 

    Lof A, Walker PU, Sedlak SM, Gruber S, Obser T, Brehm MA, Benoit M, Lipfert J. Multiplexed protein force spectroscopy reveals equilibrium protein folding dynamics and the low-force response of von Willebrand factor. Proc Natl Acad Sci U S A2019. 116: 3818798-18807

16 

    Popa I, Rivas-Pardo JA, Eckels EC, Echelman DJ, Badilla CL, Valle-Orero J, Fernandez JM. A HaloTag Anchored Ruler for Week-Long Studies of Protein Dynamics. J Am Chem Soc2016. 138: 3310546-53

17 

    Yao M, Qiu W, Liu R, Efremov AK, Cong P, Seddiki R, Payre M, Lim CT, Ladoux B, Mege RM, Yan J. Force-dependent conformational switch of alpha-catenin controls vinculin binding. Nat Commun2014. 5: 4525

18 

    Guo Z, Hong H, Yuan G, Qian H, Li B, Cao Y, Wang W, Wu CX, Chen H. Hidden Intermediate State and Second Pathway Determining Folding and Unfolding Dynamics of GB1 Protein at Low Forces. Phys Rev Lett2020. 125: 19198101

19 

    Rico-Pasto M, Pastor I, Ritort F. Force feedback effects on single molecule hopping and pulling experiments. J Chem Phys2018. 148: 12123327

20 

    Hinczewski M, Gebhardt JC, Rief M, Thirumalai D. From mechanical folding trajectories to intrinsic energy landscapes of biopolymers. Proc Natl Acad Sci U S A2013. 110: 124500-5

21 

    Berkovich R, Hermans RI, Popa I, Stirnemann G, Garcia-Manyes S, Berne BJ, Fernandez JM. Rate limit of protein elastic response is tether dependent. Proc Natl Acad Sci U S A2012. 109: 3614416-21

22 

    De Sancho D, Schonfelder J, Best RB, Perez-Jimenez R, Munoz V. Instrumental Effects in the Dynamics of an Ultrafast Folding Protein under Mechanical Force. J Phys Chem B2018. 122: 4911147-11154

23 

    Cossio P, Hummer G, Szabo A. On artifacts in single-molecule force spectroscopy. Proc Natl Acad Sci U S A2015. 112: 4614248-53

24 

    Walther KA, Brujic J, Li H, Fernandez JM. Sub-angstrom conformational changes of a single molecule captured by AFM variance analysis. Biophys J2006. 90: 103806-12

25 

    Garcia-Manyes S, Brujic J, Badilla CL, Fernandez JM. Force-clamp spectroscopy of single-protein monomers reveals the individual unfolding and folding pathways of I27 and ubiquitin. Biophys J2007. 93: 72436-46

26 

    Lei H, He C, Hu C, Li J, Hu X, Hu X, Li H. Single-Molecule Force Spectroscopy Trajectories of a Single Protein and Its Polyproteins Are Equivalent: A Direct Experimental Validation Based on A Small Protein NuG2. Angew Chem Int Ed Engl2017. 56: 226117-6121

27 

    Haining AW, von Essen M, Attwood SJ, Hytonen VP, Del Rio Hernandez A. All Subdomains of the Talin Rod Are Mechanically Vulnerable and May Contribute To Cellular Mechanosensing. ACS Nano2016. 10: 76648-58

28 

    Tapia-Rojo R, Alonso-Caballero A, Fernandez JM. Direct observation of a coil-to-helix contraction triggered by vinculin binding to talin. Sci Adv2020. 6: 21eaaz4707

29 

    Brockwell DJ, Beddard GS, Paci E, West DK, Olmsted PD, Smith DA, Radford SE. Mechanically unfolding the small, topologically simple protein L. Biophys J2005. 89: 1506-19

30 

    Tapia-Rojo R, Alonso-Caballero A, Fernandez JM. Talin folding as the tuning fork of cellular mechanotransduction. Proc Natl Acad Sci U S A2020. 117: 3521346-21353

31 

    Rajgor D, Shanahan CM. Nesprins: from the nuclear envelope and beyond. Expert Rev Mol Med2013. 15: e5

32 

    Arsenovic PT, Ramachandran I, Bathula K, Zhu R, Narang JD, Noll NA, Lemmon CA, Gundersen GG, Conway DE. Nesprin-2G, a Component of the Nuclear LINC Complex, Is Subject to Myosin-Dependent Tension. Biophys J2016. 110: 134-43

33 

    Rivas-Pardo JA, Eckels EC, Popa I, Kosuri P, Linke WA, Fernandez JM. Work Done by Titin Protein Folding Assists Muscle Contraction. Cell Rep2016. 14: 61339-1347

34 

    Ishima R, Torchia DA. Protein dynamics from NMR. Nat Struct Biol2000. 7: 9740-3

35 

    Fuertes G, Banterle N, Ruff KM, Chowdhury A, Mercadante D, Koehler C, Kachala M, Estrada Girona G, Milles S, Mishra A, Onck PR, . Decoupling of size and shape fluctuations in heteropolymeric sequences reconciles discrepancies in SAXS vs. FRET measurements. Proc Natl Acad Sci U S A2017. 114: 31E6342-E6351

36 

    Schuler B, Soranno A, Hofmann H, Nettels D. Single-Molecule FRET Spectroscopy and the Polymer Physics of Unfolded and Intrinsically Disordered Proteins. Annu Rev Biophys2016. 45: 207-31

37 

    Andrecka J, Ortega Arroyo J, Takagi Y, de Wit G, Fineberg A, MacKinnon L, Young G, Sellers JR, Kukura P. Structural dynamics of myosin 5 during processive motion revealed by interferometric scattering microscopy. Elife2015. 4: e05413

38 

    Walther KA, Grater F, Dougan L, Badilla CL, Berne BJ, Fernandez JM. Signatures of hydrophobic collapse in extended proteins captured with force spectroscopy. Proc Natl Acad Sci U S A2007. 104: 197916-21

39 

    Cecconi C, Shank EA, Bustamante C, Marqusee S. Direct observation of the three-state folding of a single protein molecule. Science2005. 309: 57432057-60

40 

    Shank EA, Cecconi C, Dill JW, Marqusee S, Bustamante C. The folding cooperativity of a protein is controlled by its chain topology. Nature2010. 465: 7298637-40

41 

    Stigler J, Ziegler F, Gieseke A, Gebhardt JC, Rief M. The complex folding network of single calmodulin molecules. Science2011. 334: 6055512-6

42 

    Fernandez JM, Li H. Force-clamp spectroscopy monitors the folding trajectory of a single protein. Science2004. 303: 56641674-8

43 

    Lannon H, Haghpanah JS, Montclare JK, Vanden-Eijnden E, Brujic J. Force-clamp experiments reveal the free-energy profile and diffusion coefficient of the collapse of protein molecules. Phys Rev Lett2013. 110: 12128301

44 

    Garcia-Manyes S, Dougan L, Badilla CL, Brujic J, Fernandez JM. Direct observation of an ensemble of stable collapsed states in the mechanical folding of ubiquitin. Proc Natl Acad Sci U S A2009. 106: 2610534-9

45 

    Brujić J, Hermans Z RI, Walther KA, Fernandez JM. Single-molecule force spectroscopy reveals signatures of glassy dynamics in the energy landscape of ubiquitin. Nature Physics2006. 2: 4282-286

46 

    Gupta AN, Neupane K, Rezajooei N, Cortez LM, Sim VL, Woodside MT. Pharmacological chaperone reshapes the energy landscape for folding and aggregation of the prion protein. Nat Commun2016. 7: 12058

47 

    Beedle AE, Lynham S, Garcia-Manyes S. Protein S-sulfenylation is a fleeting molecular switch that regulates non-enzymatic oxidative folding. Nat Commun2016. 7: 12490

48 

    Valle-Orero J, Tapia-Rojo R, Eckels EC, Rivas-Pardo JA, Popa I, Fernandez JM. Proteins Breaking Bad: A Free Energy Perspective. J Phys Chem Lett2017. 8: 153642-3647

49 

    Loescher CM, Breitkreuz M, Li Y, Nickel A, Unger A, Dietl A, Schmidt A, Mohamed BA, Kotter S, Schmitt JP, Kruger M, . Regulation of titin-based cardiac stiffness by unfolded domain oxidation (UnDOx). Proc Natl Acad Sci U S A2020. 117: 3924545-24556

50 

    Perales-Calvo J, Giganti D, Stirnemann G, Garcia-Manyes S. The force-dependent mechanism of DnaK-mediated mechanical folding. Sci Adv2018. 4: 2eaaq0243

51 

    Wiita AP, Perez-Jimenez R, Walther KA, Grater F, Berne BJ, Holmgren A, Sanchez-Ruiz JM, Fernandez JM. Probing the chemistry of thioredoxin catalysis with force. Nature2007. 450: 7166124-7

52 

    del Rio A, Perez-Jimenez R, Liu R, Roca-Cusachs P, Fernandez JM, Sheetz MP. Stretching single talin rod molecules activates vinculin binding. Science2009. 323: 5914638-41

53 

    Stirnemann G, Giganti D, Fernandez JM, Berne BJ. Elasticity, structure, and relaxation of extended proteins under force. Proc Natl Acad Sci U S A2013. 110: 103847-52

54 

    Neupane K, Foster DA, Dee DR, Yu H, Wang F, Woodside MT. Direct observation of transition paths during the folding of proteins and nucleic acids. Science2016. 352: 6282239-42

55 

    Neupane K, Manuel AP, Woodside MT. Protein folding trajectories can be described quantitatively by one-dimensional diffusion over measured energy landscapes. Nature Physics2016. 12: 7700-703

Floating objects

The effect of force on the free energy landscape of protein (un)folding.
FIG. 1

The effect of force on the free energy landscape of protein (un)folding.

Using the example of protein L, in the absence of force (grey line), there is one minimum, corresponding to the folded state (blue). When a low force is applied, a broad local minimum corresponding to an unfolded state emerges (yellow). At a certain, coexistence force (red), the protein will equally likely be found in folded and unfolded states. As force increases, the unfolded minimum becomes narrower, eventually becoming the global minimum (green).

A step-wise change in protein compliance fingerprints mechanical unfolding/refolding of a single PL monomer under force.
FIG. 2

A step-wise change in protein compliance fingerprints mechanical unfolding/refolding of a single PL monomer under force.

(A) Schematic illustration of PL (PDB: 1HZ6) construct being stretched in a magnetic tweezers set-up. (B) Raw extension-time measurement of a protein L monomer hopping in equilibrium between the folded and unfolded states. The (Ig32)2-PL-(Ig32)2 construct is initially pulled at 38 pN, which extends and unfolds the PL monomer. After 60 s, the force is quenched down to 8.1 pN for 60 minutes. During this time the PL monomer stochastically hops between the folded and unfolded states in steps of Δ<x> = 10.3 ± 0.3 nm. A final pulse at a high force (38 pN) value triggers the (re)unfolding of the PL monomer. (C) Result of flattening the protein extension data during the 8.1pN pulse to remove slow drift effects, clearly highlighting 13 individual (un)folding transitions. (D) Compliance of (C) (calculated over a 3-second moving window), showing step-wise changes (Δc = 0.86 ± 0.09 nm/pN) concomitant to the PL individual unfolding and refolding events.

The change in compliance in a PL8 polyprotein is additive and recapitulates that of an individual domain.
FIG. 3

The change in compliance in a PL8 polyprotein is additive and recapitulates that of an individual domain.

(A) Illustration of the PL8 construct being pulled in the MT set-up. (B) Typical (un)folding hopping trajectory of PL8 under force. Following a force-quench protocol, a high force pulse (49 pN) is first applied, to trigger the unfolding of each individual domain within the polyprotein chain (red asterisks), hallmarked by changes in extension of 17.2 nm. The force was subsequently dropped down to 8.1 pN (resulting in an initial decrease in extension) and held for 50 minutes. During this period of time, the PL8 polyprotein dynamically hops between well-defined states separated by Δ<x> = 10.7 ± 0.1 nm. (C) Segment flattening and colour-based assignment of the number of unfolded domains allows for clearer visualisation of the step-wise hopping dynamics. (D) The average compliance (n = 131 (un)folding events) can be calculated for a given number of unfolded domains. Linear fit (weighted by the total observation time at each level) to the data yields a gradient of 0.81 ± 0.15 (R 2 = 0.983), corresponding to the change in compliance of an individual domain.

Correlating force-dependent extension and compliance changes upon mechanical (un)folding through the FJC model of polymer elasticity.
FIG. 4

Correlating force-dependent extension and compliance changes upon mechanical (un)folding through the FJC model of polymer elasticity.

(A) Cyclic unfolding and refolding trajectories allow probing the change in protein stiffness upon PL unfolding at varying high forces. Using a force quench protocol, the protein is initially held at a seemingly low force (1.8 pN) for 30 seconds to ensure correct refolding. The force is then increased to a higher constant force value in each cycle (namely 21.9, 16.7, 12.7, and 9.7 pN) for 5 minutes, which triggers in each case the unfolding of the PL monomer (marked by a step-wise increase in length). (B) Similar cyclic force-quench protocol to capture the change in protein compliance upon mechanical refolding at low forces. Before each cycle, the force is first kept low (1.8 pN) for 30 seconds to reach the protein’s native state. Subsequently, in each cycle the protein is unfolded at high force (49.2 pN) for 40 seconds, before the force is successively quenched at low forces, namely 6.2 pN (5 min), 5.7 pN (4 min), 5.2 pN (3 min), 4.7 pN (2 min) and 4.3 pN (1 min), respectively. When held at such low force values, the protein stochastically refolds back to the native state, as evidenced by the step-wise reduction in length. While the folded state dominates at these forces, some hopping back to an unfolded state is observed. (C) Representative unfolding (left) and folding (right) events occurring at high (21.9 pN) and low (5.7 pN) forces respectively (corresponding to the green dashed box in (A) and the yellow dashed box in (B), respectively). High force unfolding is hallmarked by a large increase in average extension but a small increase in compliance, conversely low force folding is hallmarked by a small decrease in average extension but a large decrease in compliance. (D) Plot of (absolute) compliance change against (absolute) extension change combining low (yellow), hopping (red), and high (green) force experiments, clearly showing that, upon PL un/folding, an inverse correlation between its changes in extension and compliance is measured. Data points for monomers (circles) correspond to the mean ± s.e. of 9–17 events occurring at each force, data for the PL8 octamer (square) is taken from Fig. 3D. Combined, this data show very good agreement with the FJC model (grey dashed line) with Kuhn length b = 1.1 nm, k B T = 4.04 pN nm, and ΔL = 18.6 nm.

Direct observation of the change in end-to-end fluctuations upon individual talin and nesprin unfolding/refolding events.
FIG. 5

Direct observation of the change in end-to-end fluctuations upon individual talin and nesprin unfolding/refolding events.

(A) Illustration of the [(Spy0128)2-TalinR3IVVI] construct being pulled in the MT set-up (PDB: 2L7A). (B) (Left) Raw, unfiltered extension-time data corresponding to the stretching of individual talin R3 IVVI monomer at a constant force of 8.35 pN for 200 seconds. During this time, the protein monomer undergoes multiple unfolding and refolding transitions. (Right) Zoom in on the hopping section, showing 84 consecutive unfolding(red)/refolding(blue) transitions concomitant to an average change in extension of Δ<x> = 22.2 ± 0.2 nm, and an average change in compliance between the folded and unfolded states of Δc = 1.24 ± 0.14 nm/pN. Averaging over 5 such cycles gives Δ<x> = 22.2 ± 0.1 nm and Δc = 1.28 ± 0.08 nm/pN. (C) Illustration of the [Ig27-SR73]4 construct being pulled in the MT set-up (lacking crystal structure, here we display the similar SR16 (PDB: 1U4Q). (D) Typical individual (unfiltered) folding trajectory under a force quench. The force is first increased from 10.6 pN to 37.5 pN over the course of 30 s (corresponding to the linear movement of magnet position by 1.4 mm over this time), resulting in the sequential step-wise unfolding of the four SR73 domains of nesprin (green asterisks). Note that, within this force range, the mechanically rigid Ig27 domains in the construct do not unfold. The force is then quenched down to 7.4 pN for 30 seconds, triggering the collapse and refolding of nesprin. At such a low (<F 0.5) force value, refolding is relatively slow, such that the staircase-like pattern that mirrors that of unfolding can be captured. The force is subsequently ramped up again (test pulse) back to 37.5 pN, to probe that all the individual SR73 domains (green asterisks) had effectively refolded during the force-quench. (E) Zoom (yellow square) on the initial folding section of (B), where the four SR73 domains of nesprin refold sequentially, marked by a step-wise reduction of the protein extension by Δ<x> = 17.1 ± 0.5 nm. The refolding of each individual domain occurs concomitant to a reduction in the end-to-end fluctuations of the protein construct. (F) Plotting the compliance against the number of unfolded domains measured in (C) allows calculation of the change in compliance per domain, Δc/domain = 2.30 ± 0.59 nm/pN. Averaging over 10 such folding trajectories with four SR73 nesprin unfolded domains gives Δ<x> = 16.8 ± 0.2 nm and Δc/domain = 1.99 ± 0.12 nm/pN. (G) Plotting normalised compliance change (Δc/β bΔL) against normalised extension change (ΔxL) crucially reveals that the observed behaviour is protein independent, showing quantitative agreement with a master curve of the FJC model.