Spectral Graph Analyses of Water Hydrogen-Bonding Network and Osmolyte Aggregate Structures in Osmolyte-Water Solutions

Molecular Aggregation Behavior and Microscopic Heterogeneity in Binary Osmolyte-Water Solutions

Osmolytes, small organic compounds, play a key role in modulating the protein stability in aqueous solutions, but the operating mechanism of the osmolyte remains inconclusive. Here, we attempt to clarify the mode of osmolyte action by quantitatively estimating the microheterogeneity of osmolyte-water mixtures with the aid of molecular dynamics simulation, graph theoretical analysis, and spatial distribution measurement in the four osmolyte solutions of trimethylamine-N-oxide (TMAO), tetramethylurea (TMU), dimethyl sulfoxide, and urea. TMAO, acting as a protecting osmolyte, tends to remain isolated with no formation of osmolyte aggregates while preferentially interacting with water, but there is a strong aggregation propensity in the denaturant TMU solution, characterized by favored hydrophobic interactions between TMU molecules. Taken together, the mechanism of osmolyte action on protein stability is proposed as a comprehensive one that encompasses the direct interactions between osmolytes and proteins and indirect interactions through the regulation of water properties in the osmolyte-water mixtures.

Specific Proton-Donor Properties of Glycine Betaine. Metric Parameters and Enthalpy of Noncovalent Interactions in its Dimer, Water Complexes and Crystalline Hydrate

Trimethylglycine (glycine betaine, GB) is an important organic osmolyte that accumulates in various plant species in response to environmental stresses and has significant potential as a bioactive agent with low environmental impact. It is assumed that the hydration of GB is playing an important role in the protective mechanism. The hydration and aggregation properties of GB have not yet been studied in detail at the atomistic level. In this work, noncovalent interactions in the GB dimer and its complexes with water and crystalline monohydrate are studied. Depending on the object, periodic and non-periodic DFT calculations are used. Particular attention is paid to the metric parameters and enthalpies of intermolecular hydrogen bonds. The identification of noncovalent interactions is carried out by means of the Bader analysis of periodic or non-periodic electron density. The enthalpy of hydrogen bonds is estimated using the Rosenberg formula (PCCP 2 (2000) 2699). The specific proton donor properties of glycine betaine are due to its ability to form intermolecular C-H∙∙∙O bonds with the oxygen atom of a water molecule or the carboxylate group of a neighboring GB. The enthalpy of these bonds can be significantly greater than 10 kJ/mol. The water molecule that forms a hydrogen bond with the carboxylate group of GB also interacts with its CH groups through lone pairs of electrons. The C-H∙∙∙O bonds contribute up to 40% of the total entropy of the GB-water interaction, which is about 45 kJ/mol. The possibility of identifying C-H∙∙∙O bonds by the proton nuclear magnetic resonance method is discussed.

Mechanism of Osmolyte Stabilization-Destabilization of Proteins: Experimental Evidence

In this work, we investigated the influence of stabilizing (N,N,N-trimethylglycine) and destabilizing (urea) osmolytes on the hydration spheres of biomacromolecules in folded forms (trpzip-1 peptide and hen egg white lysozyme—hewl) and unfolded protein models (glycine—GLY and N-methylglycine—NMG) by means of infrared spectroscopy. GLY and NMG were clearly limited as minimal models for unfolded proteins and should be treated with caution. We isolated the spectral share of water changed simultaneously by the biomacromolecule/model molecule and the osmolyte, which allowed us to provide unambiguous experimental arguments for the mechanism of stabilization/destabilization of proteins by osmolytes. In the case of both types of osmolytes, the decisive factor determining the equilibrium folded/unfolded state of protein was the enthalpy effect exerted on the hydration spheres of proteins in both forms. In the case of stabilizing osmolytes, enthalpy was also favored by entropy, as the unfolded state of a protein was more entropically destabilized than the folded state.

Analysis of micropollutants in a marine outfall using network analysis and decision tree

The presence of micropollutants (MPs), including pharmaceutical, industrial, and pesticidal compounds, threatens both human health and the aquatic ecosystem. The development and extensive use of new chemicals have also inevitably led to the accumulation of MPs in aquatic environments. Recreational beaches are especially vulnerable to contamination, affecting humans and aquatic animals via the absorption of MPs in water during marine activities (e.g., swimming, sailing, and windsurfing). Additionally, marine outfalls in an urbanized coastal city can cause serious chemical and microbial pollution on recreational beaches, leading to an increase in adverse effects on public health and the ecological system. Therefore, the aim of this study was to, with the use of network and decision tree analyses, identify the features and factors that influence the change in MP concentrations in a marine outfall. These analyses were conducted to inspect the relationship between each MP and its hierarchical structure as well as hydrometeorological variables. Additionally, a risk analysis was conducted in this study in which the MPs were prioritized based on their optimized risk quotient values. During our monitoring of MP concentrations over time at the marine outfall, high concentrations of pharmaceutical and industrial compounds were detected when the tide level was low after rainfall. Furthermore, results of the risk analysis and the prioritization revealed that a total of 18 substances identified in our study posed a risk to the ecosystem; these include major ecotoxicologically hazardous substances such as telmisartan, mevinphos, and methiocarb. Results of the network analysis demonstrated distinct trends for pharmaceutical and industrial substances, whilst those for pesticide compounds were irregular. Additionally, the hierarchical structures for most MPs consisted of rainfall, tide level, and antecedent dry hours; this implies that these factors influence MP dynamics. These findings will be helpful for establishing chemical contamination management plans for recreational beaches in the future.

Unraveling the Influence of Osmolytes on Water Hydrogen-Bond Network: From Local Structure to Graph Theory Analysis

Water structure in aqueous osmolyte solutions, deduced from the slight alteration in the water-water radial distribution function, the decrease in water-water hydrogen bonding, and tetrahedral ordering based only on the orientation of nearest water molecules derived from the molecular dynamics simulations, appears to have been perturbed. A careful analysis, however, reveals that the hydrogen bonding and the tetrahedral ordering around a water molecule in binary solutions remain intact as in neat water when the contribution of osmolyte-water interactions is appropriately incorporated. Furthermore, the distribution of the water binding energies and the water excess chemical potential of solvation in solutions are also pretty much the same as in neat water. Osmolytes are, therefore, well integrated into the hydrogen-bond network of water. Indeed, osmolytes tend to preferentially hydrogen bond with water molecules and their interaction energies are strongly correlated to their hydrogen-bonding capability. The graph network analysis, further, illustrates that osmolytes act as hubs in the percolated hydrogen-bond network of solutions. The degree of hydrogen bonding of osmolytes predominantly determines all of the network properties. Osmolytes like ethanol that form fewer hydrogen bonds than a water molecule disrupt the water hydrogen-bond network, while other osmolytes that form more hydrogen bonds effectively increase the connectivity among water molecules. Our observation of minimal variation in the local structure and the vitality of osmolyte-water hydrogen bonds on the solution network properties clearly imply that the direct interaction between protein and osmolytes is solely responsible for the protein stability. Further, the relevance of hydrogen bonds on solution properties suggests that the hydrogen-bonding interaction among protein, water, and osmolyte could be the key determinant of the protein conformation in solutions.

Molecular heterogeneity in aqueous cosolvent systems

Aqueous cosolvent systems (ACoSs) are mixtures of small polar molecules such as amides, alcohols, dimethyl sulfoxide, or ions in water. These liquids have been the focus of fundamental studies due to their complex intermolecular interactions as well as their broad applications in chemistry, medicine, and materials science. ACoSs are fully miscible at the macroscopic level but exhibit nanometer-scale spatial heterogeneity. ACoSs have recently received renewed attention within the chemical physics community as model systems to explore the relationship between intermolecular interactions and microscopic liquid-liquid phase separation. In this perspective, we provide an overview of ACoS spatial segregation, dynamic heterogeneity, and multiscale relaxation dynamics. We describe emerging approaches to characterize liquid microstructure, H-bond networks, and dynamics using modern experimental tools combined with molecular dynamics simulations and network-based analysis techniques.

Characterization of Specific Ion Effects on PI(4,5)P2 Clustering: Molecular Dynamics Simulations and Graph-Theoretic Analysis

Numerous cellular functions mediated by phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2; PIP2) involve clustering of the lipid as well as colocalization with other lipids. Although the cation-mediated electrostatic interaction is regarded as the primary clustering mechanism, the ion-specific nature of the intermolecular network formation makes it challenging to characterize the clusters. Here we use all-atom molecular dynamics (MD) simulations of PIP2 monolayers and graph-theoretic analysis to gain insight into the phenomenon. MD simulations reveal that the intermolecular interactions preferentially occur between specific cations and phosphate groups (P1, P4, and P5) of the inositol headgroup with better-matched kosmotropic/chaotropic characters consistent with the law of matching water affinities (LMWA). Ca2+ is strongly attracted to P4/P5, while K+ preferentially binds to P1; Na+ interacts with both P4/P5 and P1. These specific interactions lead to the characteristic clustering patterns. Specificially, the size distributions and structures of PIP2 clusters generated by kosmotropic cations Ca2+ and Na+ are bimodal, with a combination of small and large clusters, while there is little clustering in the presence of only chaotropic K+; the largest clusters are obtained in systems with all three cations. The small-world network (a model with both local and long-range connections) best characterizes the clusters, followed by the random and the scale-free networks. More generally, the present results interpreted within the LMWA are consistent with the relative eukaryotic intracellular concentrations Ca2+ ≪ Na+ < Mg2+ < K+; that is, concentrations of Ca2+ and Na+ must be low to prevent damaging aggregation of lipids, DNA, RNA and phosphate-containing proteins.

GCIceNet: a graph convolutional network for accurate classification of water phases

We develop GCIceNet, which automatically generates machine-based order parameters for classifying the phases of water molecules via supervised and unsupervised learning with graph convolutional networks.

Molecular Insights into Dipole Relaxation Processes in Water-Lysine Mixtures

Dielectric spectroscopy is a robust method to investigate relaxations of molecular dipoles. It is particularly useful for studies of biological solutions because of the potential of this method to cover a broad range of dynamical time scales typical for such systems. However, this technique does not provide any information about the nature of the molecular motions, which leads to a certain underemployment of dielectric spectroscopy for gaining microscopic understanding of material properties. For such detailed understanding, computer simulations are valuable tools because they can provide information about the nature of molecular motions observed by, for example, dielectric spectroscopy and to further complement them with structural information. In this work, we acquire information about the nature of dipole relaxation, in n-lysine solutions by means of molecular dynamics simulations. Our results indicate that the experimentally observed main relaxation process of n-lysine has different origins for the single monomer and the polypeptide chains. The relaxation of 1-lysine is due to the motions of whole molecules, whereas the experimentally observed relaxation of 3-lysine and 4-lysine is due to the motions of the residues, which, in turn, are promoted by water relaxation. Furthermore, we propose a new structural model of the lysine amino acids, which can quantitatively account for the experimental dielectric relaxation data. Hydrogen bonding and the structure of water are also discussed in terms of their influence on relaxation processes.

Understanding the effect of nanoconfinement on the structure of water hydrogen bond networks

Dynamic hydrogen bond trails in water confined between two phospholipid membranes traced by the information flow model.

Structure, molecular dynamics, and interactions in aqueous xylitol solutions

Broad-band dielectric relaxation studies of xylitol-water mixture show distinctly different dynamics for distal and central –OH of xylitol molecules and indicates the presence of loose xylitol aggregates.

Phonon Spectrometric Evaluation of the Solute-Solvent Interface in Solutions of Glycine and Its N-Methylated Derivatives

From the perspective of O:H-O bond cooperativity, we analyzed the solute capability of transiting the O:H-O bond from the mode of ordinary water to the hydration state and its consequence on the solution viscosity and surface stress. Phonon spectrometric results suggest that glycine and its N-methyl derivatives strongly affect the surrounding solvent molecules through H ↔ H repulsion and dipolar polarization. The H ↔ H interproton repulsion disrupts the surface stress, and the polarization enhances the solution viscosity.

Hierarchical lattice models of hydrogen-bond networks in water

We develop a graph-based model of the hydrogen-bond network in water, with a view toward quantitatively modeling the molecular-level correlational structure of the network. The networks formed are studied by the constructing the model on two infinite-dimensional lattices. Our models are built bottom up, based on microscopic information coming from atomistic simulations, and we show that the predictions of the model are consistent with known results from ab initio simulations of liquid water. We show that simple entropic models can predict the correlations and clustering of local-coordination defects around tetrahedral waters observed in the atomistic simulations. We also find that orientational correlations between bonds are longer ranged than density correlations, determine the directional correlations within closed loops, and show that the patterns of water wires within these structures are also consistent with previous atomistic simulations. Our models show the existence of density and compressibility anomalies, as seen in the real liquid, and the phase diagram of these models is consistent with the singularity-free scenario previously proposed by Sastry and coworkers [Phys. Rev. E 53, 6144 (1996)1063-651X10.1103/PhysRevE.53.6144].

Effect of Osmolytes on the Conformational Behavior of a Macromolecule in a Cytoplasm-like Crowded Environment: A Femtosecond Mid-IR Pump-Probe Spectroscopy Study

Osmolytes found endogenously in almost all living beings play an important role in regulating cell volume under harsh environment. Here, to address the longstanding questions about the underlying mechanism of osmolyte effects, we use femtosecond mid-IR pump-probe spectroscopy with two different IR probes that are the OD stretching mode of HDO and the azido stretching mode of azido-derivatized poly(ethylene glycol) dimethyl ether (PEGDME). Our experimental results show that protecting osmolytes bind strongly with water molecules and dehydrate polymer surface, which results in promoting intramolecular interactions of the polymer. By contrast, urea behaves like water molecules without significantly disrupting water H-bonding network and favors extended and random-coil segments of the polymer chain by directly participating in solvation of the polymer. Our findings highlight the importance of direct interaction between urea and macromolecule, while protecting osmolytes indirectly affect the macromolecule through enhancing the water-osmolyte interaction in a crowded environment, which is the case that is often encountered in real biological systems.

The effect of Hofmeister anions on water structure at protein surfaces

To understand the effects of specific ions on protein-water interactions and the thermodynamic stability of proteins in salt solutions, we use a molecular dynamics (MD) simulation to examine the water structure, orientational distribution, and dynamics near the surface of ubiquitin. In particular, we consider NaCl, NaBF₄, NaSCN, and NaClO₄ salt solutions containing ubiquitin, where the anions of the latter three salts are well-known chaotropic ions in the Hofmiester anion series. The number of hydrogen bonds (H-bonds) per water molecule is found to decrease significantly at the ubiquitin-water interface, indicating a significant disruption of the water H-bonding network. The distribution of the water H-bond numbers near the protein surface is modulated by dissolved ions, and the extent of the ion effect on the H-bonding network structure follows the order of the Hofmeister anion series, while there are no specific ion effects on water properties at distances larger than 5 Å from the protein surface. From detailed analyses of the surface area, volume, and root-mean-square deviation (RMSD) of ubiquitin, we show that changes in the properties of the protein could originate from the disruption of the water H-bond network induced by ions with a higher affinity for the protein surface instead of direct protein residue-ion interactions. An interesting observation made here is that the orientational distribution of water molecules at the protein-water interface is close to random, but there is a slight preference for interfacial water molecules with a straddle structure within 2.5 Å of the protein surface, where one of the two OH groups points away from the protein surface and the other points toward the surface. In addition, comparing the MD simulation results for ubiquitin solutions with dissolved NaSCN and KSCN, we show that Na⁺ affects the water H-bonding structure at the protein surface more than K⁺. It is clear that the H-bonding network structure of water more than one water layer away from the protein surface is not distinguishably different from that of neat water. We thus anticipate that the present work will provide insights into the scale of specific ion effects on the H-bonding structure and orientational distribution of water in the vicinity of protein surfaces in aqueous solutions.

Unravelling the role of polyols with increasing carbon chain length and OH groups on the phase transition behavior of PNIPAM

In advanced applications of pharmaceutical, agricultural and biomedical research, thermoresponsive polymers (TRPs) are potential candidates which show conformational transitions at given temperatures.

Role of solvent properties of water in crowding effects induced by macromolecular agents and osmolytes

Dipolarity/polarizability of water in polymer mixtures may be additive, reduced or enhanced depending on the composition of the mixture.