Hostname: page-component-586b7cd67f-gb8f7 Total loading time: 0 Render date: 2024-11-23T17:09:30.438Z Has data issue: false hasContentIssue false

Towards the understanding of molecular motors and its relationship with local unfolding

Published online by Cambridge University Press:  08 May 2024

Zahra Alavi
Affiliation:
Department of Physics, Loyola Marymount University, Los Angeles, CA, USA
Nathalie Casanova-Morales
Affiliation:
Facultad de Artes Liberales, Universidad Adolfo Ibáñez, Santiago, Chile
Diego Quiroga-Roger
Affiliation:
Biochemistry and Molecular Biology Department, Faculty of Chemistry and Pharmaceutical Sciences, Universidad de Chile, Santiago, Chile
Christian A.M. Wilson*
Affiliation:
Biochemistry and Molecular Biology Department, Faculty of Chemistry and Pharmaceutical Sciences, Universidad de Chile, Santiago, Chile
*
Corresponding author: Christian A.M. Wilson; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Molecular motors are machines essential for life since they convert chemical energy into mechanical work. However, the precise mechanism by which nucleotide binding, catalysis, or release of products is coupled to the work performed by the molecular motor is still not entirely clear. This is due, in part, to a lack of understanding of the role of force in the mechanical–structural processes involved in enzyme catalysis. From a mechanical perspective, one promising hypothesis is the Haldane–Pauling hypothesis which considers the idea that part of the enzymatic catalysis is strain-induced. It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates. Instead, they must exert strain on the substrate upon binding, using enzyme-substrate energy interaction (binding energy) to accelerate the reaction rate. A novel idea suggests that during catalysis, significant strain energy is built up, which is then released by a local unfolding/refolding event known as ‘cracking’. Recent evidence has also shown that in catalytic reactions involving conformational changes, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect the enzyme’s integrity. Thus, it has been suggested that this released heat could promote or be linked to the cracking seen in proteins such as adenylate kinase (AK). We propose that the energy released as a consequence of ligand binding/catalysis is associated with the local unfolding/refolding events (cracking), and that this energy is capable of driving the mechanical work.

Type
Review
Copyright
© Universidad de Chile, 2024. Published by Cambridge University Press

General overview of protein structure–function relationships

In the last 70 years, many studies focused on the protein structure–function relationships have been achieved to a detailed level by high-resolution structures determined by X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, Cryo-electron microscopy and Alpha fold (Henzler-Wildman et al., Reference Henzler-Wildman, Thai, Lei, Ott, Wolf-Watz, Fenn, Pozharski, Wilson, Petsko, Karplus, Hübner and Kern2007; Klinman and Hammes-Schiffer, Reference Klinman and Hammes-Schiffer2013; Nogales, Reference Nogales2016; Jumper et al., Reference Jumper, Evans, Pritzel, Green, Figurnov, Ronneberger, Tunyasuvunakool, Bates, Žídek, Potapenko, Bridgland, Meyer, Kohl, Ballard, Cowie, Romera-Paredes, Nikolov, Jain, Adler and Hassabis2021). This is reflected in the large number of protein structures (200,069 until January 10, 2023, associated with functional studies that are now included in the protein data bank (PDB), demonstrating the importance and validity of the relationship between the protein structure and functional researches.

Folding and conformational changes in proteins

The global spatial arrangement of atoms in a protein is called protein conformation. Catalysis is not only determined by chemicals steps, but also by conformational changes in protein structure. Proteins must adopt a specific conformation, or in other words, a specific three-dimensional structure (native protein) to be catalytically active (Anfinsen et al., Reference Anfinsen, Redfield, Choate, Page and Carroll1954; Carrion-Vazquez et al., Reference Carrion-Vazquez, Oberhauser, Fowler, Marszalek, Broedel, Clarke and Fernandez1999; Klinman and Hammes-Schiffer, Reference Klinman and Hammes-Schiffer2013). In fact, proteins can be found mostly folded in physiological conditions. The stability of native proteins in folded state is only marginally higher compared to the unfolded state (Taverna and Goldstein, Reference Taverna and Goldstein2002; Magliery, Reference Magliery2015), as the difference in Gibbs free energy (ΔG) separating folded and unfolded states is only about 5 to 15 kcal/mol (Creighton, Reference Creighton1990; Voet and Voet, Reference Voet and Voet2010). Protein stability studies have usually been carried out by reversibly unfolding using denaturants (chaotropic agents such as guanidine chloride and urea), increasing or decreasing the temperature, varying the pH, applying high pressures, or cleaving disulfide bonds. These studies allow for the determination of the free energies necessary to maintain the folded state, the intermediate states if they exist, the folding energy landscape, and the population of the unfolded protein (Dill, Reference Dill1985; Shea and Brooks, Reference Shea and Brooks2001).

Many biochemical reactions proceed via large conformational changes within or between interacting molecules. In the past, protein conformational changes were viewed as rigid body movements coupled to ligand binding with just one lowest free energy conformation in each condition (with and without ligand). One of the best examples is hemoglobin, where the conformational changes are coupled to the binding of an oxygen molecule in a remote site of the active site (allosterism) (Edelstein, Reference Edelstein1975). However, the current concept of native proteins considers them as ensembles of related, interconverting, transient microstates that describe the canonical high-resolution structures observed by crystallography or NMR spectroscopy (Frauenfelder et al., Reference Frauenfelder, Parak and Young1988; Hilser et al., Reference Hilser, García-Moreno, Oas, Kapp and Whitten2006; Fenimore et al., Reference Fenimore, Frauenfelder, Magazù, McMahon, Mezei, Migliardo, Young and Stroe2013). Thus, protein movements are consequence of a collection of microstates that occur at different time scales, as seen in Figure 1.

Figure 1. The many different conformational states produce a rocky landscape.

The conformations existing under a given set of conditions are usually the ones that are thermodynamically the most stable, having similar ΔG (Creighton, Reference Creighton1990; Voet and Voet, Reference Voet and Voet2010). The existence of multiple stable conformations reflects the changes that occur in proteins as they bind to other molecules or when they catalyze reactions, validating the idea that many of the physical and functional properties of proteins (stability, solubility, binding of ligands) are influenced by the same structural fluctuations that give rise to the ensemble without a ligand (Hilser et al., Reference Hilser, García-Moreno, Oas, Kapp and Whitten2006; Klinman and Hammes-Schiffer, Reference Klinman and Hammes-Schiffer2013).

Ligand-binding and conformational changes in catalysis

Protein conformational changes play a crucial role in catalysis, enabling proteins to bind to ligands, form oligomers, and perform mechanical work (Koshland, Reference Koshland1958; Whitford et al., Reference Whitford, Miyashita, Levy and Onuchic2007; Olsson and Wolf-Watz, Reference Olsson and Wolf-Watz2010). These conformational changes can occur in time scales from milliseconds to nanoseconds (Henzler-Wildman et al., Reference Henzler-Wildman, Thai, Lei, Ott, Wolf-Watz, Fenn, Pozharski, Wilson, Petsko, Karplus, Hübner and Kern2007) and involve domain movements. According to the ‘induced fit’ theory (Koshland, Reference Koshland1958), substrate binding triggers a large conformational change in the enzyme, leading to the correct positioning of the residues involved in catalysis and the release of the product. In this view of enzyme catalysis, binding interactions between the substrate and the enzyme provoke motions in enzyme’s structure that bring the substrate into the transition state, facilitating and pushing forward the catalytic process along its reaction coordinate (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004; Voet and Voet, Reference Voet and Voet2010). Kinases are generally viewed as good examples of the induced-fit mechanism (Koshland, Reference Koshland1958; Schulz et al., Reference Schulz, Müller and Diederichs1990; Choi and Zocchi, Reference Choi and Zocchi2007).

Another view of enzyme catalysis, according to the Haldane–Pauling hypothesis (Pauling, Reference Pauling1946; Haldane, Reference Haldane1965), considers the idea that part of the enzymatic catalysis is ‘strain-induced’ (also called geometric destabilization theory (Jencks, Reference Jencks1975). It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates and they must exert strain on the substrate upon binding. Furthermore, it suggests that in order to catalyze reactions, an enzyme must be complementary to the transition state species, and not to the substrate molecule itself in its normal configuration. This hypothesis considers that the energy derived from the enzyme-substrate interaction (the binding energy) is used to accelerate the rate of the reaction. Experimental evidence of strain-induced catalysis has been indirect (Amyes and Richard, Reference Amyes and Richard2013). One example is using catalytic antibodies. They recognize and bind non-covalently to their complementary antigen, then catalyze bond-breaking or other chemical reactions involving the antigen (Hanson et al., Reference Hanson, Nishiyama and Paul2005). In those works, they used molecules that were analogs to the transition state which bind more tightly to the antibody than its actual antigen or substrate (Pollack et al., Reference Pollack, Jacobs and Schultz1986; Tramontano et al., Reference Tramontano, Janda and Lerner1986; Savinov et al., Reference Savinov, Hirschmann, Benkovic and Smith2003). Another example is the study of the strain in the carbonyl bond of the substrate of beta-lactamase (Hokenson et al., Reference Hokenson, Cope, Lewis, Oberg and Fink2000). Class A beta-lactamases hydrolyze penicillins and other beta-lactams via an acyl-enzyme catalytic mechanism. Comparison of the beta-lactam carbonyl stretch frequency in the free and enzyme-bound substrate revealed an average decrease of 13 cm−1 in its frequency, indicating substantial strain/distortion of the lactam carbonyl when bound in the enzyme-substrate (ES) complex. The substrate binding to the active site induces substantial strain and distortion that contributes significantly to the overall rate enhancement in beta-lactamase catalysis (Hokenson et al., Reference Hokenson, Cope, Lewis, Oberg and Fink2000). Although these data help us to understand strain-induced catalysis, they are not enough to quantify the contribution of strain in catalysis. It is also not clear how this event is linked with the conformational changes coupled to ligand binding.

Molecular motors

Enzymes catalyze chemical reactions that occur in life, playing a key role in almost all biological events as they accelerate the rate of chemical reactions by lowering the energy barrier of the conversion of the reactants into products (Voet and Voet, Reference Voet and Voet2010; Cornish-Bowden, Reference Cornish-Bowden2012).

Motor enzymes use the energy of nucleotide binding, hydrolysis, or product release to generate mechanical work. Thus, these motor enzymes must couple one or more of those chemical steps to mechanical transitions. There are many proteins that behave as molecular machines. For example, myosin, kinesin, and dynein families use ATP hydrolysis as a source of energy to move along a track. Also, polymerases must utilize part of the chemical energy derived from the polymerization reaction to move along the DNA or RNA in a unidirectional manner (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004). Helicases hydrolyze ATP to translocate along DNA, unwinding it into its complementary strands (Cheng et al., Reference Cheng, Dumont, Tinoco and Bustamante2007). Some viruses use ATP hydrolysis to pack the DNA inside their capsid against considerable entropic, electrostatic, and elastic forces (Smith et al., Reference Smith, Tans, Smith, Grimes, Anderson and Bustamante2001). In the case of translocation, a molecular motor uses ATP binding/hydrolysis to mechanically pull polypeptide chains (a flexible, irregular polymer) across membranes.

All these motors move along stiff tracks such as double-stranded (ds) DNA or microtubules, where the track persistence length (P; or bending stiffness of a polymer) is an order of magnitude larger than the step size of the motor (Maillard et al., Reference Maillard, Chistol, Sen, Righini, Tan, Kaiser, Hodges, Martin and Bustamante2011).

Brownian ratchet versus power stroke model

There are two general models to explain how molecular motors operate: the Brownian ratchet and the Power stroke model. In the Brownian ratchet model, the motor uses nucleotide binding and/or hydrolysis to direct and rectify its Brownian motion, with a net movement in one direction (Astumian, Reference Astumian1997). In the power stroke model, the motor utilizes the energy of ATP binding and/or hydrolysis or product release to drive the motion (Smith et al., Reference Smith, Tans, Smith, Grimes, Anderson and Bustamante2001). A typical Brownian motion ratchet mechanism molecular motor is the E. coli RNA polymerase with a work (stall force times step size) value close to 2 kBT (kB is the Boltzmann constant and T is the temperature) at stall force (force at which the velocity of the motor drops to zero); (Neuman et al., Reference Neuman, Abbondanzieri, Landick, Gelles and Block2003). In the case of the bacteriophage Phi29 (a power stroke motor), the work value is about 10 kBT (Moffitt et al., Reference Moffitt, Chemla, Aathavan, Grimes, Jardine, Anderson and Bustamante2009). It is considered a power stroke model when the work is higher than 5 kBT at the stall force (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004). How do motors convert the chemical energy into mechanical work or movement? What is the mechano-chemistry of the motor? The force can be considered as a product of the nucleotide binding/hydrolysis/product release cycle, and the velocity of a motor often depends on the external force and is related to the mechanism of the motor (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004). The rate of the motor will be force-dependent if, in the conditions of the experiments, the movement is the rate-limiting step. By systematically varying the concentration of ATP, or its hydrolysis products (ADP, Pi) at different forces, it is possible to determine the force-generating step during the chemical cycle of ATP hydrolysis and to study the mechano-chemistry of molecular motors (Wang et al., Reference Wang, Schnitzer, Yin, Landick, Gelles and Block1998; Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004). These studies have been done in many systems. For example, it has been shown that the rate constants of certain sub-steps in the mechanochemical cycle of molecular motors are force-dependent (or torque-dependent in rotary motors) as seen for Myosin and F1-ATPase (Oguchi et al., Reference Oguchi, Mikhailenko, Ohki, Olivares, De La Cruz and Ishiwata2008; Sellers and Veigel, Reference Sellers and Veigel2010; Adachi et al., Reference Adachi, Oiwa, Yoshida, Nishizaka and Kinosita2012; Watanabe et al., Reference Watanabe, Okuno, Sakakihara, Shimabukuro, Lino, Yoshida and Noji2012; Cossio et al., Reference Cossio, Hummer and Szabo2015; Volkán-Kacsó and Marcus, Reference Volkán-Kacsó and Marcus2015; Houdusse and Sweeney, Reference Houdusse and Sweeney2016). Molecular motors are enzymes. Thus to understand the mechanisms that govern their function we should focus on how enzymes work.

Importance of forces in mechanochemical processes

Different processes inside the cell, such as chromosomal segregation, DNA transcription, DNA replication, the formation and liberation of vesicles, the packing of DNA during viral replication, and the translocation of protein through channels, are mechanical processes (Schekman, Reference Schekman1994; Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004; Cecconi et al., Reference Cecconi, Shank, Bustamante and Marqusee2005; Zhang et al., Reference Zhang, Halvorsen, Zhang, Wong and Springer2009). As mentioned before, protein folding and unfolding, conformational changes, and catalysis also represent good examples of mechanical processes studied by the application of a mechanical force. One of the advantages of studying the reactions with a focus on force is that the reaction coordinate is an easily quantifiable physical parameter. The kinetics of a reaction can be studied by the application of a mechanical force (Bell, Reference Bell1978). An applied mechanical force affects the free energy, equilibrium, and rate of a reaction occurring along a mechanical reaction coordinate, so thermodynamic and kinetic parameters could be determined at zero force extrapolating the information obtained at different forces. Recently, new transition state theories stemming from the application of Kramers’s theory (Kramers, Reference Kramers1940) have emerged (Dudko et al., Reference Dudko, Hummer and Szabo2008; Cossio et al., Reference Cossio, Hummer and Szabo2016; Bullerjahn et al., Reference Bullerjahn, Sturm and Kroy2020). Another interesting and useful application of force in the study of conformational changes and catalysis, called ‘the allosteric spring probe’ (ASP), was developed by the group of G. Zocchi (Choi and Zocchi, Reference Choi and Zocchi2007), using a guanylate kinase (GK) as a model. The novelty of this technique is in the attachment of a single-stranded DNA spring to a protein, followed by using complementary strands of different sizes to apply a mechanical force, generating useful information about the mechanism of the reaction and allowing to assess quantitatively the relevance of the induced-fit theory in this enzyme.

Force is also very important in many biological signaling processes. For example, in order to start the coagulation process the Von Willebrand factor unfolds via shear forces, and exposes the binding site to start the process (Zhang et al., Reference Zhang, Halvorsen, Zhang, Wong and Springer2009). It has been shown that the folding force of SNARE proteins is enough to achieve membrane fusion in the vesicle release at the synapse (Gao et al., Reference Gao, Zorman, Gundersen, Xi, Sirinakis, Rothman and Zhang2012; Zhang et al., Reference Zhang, Rebane, Ma, Ma, Li, Jiao, Qu, Pincet, Rothman and Zhang2016). Another study has shown that the refolding of the Top7 protein is enough to release a stall sequence from the ribosome (Goldman et al., Reference Goldman, Kaiser, Milin, Righini, Tinoco and Bustamante2015). All these studies show that refolding force is crucial for biochemical processes in all cells.

Local unfolding/refolding (cracking)

In recent years, various studies have added a higher level of complexity to enzyme catalysis. These studies include computational molecular dynamics such as mixed Go model (Miyashita et al., Reference Miyashita, Onuchic and Wolynes2003; Whitford et al., Reference Whitford, Miyashita, Levy and Onuchic2007), or laboratory bench experiments such as ITC (Isothermal Titration Calorimetry) (Schrank et al., Reference Schrank, Bolen and Hilser2009). Entropy-promoting mutations (Olsson and Wolf-Watz, Reference Olsson and Wolf-Watz2010) have also provided evidence of the existence of local unfolding/refolding events (also known as ‘cracking’) during catalysis (Klinman and Hammes-Schiffer, Reference Klinman and Hammes-Schiffer2013; Schrank et al., Reference Schrank, Wrabl and Hilser2013). Adenylate kinase (AK) is an excellent example of a protein that undergoes cracking during catalysis, where significant strain energy (>20 kcal/mol) is built up during the reaction, and then released through a local unfolding event at the transition state. The protein then refolds at the downhill side of the activation energy barrier and the thermodynamic minimum is reached (Miyashita et al., Reference Miyashita, Onuchic and Wolynes2003; Whitford et al., Reference Whitford, Miyashita, Levy and Onuchic2007; Olsson and Wolf-Watz, Reference Olsson and Wolf-Watz2010). The novel Cracking idea points to make more complex the interconversion between the two simple rigid-body folded open and closed states. By means of connecting the initial and final states through a local unfolding and refolding event and adding another equilibrium step between the folded and unfolded states coupled to ligand binding or catalysis.

Another cracking event has been described for the SecA protein. SecA is a highly conserved and essential ATP-dependent motor protein of the bacterial Sec translocase machinery. This protein uses cracking to control its activation process in a precise manner and adopt alternate conformational states during the ATPase cycle (Keramisanou et al., Reference Keramisanou, Biris, Gelis, Sianidis, Karamanou, Economou and Kalodimos2006). Moreover, evidence shows that in some catalytic reactions with conformational changes coupled to catalysis, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect enzyme’s integrity. Thus, it has been suggested that this released heat could be promoting or be linked with cracking seen in AK (Riedel et al., Reference Riedel, Gabizon, Wilson, Hamadani, Tsekouras, Marqusee, Pressé and Bustamante2015). This is because enough energy is released upon catalysis, and the enthalpy change (ΔH) of the reactions that proteins catalyze is high enough to unfold the protein (Riedel et al., Reference Riedel, Gabizon, Wilson, Hamadani, Tsekouras, Marqusee, Pressé and Bustamante2015).

Interestingly the theoretical ΔH in some reactions differs from the experimental one in the presence of enzymes (Table 1), with the experimental ΔH being lower. This difference may be due to some of the energy being used to perform work in the enzyme, such as moving the center of mass of the protein or locally unfolding some region of it (Riedel et al., Reference Riedel, Gabizon, Wilson, Hamadani, Tsekouras, Marqusee, Pressé and Bustamante2015). Overall, these findings highlight the complexity of enzyme catalysis and the role that cracking events may play in the process.

Table 1. Theoretical and experimental characterization of enthalpies

This table shows the enzyme and the reaction with which the theoretical $ \boldsymbol{\Delta }\mathrm{H} $ enthalpy is calculated based on thermodynamical tables (Tinoco et al., Reference Tinoco, Sauer, Wang, Puglisi, Harbison and Rovnyak2013). Experimental enthalpy $ \boldsymbol{\Delta }{\mathrm{H}}_{\mathrm{experimental}} $ is that which is measured in an experiment under the conditions specified by the studies referenced in the last column.

Dynamic of local unfolding and the generation of work

But how fast will the protein refold? For folding of a complete protein, the experimental and theoretical approaches predict a speed limit of approximately N/100 μs for a generic N-residue single-domain protein (Kubelka et al., Reference Kubelka, Hofrichter and Eaton2004; Dill et al., Reference Dill, Ozkan, Shell and Weikl2008). Then, a small portion of the protein could fold very fast because of its size which is compatible with the catalysis timing. For example, each catalytic cycle in AK takes 25–40 ms and the size of the local unfolded region is around 20 amino acids, predicting a refolding time of 0.2 μs, one order of magnitude faster than the turnover number (Pelz et al., Reference Pelz, Žoldák, Zeller, Zacharias and Rief2016).

Despite the various hypotheses surrounding enzyme catalysis, there is still uncertainty about where the catalytic power of the enzyme resides. Most of these ideas have remained as hypotheses, partly due to the difficulty of designing experiments to demonstrate them. While the induced-fit mechanism is widely accepted (Klinman and Hammes-Schiffer, Reference Klinman and Hammes-Schiffer2013), the ‘strain-induced’ idea remains largely hypothetical as bulk chemical experiments provide only indirect support for its validity (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004). Moreover, it is challenging to correlate conventional bulk experiments with in silico dynamic simulations to demonstrate the existence of cracking. Therefore, it is crucial to find a molecular description of the mechanisms by which enzymes achieve high catalytic efficiency coupled to conformational changes, demonstrating the importance of strain in these processes. A modern approach to address this question is to consider protein folding and unfolding, conformational changes, and catalysis as mechanical processes, where basic physical concepts such as force, torque, work, energy conversion efficiency, and mechanical advantage can be determined to describe them (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004).

Based on this knowledge, we could hypothesize that the energy released by binding/catalysis in an enzyme is coupled to local unfolding and the refolding is the process that allows the mechanical work.

Another indication of local unfolding is the existence of intermediates. Intermediates of the folding pathway are very common even for very small proteins (Rivera et al., Reference Rivera, Mjaavatten, Smith, Baez and Wilson2023b). Even for typical two-states proteins, if we change more parameters to unfold the protein we could see an intermediate (Rivera et al., Reference Rivera, Mjaavatten, Smith, Baez and Wilson2023b). These intermediates could be related to the local unfolding of proteins and may be a universal characteristic. In the past unfolding regions were seen as non-functional regions of the protein. Today the view is completely contrary. There are many proteins that have disordered regions that are essential for their functions. Such proteins are known as intrinsically disordered proteins (IDPs)(Oldfield and Dunker, Reference Oldfield and Dunker2014). It has been shown that these disordered regions can fold upon binding to another protein (Malagrinò et al., Reference Malagrinò, Visconti, Pagano, Toto, Troilo and Gianni2020).

In our view of molecular motors, an unfolded region will be very important to perform work once it is re-folded. As previously mentioned, the refolding force is enough to perform work in some systems such as SNAREs, SecA, and IDPs which can fold under different conditions such as substrate binding, protein–protein interaction, and so on (Gao et al., Reference Gao, Zorman, Gundersen, Xi, Sirinakis, Rothman and Zhang2012; Goldman et al., Reference Goldman, Kaiser, Milin, Righini, Tinoco and Bustamante2015).

Enzymes as minimal molecular motors and instruments to study them

Our approach is to consider enzymes as molecular motors converting chemical energy (either in the form of binding energy, chemical bond hydrolysis, or electrochemical gradients) into mechanical work through conformational changes and displacements (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004; Zocchi, Reference Zocchi2018). The energy accumulated as a consequence of ligand binding or product release could be used to execute conformational changes during catalysis. We believe that some ATPases are minimum molecular motors, that the movement of domains coupled to catalysis can be considered as the simplest motor in nature and that with single molecule studies we can understand their behavior. By applying an external force over the moving domains, the ES complex formation, or the transition state attainment can be affected. So, if the conformational change is rate-limiting, the external force will significantly affect the rate of the reaction, providing direct evidence for the effect of strain in enzyme catalysis. Moreover, by varying the magnitude of the external force applied to the enzyme-substrate complex, it should be possible to estimate the maximum force generated within the complex. Therefore, by the optical tweezer single-molecule approach changes can be monitored as a function of velocity and force, allowing for the determination of the rate-limiting step of ligand binding and catalysis, and quantifying the force needed to stop catalysis (‘stall force’) (Bustamante et al., Reference Bustamante, Chemla, Forde and Izhaky2004; Choi and Zocchi, Reference Choi and Zocchi2007). Thus, this force can be correlated with the force that the enzyme exerts on the substrate, which in fact is the strain that the enzyme exerts on the substrate, according to the stress-induced catalysis hypothesis (Pauling, Reference Pauling1946; Haldane, Reference Haldane1965). Next, we will look at three proteins through the lense of cracking as models to study the force in unfolding events in molecular motors.

Aquifex aeolicus adenylate kinase

A suitable enzyme to study as a minimal molecular motor should involve a big conformational change upon substrate binding, it should display a correlation between turnover number and the fluctuations in the domain closure, the turnover number should not be too fast (slower than the detection system), and the molecule should be preferably a monomer (Figure 2a). In this context AK from the thermophilic organism Aquifex aeolicus and GK from Mycobacterium tuberculosis are good candidates. AK (EC 2.7.4.3) is a small (20–26 kDa) monomeric ubiquitous signal transducing enzyme that is responsible for maintaining cellular steady-state concentration of adenylate nucleotides in the cell (Schrank et al., Reference Schrank, Wrabl and Hilser2013). This enzyme has been studied in detail for decades and much is known about its structure, function, kinetics, and taxonomy. This enzyme catalyzes the reversible phosphoryl transfer reaction, using Mg2+ or Mn2+ as a cofactor, according to: $ {\mathrm{Mg}}^{2+}+\mathrm{ATP}+\mathrm{AMP}\leftrightarrow {\mathrm{Mg}}^{2+}+2\mathrm{ADP} $ (Hamada et al., Reference Hamada, Sumida, Kurokawa, Sunayashiki-Kusuzaki, Okuda, Watanabe and Kuby1985; Schulz et al., Reference Schulz, Müller and Diederichs1990). From a structural point of view, this protein is modular and is formed by three subdomains: the core subdomain that is very important for the overall stability, the ATP and AMP binding subdomains that have the ATP and AMP binding sites, and the lid subdomain that increases catalytic efficiency by covering the binding sites. The X-ray crystal structures of AK have been solved in an open (substrate-free) and a closed (inhibitor-bound) conformation (Lienhard and Secemski, Reference Lienhard and Secemski1973; Müller and Schulz, Reference Müller and Schulz1992). Analysis of the three-dimensional structures of AK, in both the bound and unbound states, reveals that the lid subdomains and the ATP and AMP subdomains undergo a significant conformational change in response to substrate binding (Schulz et al., Reference Schulz, Müller and Diederichs1990). This movement has been described as domain closure, which in the case of AK is 11 Å (Schulz et al., Reference Schulz, Müller and Diederichs1990; Sinev et al., Reference Sinev, Sineva, Ittah and Haas1996), being in the closed conformation where the phosphoryl transfer reaction occurs (Müller and Schulz, Reference Müller and Schulz1992; Ådén and Wolf-Watz, Reference Ådén and Wolf-Watz2007) and as explained previously has been shown that it undergoes cracking. Since the resolution limit of the technique is at least 10 Å, AK is a good candidate for the study of its catalysis with the optical tweezers (Wilson et al., Reference Wilson, Leachman, Cervantes, Ierokomos, Marqusee and Bustamante2013; Pelz et al., Reference Pelz, Žoldák, Zeller, Zacharias and Rief2016). AK has a catalytic constant (k cat) of 40 sec−1 (40 catalysis events per second) which is well within the instrumental capabilities, since it is possible to determine cycles of catalysis of 40 events or less with our bandwidth of 1000 Hz. The k close and k open of this enzyme were determined to be 44 and 1600 Hz, respectively (Henzler-Wildman et al., Reference Henzler-Wildman, Thai, Lei, Ott, Wolf-Watz, Fenn, Pozharski, Wilson, Petsko, Karplus, Hübner and Kern2007). Some correlations between the enzyme flexibility and catalysis have already been described using bulk techniques (Antikainen and Martin, Reference Antikainen and Martin2005). For instance, the work of Wolf-Watz et al. (Reference Wolf-Watz, Thai, Henzler-Wildman, Hadjipavlou, Eisenmesser and Kern2004) and Henzler-Wildman et al. (Reference Henzler-Wildman, Thai, Lei, Ott, Wolf-Watz, Fenn, Pozharski, Wilson, Petsko, Karplus, Hübner and Kern2007) found strong correlations between the k open and the turnover number by comparing the mesophilic and the thermophilic AK. A work from (Ben Ishay et al., Reference Ben Ishay, Rahamim, Orevi, Hazan, Amir and Haas2012) using time-resolved FRET had been shown that AK can refold locally a domain of the protein within the dead time of the stopped-flow device, showing that this local refolding is below 5 ms. This shows that the protein could unfold and refold faster than the catalytic cycle. If this intermediate corresponds to the cracked region with faster techniques and/or transition paths analysis (Cossio et al., Reference Cossio, Hummer and Szabo2018) we should be able to isolate this metastable state.

Figure 2. Single-molecule force spectroscopy techniques (open and closed conformations of the proteins to be studied). (a) Left: Open apo Aquifex AK (PDBid: 2RH5). Middle: Local unfolded or crack intermediate. Right: Close Aquifex AK in complex with the substrate analogue Zn2 + ● Ap5A (light blue) (PDBid: 2RGX) (Olsson and Wolf-Watz, Reference Olsson and Wolf-Watz2010). In yellow is the lid domain and in gray is the core of the protein. Reprinted by permission from Springer Nature. (b) Left: Open configuration GK (PDBid: 1ZNW). Right: Close configuration of GK (PDBid: 1LVG) upon substrate binding. Alpha helices are colored magenta and beta sheets are colored yellow. (Alavi, Reference Alavi2017). (c) Open configuration BiP in the presence of ADP (PDBids: 5EVZ and 5E85) and right. Close configuration of BiP in the presence of ATP (PDBid: 5E84). In blue is the NBD and orange is SBD (Pobre et al., Reference Pobre, Poet and Hendershot2019). Images were made with 3D Protein imaging software (Tomasello et al., Reference Tomasello, Armenia and Molla2020).

Another interesting feature of AK is that this enzyme has an important role in health. There exists a lot of experimental data validating the direct relationship between AK and human diseases. For example, knockout of the major human AK isoform (AK1) demonstrated a deficit in vascular and myocardial AK catalysis blunts energetic communication, AMP, and adenosine metabolic signal transduction. This compromises post-ischemic coronary reflow, providing evidence for the role of AK as a metabolic monitor supporting regulation of the reactive vascular response in the stressed heart (Dzeja et al., Reference Dzeja, Bast, Pucar, Wieringa and Terzic2007). AK1 deficiency has also been related to hemolytic anemia (Abrusci et al., Reference Abrusci, Chiarelli, Galizzi, Fermo, Bianchi, Zanella and Valentini2007). It has been observed that human AK isoform (AK4) is upregulated in lung adenocarcinoma compared with normal cells, suggesting that AK4 promotes malignant progression and recurrence by promoting metastasis in a transcription factor ATF3-dependent manner (Jan et al., Reference Jan, Tsai, Yang, Huang, Yang, Lai, Lee, Jen, Huang, Su, Chuang and Hsiao2012). Therefore, understanding the mechanisms by which AK exerts its efficient catalysis is of great importance as it could help to unravel the causes and to develop therapies involving the AK function. In this respect, we have determined the sequence identity between Aquifex aeolicus adenylate kinase (Aquifex AK) and human AK, establishing that all of them have substantial identity (over 90% between Aquifex AK and Human AK1, AK2, and AK3 isoenzymes).

Mycobacterium tuberculosis guanylate kinase

The other protein of interest is GK which catalyzes the ATP-dependent phosphorylation of GMP into GDP (ATP + GMP $ \rightleftharpoons $ ADP + GDP) (Stehle and Schulz, Reference Stehle and Schulz1992). In humans, GK is the only known enzyme responsible for GDP production and thus is essential for cellular viability and might be a potential target for cancer therapeutics (Khan et al., Reference Khan, Shah, Ban, Trigo-Mouriño, Carneiro, DeLeeuw, Dean, Trent, Beverly, Konrad, Lee and Sabo2019). This protein is around 24 kDa and 4 nm long, similar in size and conformation to AK. As seen in Figure 2b, GK’s structure consists of a clamp-like cavity with three distinct regions: a central CORE domain, a LID which closes onto the CORE domain during catalysis and a GMP binding domain (GMP-BD). Additionally, a canonical P-loop motif is involved in ATP binding (Choi and Zocchi, Reference Choi and Zocchi2007; Khan et al., Reference Khan, Shah, Ban, Trigo-Mouriño, Carneiro, DeLeeuw, Dean, Trent, Beverly, Konrad, Lee and Sabo2019).

Upon substrate binding the two lobes of the clamp close through an ∼1-nm conformational change, most of which is induced by GMP binding. This conformational change that is relatively large compared to the size of the enzyme involves several direct and indirect (via water molecules) interactions which bring the LID and the GMP-BD jointly toward the CORE region (Choi and Zocchi, Reference Choi and Zocchi2007). This substrate-induced conformational change turns the enzyme from the open state to the close state, but it is not known if this conformational change takes place via cracking. In this view, the catalytic cycle of the enzyme can be described by the kinetic cycle of conformational changes, an idea that was proposed more than 50 years ago (Bliumenfel’d, Reference Bliumenfel’d1971) and experimentally studied using NMR and Nanorheology (Eisenmesser et al., Reference Eisenmesser, Bosco, Akke and Kern2002; Qu and Zocchi, Reference Qu and Zocchi2013).

In 2007, Choi and Zocchi performed an elegant experiment on GK by inserting an externally controllable molecular spring on the protein (an allosteric spring probe, ASP). In this experiment, 60 mer DNA was used as the probe which in ss form exerts essentially no tension but significantly rigidifies upon hybridization with a complementary strand and therefore exerts a mechanical stress on the protein (Choi and Zocchi, Reference Choi and Zocchi2007).

They found that when the open structure is favored through mechanical stress the binding affinity for GMP is drastically reduced, whereas the binding affinity for ATP and the catalytic rate were unaffected. This result shows that GMP binding does not allosterically control ATP binding but must allosterically control catalysis or at least the ATP hydrolysis step, as ATP would not hydrolyze in the absence of bound GMP.

Another novel experiment in Zocchi’s lab was studying GK’s mechanical properties through a technique they call nanorheology (Figure 3). In this technique the protein tethers Gold nanoparticles, to a gold surface. The particles are driven by an AC electric field while their displacement is synchronously detected by evanescent wave scattering, yielding the mechanical response function of the protein under study in the frequency domain (Wang and Zocchi, Reference Wang and Zocchi2010). This technique can be thought of as using a rheometer (Barnes and Hutton, Reference Barnes and Hutton1989) to measure the deformation of the substance under study, in this case, the protein, in nanometer scale. Analyzing the frequency response of the protein one can determine its dynamics. For GK the results showed that the amplitude of deformation decreases as 1/ $ \omega $ (the inverse of frequency of the applied force) for low frequencies and it plateaus at higher frequencies. Therefore, it was concluded that the protein can be modeled as a viscoelastic material: its mechanics are best described by a viscous fluid at low frequencies of the applied force and by an elastic spring at high frequencies. The frequency beyond which the protein shows elastic behavior is the corner frequency. Another interpretation of this result is that the system is dissipative in low frequencies and non-dissipative in high frequencies. It has been suggested (Alavi et al., Reference Alavi, Ariyaratne and Zocchi2015) that this internal dissipation can be associated to the dynamics of the hydration layer of the enzyme. It was shown that perturbing the hydration layer of GK through small amounts of chaotropic and kosmotropic agents can have a significant effect on its internal dissipation(Casanova-Morales et al., Reference Casanova-Morales, Alavi, Wilson and Zocchi2018a). Furthermore, the closed state of the enzyme (GMP bound) was shown to behave dramatically different than the open state in the viscoelastic regime while behaving similarly in the elastic regime (Wang and Zocchi, Reference Wang and Zocchi2010; Qu et al., Reference Qu, Landy and Zocchi2012). The open state is therefore ‘softer’ as it can access the soft more dissipative modes easier through the hinge motion which is related to the viscoelastic transition. In this view, the hinge would flow which can be a result of localized melting and refolding as was previously suggested in studies using structure-based coarse-grained simulations (Miyashita et al., Reference Miyashita, Onuchic and Wolynes2003; Whitford et al., Reference Whitford, Miyashita, Levy and Onuchic2007).

Figure 3. Nanorheology setup, showing the flow chamber with enzyme-tethered GNPs (enzyme shown in green and GNPs shown as golden spheres), the parallel plates capacitor geometry used for mechanical excitation, and the evanescent wave scattering optics used for readout.

Therefore, the hinge motion can be regarded as a dissipative transition between the extended and compact conformations of the protein. When the enzyme is in its unbound state, the structure between the two conformations can be driven by applying an external oscillatory force. Then in each conformational cycle, the structure hops from one state to another on the energy landscape, this hopping is what is regarded as ‘viscoelastic transition’ in nanorheology experiments and ‘cracking’ in molecular dynamic simulations. In this view, the dissipative viscous regime is extended for the unbound state.

However, when the enzyme is in its bound state, the structure is more confined with minimal oscillations around the bound state, that is, smaller hops locally closer to the bound state on the energy landscape. One can think of this ‘viscoelastic transition’ or ‘cracking’ as local regions of the protein unfolding during conformational transitions in order to reduce the significant strain in the hinge region.

This viscoelastic model however only describes the behavior of the enzyme versus frequency and not versus force. A non-equilibrium phase diagram was obtained in the frequency force plane separating liner elasticity dynamics from softer viscoelastic dynamics (Qu et al., Reference Qu, Landy and Zocchi2012).

Therefore, enzyme’s mechanics is non-linear and displays a softening transition as a function of force as seen in experiments where the frequency of the applied force is kept constant while the amplitude of it varies. The critical force at which the transition happens depends on its frequency. Similarly, the corner frequency depends on the amplitude of the applied force. This nonlinearity can also be explained by the concept of barrier crossing in an energy landscape, similar to the problem of escape over a barrier for a particle in thermal bath which suggests local unfolding, namely cracking, can be involved in these conformational motions(Zocchi, Reference Zocchi2018). In this view, the corner frequency can be regarded as the rate of breaking the specific bond structure of the ground state confirmation of the molecule, namely the escape rate from a barrier in the energy landscape. This rate increases with the force, that is, ωc increases as the amplitude of the external force increases. In the same view, linear elasticity regime would extend at higher frequencies: at higher frequencies (on shorter timescales) a bigger force is needed to drive the structure from one regime to another. This nonlinearity is consistent with the idea of local unfolding as a relaxing mechanism in regions with high strain. The rate of this cracking would depend on the amplitude of the applied force and a bigger force is required to have local unfolding in shorter timescales.

BiP: The master regulator of the ER

One of the most important chaperones in the ER is the molecular motor BiP protein (immunoglobulin heavy-chain Binding Protein, around 72 kDa and 8 nm long, Figure 2c). BiP (also known as HSPA5 and Grp78), a monomeric ATPase conserved across species, has been referred to as the master regulator of the ER because of the broad and crucial roles that it plays in ER processes and functions (Hendershot, Reference Hendershot2004; Alfaro-Valdés et al., Reference Alfaro-Valdés, Burgos-Bravo, Casanova-Morales, Quiroga-Roger and Wilson2018), such as protein synthesis, folding, assembly, circadian cycle and as molecular motor in translocation across the ER (Zimmermann et al., Reference Zimmermann, Eyrisch, Ahmad and Helms2011; Behnke et al., Reference Behnke, Feige and Hendershot2015; Pickard et al., Reference Pickard, Chang, Alachkar, Calverley, Garva, Arvan, Meng and Kadler2019). Although at the molecular level, the study of BiP is still in its early stages, some research groups have published findings of great value. These findings suggest that this protein could be a key player in various fields, such as in the detection and treatment of serious diseases (neurodegenerative diseases, cancer, heart diseases, among others; (Shields et al., Reference Shields, Panayi and Corrigall2012; Kosakowska-Cholody et al., Reference Kosakowska-Cholody, Lin, Srideshikan, Scheffer, Tarasova and Acharya2014; Park et al., Reference Park, Eun Kim, Morales, Moda, Moreno-Gonzalez, Concha-Marambio, Lee, Hetz and Soto2017; Ichhaporia and Hendershot, Reference Ichhaporia and Hendershot2021). Also, it has been implicated recently in the binding of the coronavirus SARS-CoV-2 spike protein to its membrane receptor (Dores-Silva et al., Reference Dores-Silva, Cauvi, Coto, Kiraly, Borges and De Maio2020; Ibrahim et al., Reference Ibrahim, Abdelmalek, Elshahat and Elfiky2020; Shin et al., Reference Shin, Toyoda, Nishitani, Fukuhara, Kita, Otsuki and Shimomura2021). Structurally, BiP is formed by two domains: a nucleotide-binding domain (NBD) with ATPase activity, connected by a flexible hydrophobic linker to the substrate-binding domain (SBD). The SBD can be further divided into a compact β-sandwich domain harboring a cleft for substrate binding and an α-helical domain at its C-terminal end, the so-called ‘lid’ (Zhu et al., Reference Zhu, Zhao, Burkholder, Gragerov, Ogata, Gottesman and Hendrickson1996). Many conformational changes, such as the opening and closing movement of the lid and the variation in the distance between the SBD and NBD, have been associated with the ATPase cycle of BiP in the ER. Once BiP binds K+ and ATP, its NBD and SBD domains come into close proximity to each other and the lid of the SBD opens, adopting a form that binds substrate proteins with low affinity. Also, several BiP cofactors have been discovered that assist in controlling the substrate-binding cycle and its localization within the ER (Otero et al., Reference Otero, Lizák and Hendershot2010; Braakman and Hebert, Reference Braakman and Hebert2013; Pobre et al., Reference Pobre, Poet and Hendershot2019).

Binding of BiP to its substrate

After the Mg2+ dependent hydrolysis of ATP, BiP enters a state with low on and off rates for substrate binding/unbinding, respectively (Behnke et al., Reference Behnke, Feige and Hendershot2015). For elongated peptide substrates (such as the one that translocate), the lid closes over the bound substrate; for globular substrates, there are direct interactions between the lid and the substrate, although the lid may not be completely closed (Behnke et al., Reference Behnke, Feige and Hendershot2015). The SBD and NBD become farther apart upon substrate binding and ATP hydrolysis, being less pronounced for globular substrates. ADP must be exchanged for ATP to release the protein substrate and make BiP available for another round of client binding. Ca2+ increases the affinity for ADP, whereas Nucleotides Exchanges Factors (NEF) such as Grp170 and Sil1 facilitate the nucleotide exchange reaction (Behnke et al., Reference Behnke, Feige and Hendershot2015). Conformational changes in murine BiP during the ATPase cycle have been determined by Förster Resonance Energy Transfer (FRET) at the single molecule level, showing that NBD and SBD come into close contact with a mean distance of 58 Å–75 Å (Marcinowski et al., Reference Marcinowski, Höller, Feige, Baerend, Lamb and Buchner2011). Additionally, by using NMR residual dipolar coupling, spin labeling, and dynamic methods, it has been determined in DnaK (a bacterial Hsp70 protein similar to BiP) that the NBD and the SBD are loosely linked and can move in angles of 35° with respect to each other (Bertelsen et al., Reference Bertelsen, Chang, Gestwicki and Zuiderweg2009). Moreover, the distance between the base and the lid of the SBD domain in Hsp70 has been calculated to be 77 Å by means of FRET (Mapa et al., Reference Mapa, Sikor, Kudryavtsev, Waegemann, Kalinin, Seidel, Neupert, Lamb and Mokranjac2010). Also, a crystal structure of human BiP bound to ATP that shows similar distances has been described (Yang et al., Reference Yang, Nune, Zong, Zhou and Liu2015). The conformational changes and movements of BiP are not independent for each domain because an important allosteric communication and coupling exists between them (Marcinowski et al., Reference Marcinowski, Höller, Feige, Baerend, Lamb and Buchner2011; Casanova-Morales et al., Reference Casanova-Morales, Quiroga-Roger, Alfaro-Valdés, Alavi, Lagos-Espinoza, Zocchi and Wilson2018b). In this study, the authors determined that nucleotide binding resulted in concerted domain movements of BiP (Marcinowski et al., Reference Marcinowski, Höller, Feige, Baerend, Lamb and Buchner2011). Conformational transitions of the lid domain allowed BiP to discriminate between peptide and protein substrates. Without single molecule approaches it is very difficult to determine how BiP binds to its substrate, since the substrate of BiP is an unfolded peptide, and if we unfold the substrate, we may also unfold BiP. However, by optical tweezers manipulation, we can specifically unfold the substrate without affecting BiP. We recently developed a method to measure how BiP binds to its substrate using optical tweezers (Ramírez et al., Reference Ramírez, Rivera, Quiroga-Roger, Bustamante, Vega, Baez, Puchner and Wilson2017; Rivera et al., Reference Rivera, Burgos-Bravo, Engelberger, Asor, Lagos-Espinoza, Figueroa, Kukura, Ramírez-Sarmiento, Baez, Smith and Wilson2023a). A tethered protein is pulled and relaxed by its N-and C-terminus to mechanically unfold and refold it, respectively, while recording the force and the trap position and measuring the time during which the protein does not refold after an unfolding event. This allows us to obtain the time that BiP is bound (or the inverse k off) to its substrates while the substrate protein remains unfolded (Ramírez et al., Reference Ramírez, Rivera, Quiroga-Roger, Bustamante, Vega, Baez, Puchner and Wilson2017). Also, we developed chimera proteins in which a portion of the protein is unfolded and another is unfolded, that allows to study how BiP works in translocation (Alfaro-Valdés, Reference Alfaro-Valdés2019).

BiP’s mechanochemical mechanism

In spite of the crucial roles of BiP during translocation, it is not fully understood whether the action of BiP is through an active mechanism of pulling (as a power stroke), mediated by the binding/hydrolysis of ATP, or as a ratchet mechanism. The hypothesis of the ratchet mechanism is supported by employing antibodies against the polypeptide chains passing through the ER lumen (Matlack et al., Reference Matlack, Misselwitz, Plath and Rapoport1999) and using a protection assay with substrates that unfold at different forces (Alfaro-Valdés, Reference Alfaro-Valdés2019). Evidence for the translocation mechanism has been obtained using coarse-grained model simulations (Assenza et al., Reference Assenza, De Los Rios and Barducci2015). This study suggests that Hsp70 chaperones use an ‘entropic pulling mechanism’ applying a force of about 15 pN, and proposes that the Hsp70’s would use a combination of ratchet and power stroke mechanisms (De Los Rios et al., Reference De Los Rios, Ben-Zvi, Slutsky, Azem and Goloubinoff2006). Translocation in all eukaryotes is likely to be similar to yeast, because of the high identity of amino acids between their channels. The channel interacts with the Sec62/Sec63 complex, with BiP acting as a molecular motor to bias the passive movement of a polypeptide in the Sec61 channel. In bacterial post-translational translocation, the channel interacts with the cytoplasmic ATPase SecA. SecA moves polypeptides through the SecY channel to the periplasm by a ‘push and slide’ mechanism (Bauer et al., Reference Bauer, Shemesh, Chen and Rapoport2014). Archaea probably use both co-translational and post-translational translocation, but it is unknown how post-translational translocation occurs because these organisms lack SecA, Sec62/Sec63 complex and BiP (Zimmermann et al., Reference Zimmermann, Eyrisch, Ahmad and Helms2011; Park and Rapoport, Reference Park and Rapoport2012). Thus, these motor enzymes must couple one or more chemical steps to perform mechanical work. To understand how BiP is doing its work during post-translational translocation, it is necessary to consider that there must be a concerted mechanism engaging mechanical work with changes in the conformational arrangements (or local unfolding) of BiP, and now with single-molecule experiments is it possible establish the correlation between them. A recent paper shows a direct observation of chemo-mechanical coupling in DnaK by optical tweezers (Singh et al., Reference Singh, Rief and Žoldák2022). They observe that the SBD lid closure is strictly coupled to the chemical steps of the ATP hydrolysis cycle showing a clear domain motion dependent on this chemical step.

Viscoelastic behavior of BiP in nanorheology

The previously described nanorheology technique allows for studying the mechanical behavior of the protein in bulk and obtaining physical properties of the protein (stiffness and elasticity) which could be explained by local unfolding.

As observed for GK, BiP protein in the folded state behaves like a viscoelastic material. The protein becomes softer when bound to nucleotides ADP and ATP. In the presence of ADP, the linker separating the SBD and NBD domains is elongated which leads to a significant decrease in rigidity. In the presence of ATP, the lid is more flexible and the domains are closer (Bertelsen et al., Reference Bertelsen, Chang, Gestwicki and Zuiderweg2009; Mapa et al., Reference Mapa, Sikor, Kudryavtsev, Waegemann, Kalinin, Seidel, Neupert, Lamb and Mokranjac2010; Casanova-Morales et al., Reference Casanova-Morales, Quiroga-Roger, Alfaro-Valdés, Alavi, Lagos-Espinoza, Zocchi and Wilson2018b; Yokoyama et al., Reference Yokoyama, Fujii, Ostermann, Schrader, Nabeshima and Mizuguchi2022). However, BiP becomes more rigid when bound to the HTFPAVL peptide substrate because the lid of BiP is closed, generating a compact and thus rigid state.

Furthermore, peptide binding was reported to dramatically increase the affinity for ADP (Casanova-Morales et al., Reference Casanova-Morales, Quiroga-Roger, Alfaro-Valdés, Alavi, Lagos-Espinoza, Zocchi and Wilson2018b). This shows that the connection between the SBD and NBD domains presents an allosteric coupling (Swain et al., Reference Swain, Dinler, Sivendran, Montgomery, Stotz and Gierasch2007; Chakafana et al., Reference Chakafana, Zininga and Shonhai2019).

Nanorheology versus optical tweezers: Complimentary techniques to study local unfolding

Nanorheology and optical tweezers are complementary techniques because they operate in different force and frequency regimes. Usually in optical tweezers (for a detailed description of this instrument please see(Smith et al., Reference Smith, Cui and Bustamante2003; Bustamante et al., Reference Bustamante, Kaiser, Maillard, Goldman and Wilson2014; Sánchez et al., Reference Sánchez, Robeson, Carrasco, Figueroa, Burgos-Bravo, Wilson and Casanova-Morales2022) the force versus extension plot starts around a few picoNewtons to 67 piconewtons and in nanorheology the force is estimated to be in the low piconewton regime (Qu and Zocchi, Reference Qu and Zocchi2013). A typical tweezers pulling is around 100 nm/sec and the stiffness is 0.1 pN/nm, then the complete pulling (if you pull up to 30 pN for example) will take around 6 sec, meaning 0.1 Hz. The force versus extension plot obtained by optical tweezer can be fit to the worm-like chain interpolation formula (Bustamante et al., Reference Bustamante, Marko, Siggia and Smith1994) which does not change much when the speed of pulling is changed. Interestingly the WLC fitting works well for dsDNA in the whole range of forces, but in the case of proteins the force versus extension plot deviates from the WLC fit at low forces (Bechtluft et al., Reference Bechtluft, Van Leeuwen, Tyreman, Tomkiewicz, Nouwen, Tepper, Driessen and Tans2007; Kaiser et al., Reference Kaiser, Goldman, Chodera, Tinoco and Bustamante2011; Wilson, Reference Wilson2011; Bianco et al., Reference Bianco, Mártonfalvi, Naftz, Koszegi and Kellermayer2015; He et al., Reference He, Li, Gao, Xiao, Hu, Hu, Hu and Li2019; Li et al., Reference Li, Chen, Guo, Wang and Li2021), probably due to localized unfolding/folding. Polypeptides may exhibit non-WLC behavior at lower forces if localized structures form during relaxation, such as beads on a string caused by hydrophobic collapse, or off-pathway folding intermediates which delay the final state.

Complementary, in Nanorheology where the frequency of the applied force is 10–200 Hz, a viscoelastic behavior is seen at lower forces which may correspond to local unfolding and help explain the deviation in the force versus extension plot obtained in optical tweezers. Currently, there is no comprehensive model for this viscoelastic transition (Zocchi, Reference Zocchi2018). Developing such a model would be important for a complete description of protein’s behavior under force.

Conclusion

In the framework described here, the enzyme is viewed as a cyclic engine taking many different conformations throughout the cycle. These conformational states produce a ‘rocky’ energy landscape as seen in Figure 1. The difference in the chemical potential of substrates and products drives the cycle in one direction. The cycle is initiated when the substrate binds. The force exerted by the substrate on the different parts of the enzyme then drives the ‘open to closed’ conformation leaving the protein in the closed stressed state. This step can be thought of as the substrate pulling the structure towards the closed conformation with a constant force which results in an elastic deformation of the molecule followed by a larger viscous deformation (a.k.a flow). In the energy landscape picture, this step can be interpreted as the jump between the two conformations, a flight from one region to another in Figure 1. This jump can be due to local unfolding as a mechanism to reduce strain that resulted from the force between the substrate and the enzyme. While in the closed state, internal stresses might change due to the change in the force between the substrates and the protein, however, the overall conformational change is very small (Delalande et al., Reference Delalande, Sacquin-Mora and Baaden2011). In the energy landscape, these small conformational changes are regarded as confined diffusion within one region. The next step in the cycle involves the reaction. During this step, significant strain energy is built which could be released by local unfolding/refolding. The products are then released and the internal restoring force of the enzyme’s structure brings it back to initial open state. Typically for motor proteins, this cycle is fast-slow: the time it takes the structure to close is shorter than the time it takes the structure to open (Zocchi, Reference Zocchi2018).

The chronology above is considering only one substrate for simplicity and in general, is only an approximate model. The real system consists of ~104 atoms and thus moves in a phase space consisting of the coordinates and the momenta of all the atoms. This model however provides a qualitative representation for the ensemble-averaged trajectory of the enzyme + substrate system during the enzymatic cycle. Simply put, throughout the cycle the strength of local interactions change as the enzyme deforms, and therefore different regions of the molecule could locally unfold (for exemplification of the cycle see kinesin in Figure 4).

Figure 4. Kinesin1 as molecular machine. (a) Kinesin1 is in initial position x0–x2. (b) After ATP hydrolysis a partial unfolding of one of its dimers allows it to advance to the x3 position. (c) Finally, the folding of the dimer drags the protein to a final position x1–x3. Kinesin1 use ATP hydrolysis as a source of energy to move along a track.

Considering the gap of information in the understanding of the role of forces in the mechanical–structural processes involved in enzyme catalysis, and that the ‘strain-induced’ idea remains somewhat hypothetical, studying the forces at a single molecule level could provide a novel and feasible approach. The in singulo studies carried out by optical tweezers, coupled or not to fluorescence and nanorheology, could be used to determine the energies that govern these processes. In this context, questions such as the following can be now answered: Is AK catalysis an example of the strain-induced theory? Does GK undergo a cracking mechanism? How does the applied force affect the turnover of AK, GK, and BiP? At what levels does the applied force help or hinder catalysis, and what forces will prohibit the enzyme from binding the substrate or achieving catalytic turnover? Is BiP unfolding correlated with the translocation step and speed?

This approach could provide valuable insights into the underlying mechanisms of enzyme catalysis and help to shed light on the role of forces in these processes.

Acknowledgments

We thank Steven B. Smith, Mauricio Baez, and Carlos Bustamante for helpful discussions and to all members of the Biochemistry Laboratory of the Universidad de Chile. Finally, we would like to thank the reviewers for their thoughtful comments.

Financial support

This work was supported by FONDECYT 1181361, PCI PII20150073, and Vicerrectoría de Investigación y Desarrollo (VID) of the Universidad de Chile ENL 10/22 (C.A.M.W.).

Competing interest

All authors have no competing interests to declare.

Footnotes

Z.A. and N.C.-M. contributed equally to this paper.

References

Abrusci, P, Chiarelli, LR, Galizzi, A, Fermo, E, Bianchi, P, Zanella, A and Valentini, G (2007) Erythrocyte adenylate kinase deficiency: Characterization of recombinant mutant forms and relationship with nonspherocytic hemolytic anemia. Experimental Hematology 35(8), 11821189.CrossRefGoogle ScholarPubMed
Adachi, K, Oiwa, K, Yoshida, M, Nishizaka, T and Kinosita, K (2012) Controlled rotation of the F1-ATPase reveals differential and continuous binding changes for ATP synthesis. Nature Communications 3(1), 1022.CrossRefGoogle ScholarPubMed
Ådén, J and Wolf-Watz, M (2007) NMR identification of transient complexes critical to adenylate kinase catalysis. Journal of the American Chemical Society 129(45), 1400314012.CrossRefGoogle ScholarPubMed
Alavi, Z (2017) Probing the Surface and the Interior of an Enzyme: What Is the Origin of Dissipation at the Angstrom Scale? Los Angeles, CA: University of California Los Angeles.Google Scholar
Alavi, Z, Ariyaratne, A and Zocchi, G (2015) Nano-rheology measurements reveal that the hydration layer of enzymes partially controls conformational dynamics. Applied Physics Letters 106(20), 36.CrossRefGoogle Scholar
Alfaro-Valdés, HM (2019) Determinación Del Mecanismo Mecano Químico de la Proteína BiP en El Proceso de Translocación in Multiplo. Master degree thesis, Universidad de Chile.Google Scholar
Alfaro-Valdés, HM, Burgos-Bravo, F, Casanova-Morales, N, Quiroga-Roger, D and Wilson, CAM (2018) Mechanical properties of chaperone BiP, the master regulator of the endoplasmic reticulum. In Endoplasmic Reticulum. London: IntechOpen.Google Scholar
Amyes, TL and Richard, JP (2013) Specificity in transition state binding: The Pauling model revisited. Biochemistry 52(12), 20212035.CrossRefGoogle ScholarPubMed
Anfinsen, CB, Redfield, RR, Choate, WL, Page, J and Carroll, WR (1954) Studies on the gross structure, cross-linkages, and terminal sequences in ribonuclease. The Journal of Biological Chemistry 207(1), 201210.CrossRefGoogle ScholarPubMed
Antikainen, NM and Martin, SF (2005) Altering protein specificity: Techniques and applications. Bioorganic & Medicinal Chemistry 13(8), 27012716.CrossRefGoogle Scholar
Assenza, S, De Los Rios, P and Barducci, A (2015) Quantifying the role of chaperones in protein translocation by computational modeling. Frontiers in Molecular Biosciences 2, 8. http://doi.org/10.3389/fmolb.2015.00008.CrossRefGoogle ScholarPubMed
Astumian, RD (1997) Thermodynamics and kinetics of a Brownian motor. Science 276(5314), 917922.CrossRefGoogle ScholarPubMed
Barnes, HA and Hutton, JF (1989) An Introduction to Rheology. Amsterdam: Elsevier Science.Google Scholar
Bauer, BW, Shemesh, T, Chen, Y and Rapoport, TA (2014) A “push and slide” mechanism allows sequence-insensitive translocation of secretory proteins by the SecA ATPase. Cell 157(6), 14161429.CrossRefGoogle Scholar
Bechtluft, P, Van Leeuwen, RGH, Tyreman, M, Tomkiewicz, D, Nouwen, N, Tepper, HL, Driessen, AJM and Tans, SJ (2007) Direct observation of chaperone-induced changes in a protein folding pathway. Science 318(5855), 14581461.CrossRefGoogle Scholar
Behnke, J, Feige, MJ and Hendershot, LM (2015) BiP and its nucleotide exchange factors Grp170 and Sil1: Mechanisms of action and biological functions. Journal of Molecular Biology 427(7), 15891608.CrossRefGoogle ScholarPubMed
Bell, GI (1978) Models for the specific adhesion of cells to cells. Science 200(4342), 618627.CrossRefGoogle ScholarPubMed
Ben Ishay, E, Rahamim, G, Orevi, T, Hazan, G, Amir, D and Haas, E (2012) Fast subdomain folding prior to the global refolding transition of E. Coli adenylate kinase: A double kinetics study. Journal of Molecular Biology 423(4), 613623.CrossRefGoogle Scholar
Benini, S, Cianci, M, Mazzei, L and Ciurli, S (2014) Fluoride inhibition of Sporosarcina pasteurii urease: Structure and thermodynamics. Journal of Biological Inorganic Chemistry 19(8), 12431261.CrossRefGoogle ScholarPubMed
Bertelsen, EB, Chang, L, Gestwicki, JE and Zuiderweg, ERP (2009) Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate. Proceedings of the National Academy of Sciences 106(21), 84718476.CrossRefGoogle ScholarPubMed
Bianco, P, Mártonfalvi, Z, Naftz, K, Koszegi, D and Kellermayer, M (2015) Titin domains progressively unfolded by force are homogenously distributed along the molecule. Biophysical Journal 109(2), 340345.CrossRefGoogle ScholarPubMed
Bianconi, ML (2003) Calorimetric determination of thermodynamic parameters of reaction reveals different Enthalpic compensations of the yeast hexokinase isozymes. Journal of Biological Chemistry 278(21), 1870918713.CrossRefGoogle ScholarPubMed
Bliumenfel’d, LA (1971) [Activation parameters of enzyme reactions (applicability of the activated complex theory in enzymology)]. Biofizika 16(4), 724727.Google ScholarPubMed
Braakman, I and Hebert, DN (2013) Protein folding in the endoplasmic reticulum. Cold Spring Harbor Perspectives in Biology 5(5), a013201.CrossRefGoogle ScholarPubMed
Bullerjahn, JT, Sturm, S and Kroy, K (2020) Non-Markov bond model for dynamic force spectroscopy. The Journal of Chemical Physics 152(6), 064104.CrossRefGoogle ScholarPubMed
Bustamante, C, Chemla, YR, Forde, NR and Izhaky, D (2004) Mechanical processes in biochemistry. Annual Review of Biochemistry 73(1), 705748.CrossRefGoogle ScholarPubMed
Bustamante, CJ, Kaiser, CM, Maillard, RA, Goldman, DH and Wilson, CAM (2014) Mechanisms of cellular proteostasis: Insights from single-molecule approaches. Annual Review of Biophysics 43, 119140.CrossRefGoogle ScholarPubMed
Bustamante, C, Marko, JF, Siggia, ED and Smith, S (1994) Entropic elasticity of λ-phage DNA. Science 265(5178), 15991600.CrossRefGoogle ScholarPubMed
Carrion-Vazquez, M, Oberhauser, AF, Fowler, SB, Marszalek, PE, Broedel, SE, Clarke, J and Fernandez, JM (1999) Mechanical and chemical unfolding of a single protein: A comparison. Proceedings of the National Academy of Sciences 96(7), 36943699.CrossRefGoogle ScholarPubMed
Casanova-Morales, N, Alavi, Z, Wilson, CAM and Zocchi, G (2018a) Identifying chaotropic and kosmotropic agents by nanorheology. The Journal of Physical Chemistry B 122(14), 37543759.CrossRefGoogle ScholarPubMed
Casanova-Morales, N, Quiroga-Roger, D, Alfaro-Valdés, HM, Alavi, Z, Lagos-Espinoza, MIA, Zocchi, G and Wilson, CAM (2018b) Mechanical properties of BiP protein determined by nano-rheology. Protein Science 27(8), 14181426.CrossRefGoogle ScholarPubMed
Cecconi, C, Shank, EA, Bustamante, C and Marqusee, S (2005) Direct observation of the three-state folding of a single protein molecule. Science 309(5743), 20572060.CrossRefGoogle ScholarPubMed
Chakafana, Z, Zininga, T and Shonhai, A (2019) The link that binds: The linker of Hsp70 as a helm of the protein’s function. Biomolecules 9(10), 543.CrossRefGoogle Scholar
Cheng, W, Dumont, S, Tinoco, I and Bustamante, C (2007) NS3 helicase actively separates RNA strands and senses sequence barriers ahead of the opening fork. Proceedings of the National Academy of Sciences 104(35), 1395413959.CrossRefGoogle ScholarPubMed
Choi, B and Zocchi, G (2007) Guanylate kinase, induced fit, and the allosteric spring probe. Biophysical Journal 92(5), 16511658.CrossRefGoogle ScholarPubMed
Cornish-Bowden, A (2012) Fundamentals of Enzyme Kinetics. New York: John Wiley & Sons.Google Scholar
Cossio, P, Hummer, G and Szabo, A (2015) On artifacts in single-molecule force spectroscopy. Proceedings of the National Academy of Sciences 112(46), 1424814253.CrossRefGoogle ScholarPubMed
Cossio, P, Hummer, G and Szabo, A (2016) Kinetic ductility and force-spike resistance of proteins from single-molecule force spectroscopy. Biophysical Journal 111(4), 832840.CrossRefGoogle ScholarPubMed
Cossio, P, Hummer, G and Szabo, A (2018) Transition paths in single-molecule force spectroscopy. The Journal of Chemical Physics 148, 123309. http://doi.org/10.1063/1.5004767.CrossRefGoogle ScholarPubMed
Creighton, TE (1990) Protein folding. Biochemical Journal 270(1), 116.CrossRefGoogle ScholarPubMed
De Los Rios, P, Ben-Zvi, A, Slutsky, O, Azem, A and Goloubinoff, P (2006) Hsp70 chaperones accelerate protein translocation and the unfolding of stable protein aggregates by entropic pulling. Proceedings of the National Academy of Sciences 103(16), 61666171.CrossRefGoogle ScholarPubMed
Delalande, O, Sacquin-Mora, S and Baaden, M (2011) Enzyme closure and nucleotide binding structurally lock guanylate kinase. Biophysical Journal 101(6), 14401449.CrossRefGoogle ScholarPubMed
Dill, KA (1985) Theory for the folding and stability of globular proteins. Biochemistry 24(6), 15011509.CrossRefGoogle ScholarPubMed
Dill, KA, Ozkan, SB, Shell, MS and Weikl, TR (2008) The protein folding problem. Annual Review of Biophysics 37(1), 289316.CrossRefGoogle ScholarPubMed
Dores-Silva, PR, Cauvi, DM, Coto, ALS, Kiraly, VTR, Borges, JC and De Maio, A (2020) Interaction of HSPA5 (Grp78, BIP) with negatively charged phospholipid membranes via oligomerization involving the N-terminal end domain. Cell Stress & Chaperones 25, 979. http://doi.org/10.1007/s1s2192-020-01134-9.CrossRefGoogle ScholarPubMed
Dudko, OK, Hummer, G and Szabo, A (2008) Theory, analysis, and interpretation of single-molecule force spectroscopy experiments. Proceedings of the National Academy of Sciences of the United States of America 105(41), 1575515760.CrossRefGoogle ScholarPubMed
Dzeja, PP, Bast, P, Pucar, D, Wieringa, B and Terzic, A (2007) Defective metabolic Signaling in adenylate kinase AK1 gene knock-out hearts compromises post-ischemic coronary reflow. Journal of Biological Chemistry 282(43), 3136631372.CrossRefGoogle ScholarPubMed
Edelstein, SJ (1975) Cooperative interactions of Hemoglobin. Annual Review of Biochemistry 44(1), 209232.CrossRefGoogle ScholarPubMed
Eisenmesser, EZ, Bosco, DA, Akke, M and Kern, D (2002) Enzyme dynamics during catalysis. Science 295(5559), 15201523.CrossRefGoogle ScholarPubMed
Fenimore, PW, Frauenfelder, H, Magazù, S, McMahon, BH, Mezei, F, Migliardo, F, Young, RD and Stroe, I (2013) Concepts and problems in protein dynamics. Chemical Physics 424, 26.CrossRefGoogle Scholar
Frauenfelder, H, Parak, F and Young, RD (1988) Conformational substates in proteins. Annual Review of Biophysics and Biophysical Chemistry 17(1), 451479.CrossRefGoogle ScholarPubMed
Gao, Y, Zorman, S, Gundersen, G, Xi, Z, Sirinakis, G, Rothman, JE and Zhang, Y (2012) Single reconstituted neuronal SNARE complexes zipper in three distinct stages. Science 337(6100), 13401343.CrossRefGoogle ScholarPubMed
Goldman, DH, Kaiser, CM, Milin, A, Righini, M, Tinoco, I and Bustamante, C (2015) Mechanical force releases nascent chain–mediated ribosome arrest in vitro and in vivo. Science 348(6233), 457460.CrossRefGoogle ScholarPubMed
Haldane, JBS (1965) Enzymes. Cambridge, MA: MIT Press.Google Scholar
Hamada, M, Sumida, M, Kurokawa, Y, Sunayashiki-Kusuzaki, K, Okuda, H, Watanabe, T and Kuby, SA (1985) Studies on the adenylate kinase isozymes from the serum and erythrocyte of normal and Duchenne dystrophic patients. Isolation, physicochemical properties, and several comparisons with the Duchenne dystrophic aberrant enzyme. The Journal of Biological Chemistry 260(21), 1159511602.CrossRefGoogle ScholarPubMed
Hanson, C, Nishiyama, Y and Paul, S (2005) Catalytic antibodies and their applications. Current Opinion in Biotechnology 16(6), 631636.CrossRefGoogle ScholarPubMed
He, C, Li, S, Gao, X, Xiao, A, Hu, C, Hu, X, Hu, X and Li, H (2019) Direct observation of the fast and robust folding of a slipknotted protein by optical tweezers. Nanoscale 11(9), 39453951.CrossRefGoogle ScholarPubMed
Hendershot, LM (2004) The ER function BiP is a master regulator of ER function. Mount Sinai Journal of Medicine 71(5), 289297.Google ScholarPubMed
Henzler-Wildman, KA, Thai, V, Lei, M, Ott, M, Wolf-Watz, M, Fenn, T, Pozharski, E, Wilson, MA, Petsko, GA, Karplus, M, Hübner, CG, and Kern, D (2007) Intrinsic motions along an enzymatic reaction trajectory. Nature 450(7171), 838844.CrossRefGoogle ScholarPubMed
Hilser, VJ, García-Moreno, EB, Oas, TG, Kapp, G and Whitten, ST (2006) A statistical thermodynamic model of the protein ensemble. Chemical Reviews 106(5), 15451558.CrossRefGoogle ScholarPubMed
Hokenson, MJ, Cope, GA, Lewis, ER, Oberg, KA and Fink, AL (2000) Enzyme-induced strain/distortion in the ground-state ES complex in β-lactamase catalysis revealed by FTIR. Biochemistry 39(21), 65386545.CrossRefGoogle ScholarPubMed
Houdusse, A and Sweeney, HL (2016) How myosin generates force on actin filaments. Trends in Biochemical Sciences 41(12), 989997.CrossRefGoogle ScholarPubMed
Ibrahim, IM, Abdelmalek, DH, Elshahat, ME and Elfiky, AA (2020) COVID-19 spike-host cell receptor GRP78 binding site prediction. Journal of Infection 80(5), 554562.CrossRefGoogle ScholarPubMed
Ichhaporia, VP and Hendershot, LM (2021) Role of the hsp70 co-chaperone sil1 in health and disease. International Journal of Molecular Sciences 22(4), 123.CrossRefGoogle ScholarPubMed
Jan, Y-H, Tsai, H-Y, Yang, C-J, Huang, M-S, Yang, Y-F, Lai, T-C, Lee, C-H, Jen, Y-M, Huang, C-Y, Su, J-L, Chuang, y-J and Hsiao, M (2012) Adenylate Kinase-4 is a marker of poor clinical outcomes that promotes metastasis of lung cancer by downregulating the transcription factor ATF3. Cancer Research 72(19), 51195129.CrossRefGoogle ScholarPubMed
Jencks, WP (1975) Binding energy, specificity, and enzymic catalysis: The circe effect. Advances in Enzymology and Related Areas of Molecular Biology 43, 219410.Google ScholarPubMed
Jumper, J, Evans, R, Pritzel, A, Green, T, Figurnov, M, Ronneberger, O, Tunyasuvunakool, K, Bates, R, Žídek, A, Potapenko, A, Bridgland, A, Meyer, C, Kohl, SAA, Ballard, AJ, Cowie, A, Romera-Paredes, B, Nikolov, S, Jain, R, Adler, J and Hassabis, D (2021) Highly accurate protein structure prediction with AlphaFold. Nature 596(7873), 583589.CrossRefGoogle ScholarPubMed
Kaiser, CM, Goldman, DH, Chodera, JD, Tinoco, I and Bustamante, C (2011) The ribosome modulates nascent protein folding. Science 334(6063), 17231727.CrossRefGoogle ScholarPubMed
Kemp, RB and Guan, YH (1999) Microcalorimetric studies of animal tissues and their isolated cells. Handbook of Thermal Analysis and Calorimetry 4, 557656.CrossRefGoogle Scholar
Keramisanou, D, Biris, N, Gelis, I, Sianidis, G, Karamanou, S, Economou, A and Kalodimos, CG (2006) Disorder-order folding transitions underlie catalysis in the helicase motor of SecA. Nature Structural & Molecular Biology 13(7), 594602.CrossRefGoogle ScholarPubMed
Khan, N, Shah, PP, Ban, D, Trigo-Mouriño, P, Carneiro, MG, DeLeeuw, L, Dean, WL, Trent, JO, Beverly, LJ, Konrad, M, Lee, D and Sabo, TM (2019) Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer. Journal of Biological Chemistry 294(31), 1192011933.CrossRefGoogle Scholar
Klinman, J and Hammes-Schiffer, S (2013) Dynamics in Enzyme Catalysis. Berlin: Springer.CrossRefGoogle ScholarPubMed
Kosakowska-Cholody, T, Lin, J, Srideshikan, SM, Scheffer, L, Tarasova, NI and Acharya, JK (2014) HKH40A downregulates GRP78/BiP expression in cancer cells. Cell Death & Disease 5(5), e1240.CrossRefGoogle ScholarPubMed
Koshland, DE (1958) Application of a theory of enzyme specificity to protein synthesis. Proceedings of the National Academy of Sciences 44(2), 98104.CrossRefGoogle ScholarPubMed
Kramers, HA (1940) Brownian motion in a field of force and the diffusion model of chemical reactions. Physica 7(4), 284304.CrossRefGoogle Scholar
Kubelka, J, Hofrichter, J and Eaton, WA (2004) The protein folding ‘speed limit.’. Current Opinion in Structural Biology 14(1), 7688.CrossRefGoogle ScholarPubMed
Li, J, Chen, G, Guo, Y, Wang, H and Li, H (2021) Single molecule force spectroscopy reveals the context dependent folding pathway of the C-terminal fragment of Top7. Chemical Science 12(8), 28762884.CrossRefGoogle Scholar
Lienhard, GE, and Secemski, II (1973) P 1, P 5 -Di(adenosine-5′)pentaphosphate, a potent multisubstrate inhibitor of adenylate kinase. The Journal of Biological Chemistry 248(3), 11211123.CrossRefGoogle ScholarPubMed
Magliery, TJ (2015) Protein stability: Computation, sequence statistics, and new experimental methods. Current Opinion in Structural Biology 33, 161168.CrossRefGoogle ScholarPubMed
Maillard, RA, Chistol, G, Sen, M, Righini, M, Tan, J, Kaiser, CM, Hodges, C, Martin, A and Bustamante, C (2011) ClpX(P) generates mechanical force to unfold and translocate its protein substrates. Cell 145(3), 459469.CrossRefGoogle ScholarPubMed
Malagrinò, F, Visconti, L, Pagano, L, Toto, A, Troilo, F and Gianni, S (2020) Understanding the binding induced folding of intrinsically disordered proteins by protein engineering: Caveats and pitfalls. International Journal of Molecular Sciences 21(10), 3484.CrossRefGoogle ScholarPubMed
Mapa, K, Sikor, M, Kudryavtsev, V, Waegemann, K, Kalinin, S, Seidel, CAM, Neupert, W, Lamb, DC and Mokranjac, D (2010) The conformational dynamics of the mitochondrial Hsp70 chaperone. Molecular Cell 38(1), 89100.CrossRefGoogle ScholarPubMed
Marcinowski, M, Höller, M, Feige, MJ, Baerend, D, Lamb, DC and Buchner, J (2011) Substrate discrimination of the chaperone BiP by autonomous and cochaperone-regulated conformational transitions. Nature Structural & Molecular Biology 18(2), 150158.CrossRefGoogle ScholarPubMed
Matlack, KE, Misselwitz, B, Plath, K and Rapoport, TA (1999) BiP acts as a molecular ratchet during posttranslational transport of Prepro-α factor across the ER membrane. Cell 97(5), 553564.CrossRefGoogle ScholarPubMed
Miyashita, O, Onuchic, JN and Wolynes, PG (2003) Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins. Proceedings of the National Academy of Sciences 100(22), 1257012575.CrossRefGoogle ScholarPubMed
Moffitt, JR, Chemla, YR, Aathavan, K, Grimes, S, Jardine, PJ, Anderson, DL, and Bustamante, C (2009) Intersubunit coordination in a homomeric ring ATPase. Nature 457(7228), 446450.CrossRefGoogle Scholar
Müller, CW and Schulz, GE (1992) Structure of the complex between adenylate kinase from Escherichia coli and the inhibitor Ap5A refined at 1.9 Å resolution. Journal of Molecular Biology 224(1), 159177.CrossRefGoogle ScholarPubMed
Neuman, KC, Abbondanzieri, EA, Landick, R, Gelles, J and Block, SM (2003) Ubiquitous transcriptional pausing is independent of RNA polymerase backtracking. Cell 115(4), 437447.CrossRefGoogle ScholarPubMed
Nogales, E (2016) The development of cryo-EM into a mainstream structural biology technique. Nature Methods 13(1), 2427.CrossRefGoogle ScholarPubMed
Oguchi, Y, Mikhailenko, SV, Ohki, T, Olivares, AO, De La Cruz, EM and Ishiwata, S (2008) Load-dependent ADP binding to myosins V and VI: Implications for subunit coordination and function. Proceedings of the National Academy of Sciences 105(22), 77147719.CrossRefGoogle ScholarPubMed
Oldfield, CJ and Dunker, AK (2014) Intrinsically disordered proteins and intrinsically disordered protein regions. Annual Review of Biochemistry 83(1), 553584.CrossRefGoogle ScholarPubMed
Olsen, SN (2006) Applications of isothermal titration calorimetry to measure enzyme kinetics and activity in complex solutions. Thermochimica Acta 448(1), 1218.CrossRefGoogle Scholar
Olsson, U and Wolf-Watz, M (2010) Overlap between folding and functional energy landscapes for adenylate kinase conformational change. Nature Communications 1(1), 111.CrossRefGoogle ScholarPubMed
Otero, JH, Lizák, B and Hendershot, LM (2010) Life and death of a BiP substrate. Seminars in Cell & Developmental Biology 21(5), 472478.CrossRefGoogle ScholarPubMed
Park, KW, Eun Kim, G, Morales, R, Moda, F, Moreno-Gonzalez, I, Concha-Marambio, L, Lee, AS, Hetz, C and Soto, C (2017) The endoplasmic reticulum chaperone GRP78/BiP modulates prion propagation in vitro and in vivo. Scientific Reports 7, 113.Google ScholarPubMed
Park, E and Rapoport, TA (2012) Mechanisms of Sec61/SecY-mediated protein translocation across membranes. Annual Review of Biophysics 41(1), 2140.CrossRefGoogle ScholarPubMed
Pauling, L (1946) Molecular architecture and biological reactions. Chemical and Engineering News 24(10), 13751377.CrossRefGoogle Scholar
Pelz, B, Žoldák, G, Zeller, F, Zacharias, M and Rief, M (2016) Subnanometre enzyme mechanics probed by single-molecule force spectroscopy. Nature Communications 7(1), 10848.CrossRefGoogle ScholarPubMed
Pickard, A, Chang, J, Alachkar, N, Calverley, B, Garva, R, Arvan, P, Meng, Q-J, and Kadler, KE (2019) Preservation of circadian rhythms by the protein folding chaperone, BiP. The FASEB Journal 33(6), 74797489.CrossRefGoogle ScholarPubMed
Pobre, KFR, Poet, GJ and Hendershot, LM (2019) The endoplasmic reticulum (ER) chaperone BiP is a master regulator of ER functions: Getting by with a little help from ERdj friends. Journal of Biological Chemistry 294(6), 20982108.CrossRefGoogle ScholarPubMed
Pollack, SJ, Jacobs, JW and Schultz, PG (1986) Selective chemical catalysis by an antibody. Science 234(4783), 15701573.CrossRefGoogle ScholarPubMed
Qu, H, Landy, J and Zocchi, G (2012) Cracking phase diagram for the dynamics of an enzyme. Physical Review. E, Statistical, Nonlinear, and Soft Matter Physics 86(4), 15.Google ScholarPubMed
Qu, H and Zocchi, G (2013) How enzymes work: A look through the perspective of molecular viscoelastic properties. Physical Review X 3(1), 11009.CrossRefGoogle Scholar
Ramírez, MP, Rivera, M, Quiroga-Roger, D, Bustamante, A, Vega, M, Baez, M, Puchner, EM and Wilson, CAM (2017) Single molecule force spectroscopy reveals the effect of BiP chaperone on protein folding. Protein Science 26(7), 14041412.CrossRefGoogle ScholarPubMed
Riedel, C, Gabizon, R, Wilson, CAM, Hamadani, K, Tsekouras, K, Marqusee, S, Pressé, S and Bustamante, C (2015) The heat released during catalytic turnover enhances the diffusion of an enzyme. Nature 517(7533), 227230.CrossRefGoogle ScholarPubMed
Rivera, M, Burgos-Bravo, F, Engelberger, F, Asor, R, Lagos-Espinoza, MIA, Figueroa, M, Kukura, P, Ramírez-Sarmiento, CA, Baez, M, Smith, SB and Wilson, CAM (2023a) Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate. Protein Science 32(7), e4706. http://doi.org/10.1002/pro.4706.CrossRefGoogle ScholarPubMed
Rivera, M, Mjaavatten, A, Smith, SB, Baez, M and Wilson, CAM (2023b) Temperature dependent mechanical unfolding and refolding of a protein studied by thermo-regulated optical tweezers. Biophysical Journal 122(3), 513521.CrossRefGoogle ScholarPubMed
Sánchez, WN, Robeson, L, Carrasco, V, Figueroa, NL, Burgos-Bravo, F, Wilson, CAM and Casanova-Morales, N (2022) Determination of protein–protein interactions at the single-molecule level using optical tweezers. Quarterly Reviews of Biophysics 55, e8.CrossRefGoogle ScholarPubMed
Savinov, SN, Hirschmann, R, Benkovic, SJ and Smith, AB (2003) Investigation of an antibody-ligase. Evidence for strain-induced catalysis. Bioorganic & Medicinal Chemistry Letters 13(7), 13211324.CrossRefGoogle ScholarPubMed
Schekman, R (1994) Translocation gets a push. Cell 78(6), 911913.CrossRefGoogle ScholarPubMed
Schrank, TP, Bolen, DW and Hilser, VJ (2009) Rational modulation of conformational fluctuations in adenylate kinase reveals a local unfolding mechanism for allostery and functional adaptation in proteins. Proceedings of the National Academy of Sciences 106(40), 1698416989.CrossRefGoogle ScholarPubMed
Schrank, TP, Wrabl, JO and Hilser, VJ (2013) Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: Function follows fluctuations. Topics in Current Chemistry 337, 95121.CrossRefGoogle ScholarPubMed
Schulz, GE, Müller, CW and Diederichs, K (1990) Induced-fit movements in adenylate kinases. Journal of Molecular Biology 213(4), 627630.CrossRefGoogle ScholarPubMed
Sellers, JR and Veigel, C (2010) Direct observation of the myosin-Va power stroke and its reversal. Nature Structural & Molecular Biology 17(5), 590595.CrossRefGoogle ScholarPubMed
Shea, J-E and Brooks, CL (2001) From folding theories to folding proteins: A review and assessment of simulation studies of protein folding and unfolding. Annual Review of Physical Chemistry 52(1), 499535.CrossRefGoogle ScholarPubMed
Shields, AM, Panayi, GS and Corrigall, VM (2012) A new-age for biologic therapies: Long-term drug-free therapy with BiP? Frontiers in Immunology 3, 17. http://doi.org/10.3389/fimmu.2012.00017.CrossRefGoogle ScholarPubMed
Shin, J, Toyoda, S, Nishitani, S, Fukuhara, A, Kita, S, Otsuki, M, and Shimomura, I (2021) Possible involvement of adipose tissue in patients with older age, obesity, and diabetes with SARS-CoV-2 infection (COVID-19) via GRP78 (BIP/HSPA5): Significance of hyperinsulinemia management in COVID-19. Diabetes 70(12), 27452755.CrossRefGoogle ScholarPubMed
Sinev, MA, Sineva, EV, Ittah, V and Haas, E (1996) Domain closure in adenylate kinase. Biochemistry 35(20), 64256437.CrossRefGoogle ScholarPubMed
Singh, A, Rief, M and Žoldák, G (2022) Direct observation of chemo-mechanical coupling in DnaK by single-molecule force experiments. Biophysical Journal 121(23), 47294739.CrossRefGoogle ScholarPubMed
Sinha, AK (1972) Colorimetric assay of catalase. Analytical Biochemistry 47(2), 389394.CrossRefGoogle ScholarPubMed
Smith, SB, Cui, Y and Bustamante, C (2003) Optical-trap force transducer that operates by direct measurement of light momentum. Methods in Enzymology 361, 134162.CrossRefGoogle ScholarPubMed
Smith, DE, Tans, SJ, Smith, SB, Grimes, S, Anderson, DL and Bustamante, C (2001) The bacteriophage φ29 portal motor can package DNA against a large internal force. Nature 413(6857), 748752.CrossRefGoogle ScholarPubMed
Stehle, T and Schulz, GE (1992) Refined structure of the complex between guanylate kinase and its substrate GMP at 2·0 Å resolution. Journal of Molecular Biology 224(4), 11271141.CrossRefGoogle ScholarPubMed
Swain, JF, Dinler, G, Sivendran, R, Montgomery, DL, Stotz, M and Gierasch, LM (2007) Hsp70 chaperone ligands control domain association via an allosteric mechanism mediated by the Interdomain linker. Molecular Cell 26(1), 2739.CrossRefGoogle ScholarPubMed
Taverna, DM and Goldstein, RA (2002) Why are proteins marginally stable? Proteins: Structure, Function, and Genetics 46(1), 105109.CrossRefGoogle ScholarPubMed
Tinoco, I, Sauer, K, Wang, JC, Puglisi, JD, Harbison, G and Rovnyak, D (2013) Physical Chemistry: Principles and Applications in Biological Sciences.Google Scholar
Tomasello, G, Armenia, I and Molla, G (2020) The protein imager: A full-featured online molecular viewer interface with server-side HQ-rendering capabilities. Bioinformatics 36(9), 29092913.CrossRefGoogle ScholarPubMed
Tramontano, A, Janda, KD and Lerner, RA (1986) Catalytic antibodies. Science 234(4783), 15661570.CrossRefGoogle ScholarPubMed
Voet, D and Voet, J (2010) Biochemistry. New York: John Wiley & Sons.Google Scholar
Volkán-Kacsó, S and Marcus, RA (2015) Theory for rates, equilibrium constants, and Brønsted slopes in F 1 -ATPase single molecule imaging experiments. Proceedings of the National Academy of Sciences 112(46), 1423014235.CrossRefGoogle ScholarPubMed
Wang, MD, Schnitzer, MJ, Yin, H, Landick, R, Gelles, J and Block, SM (1998) Force and velocity measured for single molecules of RNA polymerase. Science 282(5390), 902907.CrossRefGoogle ScholarPubMed
Wang, Y and Zocchi, G (2010) Elasticity of globular proteins measured from the ac susceptibility. Physical Review Letters 105(23), 238104.CrossRefGoogle ScholarPubMed
Watanabe, R, Okuno, D, Sakakihara, S, Shimabukuro, K, Lino, R, Yoshida, M and Noji, H (2012) Mechanical modulation of catalytic power on F1-ATPase. Nature Chemical Biology 8(1), 8692.CrossRefGoogle Scholar
Whitford, PC, Miyashita, O, Levy, Y and Onuchic, JN (2007) Conformational transitions of adenylate kinase: Switching by cracking. Journal of Molecular Biology 366(5), 16611671.CrossRefGoogle ScholarPubMed
Wilson, CAM (2011) Single Molecule Studies by Optical Tweezers: Folding and Unfolding of Glucokinase from Thermococcus litoralis. PhD thesis, Universidad de Chile.Google Scholar
Wilson, CAM, Leachman, SM, Cervantes, B, Ierokomos, A, Marqusee, S and Bustamante, C (2013) A single-molecule fluorescence system for studying adenylate kinase under force. The FASEB Journal 27(S1), 996. http://doi.org/10.1096/fasebj.27.1_supplement.996.10.CrossRefGoogle Scholar
Wolf-Watz, M, Thai, V, Henzler-Wildman, K, Hadjipavlou, G, Eisenmesser, EZ and Kern, D (2004) Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair. Nature Structural & Molecular Biology 11(10), 945949.CrossRefGoogle Scholar
Yang, J, Nune, M, Zong, Y, Zhou, L and Liu, Q (2015) Close and allosteric opening of the polypeptide-binding site in a human Hsp70 chaperone BiP. Structure 23(12), 21912203.CrossRefGoogle Scholar
Yokoyama, T, Fujii, S, Ostermann, A, Schrader, TE, Nabeshima, Y and Mizuguchi, M (2022) Neutron crystallographic analysis of the nucleotide-binding domain of Hsp72 in complex with ADP. IUCrJ 9(5), 562572.CrossRefGoogle ScholarPubMed
Zhang, X, Halvorsen, K, Zhang, C-Z, Wong, WP and Springer, TA (2009) Mechanoenzymatic cleavage of the Ultralarge vascular protein von Willebrand factor. Science 324(5932), 13301334.CrossRefGoogle ScholarPubMed
Zhang, X, Rebane, AA, Ma, L, Ma, L, Li, F, Jiao, J, Qu, H, Pincet, F, Rothman, JE and Zhang, Y (2016) Stability, folding dynamics, and long-range conformational transition of the synaptic t-SNARE complex. Proceedings of the National Academy of Sciences 113(50), E8031. http://doi.org/10.1073/pnas.1605748113.CrossRefGoogle ScholarPubMed
Zhu, X, Zhao, X, Burkholder, WF, Gragerov, A, Ogata, CM, Gottesman, ME and Hendrickson, WA (1996) Structural analysis of substrate binding by the molecular chaperone DnaK. Science 272(5268), 16061614.CrossRefGoogle ScholarPubMed
Zimmermann, R, Eyrisch, S, Ahmad, M and Helms, V (2011) Protein translocation across the ER membrane. Biochimica et Biophysica Acta (BBA) – Biomembrane 1808(3), 912924.CrossRefGoogle ScholarPubMed
Zocchi, G (2018) Molecular Machines: A Materials Science Approach. Princeton, NJ: Princeton University Press.Google Scholar
Figure 0

Figure 1. The many different conformational states produce a rocky landscape.

Figure 1

Table 1. Theoretical and experimental characterization of enthalpies

Figure 2

Figure 2. Single-molecule force spectroscopy techniques (open and closed conformations of the proteins to be studied). (a) Left: Open apo Aquifex AK (PDBid: 2RH5). Middle: Local unfolded or crack intermediate. Right: Close Aquifex AK in complex with the substrate analogue Zn2 + ● Ap5A (light blue) (PDBid: 2RGX) (Olsson and Wolf-Watz, 2010). In yellow is the lid domain and in gray is the core of the protein. Reprinted by permission from Springer Nature. (b) Left: Open configuration GK (PDBid: 1ZNW). Right: Close configuration of GK (PDBid: 1LVG) upon substrate binding. Alpha helices are colored magenta and beta sheets are colored yellow. (Alavi, 2017). (c) Open configuration BiP in the presence of ADP (PDBids: 5EVZ and 5E85) and right. Close configuration of BiP in the presence of ATP (PDBid: 5E84). In blue is the NBD and orange is SBD (Pobre et al.,2019). Images were made with 3D Protein imaging software (Tomasello et al.,2020).

Figure 3

Figure 3. Nanorheology setup, showing the flow chamber with enzyme-tethered GNPs (enzyme shown in green and GNPs shown as golden spheres), the parallel plates capacitor geometry used for mechanical excitation, and the evanescent wave scattering optics used for readout.

Figure 4

Figure 4. Kinesin1 as molecular machine. (a) Kinesin1 is in initial position x0–x2. (b) After ATP hydrolysis a partial unfolding of one of its dimers allows it to advance to the x3 position. (c) Finally, the folding of the dimer drags the protein to a final position x1–x3. Kinesin1 use ATP hydrolysis as a source of energy to move along a track.