WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

Cross-correlations in the fluctuation-dissipation relation influence barrier-crossing dynamics

The Generalized Langevin Equation has been successfully used to model and understand the conformational dynamics of molecules in solution. However, recent works have demonstrated that, in these kinds of applications, the usual fluctuation-dissipation relation connecting the statistics of the random force to the memory kernel could contain a cross-correlation term. In this work, we systematically explore the origins of this cross-correlation term and argue that it plays a role, particularly in the folding dynamics of biopolymers. Finally, we propose an approximation for the cross-correlation term within the usual fluctuation-dissipation relation.

Ion-Ion Structural Correlation and Dynamics of Water in Aqueous NaCl Solutions with a Wide Range of Concentrations

The behavior of water in concentrated ionic solutions, including supersaturated conditions, is crucial for numerous material and energy conversion processes and fundamental research. All electrolytes whether they "structure-make" or "structure-break" the water structure lead to slower water motion. This study investigates the structure and dynamics of aqueous NaCl solutions across a wide range of concentrations. On the structural side, the primary focus is on ion-ion correlations. In terms of dynamics, we demonstrate that the slowing down of water dynamics continues even beyond the saturated state. We identify three distinct types of dynamics at large concentrations: ballistic, trapped, and diffusive. The van Hove correlation function exhibits no signs of relaxation within a time interval where particle motion is effectively halted. The system displays dynamical heterogeneities, confirmed by evaluating non-Gaussian parameters for the self-part of the van Hove function and identifying the mobile particles. These particles form clusters, with the largest sizes occurring when the non-Gaussian parameters are at their maximum. Additionally, we discuss the relaxation times associated with these systems using the incoherent intermediate scattering function and establish a connection with the mode-coupling theory.

Dynamic protein hydration water mediates the aggregation kinetics of amyloid β peptides at interfaces

Protein hydration water is essential for protein misfolding and amyloid formation, but how it directs the course of amyloid formation has yet to be elucidated. Here, we experimentally demonstrated that femtosecond sum frequency generation vibrational spectroscopy (SFG-VS) and the femtosecond IR pump-SFG probe technique can serve as powerful tools for addressing this issue. Using amyloid β(1-42) peptide as a model, we determined the transient misfolding intermediates by probing the amide band spectral features and the local hydration water changes by measuring the ultrafast vibrational dynamics of the amide I band. For the first time, we established a correlation between the dynamic change in protein hydration water and aggregation propensity. The aggregation propensity depends on the dynamic change in the hydration water, rather than the static hydration water content of the initial protein state. Water expulsion enhances the aggregation propensity and promotes amyloid formation, while protein hydration attenuates the aggregation propensity and inhibits amyloid formation. The suppression of water expulsion and protein hydration can prevent protein aggregation and stabilize proteins. These findings contribute to a better understanding of the underlying effect of hydration water on amyloid formation and protein structural stability and provide a strategy for maintaining long-term stabilization of biomolecules.

Impact of pH-dependent dynamics of human serum proteins on dialysis membranes: Cryptographic structure assessment, synchrotron imaging of membrane-protein adsorption, and molecular docking studies

Proteins are fundamental to biochemical processes and critical in hemodialysis. This study investigates the impact of pH on human serum albumin (HSA), fibrinogen (FB), and transferrin (TRF) interactions with polyarylethersulfone (PAES) hemodialysis membranes. A multi-method approach was utilized, including protein crystallography for structural insights, hydration layer analysis to explore solvation and interaction potentials, molecular docking using AutoDock 4.0 for binding affinity simulations, and in-situ X-ray synchrotron SR-μCT imaging to observe protein deposition dynamics. Molecular docking revealed that PAES demonstrated superior binding energies and interaction patterns with FB and TRF compared to cellulose triacetate (CTA), facilitated by specific hydrogen bonding within a water shell. CTA displayed weaker, hydration-sensitive interactions varying with pH. Imaging studies indicated that FB showed higher adsorption at pH 6 than at pH 7.2, predominantly in the middle membrane regions. Similarly, HSA and TRF exhibited increased adsorption at pH 6, suggesting a stronger affinity under acidic conditions. Mixed protein solutions also indicated higher adsorption at pH 6, emphasizing an increased risk of membrane fouling. These findings highlight the crucial role of pH in modulating protein-membrane interactions and enhancing the efficacy of hemodialysis. A deeper understanding of hydration environments and their effects on protein binding affinities provides valuable insights for optimizing membrane design and performance. Clinically, this research suggests that fine-tuning pH during hemodialysis could mitigate protein fouling on membranes, thereby improving procedural efficiency and potentially leading to better patient outcomes through enhanced dialysis effectiveness.

Computationally Efficient Polarizable MD Simulations: A Simple Water Model for the Classical Drude Oscillator Polarizable Force Field

An improvement in the computational efficiency of polarizable force field simulations is made through the development of a polarizable Drude water model, SWM3, in combination with the use of Lennard-Jones Particle Mesh Ewald (LJPME) for the treatment of long-range LJ interactions. The experimental bulk properties, density, heat of vaporization, dielectric constant, and self-diffusion constant of the SWM3 model are accurately replicated at ambient condition. The temperature dependence of the bulk properties is also captured except for the density. Microscopic properties, such as hydration free energy, dimer properties, and binding energies of clusters are well represented. The SWM3 model is ∼40% faster than the SWM4 water model traditionally used with the Drude force field, and use of LJPME offers an additional 30% speedup in heterogeneous systems. This combination makes the SWM3-based Drude model only 2.6-fold slower than the TIP3P-based additive simulations.

Water Migration through Enzyme Tunnels Is Sensitive to the Choice of Explicit Water Model

The utilization of tunnels and water transport within enzymes is crucial for their catalytic function as water molecules can stabilize bound substrates and help with unbinding processes of products and inhibitors. Since the choice of water models for molecular dynamics simulations was shown to determine the accuracy of various calculated properties of the bulk solvent and solvated proteins, we have investigated if and to what extent water transport through the enzyme tunnels depends on the selection of the water model. Here, we focused on simulating enzymes with various well-defined tunnel geometries. In a systematic investigation using haloalkane dehalogenase as a model system, we focused on the well-established TIP3P, OPC, and TIP4P-Ew water models to explore their impact on the use of tunnels for water molecule transport. The TIP3P water model showed significantly faster migration, resulting in the transport of approximately 2.5 times more water molecules compared to that of the OPC and 1.7 times greater than that of the TIP4P-Ew. Finally, the transport was 1.4-fold more pronounced in TIP4P-Ew than in OPC. The increase in migration of TIP3P water molecules was mainly due to faster transit times through dehalogenase tunnels. We observed similar behavior in two different enzymes with buried active sites and different tunnel network topologies, i.e., alditol oxidase and cytochrome P450, indicating that our findings are likely not restricted to a particular enzyme family. Overall, this study showcases the critical importance of water models in comprehending the use of enzyme tunnels for small molecule transport. Given the significant role of water availability in various stages of the catalytic cycle and the solvation of substrates, products, and drugs, choosing an appropriate water model may be crucial for accurate simulations of complex enzymatic reactions, rational enzyme design, and predicting drug residence times.

Modulating the Aggregation States of a Pd6L4 Cage for Selectivity Flipping during the Stereo-Divergent Semi-Hydrogenation of Alkynes

An enzyme-mimicking catalytic system has been established using a singular palladium-based octahedral cage as the supramolecular reactor, deftly unlocking the off-on-off selectivity in the semi-hydrogenation of alkynes. Water serves as a critical regulator, modulating the catalyst states, reaction rates, and endpoints. The choice of solvent system influences the activity of host-guest binding and the reaction types of homogeneous and heterogeneous catalysis, effectively modifying the reaction steps involved in the Z→E isomerization during the semi-hydrogenation of alkynes. Kinetic and inhibition experiments indicate that the catalyst mimics the binding and activation characteristics of enzymes towards substrates, enabling selective transformations within the confined enzyme-mimicking environment. The utility of this switchable cage-confined catalysis has been demonstrated in the synthesis and modification of complex biologically active molecules with controllable E/Z selectivity. This work sheds light on the design and control of artificial supramolecular counterparts of enzymes, offering fundamental insights into the factors influencing the activity and catalytic selectivity of biological macromolecules.

Multi-Site Red-Edge Excitation Shift Reveals the Residue-Specific Solvation Dynamics during the Native to Amyloid-like Transition of an Amyloidogenic Protein

Changes in water-protein interactions are crucial for proteins to achieve functional and nonfunctional conformations during structural transitions by modulating local stability. Amyloid-like protein aggregates in deteriorating neurons are hallmarks of neurodegenerative disorders. These aggregates form through significant structural changes, transitioning from functional native conformations to supramolecular cross-β-sheet structures via misfolded and oligomeric intermediates in a multistep process. However, the site-specific dynamics of water molecules from the native to misfolded conformations and further to oligomeric and compact amyloid structures remain poorly understood. In this study, we used the fluorescence method known as red-edge excitation shift (REES) to investigate the solvation dynamics at specific sites in various equilibrium conformations en route to the misfolding and aggregation of the functional domain of the TDP-43 protein (TDP-43tRRM). We generated three single tryptophan-single cysteine mutants of TDP-43tRRM, with the cysteines at different positions and tryptophan at a fixed position. Each sole cysteine was fluorescently labeled and used as a site-specific fluorophore along with the single tryptophan, creating four monitorable sites for REES studies. By investigating the site-specific extent of REES, we developed a residue-specific solvation dynamics map of TDP-43tRRM during its misfolding and aggregation. Our observations revealed that solvation dynamics progressively became more rigid and heterogeneous to varying extents at different sites during the transition from native to amyloid-like conformations.

Comparison of Partially Denatured Cytochrome c Structural Ensembles in Solution and Gas Phases Using Cross-Linking Mass Spectrometry

Electrospray ionization mass spectrometry (ESI-MS) can retain intact protein structures, but details about partially folded and unfolded protein structures during and after introduction to the gas phase are elusive. Here we use ESI-MS with chemical cross-linkers to compare denatured cytochrome c structures in both solution and gas phases. Solution phase cross-linking prior to ESI captures solution phase structures, while gas phase cross-linking through ion/ion reactions in the trap cell captures gas phase structures. Comparing the ECD fragmentation of the cross-linked products under both conditions shows very similar cross-linker identifications, alluding to no major structural dissimilarities between solution and gas structures. Molecular modeling of the denatured protein using the identified cross-linked sites as distant restraints allows for visualization of the denatured structures to pinpoint where unfolding begins. Our data suggest that cytochrome c likely begins to unfold due to interior hydrophobic expansion, followed by α helical unfolding. This localization of structural changes is more specific than using CCS measurements alone.

Lipid bilayer strengthens the cooperative network of membrane proteins

Although membrane proteins fold and function in a lipid bilayer constituting cell membranes, their structure and functionality can be recapitulated in diverse amphiphilic assemblies whose compositions deviate from native membranes. It remains unclear how various hydrophobic environments can stabilize membrane proteins and whether lipids play any role therein. Here, using the evolutionary unrelated α-helical and β-barrel membrane proteins of Escherichia coli , we find that the hydrophobic thickness and the strength of amphiphile- amphiphile packing are critical environmental determinants of membrane protein stability. Lipid solvation enhances stability by facilitating residue burial in the protein interior and strengthens the cooperative network by promoting the propagation of local structural perturbations. This study demonstrates that lipids not only modulate membrane proteins' stability but also their response to external stimuli.

Nuclear quantum effects on glassy water under pressure: Vitrification and pressure-induced transformations

We perform classical molecular dynamics (MD) and path-integral MD (PIMD) simulations of H2O and D2O using the q-TIP4P/F model over a wide range of temperatures and pressures to study the nuclear quantum effects (NQEs) on (i) the vitrification of liquid water upon isobaric cooling at different pressures and (ii) pressure-induced transformations at constant temperature between low-density amorphous and high-density amorphous ice (LDA and HDA) and hexagonal ice Ih and HDA. Upon isobaric cooling, classical and quantum H2O and D2O vitrify into a continuum of intermediate amorphous ices (IA), with densities in-between those of LDA and HDA (depending on pressure). Importantly, the density of the IA varies considerably if NQEs are included (similar conclusions hold for ice Ih at all pressures studied). While the structure of the IA is not very sensitive to NQE, the geometry of the hydrogen-bond (HB) is. NQE leads to longer and less linear HB in LDA, HDA, and ice Ih than found in the classical case. Interestingly, the delocalization of the H/D atoms is non-negligible and identical in LDA, HDA, and ice Ih at all pressures studied. Our isothermal compression/decompression MD/PIMD simulations show that classical and quantum H2O and D2O all exhibit LDA-HDA and ice Ih-HDA transformations, consistent with experiments. The inclusion of NQE leads to a softer HB-network, which lowers slightly the LDA/ice Ih-to-HDA transformation pressures. Interestingly, the HB in HDA is longer and less linear than in LDA, which is counterintuitive given that HDA is ≈25% denser than LDA. Overall, our results show that, while classical computer simulations provide the correct qualitative phenomenology of ice and glassy water, NQEs are necessary for a quantitative description.

Characterizing Interactions Between Small Peptides and Dimethyl Sulfoxide Using Infrared Spectroscopy and Computational Methods

This study provides a comprehensive analysis of the interactions between dimethyl sulfoxide (DMSO) and two small peptides, diglycine and N-acetyl-glycine-methylamide (NAGMA), in aqueous solutions using FTIR spectroscopy and density functional theory (DFT) calculations. ATR-FTIR spectroscopy and DFT results revealed that DMSO does not form direct bonds with the peptides, suggesting that DMSO indirectly influences both peptides by modifying the surrounding water molecules. The analysis of HDO spectra allowed for the isolation of the contribution of water molecules that were simultaneously altered by the peptide and DMSO, and it also explained the changes in the hydration shells of the peptides in the presence of DMSO. In the DMSO-diglycine system, DMSO contributes to the additional strengthening of water hydrogen bonds in the reinforced hydration sphere of diglycine. In contrast, DMSO has a more moderate effect on the water molecules surrounding NAGMA due to the similarity of their hydration shells, leading to a slight weakening of the hydrogen bonds in the NAGMA hydration sphere. DFT/ONIOM calculations confirmed these observations. These findings demonstrated that DMSO influences peptide stability differentially based on their structural characteristics.

From Gas to Solution: The Changing Neutral Structure of Proline upon Solvation

Liquid-jet photoelectron spectroscopy (LJ-PES) and electronic-structure theory were employed to investigate the chemical and structural properties of the amino acid l-proline in aqueous solution for its three ionized states (protonated, zwitterionic, and deprotonated). This is the first PES study of this amino acid in its biologically relevant environment. Proline's structure in the aqueous phase under neutral conditions is zwitterionic, distinctly different from the nonionic neutral form in the gas phase. By analyzing the carbon 1s and nitrogen 1s core levels as well as the valence spectra of aqueous-phase proline, we found that the electronic structure is dominated by the protonation state of each constituent molecular site (the carboxyl and amine groups) with small yet noticeable interference across the molecule. The site-specific nature of the core-level spectra enables the probing of individual molecular constituents. The valence photoelectron spectra are more difficult to interpret because of the overlapping signals of proline with the solvent and pH-adjusting agents (HCl and NaOH). Yet, we are able to reveal subtle effects of specific (hydrogen-bonding) interaction with the solvent on the electronic structure. We also demonstrate that the relevant conformational space is much smaller for aqueous-phase proline than for its gas-phase analogue. This study suggests that caution must be taken when comparing photoelectron spectra for gaseous- and aqueous-phase molecules, particularly if those molecules are readily protonated/deprotonated in solution.

In the Beginning: Let Hydration Be Coded in Proteins for Manifestation and Modulation by Salts and Adenosine Triphosphate

Water exists in the beginning and hydrates all matter. Life emerged in water, requiring three essential components in compartmentalized spaces: (1) universal energy sources driving biochemical reactions and processes, (2) molecules that store, encode, and transmit information, and (3) functional players carrying out biological activities and structural organization. Phosphorus has been selected to create adenosine triphosphate (ATP) as the universal energy currency, nucleic acids for genetic information storage and transmission, and phospholipids for cellular compartmentalization. Meanwhile, proteins composed of 20 α-amino acids have evolved into extremely diverse three-dimensional forms, including folded domains, intrinsically disordered regions (IDRs), and membrane-bound forms, to fulfill functional and structural roles. This review examines several unique findings: (1) insoluble proteins, including membrane proteins, can become solubilized in unsalted water, while folded cytosolic proteins can acquire membrane-inserting capacity; (2) Hofmeister salts affect protein stability by targeting hydration; (3) ATP biphasically modulates liquid-liquid phase separation (LLPS) of IDRs; (4) ATP antagonizes crowding-induced protein destabilization; and (5) ATP and triphosphates have the highest efficiency in inducing protein folding. These findings imply the following: (1) hydration might be encoded in protein sequences, central to manifestation and modulation of protein structures, dynamics, and functionalities; (2) phosphate anions have a unique capacity in enhancing μs-ms protein dynamics, likely through ionic state exchanges in the hydration shell, underpinning ATP, polyphosphate, and nucleic acids as molecular chaperones for protein folding; and (3) ATP, by linking triphosphate with adenosine, has acquired the capacity to spacetime-specifically release energy and modulate protein hydration, thus possessing myriad energy-dependent and -independent functions. In light of the success of AlphaFolds in accurately predicting protein structures by neural networks that store information as distributed patterns across nodes, a fundamental question arises: Could cellular networks also handle information similarly but with more intricate coding, diverse topological architectures, and spacetime-specific ATP energy supply in membrane-compartmentalized aqueous environments?

Comparative Bindings of Glycosaminoglycans with CXCL8 Monomer and Dimer: Insights from Conformational Dynamics and Kinetics of Hydrogen Bonds

GAGs bind to both the monomeric and dimeric forms of CXCL8, helping to form a concentration gradient of the chemokine that facilitates the recruitment of neutrophils to an injury site and supports other biological functions. In this study, atomistic molecular dynamics simulations were conducted to investigate the binding behavior of two hexameric GAGs sulfated at two different positions, chondroitin sulfate (CS) and heparan sulfate (HS), with the monomer (SIL8) and dimer (DIL8) forms of the CXCL8 protein. The results support that the conformational diversity of CS and HS appeared to be more when binding with monomer SIL8 than dimer DIL8. CS gained more configurational entropy from glycosidic linkage flexibility when bound to SIL8 than DIL8, with a higher energy barrier, whereas HS exhibited a lower energy barrier for configurational entropy when bound to SIL8 and DIL8. The monomer SIL8 exhibited more favorable and preferential binding with GAGs compared to DIL8. Formation of hydrogen bonds with the basic amino acids of SIL8 and GAG was more rigid and required higher activation energy to break than the other identified hydrogen bondings. Water molecules involved in hydrogen bonding with GAGs, excluding those with basic amino acids of DIL8, showed longer lifetimes and slower relaxation compared to SIL8. This suggests that water-mediated interactions also favor binding of DIL8 with GAGs. Despite having more basic amino acids, DIL8 did not display stronger binding than SIL8, indicating the significant role of basic residues in stabilizing the GAG-protein interactions in the monomers. The reason could be that the greater number of nonbasic amino acids in dimeric CXCL8 stabilizes the complex by forming water-mediated hydrogen bonds, reducing the conformational preferences for binding with GAGs. In contrast, the monomeric form of CXCL8 exhibits a higher conformational preference for protein-GAG interactions.

The low-entropy hydration shell mediated ice-binding mechanism of antifreeze proteins

Antifreeze proteins (AFPs) can inhibit ice crystal growth. The ice-binding mechanism of AFPs remains unclear, yet the hydration shells of AFPs are thought to play an important role in modulating the binding of AFPs and ice. Here, we performed all-atom molecular dynamics simulations of an AFP from Choristoneura fumiferana (CfAFP) at four different temperatures, with a focus on analysis at 240 and 300 K, to investigate the dynamic and thermodynamic characteristics of hydration shells around ice-binding surfaces (IBS) and non-ice-binding surfaces (NIBS). Our results revealed that the dynamics of CfAFP hydration shells were highly heterogeneous, with its IBS favoring a less dense and more tetrahedral solvation shell, and NIBS hydration shells having opposite features to those of the IBS. The IBS of nine typical hyperactive AFPs were found to be in pure low-entropy hydration shell region, indicating that low-entropy hydration shell region of IBS and the tetrahedral arrangements of water molecules around them mediate the ice-binding mechanism of AFPs. It is because the entropy increase of the low-entropy hydration shell around IBS, while the higher entropy water molecules at NIBS most likely prevent ice crystal growth. These findings provide new mechanistic insights into the ice-binding of AFPs.

Osmolyte-induced conformational stabilization of a hydrophobic polymer

Elucidating the mechanistic role of osmolytes on conformations of hydrophobic prototypical macromolecules in principle is the stepping stone towards understanding the effect of osmolytes on proteins. Motivated by this, we use equilibrium simulations and umbrella sampling techniques to dissect the underlying mechanism of osmolyte-induced conformational stability of a hydrophobic polymer. Our results unveil a remarkable osmolyte-dependent conformational stabilization of the polymer. In an aqueous solution of 4 M choline chloride (ChCl), the polymer has an even more compact structure than in water. On the other hand, an aqueous solution of 8 M urea stabilizes the extended state of the polymer. Interestingly, the polymer adopts an intermediate hairpin conformation in a mixed osmolyte solution of 4 M ChCl and 8 M urea in water due to the interplay of ChCl and urea. Our simulations identify the relative accumulation of water and the hydrophilic part of choline or preferential binding of urea near the collapsed and the extended states, respectively. Analyses split out the enthalpic and entropic contributions to the overall free energy. This decides the stabilization of the preferred conformation in the chosen osmolyte solution. Our simulations show that in an aqueous solution of ChCl, the hairpin state is stabilized by entropy gain. In contrast, the enthalpic contribution stabilizes the hairpin state in mixed environments. However, a collapsed state is energetically not favored in the presence of urea. In brief, via employing an in silico approach, the current findings indicate the importance of osmolytes in stabilizing the conformational states of hydrophobic polymers.

Graph-based analysis of H-bond networks and unsupervised learning reveal conformational coupling in prion peptide segments

In this study, we employed a comprehensive computational approach to investigate the physical chemistry of the water networks surrounding hydrated peptide segments, as derived from molecular dynamics simulations. Our analysis uncovers a complex interplay of direct and water-mediated hydrogen bonds that intricately weave through the peptides. We demonstrate that these hydrogen bond networks encode critical information about the peptides' conformational behavior, with the dimensionality of these networks showing sensitivity to the peptides' conformations. Additionally, we estimated the free-energy landscape of the peptides across various conformations, revealing that their structures are predominantly characterized by unfolded, partially folded, and folded configurations, resulting in broad and rugged free-energy surfaces due to the numerous degrees of freedom contributed by the surrounding solvent. Importantly, the structured nature of this free-energy landscape becomes obscured when conventional collective variables, such as the number of hydrogen bonds, are used. Our findings provide new insights into the molecular mechanisms that couple protein and solvent degrees of freedom, highlighting their significance in the functioning of biological systems.

Interfacial dynamics mediate surface binding events on supramolecular nanostructures

The dynamic behavior of biological materials is central to their functionality, suggesting that interfacial dynamics could also mediate the activity of chemical events at the surfaces of synthetic materials. Here, we investigate the influence of surface flexibility and hydration on heavy metal remediation by nanostructures self-assembled from small molecules that are decorated with surface-bound chelators in water. We find that incorporating short oligo(ethylene glycol) spacers between the surface and interior domain of self-assembled nanostructures can drastically increase the conformational mobility of surface-bound lead-chelating moieties and promote interaction with surrounding water. In turn, we find the binding affinities of chelators tethered to the most flexible surfaces are more than ten times greater than the least flexible surfaces. Accordingly, nanostructures composed of amphiphiles that give rise to the most dynamic surfaces are capable of remediating thousands of liters of 50 ppb Pb2+-contaminated water with single grams of material. These findings establish interfacial dynamics as a critical design parameter for functional self-assembled nanostructures.

High-Speed Three-Dimensional Scanning Force Microscopy Visualization of Subnanoscale Hydration Structures on Dissolving Calcite Step Edges

Hydration at solid-liquid interfaces plays an essential role in a wide range of phenomena in biology and in materials and Earth sciences. However, the atomic-scale dynamics of hydration have remained elusive because of difficulties associated with their direct visualization. In this work, a high-speed three-dimensional (3D) scanning force microscopy technique that produces 3D images of solid-liquid interfaces with subnanoscale resolution at a rate of 1.6 s per 3D image was developed. Using this technique, direct 3D images of moving step edges were acquired during calcite dissolution in water, and hydration structures on transition regions were visualized. A Ca(OH)2 monolayer was found to form along the step edge as an intermediate state during dissolution. This imaging process also showed that hydration layers extended from the upper terraces to the transition regions to stabilize adsorbed Ca(OH)2. This technique provides information that cannot be obtained via conventional 1D/2D measurement methods.

De novo design of mini-binder proteins against IL-2 receptor β chain

IL-2 regulates the immune response by interacting with different IL-2 receptor (IL-2R) subunits. High dose of IL-2 binds to IL-2Rβγc heterodimer, which induce various side effects while activating immune function. Disrupting IL-2 and IL-2R interactions can block IL-2 mediated immune response. Here, we used a computational approach to de novo design mini-binder proteins against IL-2R β chain (IL-2Rβ) to block IL-2 signaling. The hydrophobic region where IL-2 binds to IL-2Rβ was selected and the promising binding mode was broadly explored. Three mini-binders with amino acid numbers ranging from 55 to 65 were obtained and binder 1 showed the best effects in inhibiting CTLL-2 cells proliferation and STAT5 phosphorylation. Molecular dynamics simulation showed that the binding of binder 1 to IL-2Rβ was stable; the free energy of binder1/IL-2Rβ complex was lower, indicating that the affinity of binder 1 to IL-2Rβ was higher than that of IL-2. Free energy decomposition suggested that the ARG35 and ARG131 of IL-2Rβ might be the key to improve the affinity of binder. Our efforts provided new insights in developing of IL-2R blocker, offering a potential strategy for ameliorating the side effects of IL-2 treatment.

NMR Dynamic View of the Stabilization of the WW4 Domain by Neutral NaCl and Kosmotropic Na2SO4 and NaH2PO4

The Hofmeister series categorizes ions based on their effects on protein stability, yet the microscopic mechanism remains a mystery. In this series, NaCl is neutral, Na2SO4 and Na2HPO4 are kosmotropic, while GdmCl and NaSCN are chaotropic. This study employs CD and NMR to investigate the effects of NaCl, Na2SO4, and Na2HPO4 on the conformation, stability, binding, and backbone dynamics (ps-ns and µs-ms time scales) of the WW4 domain with a high stability and accessible side chains at concentrations ≤ 200 mM. The results indicated that none of the three salts altered the conformation of WW4 or showed significant binding to the four aliphatic hydrophobic side chains. NaCl had no effect on its thermal stability, while Na2SO4 and Na2HPO4 enhanced the stability by ~5 °C. Interestingly, NaCl only weakly interacted with the Arg27 amide proton, whereas Na2SO4 bound to Arg27 and Phe31 amide protons with Kd of 32.7 and 41.6 mM, respectively. Na2HPO4, however, bound in a non-saturable manner to Trp9, His24, and Asn36 amide protons. While the three salts had negligible effects on ps-ns backbone dynamics, NaCl and Na2SO4 displayed no effect while Na2HPO4 significantly increased the µs-ms backbone dynamics. These findings, combined with our recent results with GdmCl and NaSCN, suggest a microscopic mechanism for the Hofmeister series. Additionally, the data revealed a lack of simple correlation between thermodynamic stability and backbone dynamics, most likely due to enthalpy-entropy compensation. Our study rationalizes the selection of chloride and phosphate as the primary anions in extracellular and intracellular spaces, as well as polyphosphate as a primitive chaperone in certain single-cell organisms.

Water will Find Its Way: Transport through Narrow Tunnels in Hydrolases

An aqueous environment is vital for life as we know it, and water is essential for nearly all biochemical processes at the molecular level. Proteins utilize water molecules in various ways. Consequently, proteins must transport water molecules across their internal network of tunnels to reach the desired action sites, either within them or by functioning as molecular pipes to control cellular osmotic pressure. Despite water playing a crucial role in enzymatic activity and stability, its transport has been largely overlooked, with studies primarily focusing on water transport across membrane proteins. The transport of molecules through a protein's tunnel network is challenging to study experimentally, making molecular dynamics simulations the most popular approach for investigating such events. In this study, we focused on the transport of water molecules across three different α/β-hydrolases: haloalkane dehalogenase, epoxide hydrolase, and lipase. Using a 5 μs adaptive simulation per system, we observed that only a few tunnels were responsible for the majority of water transport in dehalogenase, in contrast to a higher diversity of tunnels in other enzymes. Interestingly, water molecules could traverse narrow tunnels with subangstrom bottlenecks, which is surprising given the commonly accepted water molecule radius of 1.4 Å. Our analysis of the transport events in such narrow tunnels revealed a markedly increased number of hydrogen bonds formed between the water molecules and protein, likely compensating for the steric penalty of the process. Overall, these commonly disregarded narrow tunnels accounted for ∼20% of the total water transport observed, emphasizing the need to surpass the standard geometrical limits on the functional tunnels to properly account for the relevant transport processes. Finally, we demonstrated how the obtained insights could be applied to explain the differences in a mutant of the human soluble epoxide hydrolase associated with a higher incidence of ischemic stroke.

An integrative structural study of the human full-length RAD52 at 2.2 Å resolution

Human RAD52 (RAD52) is a DNA-binding protein involved in many DNA repair mechanisms and genomic stability maintenance. In the last few years, this protein was discovered to be a promising novel pharmacological target for anticancer strategies. Although the interest in RAD52 has exponentially grown in the previous decade, most information about its structure and mechanism still needs to be elucidated. Here, we report the 2.2 Å resolution cryo-EM reconstruction of the full-length RAD52 (FL-RAD52) protein. This allows us to describe the hydration shell of the N-terminal region of FL-RAD52, which is structured in an undecamer ring. Water molecules coordinate with protein residues to promote stabilization inside and among the protomers and within the inner DNA binding cleft to drive protein-DNA recognition. Additionally, through a multidisciplinary approach involving SEC-SAXS and computational methods, we comprehensively describe the highly flexible and dynamic organization of the C-terminal portion of FL-RAD52. This work discloses unprecedented structural details on the FL-RAD52, which will be critical for characterizing its mechanism of action and inhibitor development, particularly in the context of novel approaches to synthetic lethality and anticancer drug discovery.

Under the Influence of Water: Molecular Recognition of Organic Hydrophilic Molecules in Water with a Prismarene Host Driven by Hydration Effects

A water-soluble prism[5]arene host can form endo-cavity complexes with hydrophilic organic substances in water by displacing frustrated water molecules from its deep cavity. Water molecules structured at both rims of the prismarene host can mediate hydrogen bonding interactions with the guest. Water-mediated hydrogen bonding interactions were invoked here to elucidate the elevated binding affinities and selectivity of the prismarene host toward hydrophilic organic guests. We show that water at the interface of a host-guest complex can act as an extension of the host structure, facilitating the accommodation of neutral guests within the binding site. This study highlights the crucial role of water in facilitating supramolecular interactions between a deep-cavity prismarene host and organic hydrophilic guests in aqueous medium.

Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.

How Molecular Orientation Affects the Static Permittivity Profile of the Polar and Nonpolar Liquid-Liquid Interface

The dielectric permittivity across the liquid-liquid interface presents an intrinsic response, with respect to the instantaneous interface reference. We hypothesize that dielectric responses across the nonpolar and polar liquid-liquid interfaces have different behaviors and underlying mechanisms. Molecular dynamics simulations were used to compare and contrast the dielectric response of a nonpolar (1,2-dichloroethane/water) and a polar (1-octanol/water) liquid-liquid interface system. We found that the enhanced dielectric permittivity at the nonpolar interface is attributed to the increased water dipole orientation and polarization density. In the case of the polar interface, strong association of the immiscible solvents inhibits the molecular dipole orientation, counteracting the effect from the enhanced surface water polarization density and resulting in a standard dielectric response. Detailed knowledge of the hydrogen bond networks and molecular dipole orientation with respect to the specific instantaneous interfacial and bulk regions reveals the effect of molecular proximity and the interaction with the opposing interfacial molecules on the mechanism of the dielectric permittivity response across the liquid-liquid interface phase boundary.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

NMR Dynamic View of the Destabilization of WW4 Domain by Chaotropic GdmCl and NaSCN

GdmCl and NaSCN are two strong chaotropic salts commonly used in protein folding and stability studies, but their microscopic mechanisms remain enigmatic. Here, by CD and NMR, we investigated their effects on conformations, stability, binding and backbone dynamics on ps-ns and µs-ms time scales of a 39-residue but well-folded WW4 domain at salt concentrations ≤200 mM. Up to 200 mM, both denaturants did not alter the tertiary packing of WW4, but GdmCl exerted more severe destabilization than NaSCN. Intriguingly, GdmCl had only weak binding to amide protons, while NaSCN showed extensive binding to both hydrophobic side chains and amide protons. Neither denaturant significantly affected the overall ps-ns backbone dynamics, but they distinctively altered µs-ms backbone dynamics. This study unveils that GdmCl and NaSCN destabilize a protein before the global unfolding occurs with differential binding properties and µs-ms backbone dynamics, implying the absence of a simple correlation between thermodynamic stability and backbone dynamics of WW4 at both ps-ns and µs-ms time scales.

Exploring the Application of SiteMap and Site Finder for Focused Cryptic Pocket Identification

The characterization of cryptic pockets has been elusive, despite substantial efforts. Computational modeling approaches, such as molecular dynamics (MD) simulations, can provide atomic-level details of binding site motions and binding pathways. However, the time scale that MD can achieve at a reasonable cost often limits its application for cryptic pocket identification. Enhanced sampling techniques can improve the efficiency of MD simulations by focused sampling of important regions of the protein, but prior knowledge of the simulated system is required to define the appropriate coordinates. In the case of a novel, unknown cryptic pocket, such information is not available, limiting the application of enhanced sampling techniques for cryptic pocket identification. In this work, we explore the ability of SiteMap and Site Finder, widely used commercial packages for pocket identification, to detect focus points on the protein and further apply other advanced computational methods. The information gained from this analysis enables the use of computational modeling, including enhanced MD sampling techniques, to explore potential cryptic binding pockets suggested by SiteMap and Site Finder. Here, we examined SiteMap and Site Finder results on 136 known cryptic pockets from a combination of the PocketMiner dataset (a recently curated set of cryptic pockets), the Cryptosite Set (a classic set of cryptic pockets), and Natural killer group 2D (NKG2D, a protein target where a cryptic pocket is confirmed). Our findings demonstrate the application of existing, well-studied tools in efficiently mapping potential regions harboring cryptic pockets.

Revealing the synergistic effect of hydration and pulsed ultrasound on the emulsifying properties of silkworm pupa protein and its stabilized emulsion

BACKGROUND

Silkworm (Bombyx moil L.) Pupa protein (SPP) is a high-quality insect protein and is considered a sustainable alternative source for traditional animal food protein. However, the utilization of SPP is limited because of its low solubility and emulsifying ability. In the present study, the synergistic effect of hydration and pulsed ultrasound on the physicochemical properties of SPP and SPP-stabilized Pickering emulsions was evaluated.

RESULTS

Pulsed ultrasound changed the particle size of SPP and its conformation. As the pulsed ultrasound increased from 0 s to 5 s, the α-helix and SS contents of SPP decreased, whereas the β-sheet and SH contents increased, which in turn improved its solubility and amphiphilicity. As a result, the SPP treated by a combination of 12 h of hydration and 3 s of ultrasound exhibited a contact angle of 74.95°, hydrophobicity of 904.83, EAI of 6.66 m2 g-1 and ESI of 190.69 min. Compared with the combination of 1 h of hydration and 5 s of ultrasound, the combination of 12 h of hydration and 3 s of ultrasound exerted more soluble and hydrophobic SPP, whereas the EAI and ESI of the samples were higher. Notably, the ultrasound-treated SPP can form a stable gel-like emulsion (oil fraction ranging from 70% to 80%).

CONCLUSION

The combination of hydration and ultrasound can effectively improve the physicochemical characteristics of SPP as well as its emulsion stability. Sufficient hydration is a cost-effective method for facilitating the modification of proteins by ultrasound treatment. © 2024 Society of Chemical Industry.

Probing the structure of water in individual living cells

Water regulates or even governs a wide range of biological processes. Despite its fundamental importance, surprisingly little is known about the structure of intracellular water. Herein we employ a Raman micro-spectroscopy technique to uncover the composition, abundance and vibrational spectra of intracellular water in individual living cells. In three different cell types, we show a small but consistent population (~3%) of non-bulk-like water. It exhibits a weakened hydrogen-bonded network and a more disordered tetrahedral structure. We attribute this population to biointerfacial water located in the vicinity of biomolecules. Moreover, our whole-cell modeling suggests that all soluble (globular) proteins inside cells are surrounded by, on average, one full molecular layer (about 2.6 Angstrom) of biointerfacial water. Furthermore, relative invariance of biointerfacial water is observed among different single cells. Overall, our study not only opens up experimental possibilities of interrogating water structure in vivo but also provides insights into water in life.

Structural and dynamic insights into the activation of the μ-opioid receptor by an allosteric modulator

G-protein-coupled receptors (GPCRs) play pivotal roles in various physiological processes. These receptors are activated to different extents by diverse orthosteric ligands and allosteric modulators. However, the mechanisms underlying these variations in signaling activity by allosteric modulators remain largely elusive. Here, we determine the three-dimensional structure of the μ-opioid receptor (MOR), a class A GPCR, in complex with the Gi protein and an allosteric modulator, BMS-986122, using cryogenic electron microscopy. Our results reveal that BMS-986122 binding induces changes in the map densities corresponding to R1673.50 and Y2545.58, key residues in the structural motifs conserved among class A GPCRs. Nuclear magnetic resonance analyses of MOR in the absence of the Gi protein reveal that BMS-986122 binding enhances the formation of the interaction between R1673.50 and Y2545.58, thus stabilizing the fully-activated conformation, where the intracellular half of TM6 is outward-shifted to allow for interaction with the Gi protein. These findings illuminate that allosteric modulators like BMS-986122 can potentiate receptor activation through alterations in the conformational dynamics in the core region of GPCRs. Together, our results demonstrate the regulatory mechanisms of GPCRs, providing insights into the rational development of therapeutics targeting GPCRs.

Evaporation-driven interfacial restructuring induces highly efficient methanogenesis of waste biomass

Methanogenesis of waste biomass (WB) is a promising method for global sustainable development, reduction of pollution and carbon emission levels, and recovering bioenergy. Unlike in the methanogenesis of organic wastewater, in which microbial cells come into direct contact with the dissolved substrate, the 'solid-liquid-solid' modes in WB and between WB and microbial cells, which involve numerous solid-liquid interfaces, greatly hinder the methanogenesis efficiency of WB. Amongst all WB, waste activated sludge is the most complex, poorly biodegradable and representative. Herein, we highlight the role of water evaporation-driven solid-liquid interfacial restructuring of sludge in determining its methanogenesis efficiency. Non-free water evaporation increased surface roughness and adhesion, and compressed pore structure with numerous capillaries in sludge, resulting in a new solid-liquid interface of sludge with great capillary force and highly ordered interfacial water molecules, which provides an extremely favourable condition for high mass transfer and proton-coupled electron transfer (PCET) in sludge. This restructuring was confirmed to induce the enhancement of solid-liquid interfacial noncovalent interactions and electron transfer efficiency in the subsequent methanogenesis process (P < 0.05), promoting the effective contact between the sludge substrate and microbial cells, thereby enriching the methanogenic consortia (i.e., Clostridia and Methanosarcina were increased by 290.0 % and 239.7 %, respectively) and improving the activities of key enzymes. Stable isotope tracing and metagenomic analysis further reveal that this restructuring promoted the participation of water molecules in the methane formation by PCET-driven release of protons from water, and enhanced main methanogenesis metabolic pathways, especially the metabolic pathway of CO2-reduction methanogenesis (+65.2 %), thereby resulting in a great advance in methane generation (+147 %, P < 0.001). The findings can provide a reference for regulating directional anaerobic biotransformation of water-rich multiphase complex substrates by interfacial restructuring inducement.

Adenosine Triphosphate: The Primordial Molecule That Controls Protein Homeostasis and Shapes the Genome-Proteome Interface

Adenosine triphosphate (ATP) acts as the universal energy currency that drives various biological processes, while nucleic acids function to store and transmit genetic information for all living organisms. Liquid-liquid phase separation (LLPS) represents the common principle for the formation of membrane-less organelles (MLOs) composed of proteins rich in intrinsically disordered regions (IDRs) and nucleic acids. Currently, while IDRs are well recognized to facilitate LLPS through dynamic and multivalent interactions, the precise mechanisms by which ATP and nucleic acids affect LLPS still remain elusive. This review summarizes recent NMR results on the LLPS of human FUS, TDP-43, and the viral nucleocapsid (N) protein of SARS-CoV-2, as modulated by ATP and nucleic acids, revealing the following: (1) ATP binds to folded domains overlapping with nucleic-acid-binding interfaces; (2) ATP and nucleic acids interplay to biphasically modulate LLPS by competitively binding to overlapping pockets of folded domains and Arg/Lys within IDRs; (3) ATP energy-independently induces protein folding with the highest efficiency known so far. As ATP likely emerged in the prebiotic monomeric world, while LLPS represents a pivotal mechanism to concentrate and compartmentalize rare molecules for forming primordial cells, ATP appears to control protein homeostasis and shape genome-proteome interfaces throughout the evolutionary trajectory, from prebiotic origins to modern cells.

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Free Energy Density of a Fluid and Its Role in Solvation and Binding

The concept that a fluid has a position-dependent free energy density appears in the literature but has not been fully developed or accepted. We set this concept on an unambiguous theoretical footing via the following strategy. First, we set forth four desiderata that should be satisfied by any definition of the position-dependent free energy density, f(R), in a system comprising only a fluid and a rigid solute: its volume integral, plus the fixed internal energy of the solute, should be the system free energy; it deviates from its bulk value, fbulk, near a solute but should asymptotically approach fbulk with increasing distance from the solute; it should go to zero where the solvent density goes to zero; and it should be well-defined in the most general case of a fluid made up of flexible molecules with an arbitrary interaction potential. Second, we use statistical thermodynamics to formulate a definition of the free energy density that satisfies these desiderata. Third, we show how any free energy density satisfying the desiderata may be used to analyze molecular processes in solution. In particular, because the spatial integral of f(R) equals the free energy of the system, it can be used to compute free energy changes that result from the rearrangement of solutes as well as the forces exerted on the solutes by the solvent. This enables the use of a thermodynamic analysis of water in protein binding sites to inform ligand design. Finally, we discuss related literature and address published concerns regarding the thermodynamic plausibility of a position-dependent free energy density. The theory presented here has applications in theoretical and computational chemistry and may be further generalizable beyond fluids, such as to solids and macromolecules.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

Asparagine-85 Stabilizes a Structural Active Site Water Network in CYP121A1 of Mycobacterium tuberculosis

The cytochrome P450 enzyme CYP121A1 endogenously catalyzes the formation of a carbon-carbon bond between the two phenol groups of dicyclotyrosine (cYY) in Mycobacterium tuberculosis (Mtb). One of 20 CYP enzymes in Mtb, CYP121A1 continues to garner significant interest as a potential drug target. The accompanying reports the use of 19F NMR spectroscopy, reconstituted activity assays, and molecular dynamics simulations to investigate the significance of hydrogen bonding interactions that were theorized to stabilize a static active site water network. The active site residue Asn-85, whose hydrogen bonds with the diketopiperazine ring of cYY contributes to a contiguous active site water network in the absence of cYY, was mutated to a serine (N85S) and to a glutamine (N85Q). These conservative changes in the hydrogen bond donor side chain result in inactivation of the enzyme. Moreover, the N85S mutation induces reverse type-I binding as measured by absorbance difference spectra. NMR spectra monitoring the ligand-adaptive FG-loop and the active site Trp-182 side chain confirm that disruption of the active site water network also significantly alters the structure of the active site. These data were consistent with dynamics simulations of N85S and N85Q that demonstrate that a compromised water network is responsible for remodeling of the active site B-helix and a repositioning of cYY toward the heme. These findings implicate a slowly exchanging water network as a critical factor in CYP121A1 function and a likely contributor to the unusual rigidity of the structure.

Impact of Confinement within a Hydrogel Mesh on Protein Thermodynamic Stability and Aggregation Kinetics

Though protein stability and aggregation have been well characterized in dilute solutions, the influence of a confining environment that exists (e.g., in intercellular and tissue spaces and therapeutic formulations) on the protein structure is largely unknown. Herein, the effects of confinement on stability and aggregation were explored for proteins of different sizes, stability, and hydrophobicity when encapsulated in hydrophilic poly(ethylene glycol) hydrogels. Denaturation curves show linear correlations between confinement size (mesh size) and thermodynamic stability, i.e., unfolding free energy and surface area accessible for solvation (represented by m-value). Two clusters of protein types are identifiable from these correlations; the clusters are defined by differences in protein stability, surface area, and aggregation propensity. Proteins with higher stability, larger surface area, and lower aggregation propensity (e.g., lysozyme and myoglobin) are less affected by the confinement imposed by mesh size than proteins with lower stability, smaller surface area, and higher aggregation propensity (e.g., growth hormone and aldehyde dehydrogenase). According to aggregation kinetics measured by thioflavin T fluorescence, confinement in smaller mesh sizes resulted in slower aggregation rates than that in larger mesh sizes. Compared to that in buffer solution, the confinement of a hydrophobic protein (e.g., human insulin) in the hydrogels accelerates protein aggregation. Conversely, the confinement of a hydrophilic protein (e.g., human amylin) in the hydrogels decelerates or prevents aggregation, with the rates of aggregation inversely proportional to mesh size. These findings provide new insights into protein conformational stability in confined microenvironments relevant to various cellular, tissue, and therapeutics scenarios.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Role of tungsten disulfide quantum dots in specific protein-protein interactions at air-water interface

The intriguing network of antibody-antigen (Ab-Ag) interactions is highly governed by environmental perturbations and the nature of biomolecular interaction. Protein-protein interactions (PPIs) have potential applications in developing protein-adsorption-based sensors and nano-scale materials. Therefore, characterizing PPIs in the presence of a nanomaterial at the molecular level becomes imperative. The present work involves the investigation of antiferritin-ferritin (Ab-Ag) protein interactions under the influence of tungsten disulfide quantum dots (WS2 QDs). Isothermal calorimetry and contact angle measurements validated the strong influence of WS2 QDs on Ab-Ag interactions. The interfacial signatures of nano-bio-interactions were evaluated using sum frequency generation vibration spectroscopy (SFG-VS) at the air-water interface. Our SFG results reveal a variation in the tilt angle of methyl groups by ∼12° ± 2° for the Ab-Ag system in the presence of WS2 QDs. The results illustrated an enhanced ordering of water molecules in the presence of QDs, which underpins the active role of interfacial water molecules during nano-bio-interactions. We have also witnessed a differential impact of QDs on Ab-Ag by raising the concentration of the Ab-Ag combination, which showcased an increased inter-molecular interaction among the Ab and Ag molecules and a minimal influence on the methyl tilt angle. These findings suggest the formation of stronger and ordered Ab-Ag complexes upon introducing WS2 QDs in the aqueous medium and signify the potentiality of WS2 QDs relevant to protein-based sensing assays.

Spontaneous Dipole Reorientation in Confined Water and Its Effect on Wetting/Dewetting of Hydrophobic Nanopores

The properties of nanoconfined fluids are important for a broad range of natural and engineering systems. In particular, wetting/dewetting of hydrophobic nanoporous materials is crucial due to their broad applicability for molecular separation and liquid purification; energy storage, conversion, recuperation, and dissipation; for catalysis, chromatography, and so on. In this work, a rapid, orchestrated, and spontaneous dipole reorientation was observed in hydrophobic nanotubes of various pore sizes d (7.9-16.5 Å) via simulations. This phenomenon leads to the fragmentation of water clusters in the narrow nanopores (d = 7.9, 10 Å) and strongly affects dewetting through cluster repulsion. The cavitation in these pores has an electrostatic origin. The dependence of hydrogen-bonded network properties on the tube aperture is obtained and is used to explain wetting (intrusion)-dewetting (extrusion) hysteresis. Computer simulations and experimental data demonstrate that d equals ca. 12.5 Å is a threshold between a nonhysteretic (spring) behavior, where intrusion-extrusion is reversible, and a hysteretic one (shock absorber), where hysteresis is prominent. This work suggests that water clustering and the electrostatic nature of cavitation are important factors that can be effectively exploited for controlling the wetting/dewetting of nanoporous materials.

Phase separation of multicomponent peptide mixtures into dehydrated clusters with hydrophilic cores

Phase separation of biomolecules underlies the formation and regulation of various membraneless condensates in cells. How condensates function reliably while surrounded by heterogeneous and dynamic mixtures of biomolecular components with specific and nonspecific interactions is yet to be understood. Studying multicomponent biomolecular mixtures with designer peptides has recently become an attractive avenue for learning about physicochemical principles governing cellular condensates. In this work, we employed long-timescale atomistic simulations of multicomponent tripeptide mixtures with all residue substitutions to illuminate the nature of direct and water-mediated interactions in a prototypical cellular condensate environment. We find that peptide mixtures form clusters with inverse hydrophobic order. Most multivalent and charged residues are localized in the cluster's core, with a large fraction of nonaromatic hydrophobic residues remaining on the surface. This inverse hydrophobic order in peptide clusters is partly driven by the expulsion of nonspecifically bound water molecules following peptide cluster growth. The growth of clusters is also accompanied by the formation of increasing numbers of specific water-mediated interactions between polar and charged residues. While the present study focused on the condensation of short peptide motifs, the general findings and analysis techniques should be helpful for future studies on larger peptides and protein condensation.

De novo design of a protein binder against Staphylococcus enterotoxin B

Staphylococcus enterotoxin B (SEB) interacts with MHC-II molecules to overactivate immune cells and thereby to produce excessive pro-inflammatory cytokines. Disrupting the interactions between SEB and MHC-II helps eliminate the lethal threat posed by SEB. In this study, a de novo computational approach was used to design protein binders targeting SEB. The MHC-II binding domain of SEB was selected as the target, and the possible promising binding mode was broadly explored. The obtained original binder was folded into triple-helix bundles and contained 56 amino acids with molecular weight 5.9 kDa. The interface of SEB and the binder was highly hydrophobic. ProteinMPNN optimization further enlarged the hydrophobic region of the binder and improved the stability of the binder-SEB complex. In vitro study demonstrated that the optimized binder significantly inhibited the inflammatory response induced by SEB. Overall, our research demonstrated the applicability of this approach in de novo designing protein binders against SEB, and thereby providing potential therapeutics for SEB induced diseases.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes

The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.

Determination of the Thickness of Interfacial Water by Time-Resolved Sum-Frequency Generation Vibrational Spectroscopy

The physics and chemistry of a charged interface are governed by the structure of the electrical double layer (EDL). Determination of the interfacial water thickness (diw) of the charged interface is crucial to quantitatively describe the EDL structure, but it can be utilized with very scarce experimental methods. Here, we propose and verify that the vibrational relaxation time (T1) of the OH stretching mode at 3200 cm-1, obtained by time-resolved sum frequency generation vibrational spectroscopy with ssp polarizations, provides an effective tool to determine diw. By investigating the T1 values at the SiO2/NaCl solution interface, we established a time-space (T1-diw) relationship. We find that water has a T1 lifetime of ≥0.5 ps for diw ≤ 3 Å, while it displays bulk-like dynamics with T1 ≤ 0.2 ps for diw ≥ 9 Å. T1 decreases as diw increases from ∼3 Å to 9 Å. The hydration water at the DPPG lipid bilayer and LK15β protein interfaces has a thickness of ≥9 Å and shows a bulk-like feature. The time-space relationship will provide a novel tool to pattern the interfacial topography and heterogeneity in Ångstrom-depth resolution by imaging the T1 values.

Scrutinizing the protein hydration shell from molecular dynamics simulations against consensus small-angle scattering data

Biological macromolecules in solution are surrounded by a hydration shell, whose structure differs from the structure of bulk solvent. While the importance of the hydration shell for numerous biological functions is widely acknowledged, it remains unknown how the hydration shell is regulated by macromolecular shape and surface composition, mainly because a quantitative probe of the hydration shell structure has been missing. We show that small-angle scattering in solution using X-rays (SAXS) or neutrons (SANS) provide a protein-specific probe of the protein hydration shell that enables quantitative comparison with molecular simulations. Using explicit-solvent SAXS/SANS predictions, we derived the effect of the hydration shell on the radii of gyration Rg of five proteins using 18 combinations of protein force field and water model. By comparing computed Rg values from SAXS relative to SANS in D2O with consensus SAXS/SANS data from a recent worldwide community effort, we found that several but not all force fields yield a hydration shell contrast in remarkable agreement with experiments. The hydration shell contrast captured by Rg values depends strongly on protein charge and geometric shape, thus providing a protein-specific footprint of protein-water interactions and a novel observable for scrutinizing atomistic hydration shell models against experimental data.

Osmolyte induced protein stabilization: modulation of associated water dynamics might be a key factor

The mechanism of protein stabilization by osmolytes remains one of the most important and long-standing puzzles. The traditional explanation of osmolyte-induced stability through the preferential exclusion of osmolytes from the protein surface has been seriously challenged by the observations like the concentration-dependent reversal of osmolyte-induced stabilization/destabilization. The more modern explanation of protein stabilization/destabilization by osmolytes considers an indirect effect due to osmolyte-induced distortion of the water structure. It provides a general mechanism, but there are numerous examples of protein-specific effects, i.e., a particular osmolyte might stabilize one protein, but destabilize the other, that could not be rationalized through such an explanation. Herein, we hypothesized that osmolyte-induced modulation of associated water might be a critical factor in controlling protein stability in such a medium. Taking different osmolytes and papain as a protein, we proved that our proposal could explain protein stability in osmolyte media. Stabilizing osmolytes rigidify associated water structures around the protein, whereas destabilizing osmolytes make them flexible. The strong correlation between the stability and the associated water dynamics, and the fact that such dynamics are very much protein specific, established the importance of considering the modulation of associated water structures in explaining the osmolyte-induced stabilization/destabilization of proteins. More interestingly, we took another protein, bromelain, for which a traditionally stabilizing osmolyte, sucrose, acts as a stabilizer at higher concentrations but as a destabilizer at lower concentrations. Our proposal successfully explains such observations, which is probably impossible by any known mechanisms. We believe this report will trigger much research in this area.