Protein Folding Under Mechanical Forces: A Physiological View

On the free energy of protein folding in optical tweezers experiments

Free energy is a critical parameter in understanding the equilibrium in chemical reactions. It enables us to determine the equilibrium proportion between the different species in the reaction and to predict in which direction the reaction will proceed if a change is performed in the system. Historically, to calculate this value, bulk experiments were performed where a parameter was altered at a gradual rate to change the population until a new equilibrium was established. In protein folding studies, it is common to vary the temperature or chaotropic agents in order to change the population and then to extrapolate to physiological conditions. Such experiments were time-consuming due to the necessity of ensuring equilibrium and reversibility. Techniques of single-molecule manipulation, such as optical/magnetic tweezers and atomic force microscopy, permit the direct measurement of the work performed by a protein undergoing unfolding/refolding at particular forces. Also, with the development of non-equilibrium free energy theorems (Jarzynski equality, Crooks fluctuation theorem, Bennett acceptance ratio, and overlapping method), it is possible to obtain free energy values in experiments far from equilibrium. This review compares different methodologies and their application in optical tweezers. Interestingly, in many proteins, discrepancies in free energy values obtained through different methods suggest additional complexities in the folding pathway, possibly involving intermediate states such as the molten globule. Further studies are needed to confirm their presence and significance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-025-01310-0.

Structurally different chemical chaperones show similar mechanical roles with independent molecular mechanisms

Osmolytes are well known to protect the protein structure against different chemical and physical denaturants. Since their actions with protein surfaces are mechanistically complicated and context dependent, the underlying molecular mechanism is not fully understood. Here, we combined single-molecule magnetic tweezers and molecular dynamics (MD) simulation to explore the mechanical role of osmolytes from two different classes, trimethylamine N-oxide (TMAO) and trehalose, as mechanical stabilizers of protein structure. We observed that these osmolytes increase the protein L mechanical stability by decreasing unfolding kinetics while accelerating the refolding kinetics under force, eventually shifting the energy landscape toward the folded state. These osmolytes mechanically stabilize the protein L and plausibly guide them to more thermodynamically robust states. Finally, we observed that osmolyte-modulated protein folding increases mechanical work output up to twofold, allowing the protein to fold under a higher force regime and providing a significant implication for folding-induced structural stability in proteins.

Unfolding and Refolding Proteins Using Single-Molecule AFM

Single-molecule atomic force microscopy (AFM) allows capturing the conformational dynamics of an individual molecule under force. In this chapter, we describe a protocol for conducting a protein nanomechanical experiment using AFM, covering both the force-extension and force-clamp modes. Combined, these experiments provide an integrated vista of the molecular mechanisms-and their associated kinetics-underpinning the mechanical unfolding and refolding of individual proteins when exposed to mechanical load.

Multifaceted folding-unfolding landscape of the TrpZip2 β-hairpin and the role of external sub-piconewton mechanical tensions

Proteins can experience uneven tensions of the order of tens of piconewtons when exposed to different solvent environment due to the thermal motion of the solvent. It is also true that biomolecules, especially proteins, are subjected to a variety of mechanical tensions generated by several factors, including mechanically assisted translocation and pressure gradients within living systems. Here, we use metadynamics simulations to revisit the folding-unfolding of the TrpZip2 β-hairpin and redefine it from the perspective of an external force of a sub-piconewton magnitude acting on the ends of the hairpin. The chosen forces, while preserving the morphology of the β-hairpin chain when it is pulled, are capable of influencing the conformational behavior of the chain during folding and unfolding. Our investigations confirm that the TrpZip2 β-hairpin exhibits a zipper (zip-out) mechanism for folding-unfolding in both mechanically unbiased and biased (with a 30 pN end force) situations. However, it is important to note that they present marked differences in their folding and unfolding paths, with the mechanically biased system capable of becoming trapped in various intermediate states. Both unbiased and biased scenarios of the hairpin indicate that the hairpin turn is highly stable during the folding-unfolding event and initiates folding. More importantly we confirm that the existing heterogeneity in the TrpZip2 β-hairpin folding-unfolding is a consequence of the wide range of conformations observed, owing to the different trapped intermediates caused by the uneven forces it may experience in solution.

Single Molecule Force Spectroscopy Studies on Metalloproteins: Opportunities and Challenges

Metalloproteins play important roles in a wide range of biological processes. Elucidating the mechanisms via which metalloproteins fold and constitute their metal centers is critical to the understanding of the functions and dynamics of metalloproteins. Owing to its superior force and length resolution, single-molecule force spectroscopy (SMFS) has evolved into a powerful tool to probe the unfolding and folding mechanisms of metalloproteins at the single level by forcing metalloproteins to unfold and then refold along a reaction coordinate defined by the applied stretching force. The folding of metalloproteins is complex and involves two interwound processes, the folding of the polypeptide chain and the constitution of the metal center. Experimental studies of the folding of metalloproteins are challenging. SMFS studies have allowed researchers to directly probe the folding and unfolding of metalloproteins at the single-molecule level and the effect of metal centers on the folding-unfolding energy landscape of metalloproteins. New mechanistic insights on the folding and unfolding of some metalloproteins have been obtained, demonstrating the power and unique advantages that SMFS techniques may offer. In this Perspective, using calcium-binding proteins and small iron-sulfur proteins as examples, I provide a concise overview of the information and insights that SMFS studies have provided to understand the folding and unfolding of metalloproteins. I also discuss the opportunities and challenges that are present in this fast-progressing area of research.

Measuring (biological) materials mechanics with atomic force microscopy. 3. Mechanical unfolding of biopolymers

Biopolymers, such as polynucleotides, polypeptides and polysaccharides, are macromolecules that direct most of the functions in living beings. Studying the mechanical unfolding of biopolymers provides important information about their molecular elasticity and mechanical stability, as well as their energy landscape, which is especially important in proteins, since their three‐dimensional structure is essential for their correct activity. In this primer, we present how to study the mechanical properties of proteins with atomic force microscopy and how to obtain information about their stability and energetic landscape. In particular, we discuss the preparation of polyprotein constructs suitable for AFM single molecule force spectroscopy (SMFS), describe the parameters used in our force‐extension SMFS experiments and the models and equations employed in the analysis of the data. As a practical example, we show the effect of the temperature on the unfolding force, the distance to the transition state, the unfolding rate at zero force, the height of the transition state barrier, and the spring constant of the protein for a construct containing nine repeats of the I27 domain from the muscle protein titin.

Modeling multiple duplex DNA attachments in a force-extension experiment

Optical tweezers-based DNA stretching often relies on tethering a single end-activated DNA molecule between optically manipulated end-binding beads. Measurement success can depend on DNA concentration. At lower DNA concentrations tethering is less common, and many trials may be required to observe a single-molecule stretch. At higher DNA concentrations tethering is more common; however, the resulting force-extensions observed are more complex and may vary from measurement to measurement. Typically these more complex results are attributed to the formation of multiple tethers between the beads; however, to date there does not appear to have been a critical examination of this hypothesis or the potential usefulness of such data. Here we examine stretches at a higher DNA concentration and use analysis and simulation to show how the more complex force-extensions observed can be understood in terms of multiple DNA attachments.

New insights into the folding–unfolding mechanism and conformations of cytochrome C

Optical trapping experiments offer new insights into the folding and unfolding of cytochrome C.

Out-of-Equilibrium Biophysical Chemistry: The Case for Multidimensional, Integrated Single-Molecule Approaches

Out-of-equilibrium processes are ubiquitous across living organisms and all structural hierarchies of life. At the molecular scale, out-of-equilibrium processes (for example, enzyme catalysis, gene regulation, and motor protein functions) cause biological macromolecules to sample an ensemble of conformations over a wide range of time scales. Quantifying and conceptualizing the structure–dynamics to function relationship is challenging because continuously evolving multidimensional energy landscapes are necessary to describe nonequilibrium biological processes in biological macromolecules. In this perspective, we explore the challenges associated with state-of-the-art experimental techniques to understanding biological macromolecular function. We argue that it is time to revisit how we probe and model functional out-of-equilibrium biomolecular dynamics. We suggest that developing integrated single-molecule multiparametric force–fluorescence instruments and using advanced molecular dynamics simulations to study out-of-equilibrium biomolecules will provide a path towards understanding the principles of and mechanisms behind the structure–dynamics to function paradigm in biological macromolecules.

Single molecule force spectroscopy reveals the context dependent folding pathway of the C-terminal fragment of Top7

Top7 is a de novo designed protein with atomic level accuracy and shows a folded structure not found in nature. Previous studies showed that the folding of Top7 is not cooperative and involves various folding intermediate states. In addition, various fragments of Top7 were found to fold on their own in isolation. These features displayed by Top7 are distinct from those of naturally occurring proteins of a similar size and suggest a rough folding energy landscape. However, it remains unknown if and how the intra-polypeptide chain interactions among the neighboring sequences of Top7 affect the folding of these Top7 fragments. Here we used single-molecule optical tweezers to investigate the folding–unfolding pathways of full length Top7 as well as its C-terminal fragment (CFr) in different sequence environments. Our results showed that the mechanical folding of Top7 involves an intermediate state that likely involves non-native interactions/structure. More importantly, we found that the folding of CFr is entirely dependent upon its sequence context in which it is located. When in isolation, CFr indeed folds into a cooperative structure showing near-equilibrium unfolding–folding transitions at ∼6.5 pN in OT experiments. However, CFr loses its autonomous cooperative folding ability and displays a folding pathway that is dependent on its interactions with its neighboring sequence/structure. This context-dependent folding dynamics and pathway of CFr are distinct from those of naturally occurring proteins and highlight the critical importance of intra-chain interactions in shaping the overall energy landscape and the folding pathway of Top7. These new insights may have important implications on the de novo design of proteins.

Optical tweezers experiments reveal that the folding of the C-terminal fragment of Top7 (cFr) is context-dependent. Depending on its neighboring sequence, cFr shows very different folding pathways and folding kinetics.

Single molecule force spectroscopy reveals that a two-coordinate ferric site is critical for the folding of holo-rubredoxin

Metalloproteins play important roles in a wide range of biological processes. The folding process of metalloproteins is complex due to the synergistic effects of the folding of their polypeptide chains and the incorporation of metal cofactors. The folding mechanism of the simplest iron-sulfur protein rubredoxin, which contains one ferric ion coordinated by four cysteinyl sulfurs, is revealed using optical tweezers for the first time. The folding of the rubredoxin polypeptide chain is rapid and robust, while the reconstitution of the iron-sulfur center is greatly dependent upon the coordination state of the ferric ion on the unfolded polypeptide chain. If the ferric ion is coordinated by two neighboring cysteines, rubredoxin can readily fold with the iron-sulfur center fully reconstituted. However, if the ferric ion is only mono-coordinated, rubredoxin can fold but the iron-sulfur center is not reconstituted. Our results suggested that the folding of holo-rubredoxin follows a novel binding-folding-reconstitution mechanism, which is distinct from the folding mechanisms proposed for the folding of metalloproteins. Our study highlights the critical importance of the two-coordinate ferric site in the folding of holo-rubredoxin, which may have some important implications to our understanding of the folding mechanism of more complex metalloproteins in vivo.

The nanomechanics of individual proteins

This tutorial review provides an overview of the single protein force spectroscopy field, including the main techniques and the basic tools for analysing the data obtained from the single molecule experiments.

Dihedral Angle Calculations To Elucidate the Folding of Peptides through Its Main Mechanical Forces

This study reports a general method to calculate dihedral angles (φ and ψ) of a given amino acid sequence, focusing on potential energy and torque moment concepts. By defining these physical measures in relation to the chemical interactions that occur on each single amino acid residue within a peptide, we analyze the folding process as the result of main mechanical forces (MMFs) exerted in the specific amino acid chain of interest. As a proof of concept, Leu-enkephalin was initially used as a model peptide to carry out the theoretical study. Our data show agreement between calculated Leu-enkephalin backbone dihedral angles and the corresponding experimentally determined X-ray values. Hence, we used calcitonin to validate our MMF-based method on a larger peptide, i.e., 32 amino acid residues forming an α-helix. Through a similar approach (although simplified with regard to electrostatic interactions), the calculations for calcitonin also demonstrate a good agreement with experimental values. This study offers new opportunities to analyze peptides' amino acid sequences and to help in the prediction of how they must fold, assisting in the development of new computational techniques in the field.

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.

Mechanical Unfolding Pathway of the High-Potential Iron-Sulfur Protein Revealed by Single-Molecule Atomic Force Microscopy: Toward a General Unfolding Mechanism for Iron-sulfur Proteins

High-potential iron-sulfur proteins (HiPIPs) are an important class of metalloproteins with a [4Fe-4S] cluster coordinated by four cysteine residues. Distinct from other iron-sulfur proteins, the cluster in HiPIP has a high reduction potential, making it an essential electron carrier in bacterial photosynthesis. Here, we combined single-molecule atomic force microscopy and protein engineering techniques to investigate the mechanical unfolding mechanism of HiPIP from Chromatium tepidum (cHiPIP). We found that cHiPIP unfolds in a two-step fashion with the protein sequence sequestered by the iron-sulfur center as a stable unfolding intermediate state. The rupture of the iron-sulfur center of cHiPIP proceeds in two distinct parallel pathways; one pathway involves the concurrent rupture of multiple iron-thiolate bonds, and the other one involves the sequential rupture of the iron-thiolate bonds. This mechanistic information was further confirmed by mutational studies. We found that the rupture of the iron-thiolate bonds in reduced and oxidized cHiPIP occurred in the range of 150-180 pN at a pulling speed of 400 nm/s, similar to that measured for iron-thiolate bonds in rubredoxin and ferredoxin. Our results may have important implications for understanding the general unfolding mechanism governing iron-sulfur proteins, as well as the mechanism governing the mechanical rupture of the iron-sulfur center.

The life of proteins under mechanical force

Although much of our understanding of protein folding comes from studies of isolated protein domains in bulk, in the cellular environment the intervention of external molecular machines is essential during the protein life cycle. During the past decade single molecule force spectroscopy techniques have been extremely useful to deepen our understanding of these interventional molecular processes, as they allow for monitoring and manipulating mechanochemical events in individual protein molecules. Here, we review some of the critical steps in the protein life cycle, starting with the biosynthesis of the nascent polypeptide chain in the ribosome, continuing with the folding supported by chaperones and the translocation into different cell compartments, and ending with proteolysis in the proteasome. Along these steps, proteins experience molecular forces often combined with chemical transformations, affecting their folding and structure, which are measured or mimicked in the laboratory by the application of force with a single molecule apparatus. These mechanochemical reactions can potentially be used as targets for fighting against diseases. Inspired by these insightful experiments, we devise an outlook on the emerging field of mechanopharmacology, which reflects an alternative paradigm for drug design.

Differences in the mechanical unfolding pathways of apo- and copper-bound azurins

Metalloproteins carry out diverse biological functions including metal transport, electron transfer, and catalysis. At present, the influence of metal cofactors on metalloprotein stability is not well understood. Here, we report the mechanical stability and unfolding pathway of azurin, a cupredoxin family protein with β-barrel topology and type I copper-binding centre. Single-molecule force spectroscopy (SMFS) experiments reveal 2-state and 3-state unfolding pathways for apo-azurin. The intermediate in the 3-state pathway occurs at an unfolding contour length of 7.5 nm from the native state. Steered molecular dynamics (SMD) simulations show that apo-azurin unfolds via a first transition state (TS) where β2Β–β8 and β7–β8 strand pairs rupture to form the intermediate, which subsequently unfolds by the collective rupture of remaining strands. SMFS experiments on holo-azurin exhibit an additional 4-state pathway besides the 2-state and 3-state pathways. The unfolding contour length leading to the first intermediate is 6.7 nm suggesting a sequestration of ~1 nm polypeptide chain length by the copper. SMD simulations reveal atomistic details of the copper sequestration and predict a combined β4–β7 pair and copper coordination sphere rupture to create the third TS in the 4-state pathway. Our systematic studies provide detailed mechanistic insights on modulation of protein mechanical properties by metal-cofactors.

AFM-Based Single-Molecule Force Spectroscopy of Proteins

Single-molecule force spectroscopy by AFM (AFM-SMFS) is an experimental methodology that allows unequivocal sensitivity and control for investigating and manipulating the mechanical properties of single molecules. The past 20 years of AFM-SMFS has provided numerous breakthroughs in the understanding of the mechanical properties and force-induced structural rearrangements of sugars, DNA, and proteins. Here, we focus on the application of AFM-SMFS to study proteins, since AFM-SMFS has succeeded in providing abundant information about protein folding pathways, kinetics, interactions, and misfolding. In this chapter we describe the experimental procedures for conducting a SMFS-AFM experiment-including purification of protein samples, setup and calibration of the AFM instrumentation, and the thorough and unbiased analysis of resulting AFM data.

The Power of Force: Insights into the Protein Folding Process Using Single-Molecule Force Spectroscopy

One of the major challenges in modern biophysics is observing and understanding conformational changes during complex molecular processes, from the fundamental protein folding to the function of molecular machines. Single-molecule techniques have been one of the major driving forces of the huge progress attained in the last few years. Recent advances in resolution of the experimental setups, aided by theoretical developments and molecular dynamics simulations, have revealed a much higher degree of complexity inside these molecular processes than previously reported using traditional ensemble measurements. This review sums up the evolution of these developments and gives an outlook on prospective discoveries.

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

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

Steered molecular dynamics study reveals insights into the function of the repetitive B region of collagen- and fibrinogen-binding MSCRAMMs

MSCRAMMs (microbial surface components recognizing adhesive matrix molecules) are modular proteins covalently anchored in the bacterial cell wall of many Gram-positive bacteria. The N-terminal region of most MSCRAMMs carries the ligand-binding domains (A region) which specifically target the host extracellular matrix (ECM) proteins such as collagen, fibrinogen and fibronectin. In Staphylococcus aureus Cna, the prototype collagen-binding MSCRAMM, the A region is followed by a repetitive B region which is found to be conserved among many Gram-positive bacteria. This conservation signifies an important functional role for the B region which is made of repetitive domains. It was suggested that this region could act as a 'stalk' as well as a 'spring' to present the ligand-binding A region, away from the bacterial surface. But there is no clear functional implication of this region available till date. Each repetitive domain in the B region possesses a variant of the Ig fold called the CnaB fold. Additionally, the B repeats are also paired and the pairs are clustered together. To investigate if the B domains have a function similar to the Ig domains in the I-band region of the giant muscle protein, titin, steered molecular dynamics simulations of one, two and four B repeats of Cna were carried out. The results of the simulations suggest that the B region could provide mechanical stability, extensibility and elasticity to Cna due to the CnaB fold as well as the clustered arrangement of their domains. This study thus provided further insights into the biological underpinnings of adhesin-host interaction.

Quantitative interpretation of FRET experiments via molecular simulation: force field and validation

Molecular simulation is a valuable and complementary tool that may assist with the interpretation of single-molecule Förster resonance energy transfer (FRET) experiments, if the energy function is of sufficiently high quality. Here we present force-field parameters for one of the most common pairs of chromophores used in experiments, AlexaFluor 488 and 594. From microsecond molecular-dynamics simulations, we are able to recover both experimentally determined equilibrium constants and association/dissociation rates of the chromophores with free tryptophan, as well as the decay of fluorescence anisotropy of a labeled protein. We find that it is particularly important to obtain a correct balance of solute-water interactions in the simulations in order to faithfully capture the experimental anisotropy decays, which provide a sensitive benchmark for fluorophore mobility. Lastly, by a combination of experiment and simulation, we address a potential complication in the interpretation of experiments on polyproline, used as a molecular ruler for FRET experiments, namely the potential association of one of the chromophores with the polyproline helix. Under conditions where simulations accurately capture the fluorescence anisotropy decay, we find at most a modest, transient population of conformations in which the chromophores associate with the polyproline. Explicit calculation of FRET transfer efficiencies for short polyprolines yields results in good agreement with experiment. These results illustrate the potential power of a combination of molecular simulation and experiment in quantifying biomolecular dynamics.

Higher impact energy in traumatic brain injury interferes with noncovalent and covalent bonds resulting in cytotoxic brain tissue edema as measured with computational simulation

BACKGROUND

Cytotoxic brain tissue edema is a complicated secondary consequence of ischemic injury following cerebral diseases such as traumatic brain injury and stroke. To some extent the pathophysiological mechanisms are known, but far from completely. In this study, a hypothesis is proposed in which protein unfolding and perturbation of nucleotide structures participate in the development of cytotoxic edema following traumatic brain injury (TBI).

METHODS

An advanced computational simulation model of the human head was used to simulate TBI. The consequences of kinetic energy transfer following an external dynamic impact were analyzed including the intracranial pressure (ICP), strain level, and their potential influences on the noncovalent and covalent bonds in folded protein structures.

RESULTS

The result shows that although most of the transferred kinetic energy is absorbed in the skin and three bone layers, there is a substantial amount of energy reaching the gray and white matter. The kinetic energy from an external dynamic impact has the theoretical potential to interfere not only with noncovalent but also covalent bonds when high enough. The induced mechanical strain and pressure may further interfere with the proteins, which accumulate water molecules into the interior of the hydrophobic structures of unfolded proteins. Simultaneously, the noncovalent energy-rich bonds in nucleotide adenosine-triphosphates may be perturbed as well.

CONCLUSIONS

Based on the analysis of the numerical simulation data, the kinetic energy from an external dynamic impact has the theoretical potential to interfere not only with noncovalent, but also with covalent bonds when high enough. The subsequent attraction of increased water molecules into the unfolded protein structures and disruption of adenosine-triphosphate bonds could to some extent explain the etiology to cytotoxic edema.

Single-molecule force spectroscopy reveals force-enhanced binding of calcium ions by gelsolin

Force is increasingly recognized as an important element in controlling biological processes. Forces can deform native protein conformations leading to protein-specific effects. Protein–protein binding affinities may be decreased, or novel protein–protein interaction sites may be revealed, on mechanically stressing one or more components. Here we demonstrate that the calcium-binding affinity of the sixth domain of the actin-binding protein gelsolin (G6) can be enhanced by mechanical force. Our kinetic model suggests that the calcium-binding affinity of G6 increases exponentially with force, up to the point of G6 unfolding. This implies that gelsolin may be activated at lower calcium ion levels when subjected to tensile forces. The demonstration that cation–protein binding affinities can be force-dependent provides a new understanding of the complex behaviour of cation-regulated proteins in stressful cellular environments, such as those found in the cytoskeleton-rich leading edge and at cell adhesions.

The application of force can influence biological processes such as ligand and protein–protein binding, with mechanical stress typically hindering such interactions. Here, the authors use atomic force microscopy to show that the binding of calcium to gelsolin can be improved under stress.