An official website of the United States government
Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane–protein interactions had only moderate influences, except for the extraction of the last helix where the membrane–protein interactions had a stronger influence. The up–down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.
more »
« less
Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification
Protein–protein and protein–water hydrogen bonding interactions play essential roles in the way a protein passes through the transition state during folding or unfolding, but the large number of these interactions in molecular dynamics (MD) simulations makes them difficult to analyze. Here, we introduce a state space representation and associated “rarity” measure to identify and quantify transition state passage (transit) events. Applying this representation to a long MD simulation trajectory that captured multiple folding and unfolding events of the GTT WW domain, a small protein often used as a model for the folding process, we identified three transition categories: Highway (faster), Meander (slower), and Ambiguous (intermediate). We developed data sonification and visualization tools to analyze hydrogen bond dynamics before, during, and after these transition events. By means of these tools, we were able to identify characteristic hydrogen bonding patterns associated with “Highway” versus “Meander” versus “Ambiguous” transitions and to design algorithms that can identify these same folding pathways and critical protein–water interactions directly from the data. Highly cooperative hydrogen bonding can either slow down or speed up transit. Furthermore, an analysis of protein–water hydrogen bond dynamics at the surface of WW domain shows an increase in hydrogen bond lifetime from folded to unfolded conformations with Ambiguous transitions as an outlier. In summary, hydrogen bond dynamics provide a direct window into the heterogeneity of transits, which can vary widely in duration (by a factor of 10) due to a complex energy landscape.
more »
« less
- Award ID(s):
- 2205665
- PAR ID:
- 10575350
- Publisher / Repository:
- NSF DAR
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 121
- Issue:
- 22
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Single-molecule force spectroscopy is a powerful tool for studying protein folding. Over the last decade, a key question has emerged: how are changes in intrinsic biomolecular dynamics altered by attachment to μm-scale force probes via flexible linkers? Here, we studied the folding/unfolding of α3D using atomic force microscopy (AFM)–based force spectroscopy. α3D offers an unusual opportunity as a prior single-molecule fluorescence resonance energy transfer (smFRET) study showed α3D’s configurational diffusion constant within the context of Kramers theory varies with pH. The resulting pH dependence provides a test for AFM-based force spectroscopy’s ability to track intrinsic changes in protein folding dynamics. Experimentally, however, α3D is challenging. It unfolds at low force (<15 pN) and exhibits fast-folding kinetics. We therefore used focused ion beam–modified cantilevers that combine exceptional force precision, stability, and temporal resolution to detect state occupancies as brief as 1 ms. Notably, equilibrium and nonequilibrium force spectroscopy data recapitulated the pH dependence measured using smFRET, despite differences in destabilization mechanism. We reconstructed a one-dimensional free-energy landscape from dynamic data via an inverse Weierstrass transform. At both neutral and low pH, the resulting constant-force landscapes showed minimal differences (∼0.2 to 0.5kBT) in transition state height. These landscapes were essentially equal to the predicted entropic barrier and symmetric. In contrast, force-dependent rates showed that the distance to the unfolding transition state increased as pH decreased and thereby contributed to the accelerated kinetics at low pH. More broadly, this precise characterization of a fast-folding, mechanically labile protein enables future AFM-based studies of subtle transitions in mechanoresponsive proteins.more » « less
-
Abstract The forces that stabilize membrane proteins remain elusive to precise quantification. Particularly important, but poorly resolved, are the forces present during the initial unfolding of a membrane protein, where the most native set of interactions is present. A high‐precision, atomic force microscopy assay was developed to study the initial unfolding of bacteriorhodopsin. A rapid near‐equilibrium folding between the first three unfolding states was discovered, the two transitions corresponded to the unfolding of five and three amino acids, respectively, when using a cantilever optimized for 2 μs resolution. The third of these states was retinal‐stabilized and previously undetected, despite being the most mechanically stable state in the whole unfolding pathway, supporting 150 pN for more than 1 min. This ability to measure the dynamics of the initial unfolding of bacteriorhodopsin provides a platform for quantifying the energetics of membrane proteins under native‐like conditions.more » « less
-
Abstract The influenza B M2 protein forms a water-filled tetrameric channel to conduct protons across the lipid membrane. To understand how channel water mediates proton transport, we have investigated the water orientation and dynamics using solid-state NMR spectroscopy and molecular dynamics (MD) simulations. 13 C-detected water 1 H NMR relaxation times indicate that water has faster rotational motion in the low-pH open channel than in the high-pH closed channel. Despite this faster dynamics, the open-channel water shows higher orientational order, as manifested by larger motionally-averaged 1 H chemical shift anisotropies. MD simulations indicate that this order is induced by the cationic proton-selective histidine at low pH. Furthermore, the water network has fewer hydrogen-bonding bottlenecks in the open state than in the closed state. Thus, faster dynamics and higher orientational order of water molecules in the open channel establish the water network structure that is necessary for proton hopping.more » « less
-
Proteins are constantly undergoing folding and unfolding transitions, with rates that determine their homeostasis in vivo and modulate their biological function. The ability to optimize these rates without affecting overall native stability is hence highly desirable for protein engineering and design. The great challenge is, however, that mutations generally affect folding and unfolding rates with inversely complementary fractions of the net free energy change they inflict on the native state. Here we address this challenge by targeting the folding transition state (FTS) of chymotrypsin inhibitor 2 (CI2), a very slow and stable two‐state folding protein with an FTS known to be refractory to change by mutation. We first discovered that the CI2's FTS is energetically taxed by the desolvation of several, highly conserved, charges that form a buried salt bridge network in the native structure. Based on these findings, we designed a CI2 variant that bears just four mutations and aims to selectively stabilize the FTS. This variant has >250‐fold faster rates in both directions and hence identical native stability, demonstrating the success of our FTS‐centric design strategy. With an optimized FTS, CI2 also becomes 250‐fold more sensitive to proteolytic degradation by its natural substrate chymotrypsin, and completely loses its activity as inhibitor. These results indicate that CI2 has been selected through evolution to have a very unstable FTS in order to attain the kinetic stability needed to effectively function as protease inhibitor. Moreover, the CI2 case showcases that protein (un)folding rates can critically pivot around a few key residues‐interactions, which can strongly modify the general effects of known structural factors such as domain size and fold topology. From a practical standpoint, our results suggest that future efforts should perhaps focus on identifying such critical residues‐interactions in proteins as best strategy to significantly improve our ability to predict and engineer protein (un)folding rates.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2319094121
