Water structure in vitro and within Saccharomyces cerevisiae yeast cells under conditions of heat shock

The OH stretch mode from water and organic hydroxyl groups have strong infrared absorption, the position of the band going to lower frequency with increased H-bonding. This band was used to study water in trehalose and glycerol solutions and in genetically modified yeast cells containing varying amounts of trehalose. Concentration-dependent changes in water structure induced by trehalose and glycerol in solution were detected, consistent with an increase of lower-energy H-bonds and interactions at the expense of higher-energy interactions. This result suggests that these molecules disrupt the water H-bond network in such a way as to strengthen molecule-water interactions while perturbing water-water interactions. The molecule-induced changes in the water H-bond network seen in solution do not translate to observable differences in yeast cells that are trehalose-deficient and trehalose-rich. Although comparison of yeast with low and high trehalose showed no observable effect on intracellular water structure, the structure of water in cells is different from that in bulk water. Cellular water exhibits a larger preference for lower-energy H-bonds or interactions over higher-energy interactions relative to that shown in bulk water. This effect is likely the result of the high concentration of biological molecules present in the cell. The ability of water to interact directly with polar groups on biological molecules may cause the preference seen for lower-energy interactions.     

Water structure in vitro and within Saccharomyces cerevisiae yeast cells under conditions of heat shock

The OH stretch mode from water and organic hydroxyl groups have strong infrared absorption, the position of the band going to lower frequency with increased H-bonding. This band was used to study water in trehalose and glycerol solutions and in genetically modified yeast cells containing varying amounts of trehalose. Concentration-dependent changes in water structure induced by trehalose and glycerol in solution were detected, consistent with an increase of lower-energy H-bonds and interactions at the expense of higher-energy interactions. This result suggests that these molecules disrupt the water H-bond network in such a way as to strengthen molecule–water interactions while perturbing water–water interactions. The molecule-induced changes in the water H-bond network seen in solution do not translate to observable differences in yeast cells that are trehalose-deficient and trehalose-rich. Although comparison of yeast with low and high trehalose showed no observable effect on intracellular water structure, the structure of water in cells is different from that in bulk water. Cellular water exhibits a larger preference for lower-energy H-bonds or interactions over higher-energy interactions relative to that shown in bulk water. This effect is likely the result of the high concentration of biological molecules present in the cell. The ability of water to interact directly with polar groups on biological molecules may cause the preference seen for lower-energy interactions.     

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.     

Water structure and interactions with protein surfaces

The structure of liquid water and its interaction with biological molecules is a very active area of experimental and theoretical research. The chemically complex surfaces of protein molecules alter the structure of the surrounding layer of hydrating water molecules. The dynamics of hydration water can be detected by a series of experimental techniques, which show that hydration waters typically have slower correlation times than water in bulk. Specific water-mediated interactions in protein complexes have been studied in detail, and these interactions have been incorporated into potential energy functions for protein folding and design. The subtle changes in the structure of hydration water have been investigated by theoretical studies.     

Some ways of looking at compensatory kosmotropes and different water environments

Hydration of macromolecular structures determines biological activity. Stabilizing solutes are kosmotropic (increase order of water) rather than chaotropic (decrease order). Preferential hydration of surfaces is a thermodynamic consequence of the solution behavior of kosmotropic solutes, but inconsistencies imply interactions such as the hydration of specific sites within macromolecules. Thermodynamic measures require bulk pure solutes; here simpler measures of the effects on bulk water, water at surfaces and hydration water of probes have been applied to solutes including natural stabilizers, analogues and example chaotropes. Changes in the near-infrared spectra, water proton NMR chemical shifts and relaxation times measure changes in the bulk liquid; HPLC-column retention of solutes indicate interactions with hydration water at different surfaces, and fluorescence probes detect effects on functional group hydration water. Ab initio calculations and Monte-Carlo simulations of the solutes in water measure the energetics of the solute-water interactions, the dipole moments of these molecules, their charge distributions and the effect of the solute molecules on the structure of water. The rankings of the test solutes by these measures are not consistent. Thus, stabilizing solutes are not interchangeable in biological systems and the intracellular replacement of one by another could affect the integration of cell metabolism.     

Atomic and residue hydrophilicity in the context of folded protein structures

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes. © 1995 Wiley-Liss, Inc.     

A model for the simulation of an aqueous dipeptide solution

A model for the simulation of a solution of an alanine dipeptide in water is presented that combines a previous model for bulk water (ST2) with that used in conformational energy studies on small molecules and proteins. The results of a pilot molecular dynamics study indicate that the model leads to reasonable solvent–solute interactions. No evidence is found for substantial changes in the structure or dynamics of the dipeptide in solution as compared to in vacuo. Furthermore, at the elevated temperature examined, there appear to be no significant effects on the dynamics or intermolecular bonding of the water molecules in contact with the solute.     

The Hofmeister effect and the behaviour of water at interfaces

Starting from known properties of non-specific salt effects on the surface tension at an air-water interface, we propose the first general, detailed qualitative molecular mechanism for the origins of ion-specific (Hofmeister) effects on the surface potential difference at an air-water interface; this mechanism suggests a simple model for the behaviour of water at all interfaces (including water-solute interfaces), regardless of whether the non-aqueous component is neutral or charged, polar or non-polar. Specifically, water near an isolated interface is conceptually divided into three layers, each layer being I water-molecule thick. We propose that the solute determines the behaviour of the adjacent first interfacial water layer (I1); that the bulk solution determines the behaviour of the third interfacial water layer (I3), and that both I1 and I3 compete for hydrogen-bonding interactions with the intervening water layer (I2), which can be thought of as a transition layer. The model requires that a polar kosmotrope (polar water-structure maker) interact with I1 more strongly than would bulk water in its place; that a chaotrope (water-structure breaker) interact with I1 somewhat less strongly than would bulk water in its place; and that a non-polar kosmotrope (non-polar water-structure maker) interact with I1 much less strongly than would bulk water in its place. We introduce two simple new postulates to describe the behaviour of I1 water molecules in aqueous solution. The first, the 'relative competition' postulate, states that an I1 water molecule, in maximizing its free energy (--delta G), will favour those of its highly directional polar (hydrogen-bonding) interactions with its immediate neighbours for which the maximum pairwise enthalpy of interaction (--delta H) is greatest; that is, it will favour the strongest interactions. We describe such behaviour as 'compliant', since an I1 water molecule will continually adjust its position to maximize these strong interactions. Its behaviour towards its remaining immediate neighbours, with whom it interacts relatively weakly (but still favourably), we describe as 'recalcitrant', since it will be unable to adjust its position to maximize simultaneously these interactions.(ABSTRACT TRUNCATED AT 400 WORDS)     

Interfacial interactions of glutamate, water and ions with carbon nanopore evaluated by molecular dynamics simulations

Molecular dynamics simulations were used to assess the transport of glutamate, water and ions (Na(+) and Cl(-)) in a single wall carbon nanopore. The spatial profiles of Na(+) and Cl(-) ions are largely determined by the pore wall charges. Co-ions are repelled whereas the counter-ions are attracted by the pore charges, but this 'rule' breaks down when the water concentration is set to a level significantly below that in the physiological bulk solution. In such cases water is less able to counteract the ion-wall interactions (electrostatic or non-electrostatic), co-ions are layered near the counter-ions attracted by the wall charges and are thus layered as counter-ions. Glutamate is concentrated near the pore wall even at physiological water concentration, and irrespective of whether the pore wall is neutral or charged (positively or negatively), and its peak levels are up to 40 times above mean values. The glutamate is thus always layered as a counter-ion. Layering of water near the wall is independent of charges on the pore wall, but its peak levels near the wall are 'only' 6-8 times above the pore mean values. However, if the mean concentration of water is significantly below the level in the physiological bulk solution, its layering is enhanced, whereas its concentration in the pore center diminishes to very low levels. Reasons for such a 'paradoxical' behavior of molecules (glutamate and water) are that the non-electrostatic interactions are (except at very short distances) attractive, and electrostatic interactions (between the charged atoms of the glutamate or water and the pore wall) are also attractive overall. Repulsive interactions (between equally charged atoms) exist, and they order the molecules near the wall, whereas in the pore center the glutamate (and water) angles are largely randomly distributed, except in the presence of an external electric field. Diffusion of molecules and ions is complex. The translational diffusion is in general both inhomogeneous and anisotropic. Non-electrostatic interactions (ion-wall, glutamate-wall or water-wall) powerfully influence diffusion. In the neutral nanopore the effective axial diffusion constants of glutamate, water and Na(+) and Cl(-) ions are all <10% of their values in the bulk, and the electrostatic interactions can reduce them further. Diffusion of molecules and ions is further reduced if the water concentration in the pore is low. Glutamate(-) is slowed more than water, and ions are reduced the most especially co-ions. In conclusion the interfacial interactions influence the spatial distribution of glutamate, water and ions, and regulate powerfully, in a complex manner and over a very wide range their transport through nanosize pores.     

Water in channel-like cavities: structure and dynamics.

Ion channels contain narrow columns of water molecules. It is of interest to compare the structure and dynamics of such intrapore water with those of the bulk solvent. Molecular dynamics simulations of modified TIP3P water molecules confined within channel-like cavities have been performed and the orientation and dynamics of the water molecules analyzed. Channels were modeled as cylindrical cavities with lengths ranging from 15 to 60 A and radii from 3 to 12 A. At the end of the molecular dynamics simulations water molecules were observed to be ordered into approximately concentric cylindrical shells. The waters of the outermost shell were oriented such that their dipoles were on average perpendicular to the normal of the wall of the cavity. Water dynamics were analyzed in terms of self-diffusion coefficients and rotational reorientation rates. For cavities of radii 3 and 6 A, water mobility was reduced relative to that of simulated bulk water. For 9- and 12-A radii confined water molecules exhibited mobilities comparable with that of the bulk solvent. If water molecules were confined within an hourglass-shaped cavity (with a central radius of 3 A increasing to 12 A at either end) a gradient of water mobility was observed along the cavity axis. Thus, water within simple models of transbilayer channels exhibits perturbations of structure and dynamics relative to bulk water. In particular the reduction of rotational reorientation rate is expected to alter the local dielectric constant within a transbilayer pore.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Molecular dynamics simulations of interaction between sub-bituminous coal and water

The high moisture content of sub-bituminous coal is associated with the interactions between coal and water. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of sub-bituminous coal according to XPS results. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxy and carbonyl is similar.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Interfacial interactions of glutamate, water and ions with carbon nanopore evaluated by molecular dynamics simulations

Molecular dynamics simulations were used to assess the transport of glutamate, water and ions (Na⁺ and Cl⁻) in a single wall carbon nanopore. The spatial profiles of Na⁺ and Cl⁻ ions are largely determined by the pore wall charges. Co-ions are repelled whereas the counter-ions are attracted by the pore charges, but this ‘rule’ breaks down when the water concentration is set to a level significantly below that in the physiological bulk solution. In such cases water is less able to counteract the ion-wall interactions (electrostatic or non-electrostatic), co-ions are layered near the counter-ions attracted by the wall charges and are thus layered as counter-ions. Glutamate is concentrated near the pore wall even at physiological water concentration, and irrespective of whether the pore wall is neutral or charged (positively or negatively), and its peak levels are up to 40 times above mean values. The glutamate is thus always layered as a counter-ion. Layering of water near the wall is independent of charges on the pore wall, but its peak levels near the wall are ‘only’ 6-8 times above the pore mean values. However, if the mean concentration of water is significantly below the level in the physiological bulk solution, its layering is enhanced, whereas its concentration in the pore center diminishes to very low levels. Reasons for such a ‘paradoxical’ behavior of molecules (glutamate and water) are that the non-electrostatic interactions are (except at very short distances) attractive, and electrostatic interactions (between the charged atoms of the glutamate or water and the pore wall) are also attractive overall. Repulsive interactions (between equally charged atoms) exist, and they order the molecules near the wall, whereas in the pore center the glutamate (and water) angles are largely randomly distributed, except in the presence of an external electric field. Diffusion of molecules and ions is complex. The translational diffusion is in general both inhomogeneous and anisotropic. Non-electrostatic interactions (ion-wall, glutamate-wall or water-wall) powerfully influence diffusion. In the neutral nanopore the effective axial diffusion constants of glutamate, water and Na⁺ and Cl⁻ ions are all < 10% of their values in the bulk, and the electrostatic interactions can reduce them further. Diffusion of molecules and ions is further reduced if the water concentration in the pore is low. Glutamate⁻ is slowed more than water, and ions are reduced the most especially co-ions. In conclusion the interfacial interactions influence the spatial distribution of glutamate, water and ions, and regulate powerfully, in a complex manner and over a very wide range their transport through nanosize pores.     

Aspects of native proteins are retained in vacuum

It is now well established that protein complexes can be maintained in the gas phase, under the reduced pressure conditions experienced within mass spectrometers. This remains an intriguing phenomenon, given the importance placed on water molecules in stabilizing protein structure. As a consequence, the overall structure of protein molecules and their complexes in the gas phase continues to stimulate much research. The eventual outcomes of this research will not only increase the number of applications of mass spectrometry to structural biology, but will also enhance our understanding of the interactions important for maintaining protein complexes in the absence of bulk solvent.     

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein–protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non‐interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water‐mediated structures may participate in the shaping of the PPI interface. The pairwise water‐mediated interaction was then investigated, and the water‐mediated residue–residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water–protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function. Proteins 2016; 84:43–51. © 2015 Wiley Periodicals, Inc.     

Structural details of urea binding to barnase: a molecular dynamics analysis.

BACKGROUND: The molecular mechanism of urea-induced protein unfolding has not been established. It is generally thought that denaturation results from the stabilizing interactions of urea with portions of the protein that are buried in the native state and become exposed upon unfolding of the protein. RESULTS: We have performed molecular dynamics simulations of barnase (a 110 amino acid RNase from Bacillus amyloliquefaciens) with explicit water and urea molecules at 300 K and 360 K. The native conformation was unaffected in the 300 K simulations at neutral and low pH. Two of the three runs at 360 K and low pH showed some denaturation, with partial unfolding of the hydrophobic core 2. The first solvation shell has a much higher density of urea molecules (water/urea ratio ranging from 2.07 to 2.73) than the bulk (water/urea ratio of 4.56). About one half of the first-shell urea molecules are involved in hydrogen bonds with polar or charged groups on the barnase surface, and between 15% and 18% of the first-shell urea molecules participate in multiple hydrogen bonds with barnase. The more stably bound urea molecules tend to be in crevices or pockets on the barnase surface. CONCLUSIONS: The simulation results indicate that an aqueous urea solution solvates the surface of a polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the protein than water does, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the protein. The results suggest that urea denaturation involves effects on both nonpolar and polar groups of proteins.     

Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme

The stochastic boundary molecular dynamics methodology is applied to the active site of the enzyme lysozyme. A comparison is made of in vacuo dynamics results from the stochastic boundary method and a full conventional molecular dynamics simulation of lysozyme. Excellent agreement between the two approaches is obtained. The influence of solvent on the residues in the active site region is explored and it is shown that both the structure and dynamics are affected. Of particular importance for the structure of the protein is the solvation of polar residues and the stabilization of like-charged ion pairs. The magnitude of the fluctuations is only slightly altered by the solvent; the overall increase in the root-mean-square fluctuations, relative to the vacuum run, is 11%. The solvent effect on dynamical properties is found not to be simply related to the solvent viscosity. Both the solvent exposure and dynamic aspects of protein-solvent interactions, including the relative time scales of the motions, are shown to play a role. The effects of the protein on solvent dynamics and structure are also observed to be significant. The solvent molecules around atoms in charged, polar and apolar side-chains show markedly different diffusion coefficients as well as exhibiting different solvation structures. One key example is the water around apolar groups, which is much less mobile than bulk water, or water solvating polar groups.     

Inspecting fluctuation and coordination around chromophore inside green fluorescent protein from water to nonpolar solvent

Green fluorescent protein (GFP) is a widely used biomarker that demands systematical rational approaches to its structure function redesign. In this work, we mainly utilized atomistic molecular dynamics simulations to inspect and visualize internal fluctuation and coordination around chromophore inside GFP, from water to nonpolar octane solvent. We found that GFP not only maintains its β-barrel structure well into the octane, but also sustains internal residue and water coordination to position the chromophore stably while suppress dihedral fluctuations of the chromophore, so that functional robustness of GFP is achieved. Our accompanied fluorescence microscope measurements accordingly confirmed the GFP functioning into the octane. Furthermore, we identified that crucial water sites inside GFP along with permeable pores on the β-barrel of the protein are largely preserved from the water to the octane solvent, which allows sufficiently fast exchanges of internal water with the bulk or with the water layer kept on the surface of the protein. By additionally pulling GFP from bulk water to octane, we suggest that the GFP function can be well maintained into the nonpolar solvent as long as, first, the protein does not denature in the nonpolar solvent nor across the polar-nonpolar solvent interface; second, a minimal set of water molecules are in accompany with the protein; third, the nonpolar solvent molecules may need to be large enough to be nonpermeable via the water pores on the β-barrel.     

SANS/SAXS study of the BSA solvation properties in aqueous urea solutions via a global fit approach

We report on the solvation properties and intermolecular interactions of a model protein (bovine serum albumine, BSA) in urea aqueous solutions, as obtained by combining small-angle neutron and X-ray scattering experiments. According to a global fit strategy, all the whole set of scattering curves are analysed by considering a unique model which includes the BSA structure, the protein-protein interactions and the thermodynamic exchange process of water/urea molecules at the protein solvent interface. As a main result, the equilibrium constant that accounts for the difference in composition between the bulk solvent and the protein solvation layer is derived. Results confirm that urea preferentially sticks to the protein surface, inducing a noticeable change in both the repulsive and the attractive interaction potentials.