WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

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.

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 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.

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.

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.

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.

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.

Water molecule ordering on the surface of an intrinsically disordered protein

The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.

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 lineshapes 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 non-adjacent 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 were observed on a picosecond timescale. These dynamics are in stark contrast with liquid water that remains highly dynamic even in the presence of ions at high concentrations.

Using Implicit-Solvent Potentials to Extract Water Contributions to Enthalpy-Entropy Compensation in Biomolecular Associations

Biomolecular assembly typically exhibits enthalpy-entropy compensation (EEC) behavior whose molecular origin remains a long-standing puzzle. While water restructuring is believed to play an important role in EEC, its contribution to the entropy and enthalpy changes, and how these changes relate to EEC, remains poorly understood. Here, we show that water reorganization entropy/enthalpy can be obtained by exploiting the temperature dependence in effective, implicit-solvent potentials. We find that the different temperature dependencies in the hydrophobic interaction, rooted in water reorganization, result in substantial variations in the entropy/enthalpy change, which are responsible for EEC. For lower-critical-solution-temperature association, water reorganization entropy dominates the free-energy change at the expense of enthalpy; for upper-critical-solution-temperature association, water reorganization enthalpy drives the process at the cost of entropy. Other effects, such as electrostatic interaction and conformation change of the macromolecules, contribute much less to the variations in entropy/enthalpy.

LAWS: Local alignment for water sites-Tracking ordered water in simulations

Accurate modeling of protein-water interactions in molecular dynamics (MD) simulations is important for understanding the molecular basis of protein function. Data from x-ray crystallography can be useful in assessing the accuracy of MD simulations, in particular, the locations of crystallographic water sites (CWS) coordinated by the protein. Such a comparison requires special methodological considerations that take into account the dynamic nature of proteins. However, existing methods for analyzing CWS in MD simulations rely on global alignment of the protein onto the crystal structure, which introduces substantial errors in the case of significant structural deviations. Here, we propose a method called local alignment for water sites (LAWS), which is based on multilateration-an algorithm widely used in GPS tracking. LAWS considers the contacts formed by CWS and protein atoms in the crystal structure and uses these interaction distances to track CWS in a simulation. We apply our method to simulations of a protein crystal and to simulations of the same protein in solution. Compared with existing methods, LAWS defines CWS characterized by more prominent water density peaks and a less-perturbed protein environment. In the crystal, we find that all high-confidence crystallographic waters are preserved. Using LAWS, we demonstrate the importance of crystal packing for the stability of CWS in the unit cell. Simulations of the protein in solution and in the crystal share a common set of preserved CWS that are located in pockets and coordinated by residues of the same domain, which suggests that the LAWS algorithm will also be useful in studying ordered waters and water networks in general.

Correlation between chemical denaturation and the unfolding energetics of Acanthamoeba actophorin

The actin filament network is in part remodeled by the action of a family of filament severing proteins that are responsible for modulating the ratio between monomeric and filamentous actin. Recent work on the protein actophorin from the amoeba Acanthamoeba castellani identified a series of site-directed mutations that increase the thermal stability of the protein by 22°C. Here, we expand this observation by showing that the mutant protein is also significantly stable to both equilibrium and kinetic chemical denaturation, and employ computer simulations to account for the increase in thermal or chemical stability through an accounting of atomic-level interactions. Specifically, the potential of mean force (PMF) can be obtained from steered molecular dynamics (SMD) simulations in which a protein is unfolded. However, SMD can be inefficient for large proteins as they require large solvent boxes, and computationally expensive as they require increasingly many SMD trajectories to converge the PMF. Adaptive steered molecular dynamics (ASMD) overcomes the second of these limitations by steering the particle in stages, which allows for convergence of the PMF using fewer trajectories compared with SMD. Use of the telescoping water scheme within ASMD partially overcomes the first of these limitations by reducing the number of waters at each stage to only those needed to solvate the structure within a given stage. In the PMFs obtained from ASMD, the work of unfolding Acto-2 was found to be higher than the Acto-WT by approximately 120 kCal/mol and reflects the increased stability seen in the chemical denaturation experiments. The evolution of the average number of hydrogen bonds and number of salt bridges during the pulling process provides a mechanistic view of the structural changes of the actophorin protein as it is unfolded, and how it is affected by the mutation in concert with the energetics reported through the PMF.

Spatial Layouts of Low-Entropy Hydration Shells Guide Protein Binding

Protein-protein binding enables orderly biological self-organization and is therefore considered a miracle of nature. Protein‒protein binding is driven by electrostatic forces, hydrogen bonding, van der Waals force, and hydrophobic interactions. Among these physical forces, only hydrophobic interactions can be considered long-range intermolecular attractions between proteins due to the electrostatic shielding of surrounding water molecules. Low-entropy hydration shells around proteins drive hydrophobic attraction among them that essentially coordinate protein‒protein binding. Here, an innovative method is developed for identifying low-entropy regions of hydration shells of proteins by screening off pseudohydrophilic groups on protein surfaces and revealing that large low-entropy regions of the hydration shells typically cover the binding sites of individual proteins. According to an analysis of determined protein complex structures, shape matching between a large low-entropy hydration shell region of a protein and that of its partner at the binding sites is revealed as a universal law. Protein‒protein binding is thus found to be mainly guided by hydrophobic collapse between the shape-matched low-entropy hydration shells that is verified by bioinformatics analyses of hundreds of structures of protein complexes, which cover four test systems. A simple algorithm is proposed to accurately predict protein binding sites.

Detecting Interplay of Chirality, Water, and Interfaces for Elucidating Biological Functions

Chemists have long been fascinated by chirality, water, and interfaces, making tremendous progress in each research area. However, the chemistry emerging from the interplay of chirality, water, and interfaces has been difficult to study due to technical challenges, creating a barrier to elucidating biological functions at interfaces. Most biopolymers (proteins, DNA, and RNA) fold into macroscopic chiral structures to perform biological functions. Their folding requires water, but water behaves differently at interfaces where the bulk water hydrogen-bonding network terminates. A question arises as to how water molecules rearrange to minimize free energy at interfaces while stabilizing the macroscopic folding of biopolymers to support biological function. This question is central to solving many research challenges, including the molecular origin of biological homochirality, folding and insertion of proteins into cell membranes, and the design of heterogeneous biocatalysts. Researchers can resolve these challenges if they have the theoretical tools to accurately predict molecular behaviors of water and biopolymers at various interfaces. However, developing such tools requires validation by the experimental data. These experimental data are scarce because few physical methods can simultaneously distinguish chiral folding of the biopolymers, separate signals of interfaces from the overwhelming background of bulk solvent, and differentiate water in hydration shells of the polymers from water elsewhere.We recently illustrated these very capacities of chirality-sensitive vibrational sum frequency generation spectroscopy (chiral SFG). While chiral SFG theory dictates that the method is surface-specific under the condition of electronic nonresonance, we show the method can distinguish chiral folding of proteins and DNA and probe water structures in the first hydration shell of proteins at interfaces. Using amide I signals, we observe protein folding into β-sheets without background signals from α-helices and disordered structures at interfaces, thereby demonstrating the effect of 2D crowding on protein folding. Also, chiral SFG signals of C-H stretches are silent from single-stranded DNA, but prominent for canonical antiparallel duplexes as well as noncanonical parallel duplexes at interfaces, allowing for sensing DNA secondary structures and hybridization. In establishing chiral SFG for detecting protein hydration structures, we observe an H218O isotopic shift that reveals water contribution to the chiral SFG spectra. Additionally, the phase of the O-H stretching bands flips when the protein chirality is switched from L to D. These experimental results agree with our simulated chiral SFG spectra of water hydrating the β-sheet protein at the vacuum-water interface. The simulations further reveal that over 90% of the total chiral SFG signal comes from water in the first hydration shell. We conclude that the chiral SFG signals originate from achiral water molecules that assemble around the protein into a chiral supramolecular structure with chirality transferred from the protein. As water O-H stretches can reveal hydrogen-bonding interactions, chiral SFG shows promise in probing the structures and dynamics of water-biopolymer interactions at interfaces. Altogether, our work has created an experimental and computational framework for chiral SFG to elucidate biological functions at interfaces, setting the stage for probing the intricate chemical interplay of chirality, water, and interfaces.

Lipid Bilayer Strengthens the Cooperative Network of a Membrane-Integral Enzyme

Lipid bilayer provides a two-dimensional hydrophobic solvent milieu for membrane proteins in cells. Although the native bilayer is widely recognized as an optimal environment for folding and function of membrane proteins, the underlying physical basis remains elusive. Here, employing the intramembrane protease GlpG of Escherichia coli as a model, we elucidate how the bilayer stabilizes a membrane protein and engages the protein's residue interaction network compared to the nonnative hydrophobic medium, micelles. We find that the bilayer enhances GlpG stability by promoting residue burial in the protein interior compared to micelles. Strikingly, while the cooperative residue interactions cluster into multiple distinct regions in micelles, the whole packed regions of the protein act as a single cooperative unit in the bilayer. Molecular dynamics (MD) simulation indicates that lipids less efficiently solvate GlpG than detergents. Thus, the bilayerinduced enhancement of stability and cooperativity likely stems from the dominant intraprotein interactions outcompeting the weak lipid solvation. Our findings reveal a foundational mechanism in the folding, function, and quality control of membrane proteins. The enhanced cooperativity benefits function facilitating propagation of local structural perturbation across the membrane. However, the same phenomenon can render the proteins' conformational integrity vulnerable to missense mutations causing conformational diseases1,2.

The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice

We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (I_h) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H₂O and D₂O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice I_h. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice I_h increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice I_h. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient α_P(T) and bulk modulus B(T) of both LDA and ice I_h. Remarkably, the ρ(T), α_P(T), and B(T) of LDA are practically identical to those of ice I_h. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice I_h. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.

Water potential governs the effector specificity of the transcriptional regulator XylR of Pseudomonas putida

The biodegradative capacity of bacteria in their natural habitats is affected by water availability. In this work, we have examined the activity and effector specificity of the transcriptional regulator XylR of the TOL plasmid pWW0 of Pseudomonas putida mt-2 for biodegradation of m-xylene when external water potential was manipulated with polyethylene glycol PEG8000. By using non-disruptive luxCDEAB reporter technology, we noticed that the promoter activated by XylR (Pu) restricted its activity and the regulator became more effector-specific towards head TOL substrates when cells were grown under water subsaturation. Such a tight specificity brought about by water limitation was relaxed when intracellular osmotic stress was counteracted by the external addition of the compatible solute glycine betaine. With these facts in hand, XylR variants isolated earlier as effector-specificity responders to the non-substrate 1,2,4-trichlorobenzene under high matric stress were re-examined and found to be unaffected by water potential in vivo. All these phenomena could be ultimately explained as the result of water potential-dependent conformational changes in the A domain of XylR and its effector-binding pocket, as suggested by AlphaFold prediction of protein structures. The consequences of this scenario for the evolution of specificities in regulators and the emergence of catabolic pathways are discussed.

Light, Water, and Melatonin: The Synergistic Regulation of Phase Separation in Dementia

The swift rise in acceptance of molecular principles defining phase separation by a broad array of scientific disciplines is shadowed by increasing discoveries linking phase separation to pathological aggregations associated with numerous neurodegenerative disorders, including Alzheimer’s disease, that contribute to dementia. Phase separation is powered by multivalent macromolecular interactions. Importantly, the release of water molecules from protein hydration shells into bulk creates entropic gains that promote phase separation and the subsequent generation of insoluble cytotoxic aggregates that drive healthy brain cells into diseased states. Higher viscosity in interfacial waters and limited hydration in interiors of biomolecular condensates facilitate phase separation. Light, water, and melatonin constitute an ancient synergy that ensures adequate protein hydration to prevent aberrant phase separation. The 670 nm visible red wavelength found in sunlight and employed in photobiomodulation reduces interfacial and mitochondrial matrix viscosity to enhance ATP production via increasing ATP synthase motor efficiency. Melatonin is a potent antioxidant that lowers viscosity to increase ATP by scavenging excess reactive oxygen species and free radicals. Reduced viscosity by light and melatonin elevates the availability of free water molecules that allow melatonin to adopt favorable conformations that enhance intrinsic features, including binding interactions with adenosine that reinforces the adenosine moiety effect of ATP responsible for preventing water removal that causes hydrophobic collapse and aggregation in phase separation. Precise recalibration of interspecies melatonin dosages that account for differences in metabolic rates and bioavailability will ensure the efficacious reinstatement of the once-powerful ancient synergy between light, water, and melatonin in a modern world.

Quantifying the Molecular Polarization Response of Liquid Water Interfaces at Heterogeneously Charged Surfaces

The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water's hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water's collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein.

Enhanced Grand Canonical Sampling of Occluded Water Sites Using Nonequilibrium Candidate Monte Carlo

Water molecules play a key role in many biomolecular systems, particularly when bound at protein-ligand interfaces. However, molecular simulation studies on such systems are hampered by the relatively long time scales over which water exchange between a protein and solvent takes place. Grand canonical Monte Carlo (GCMC) is a simulation technique that avoids this issue by attempting the insertion and deletion of water molecules within a given structure. The approach is constrained by low acceptance probabilities for insertions in congested systems, however. To address this issue, here, we combine GCMC with nonequilibium candidate Monte Carlo (NCMC) to yield a method that we refer to as grand canonical nonequilibrium candidate Monte Carlo (GCNCMC), in which the water insertions and deletions are carried out in a gradual, nonequilibrium fashion. We validate this new approach by comparing GCNCMC and GCMC simulations of bulk water and three protein binding sites. We find that not only is the efficiency of the water sampling improved by GCNCMC but that it also results in increased sampling of ligand conformations in a protein binding site, revealing new water-mediated ligand-binding geometries that are not observed using alternative enhanced sampling techniques.

Elucidating the Sluggish Water Dynamics at the Ice-Binding Surface of the Hyperactive Tenebrio molitor Antifreeze Protein

Quasi-ice-like hydration waters on the ice-binding surface (IBS) of an antifreeze protein (AFP) commonly exhibit sluggish dynamics especially at low temperatures. In this work, we have analyzed molecular dynamics (MD) simulation trajectories at two different temperatures for Tenebrio molitor antifreeze protein (TmAFP) to explore whether the unique quasi-ice-like structuring of hydration water has any impact on making their dynamics slower on the IBS of the protein. Our calculation reveals that, as translational dynamics is coupled with the conformational fluctuations, hydration water on the IBS exhibits sluggish translational motion due to reduced flexibility of the IBS compared to that on the non-ice-binding surface (NIBS) of the protein. Interestingly, it is noticed that rotational motion of hydration water is not coupled with the conformational fluctuations of the surfaces. In that case, structural relaxations of the protein-water (PW) and water-water (WW) hydrogen bonds compete with each other to make the rotational dynamics of hydration water around the IBS either faster or slower with respect to those around the NIBS. At low temperature, the slower structural relaxation of water-water hydrogen bonds dominates and imparts sluggish rotational motion of the hydration water on the IBS of the protein. The slower structural relaxation of water-water hydrogen bonds and hence the retarded rotational dynamics, despite the weak short-lived PW hydrogen bonds on the IBS, is clearly a manifestation of the rigid quasi-ice-like structure of the hydration shell on that surface.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Inconsistent hydrogen bond-mediated vibrational coupling of amide I

Using infrared spectroscopy and density functional theory (DFT) calculations, we scrutinized an amide (dimethylformamide) as a "model" compound to interpret the interactions of amide 1 with different phenol derivatives (para-chlorophenol (PCP) and para-cresol (CP)) as "model guest molecules". We established the involvement of amide I in vibrational coupling with symmetric and asymmetric C[double bond, length as m-dash]C modes of different phenolic derivatives and how their coupling was dependent upon different guest aromatic phenolic compounds. Interestingly, substitution of phenol perturbed the pattern of vibrational coupling with amide I. The symmetric and asymmetric C[double bond, length as m-dash]C modes of PC were coupled significantly with amide 1. For PCP, the symmetric C[double bond, length as m-dash]C mode coupled significantly, but the asymmetric mode coupled negligibly, with amide I. Here, we reveal the nature of vibrational coupling based on the structure of a guest molecule hydrogen-bonded with amide I. Our conclusions could be valuable for depiction of the unusual dynamics of coupled amide-I modes as well as the dependency of vibrational coupling on altered factors.

The significance of bioelectricity on all levels of organization of an organism. Part 1: From the subcellular level to cells

Bioelectricity plays an essential role in the structural and functional organization of biological organisms. In this first article of our three-part series, we summarize the importance of bioelectricity for the basic structural level of biological organization, i.e. from the subcellular level (charges, ion channels, molecules and cell organelles) to cells.

Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS

The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.

Antibody CDR amino acids underlying the functionality of antibody repertoires in recognizing diverse protein antigens

Antibodies recognize protein antigens with exquisite specificity in a complex aqueous environment, where interfacial waters are an integral part of the antibody–protein complex interfaces. In this work, we elucidate, with computational analyses, the principles governing the antibodies’ specificity and affinity towards their cognate protein antigens in the presence of explicit interfacial waters. Experimentally, in four model antibody–protein complexes, we compared the contributions of the interaction types in antibody–protein antigen complex interfaces with the antibody variants selected from phage-displayed synthetic antibody libraries. Evidently, the specific interactions involving a subset of aromatic CDR (complementarity determining region) residues largely form the predominant determinant underlying the specificity of the antibody–protein complexes in nature. The interfacial direct/water-mediated hydrogen bonds accompanying the CDR aromatic interactions are optimized locally but contribute little in determining the epitope location. The results provide insights into the phenomenon that natural antibodies with limited sequence and structural variations in an antibody repertoire can recognize seemingly unlimited protein antigens. Our work suggests guidelines in designing functional artificial antibody repertoires with practical applications in developing novel antibody-based therapeutics and diagnostics for treating and preventing human diseases.

Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor

Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.

Experimental Detection and Simulation of Terahertz Spectra of Aqueous L-Arginine

Terahertz (THz) wave is a good candidate for biological sample detection, because vibration and rotation energy levels of biomolecule are in THz band. However, the strong absorption of THz wave by water in biological samples hinders its development. In this paper, a method for direct detection of THz absorption spectra of L-arginine suspension was proposed by using a strong field THz radiation source combined with a polyethylene cell with micrometer thickness in a THz time-domain spectroscopy system. And the THz absorption spectrum of L-arginine solution was simulated by the density functional theory and the simulation result is in good agreement with the experimental results. Finally, the types of chemical bond interaction that cause the absorption peak are identified based on the experimental and simulation results. This work paves a way to investigate the THz absorption spectra and intramolecular interactions of aqueous biological samples.