An official website of the United States government
This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(Nisopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.
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 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
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2319094121
